samp5_12.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 12 行

%Samp5_12
b=[0.2 0.3 1];   %滤波器传递函数分子多项式系数
a=[1 0.4 1]; %滤波器传递函数分母多项式系数
figure(1),freqs(b,a);  %第一种输出方法
[h,w]=freqs(b,a); %计算滤波器的复数频率响应
mag=abs(h);pha=angle(h);  %得到滤波器的幅频和相频响应
figure(2),subplot(2,1,1),loglog(w,mag);  %运用双对数坐标绘制幅频响应
grid on;xlabel('角频率/rad\cdots^-^1');ylabel('振幅');
subplot(2,1,2),semilogx(w,pha*180/pi) %运用半对数坐标绘相频响应
grid on;xlabel('角频率/rad\cdots^-^1');
ylabel('相位/^o');