数字信号处理实验的代码,仅参考

B

N=50; Nx=10; n=0:N-1; h=(0.9).^n;

R=sign(sign(Nx-1-n)+1); w=conv(h,R); stem(w);

title('求h(n)与Rn的卷积'

);

c.

N=20;

b=[1,-1,0.9]; a=[1,0,0];

h1=[1 zeros(1,N-1)]; h2=ones(1,N); n=0:N-1;

y1=filter(a,b,h1);%单位冲激响应 y2=filter(a,b,h2);%单位阶跃响应 subplot(1,2,1); stem(n,y1);

legend('单位冲激响应'); subplot(1,2,2); stem(n,y2);

legend('单位阶跃响应');