数字信号处理课设.docx

数字信号处理课设.docx

《数字信号处理课设.docx》由会员分享，可在线阅读，更多相关《数字信号处理课设.docx（10页珍藏版）》请在冰点文库上搜索。

数字信号处理课设

【2】

clear;closeall;

wp=0.2;ws=0.3;Rp=1;Rs=25;

[n,Wn]=buttord（wp,ws,Rp,Rs）;

[b,a]=butter（n,Wn）;

Freqz（b,a,600,1）;

双线性:

clear;closeall;

wp=0.2;ws=0.3;Rp=1;Rs=25;

[n,Wn]=buttord（wp,ws,Rp,Rs）;

[b,a]=butter（n,Wn,'s'）;

freqs（b,a）

[bz,az]=bilinear（b,a,1）;

freqz（bz,az,600,1）

脉冲：

wp=0.2;ws=0.3;Rp=1;Rs=25;

[n,Wn]=buttord（wp,ws,Rp,Rs,'s'）;

[b,a]=butter（n,Wn,'s'）;

freqs（b,a）

[bz,az]=impinvar（b,a,1）;

freqz（bz,az,600,1）;

【3】

带同：

T=10us;fp1=20khz;fp2=30khz;fs1=15khz;fs2=35khz;ap=1db;as=40db;

T=10*10^（-6）;fp1=20000;fp2=30000;

fs1=15000;fs2=35000;ap=1;as=40;

Wp1=fp1*T;Wp2=fp2*T;Ws1=fs1*T;Ws2=fs2*T;

fp=[0.2,0.3];fs=[0.15,0.35];

[n,Wn]=buttord（fp,fs,1,40）;

[b,a]=butter（n,Wn）;

Freqz（b,a,600,1）;

>>T=10*10^（-6）;fp1=20000;fp2=30000;

fs1=15000;fs2=35000;ap=1;as=40;

Wp1=fp1*T;Wp2=fp2*T;Ws1=fs1*T;Ws2=fs2*T;

fp=[0.2,0.3];fs=[0.15,0.35];

>>[n,Wn]=buttord（fp,fs,1,40,'s'）;

>>[b,a]=butter（n,Wn,'s'）;

>>freqs（b,a）

[bz,az]=impinvar（b,a,1）;freqz（bz,az,600,1）

>>T=10*10^（-6）;fp1=20000;fp2=30000;

fs1=15000;fs2=35000;ap=1;as=40;

Wp1=fp1*T;Wp2=fp2*T;Ws1=fs1*T;Ws2=fs2*T;

fp=[0.2,0.3];fs=[0.15,0.35];

[n,Wn]=buttord（fp,fs,1,40,'s'）;

[b,a]=butter（n,Wn,'s'）;freqs（b,a）

>>[bz,az]=bilinear（b,a,1）;freqz（bz,az,600,1）