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

H．

clc; b=[1 0]; a=[1 -0.9]; figure(1); zplane(b,a);

title('零极点图'); figure(2); H=freqz(b,a); magH=abs(H); angH=angle(H); figure(2); stem(magH);

title('频率响应的幅值谱'); figure(3); stem(angH);

title('频率响应的相位谱'); x=[1,zeros(1,100)]; y=filter(b,a,x); figure(4); stem(y); xlabel('n'); ylabel('h(n)');

title('单位脉冲响应');

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

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

I.

n=[];

b=[2.2403 2.4908 2.2403]; a=[1 -0.4 0.75]; x1=[1,2,3,1,0]; x2=[0,1,2,3,1,0]; y1=filter(b,a,x1)

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

y2=filter(b,a,x2) y=[0,y1] figure(1); stem(y1) xlabel('n'); ylabel('y1(n)'); figure(2); stem(y2) xlabel('n'); ylabel('y2(n)'); figure(3); stem(y) xlabel('n');

ylabel('y(n)');