matlab 高手请进来看看 关于matlab 巴特沃斯带通及低通滤波器的问题w1=(fc-600000)*2*pi*0.00000007;w2=(fc+600000)*2*pi*0.00000007;[a,b]=butter(5,[w1 w2],'bandpass');h2=filter(a,b,h1); y=h2.*c; %完成相干解调wp=10100*2*pi/50000

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/29 20:35:42

matlab 高手请进来看看 关于matlab 巴特沃斯带通及低通滤波器的问题w1=(fc-600000)*2*pi*0.00000007;w2=(fc+600000)*2*pi*0.00000007;[a,b]=butter(5,[w1 w2],'bandpass');h2=filter(a,b,h1); y=h2.*c; %完成相干解调wp=10100*2*pi/50000
matlab 高手请进来看看 关于matlab 巴特沃斯带通及低通滤波器的问题
w1=(fc-600000)*2*pi*0.00000007;
w2=(fc+600000)*2*pi*0.00000007;
[a,b]=butter(5,[w1 w2],'bandpass');
h2=filter(a,b,h1);
y=h2.*c; %完成相干解调
wp=10100*2*pi/5000000; %Butterworth低通滤波器
ws=20500*2*pi/5000000;
ap=0.2;
as=1;
[N,wn]=buttord(wp,ws,ap,as);
[a,b]=butter(N,wn);
s=filter(a,b,y);
请问这两段程序的意思是什么啊 有点看不懂,

matlab 高手请进来看看 关于matlab 巴特沃斯带通及低通滤波器的问题w1=(fc-600000)*2*pi*0.00000007;w2=(fc+600000)*2*pi*0.00000007;[a,b]=butter(5,[w1 w2],'bandpass');h2=filter(a,b,h1); y=h2.*c; %完成相干解调wp=10100*2*pi/50000
wp=2*pi*6000;ws=2*pi*12000; rp=3; as=25;
dispf=14; dispa=-90;
grid on;
axis([0,dispf,dispa,5]);
[n,wc]=buttord(wp,ws,rp,as,'s');
[B,A]=butter(n,wc,'s');
k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;
Hk=freqs(B,A,wk);
[n1,wp1]=cheb1ord(wp,ws,rp,as,'s');
[B1,A1]=cheby1(n1,rp,wp1,'s');
k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;
Hk1=freqs(B1,A1,wk);
[n2,wp2]=cheb2ord(wp,ws,rp,as,'s');
[B2,A2]=cheby2(n2,as,wp2,'s');
k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;
Hk2=freqs(B2,A2,wk);
[n3,wp3]=ellipord(wp,ws,rp,as,'s');
[B3,A3]=ellip(n3,rp,as,wp3,'s');
k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;
Hk3=freqs(B3,A3,wk);
plot(fk/1000,20*log10(abs(Hk)),fk/1000,20*log10(abs(Hk1)),':',fk/1000,20*log10(abs(Hk2)),'--',fk/1000,20*log10(abs(Hk3)));
xlabel('频率(kHz)','FontSize',16,'FontWeight','bold');
ylabel('幅度(dB)','FontSize',16,'FontWeight','bold');
title('4种滤波器性能比较','FontSize',16,'FontWeight','bold');
legend('巴特沃兹滤波','切比雪夫I滤波','切比雪夫II滤波','椭圆滤波','Location', 'East')
grid;