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

E.(稳态响应不会求)

clc;

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

[H,w]=freqz(b,a,100); magH=abs(H); phaH=angle(H); subplot(2,1,1); plot(w,magH);

title('频率响应的幅值'); subplot(2,1,2); plot(w,phaH);

title('频率响应的相位

');

F．

clc;

Dt=0.00005;

t=-0.005:Dt:0.005; xa=exp(-1000*abs(t)); Wmax=2*pi*2000; K=500; k=0:1:K;

W=k*Wmax/K;

Xa=xa*exp(-j*t'*W)*Dt; Xa=real(Xa);

W=[-fliplr(W),W(2:501)]; Xa=[fliplr(Xa),Xa(2:501)]; subplot(4,2,1); plot(t*1000,xa); ylabel('xa(t)');