加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

matlab设计切比雪夫低通滤波器

(2010-12-13 21:47:07)
标签:

matlab

滤波器

切比雪夫

chebyshev

数字信号处理

dsp

杂谈

 

绘制Chebshev 1型模拟低通滤波器的平方幅频响应曲线,阶数为2、4、6、8:

n=0:0.01:2;
for i=1:4
switch i
    case 1 ,N=2;
    case 2; N=4;
    case 3;N=6;
    case 4; N=8;
end
Rp = 1;                         %设置通滤波纹为1db
[z,p,k] =cheb1ap(N,Rp);         %设计chebyshev 1 滤波器
[b,a] = zp2tf(z,p,k);           %将零点、级点、增益形式转换为传递函数的形式
[H,w] =freqs(b,a,n);            %按n指定的频率点给出频率响应
magH2 = (abs(H)).^2;            %给出传递函数的幅度平方
posplot=['2,2',num2str(i)];     %将数字i转换为字符转,与‘2,2’合并并赋给posplot
subplot(posplot);
plot(w,magH2);
xlabel ('w/wc');

ylabel('Chebyshev 1 | H(jw)|^2');
grid on;
end

http://s2/middle/65de97a9g974b22e17731&690
Chebyshev 1行滤波器的特点:通带内具有等波纹起伏特性,而在阻带内单调下降,具有更大的衰减特性;阶数愈高,特性愈接近矩形;传递函数是没有零点的(这个就要有数学的推到了)。

绘制Chebshev 2型模拟低通滤波器的平方幅频响应曲线,阶数为2、4、6、8:

n=0:0.01:2;
for i=1:4
switch i
    case 1 ,N=2;
    case 2; N=4;
    case 3;N=6;
    case 4; N=8;
end
Rs = 16;                           %设置通滤波纹为1db
[z,p,k] =cheb2ap(N,Rs);            %设计chebyshev 1 滤波器
[b,a] = zp2tf(z,p,k);              %将零点级点增益转换为传递函数的形式
[H,w] =freqs(b,a,n);               %按n指定的频率点给出频率响应
magH2 = (abs(H)).^2;               %给出传递函数幅度平方
posplot=['2,2',num2str(i)];        %将数字i转换为字符串‘N=’合并作为标题
subplot(posplot);
plot(w,magH2);
xlabel ('w/wc');

ylabel('Chebyshev 1 | H(jw)|^2');
grid on;
end

http://s7/middle/65de97a9g974b2515da86&690

Chebyshev 1行滤波器的特点:通带内具单调的、平滑的,而在阻带内有等纹波起伏特性;阶数愈高,特性愈接近矩形;传递函数是既有极点又有零点(这个就要有数学的推到了)。

 


 

 

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有