克娄巴特拉计划 动漫:stormslowly :请问,在MATLAB环境下,怎么画出频率响应图啊?请帮帮我!!

来源:百度文库 编辑:杭州交通信息网 时间:2024/05/07 12:09:01
%数字低通指标:
wp=0.2*pi;
ws=0.3*pi;
Rp=1;
Rs=15;
T=1;
fs=1/T;
%模拟低通指标:
mp=0.2*pi;
ms=0.3*pi;
Rp=1;
Rs=15
[N,mc]=buttord(mp,ms,Rp,Rs,'s');
N=6;
mc=0.7081;
%设计阶数为N和截止频率为mcde的Butterworth模拟原形低通滤波器
[z,p,k]=buttap(N);
%计算单位截止频率滤波器的零点和增益;
[b,a]=zp2tf(z,p,k);
%计算单位截止频率的模拟低通滤波器的系数
[bt,at]=lp2lp(b,a,mc);
%计算截止频率为mc的模拟低通滤波器的系数
%将模拟滤波器映射位数字滤波器
[Bd,Ad]=impinvar(bt,at,fs);
%作出频率响应图

%模拟低通指标:
mp=0.2*pi;
ms=0.3*pi;
Rp=1;
Rs=15
[N,mc]=buttord(mp,ms,Rp,Rs,'s');
N=6;
mc=0.7081;

%模拟低通指标:
mp=0.2*pi;
ms=0.3*pi;
Rp=1;
Rs=15
[N,mc]=buttord(mp,ms,Rp,Rs,'s');
N=6;
mc=0.7081;

上面的ikanchi说的对,用freqz命令,具体代码如下:
[h,omega]=freqz(Bd,Ad,256); %求频率响应
plot(omega/pi,20*log10(abs(h)));grid; %作图

请参考 freqz 函数。

给你一本参考书,我现在手头上正是这本,有比较详细的例子,你先去图书馆看看能不能借到这本书,不能的话我再给你打出来吧。
《MATLAB及在电子信息课程中的应用》电子冲业出版社 陈怀琛等人编的
第114页就是频率响应的内容