hc233.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——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])