%《數字信號處理教程——MATLAB釋義與實現》第二章例2.2.2程序hc222
% 離散序列的合成
% 電子工業出版社出版  陳懷琛編著 2004年9月
% 
disp('例題2.2.2：畫出以下各序列在給定區間的波形圖')
disp('a. x(n)=1.5*delta(n+1)-delta(n-4).  -4<=n<=5')
disp('b. x(n)=n*(u(n)-u(n-8))-10*exp(-0.3*(n-10))*(u(n-10)-u(n-16)). 0<=n<=20')
disp('c. x(n)=cos(0.07*pi*n)+0.2*w(n)). 0<=n<=30')
disp('d. b中函數x(n)的能量E(n)')
clf
% 對題(a)
n1=[-4:5]; x1=1.5*impseq(-1,-4,5)-impseq(3,-4,5);           % 列出x1序列
subplot(2,2,1); stem(n1,x1,'.'); title('例2.2.2a 的序列圖')
ylabel('x1(n)'); axis([-5,5,-2,3]);text(5.5,-2,'n')
% 對題(b)
n2=[0:20];x21 = n2.*(stepseq(0,0,20)-stepseq(8,0,20));          % 列出x21序列
x22 = 10*exp(-0.3*(n2-10)).*(stepseq(10,0,20)-stepseq(16,0,20));% 列出x22序列
x2 = x21-x22;                                           % x2序列是x21和x22之和
subplot(2,2,2);stem(n2,x2,'.');
title('例2.2.2b 的序列圖')
ylabel('x2(n)');
% 對題(c)
n3 = [0:30]; x3 = cos(0.07*pi*n3)+0.2*randn(size(n3));     % 列出x3序列
subplot(2,2,3);stem(n3,x3,'.');title('例2.2.2c 的序列圖')
ylabel('x3(n)');axis([0,30,-1.4,1.4]);text(53,-1.4,'n')
% 對題(d)
n4=n2; E=x2.^2; 
subplot(2,2,4);stem(n4,E,'.');title('例2.2.2d 的序列圖')   % 列出x4序列
ylabel('E(n)');         %
set(gcf,'color','w')                        % 置圖形背景色為白