有哪些信誉好的足球投注网站- 1、本文档共7页，可阅读全部内容。
- 2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
源程序 数字信号处理课程设计
1 IIR数字巴特沃思滤波器设计（冲激不变法）
function IIR_Imp
%
% Impulse Invariance Transformation
% Butterworth Lowpass Filter Design
%
wp=0.2*pi;
ws=0.3*pi;
Rp=1;
As=15;
T=1;
OmegaP=wp*T;
OmegaS=ws*T;
ep=sqrt(10^(Rp/10)-1);
Ripple=sqrt(1/(1+ep*ep));
Attn=1/(10^(As/20));
%
% Analog Butterworth Prototype Filter Calculation
[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);
[b,a]=imp_invr(cs,ds,T);
[C,B,A]=dir2par(b,a);
%
figure(1);
[db,mag,pha,grd,w]=freqz_m(b,a);
subplot(2,2,1);plot(w/pi,mag);title(Magnitude Response);
xlabel(1frequency in pi units);ylabel(/H/);axis([0,1,0,1.1]);
set(gca,XTickMode,manual,XTick,[0,0.2,0.3,1]);
set(gca,YTickMode,manual,YTick,[0,Attn,Ripple,1]);grid;
subplot(2,2,3);plot(w/pi,db);title(Magnitude in dB);
xlabel(2frequency in pi units);ylabel(decibels);
axis([0,1,-40,5]);
set(gca,XTickMode,manual,XTick,[0,0.2,0.3,1]);
set(gca,YTickMode,manual,YTick,[-50,-15,-1,0]);grid;
set(gca,YTickLabelMode,manual,YTickLabels,[50,15,1,0]);
subplot(2,2,2);plot(w/pi,pha/pi);title(Phase Response);
xlabel(3frequency in pi units);ylabel(pi units);
axis([0,1,-1,1]);
set(gca,XTickMode,manual,XTick,[0,0.2,0.3,1]);
set(gca,YTickmode,manual,YTick,[-1,0,1]);grid;
subplot(2,2,4);plot(w/pi,grd);title(Group Delay);
xlabel(4frequency in pi units);ylabel(Samples);
axis([0,1,0,10]);
set(gca,XTickMode,manual,XTick,[0,0.2,0.3,1]);
set(gca,YTickMode,manual,YTick,[0:2:10]);
grid;
2 IIR数字巴特沃思滤波器设计（双线性变换法）
function IIR_BIL
%
%利用双线性法设计Butterworth数字滤波器
%
wp=0.2*pi;
ws=0.3*pi;
Rp=1;
As=15;
T=1;
Fs=1/T;
OmegaP=(2/T)*tan(wp/2);
OmegaS=(2/T)*tan(ws/2);
ep=sqrt(10^(Rp/10)-1);
Ripple=sqrt(1/(1+ep*ep));
Attn=1/(10^(As/20));
%
[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);
[b,a]=bilinear(cs,ds,T);
[C,B,A]=dir2cas(b,a);
figure(1);
[db,mag,pha,grd,w]=freqz_m(b,a);
subplot(2,2,1);plot(w/pi,mag);title(Magnitude Response)
xlabel(frequency in pi units);ylabel(|H|);
axis([0,1,0,1.1])
set(gca,XTickMode,manual,XTick,[0,0.2,0.3,1]);
set(gca,YTickMode,manual,YTick,[0,Attn,Ripple,1]);grid
subplot(2,2,3);plot(w/pi,db);title(Magnitud
文档评论（0）