%《數字信號處理教程——MATLAB釋義與實現》第二章例2.3.3程序hc233
% 信號合成
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
disp('例2.3.3 求下列序列的結果及圖形。')
disp('(a)	n1=2:8,x1=sin(n1);n2=-1:3;x2=ones(1,5);求x1與x2之和')
disp('(b)	求x1與x2之和')
disp('(c)	n3=0:4;x3=[6,5,4,3,2]; 求它對n=0點折疊后的序列y3。')
disp('(d)	求把(b)中的原序列及折疊后的序列y3合在一起，構成對稱序列y4。')
disp('(e)	求按序列y3構成的反對稱序列y5。')
clear
n1=2:8;x1=sin(n1);n2=-1:3;x2=ones(1,5);             % 給出序列x1和x2
[y1,ny1]=seqadd(x1,n1,x2,n2);                       % 序列x1,x2相加，得出y1
[y2,ny2]=seqmult(x1,n1,x2,n2);                      % 序列x1,x2相乘，得出y2
n3=0:4;x3=[6,5,4,3,2]; [y3,ny3]=seqfold(x3,n3);     % 序列x3折疊,得出y3
ny4=-4:4; y4=[y3(1:4),x3];                          % 構成y4
subplot(2,2,1); stem(ny1,y1,'.');                   % 繪制y1
title('例 2.3.3(a) 的結果圖')
ylabel('y1(n)');line([-5,10],[0,0])
subplot(2,2,2);stem(ny2,y2,'.');                    % 繪制y2
title('(b) 的結果圖')
ylabel('y2(n)');
set(gcf,'color','w')                                % 置圖形背景色為白
subplot(2,2,3);stem(ny3,y3,'.');title('(c) 的結果圖')       % 繪制y3
ylabel('y3(n)');%
subplot(2,2,4);stem(ny4,y4);title('(d) 的結果圖'),hold on   % 繪制y4
y5=[-y3(1:4),0,x3(2:5)];stem(ny4,y5,'*:')
ylabel('y4(n)');line([-5,5],[0,0])