ex1_5.m

来自「数字信号处理 实验 的matlab 源码」· M 代码 · 共 12 行

%1-5
clear
clf
R=51;d=0.8*(rand(1,R)-0.5);m=0:R-1;
s=2*m.*(0.9.^m);x=s+d;
subplot(211),plot(m,d,'r-',m,s,'g--',m,x,'b-.')
legend('d[n] ','s[n] ','x[n] '),xlabel(' 时间序号 n'),ylabel('振幅')
x1=[0 0 x];x2=[0 x 0];x3=[x 0 0];
y=(x1+x2+x3)/3;
subplot(212),plot(m,y(2:R+1),'-r',m,s,'g--')
legend('y[n]','s[n]'),xlabel(' 时间序号 n'),ylabel('振幅')