Matlab论坛
400万+工程师在用
华为|鸿蒙开发者日
直播报名
400万+工程师在用
华为|鸿蒙开发者日
直播报名

送你一把好柴刀

4年用户 71经验值
私信 关注
[问答]

请问一下这是切比雪夫带通滤波器的图像吗

2016-12-17 22:45

程序应该是对的,我修改一下技术指标, 发现出来的 带通滤波器图像很好,但用下面的指标出来的图像感觉不像带通滤波
wp1=0.25*pi;wp2=0.65*pi;wp0=0.45*pi;ws1=0.2*pi;ws2=0.8*pi;Rp=1;As=60;
T=0.001;
Wp1=(2/T)*tan(wp1/2);Wp2=(2/T)*tan(wp2/2);
Ws1=(2/T)*tan(ws1/2);Ws2=(2/T)*tan(ws2/2);Wp=(2/T)*tan(wp0/2);
BW=Wp2-Wp1;
W0=Wp1*Wp2;
WP=1;
WS=WP*(W0^2-Ws1^2)/(Ws1*BW);
[N,Wn]=cheb1ord(WP,WS,Rp,As,'s');
[B1,A1]=cheby1(N,Rp,Wn,'s');
[h1,w1]=freqs(B1,A1);
subplot(3,2,1);
plot(w1/pi,20*log(abs(h1)));grid on;
xlabel('w(rad)');
ylabel('|H(jw)|.dB');
title('模拟低通滤波器幅频特性曲线');
[B2,A2]=lp2bp(B1,A1,Wp,BW);
[h2,w2]= freqs(B2,A2);
subplot(3,2,2);
plot(w2,20*log(abs(h2)));
axis([0,3,-400,50]);grid on;
xlabel('w(rad)');
ylabel('|H(jw)|.dB');
title('模拟带通滤波器幅频特性曲线');
[B3,A3]=bilinear(B2,A2,0.5);
[h,w]=freqz(B3,A3,64);
phz=unwrap(angle(h));
subplot(3,2,3);
plot(w/pi,abs(h));
xlabel('w(rad)');ylabel('|H(z)|');
title('数字带通滤波器幅频特性曲线');
subplot(3,2,4);
plot(w/pi,20*log10(abs(h)));
axis([-1,2,-250,50]);
xlabel('w(rad)');ylabel('|H(z)|.dB');
title('数字带通滤波器幅频特性曲线');
subplot(3,2,5);
plot(w/pi,phz);grid on;
xlabel('w(rad)');ylabel('H(z)');
title('数字带通滤波器相频特性曲线');

  • 2.PNG
已退回5积分

回帖(8)

王栋春

2016-12-17 23:01:11
感觉是的  这个资料不错

送你一把好柴刀

2016-12-18 09:19:00
原来是T取得太小的缘故

王栋春

2016-12-18 19:54:49
引用: 送你一把好柴刀 发表于 2016-12-18 09:19
原来是T取得太小的缘故

楼主能讲一下这个滤波器的用处吗  对这一块的知识不了解  谢谢

送你一把好柴刀

2017-1-21 21:00:10
引用: 王栋春 发表于 2016-12-18 19:54
楼主能讲一下这个滤波器的用处吗  对这一块的知识不了解  谢谢

这个程序还需要改动一些地方,这就是把一个频率范围内的信号以增益为1的方式通过,而对不在这个范围的信号进行抑制

王栋春

2017-1-21 21:08:58
引用: 送你一把好柴刀 发表于 2017-1-21 21:00
这个程序还需要改动一些地方,这就是把一个频率范围内的信号以增益为1的方式通过,而对不在这个范围的信号进行抑制

楼主说的是带通滤波器吧

送你一把好柴刀

2017-1-24 20:09:35
引用: 王栋春 发表于 2017-1-21 21:08
楼主说的是带通滤波器吧

没错,做的就是一个带通滤波器

王栋春

2017-1-24 20:29:56
引用: 送你一把好柴刀 发表于 2017-1-24 20:09
没错,做的就是一个带通滤波器

好的  你所说的那个外国名字我第一次听说 见笑了

丁飞

2017-9-27 10:50:58
想问下楼主这个滤波器物理实现了没?

更多回帖

打开APP