maltab文件.rar_基于matlab函数稳定性判别

%使用双线性法设计IIR数字低通滤波器 %脉冲响应不变法的主要缺点是会产生频谱混叠 %设计步骤： %先要确定数字低通滤波器的技术指标便带频率wp，最呆最大衰减ap,阻带截止频率ws，阻带最小衰减as %将数字低通滤波器的技术指标转化为相应模拟低通滤波器的技术指标，这里主要是边界频率wp和ws的转化，ap和as指标不变 %按照双线性变化法公式进行wp,ws的转变 %按照模拟低通滤波器的技术指标设计过度模拟低通滤波器 %最后用转换的方法,将模拟滤波器H_a(s)转化为数字滤波器H(Z) %要求：设计低通数字滤波器，要求频率低于 0.2pi（rad）时，容许幅度误差在1dB以内； %在频率为0.3pi到pi之间的阻带衰减大于15dB.用双线性法设计巴特沃斯低通滤波器 T=1;Fs=1/T; rp=1;rs=1; wpz=0.2;wsz=0.3; %设置技术指标 wp=2*tan(wpz*pi/2); ws=2*tan(wsz*pi/2); %利用双线性法公式（预畸变校正转换指标） [N,wc]=buttord(wp,ws,rp,rs,'s'); %设计过渡模拟滤波器 [B,A]=butter(N,wc,'s'); [Bz,Az]=bilinear(B,A,Fs); %用双线性法转化为数字低通滤波器 fs=ws/(2*pi); fk=0:fs/512:0.9; wk=2*pi*fk Hk=freqs(B,A,Fs); subplot(2,1,1); plot(fk,20*log10(abs(Hk)));