信号与线性系统实验报告(9)

magh=20*log10(magh); % 以分贝 magh(zerosIndx)=-inf; angh=angle(Hz);

angh=unwrap(angh)*180/pi; % 角度换算 figure

subplot(1,2,1) plot(w,magh); grid on

%set(H(2,2),'xlim',[0,1])

xlabel('特征角频率(\times\pi rads/sample)') title('幅频特性曲线 |H(w)| (dB)'); subplot(1,2,2) plot(w,angh); grid on

xlabel('特征角频率 (\times\pi rads/sample)') title('相频特性曲线 \theta(w) (degrees)');

这里H s

1

。其中m=0.3 2

1 m ms

仿真结果为：