MATLAB信号处理工具箱幻灯片.ppt

例7.16 二阶系统 MATLAB程序如下： a=[-0.55, -0.78;0.78, 0];b=[1;0]; c=[5.96 6.45];d=[0]; x0=[1;0];t0=0:0.1:20; initial(a, b, c, d, x0, t0); title(LTI系统的零输入响应曲线) 运行结果如图7.15所示。 3．滤波函数filter MATLAB信号处理工具箱中提供了滤波函数filter。 格式：y=filter(B, A, x) 功能：对向量x中的数据进行滤波处理 例7.18 设系统差分方程为 MATLAB程序如下： B=1; A=[1, -0.8]; n=0:31; x=0.8.^n; y=filter(B, A, x); subplot(2, 1, 1);stem(x) subplot(2, 1, 2);stem(y) 程序运行结果如图7.17所示。 例7.20 利用函数impinvar将一个低通模拟滤波器转变成数字滤波器，采样频率为10 Hz。 程序运行结果如图7.20所示。 例7.21 一个三阶模拟Butterworth低通滤波器的传递函数为 MATLAB程序如下： num=1; ％ 模拟滤波器系统函数的分子 den=[1, sqrt(3), sqrt(2), 1]; ％ 模拟滤波器系统函数的分母 [num1, den1]=bilinear(num, den, 1) ％ 求数字滤波器的传递函数 程序运行结果为 num1 = 0.0533 0.1599 0.1599 0.0533 den1 = 1.0000 -1.3382 0.9193 -0.1546 例7.24 设计一个在通带内的最大衰减为3 dB，在阻带内的最小衰减为40 dB的4阶低通模拟椭圆滤波器原型。 MATLAB程序如下： n=4; rp=3; rs=40; [z, p, k]=ellipap(n, rp, rs); [b, a]=zp2tf(z, p, k); w=logspace(-1, 1); freqs(b, a) 运行结果如图7.24所示。 例7.27 用矩形窗、三角窗、汉宁窗、哈明窗分别设计低通数字滤波器。信号的采样频率为1000 Hz，数字滤波器的截止频率为100 Hz，滤波器的阶数为80。 MATLAB程序如下： passrad=0.2*pi; w1=boxcar(81); w2=triang(81); w3=hanning(81); w4=hamming(81); n=1:1:81; hd=sin(passrad*(n-41))./(pi*(n-41)); hd(41)=passrad/pi; h1=hd.*rot90(w1); h2=hd.*rot90(w2); h3=hd.*rot90(w3); h4=hd.*rot90(w4); [MAG1, RAD]=freqz(h1); [MAG2, RAD]=freqz(h2); [MAG3, RAD]=freqz(h3); [MAG4, RAD]=freqz(h4); subplot(2, 2, 1); plot(RAD, 20*log10(abs(MAG1))); grid on subplot(2, 2, 2); plot(RAD, 20*log10(abs(MAG2))); grid on subplot(2, 2, 3); plot(RAD, 20*log10(abs(MAG3))); grid on subplot(2, 2, 4); plot(RAD, 20*log10(abs(MAG4))); grid on 程序运行结果如图7.28所示。 利用在MATLAB设计IIR数字滤波器可分以下几步来实现 (1)按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标； (2)根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率Wc； (3)利用最小阶数N产生模拟低通滤波原型； (4)利用截止频率Wc把模拟低通滤波器原型转换成模拟低通、高通、带通或带阻滤波器； (5)利用冲激响应不变法或双线性不变法把模拟滤波器转换成数字滤波器。 2．MATLAB的直接设计 图7.25 IIR数字滤波器MATLAB设计步骤流程图 [例7-25] 试设计一个带阻IIR数字滤波器，其具体的要求是：通带的截止频率：wp1＝650Hz、wp2＝850Hz；阻带的截止频率：ws1＝700Hz、ws2＝800Hz；通带内的最大衰减为rp＝0.1dB；阻带内的最小衰减为rs＝50dB；采样频率为Fs＝2000Hz。MATLAB源程序设计如下： w