[n1,wn1]=buttord(wp,ws,Ap,As);%估计巴特沃什滤波器的阶数 %wn1为归一化截止频率
[num,den]=butter(n1,wn1);%设计巴特沃什滤波器 %计算单位为dB的gain函数
%在单位圆上将函数转化为256个等份 w=0:pi/255:pi;%计算相位 h=freqz(num,den,w);
g=20*log10(abs(h));%计算幅度 pha=angle(h); %画图
subplot(1,2,1);plot(w/pi,g);grid%画带通滤波器的幅度响应曲线图 axis([0 1 -60 5]);
xlabel('\omega/\pi');ylabel('gain in db');%横坐标表示频率,纵坐标表示分贝数 title('巴特沃什带通滤波器幅频响应');%写标题
subplot(1,2,2);plot(w/pi,pha);grid%画带通滤波器的相位响应曲线图 xlabel('\omega/\pi');ylabel('phase'); title('巴特沃什带通滤波器相位响应');
3.1.3IIR数字滤波器MATLAB仿真结果与分析
可以得出:
IIR数字滤波器的阶数为:n1=7
3dB截止频率为:wn1=0.4284π,0.6693π 系统函数H(z)的分子系数为:
num=0.0003,0,-0.0019,0,0.0057,0,-0.0095,0.0095,0,-0.0057,0,0.0019,
0,-0.0003
系统函数H(z)的分母系数为:
den=1.0000,1.7451,4.9282,6.1195,9.8134,9.2245,10.4323,7.5154,6.4091,
3.4595,2.2601,0.8470,0.4167,0.0856,0.0299 根据运行结果可写出IIR滤波器的系统函数为
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说综合文库基于MATLAB的数字滤波器设计与仿真(很全--免费下载)(9)在线全文阅读。
相关推荐: