数字信号处理实验3数字滤波器的应用李永杰
—— 数字信号处理 —— 李 永 杰 实验3 数字滤波器的应用
·目的：给定IIR数字滤波器的系统函数从频谱图上观察滤波器的特性；进行数字滤波；·内容：(1)编写程序观给定的滤波器的频谱(2)根据给定的滤波器的频率特性，对构造的信号进行滤波
• 目的：给定IIR数字滤波器的系统函数， 从频谱图上观察滤波器的特性；进行 数字滤波； • 内容： （1）编写程序观给定的滤波器的频谱。 （2）根据给定的滤波器的频率特性，对 构造的信号进行滤波
(一)滤波器1已知低通DF的3dB带宽频率。=0.2元，阻带起始频率の，=0.4元，阻带衰减20lgH(ews)=-15dBH(elo1-a1通带阻带过:渡a2十带0West元
（一）滤波器1 已知低通DF的3dB带宽频率 ，阻带起始 频率 ，阻带衰减 c = 0.2 s = 0.4 H e dB s j w 20lg ( ) = −15
0.0343 +0.10289z-1 +0.10289z-2 + 0.0343z-3H(z) =1.8952 - 3.3358z-1 + 2.242z-2 - 0.527z0.102890.102890.034301;b=[0.034302.2420a=[1.8952-3.3358-0.52701;; %计算频谱,w=2*pi/100[h,w]=freqz(b,a,100);plot(w/pi, abs(h));0.90.80.70.60.50.40.30.20.100.20.10.30.40.50.60.70.80.90
1 2 3 1 2 3 1.8952 3.3358 2.242 0.527 0.0343 0.10289 0.10289 0.0343 ( ) − − − − − − − + − + + + = z z z z z z H z b=[0.03430 0.10289 0.10289 0.03430]; a=[1.8952 -3.3358 2.2420 -0.5270]; [h,w]=freqz(b,a,100); %计算频谱, w=2*pi/100 plot(w/pi, abs(h) ); 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
n=0:99;x signal = cos(n*pi*0.1):x noise = 2*cos(n*pi*0.5);x = x signal + x noise;figure(1); plot(n,x signal,'b',n,x,'r'):title(原始信号为蓝色，含噪声的信号为红色");0.10289b=[0.034300.102890.03430];-3.33582.2420a=[1.8952-0.5270];y=filtfilt(b, a, x); %滤波figure(2); stem(abs(fft(x);title(滤波前信号的频谱");figure(3);stem(abs(fft(y));title(滤波后信号的频谱);figure(4);plot(n,x_signal,'b',n,y,'r');title(原始信号为蓝色,滤波得到的信号为红色');
n=0:99; x_signal = cos(n*pi*0.1); x_noise = 2*cos(n*pi*0.5); x = x_signal + x_noise; figure(1); plot(n,x_signal,'b',n,x,'r'); title('原始信号为蓝色，含噪声的信号为红色'); b=[0.03430 0.10289 0.10289 0.03430]; a=[1.8952 -3.3358 2.2420 -0.5270]; y=filtfilt(b, a, x); %滤波 figure(2); stem(abs(fft(x)));title('滤波前信号的频谱'); figure(3);stem(abs(fft(y)));title('滤波后信号的频谱'); figure(4);plot(n,x_signal,'b',n,y,'r'); title('原始信号为蓝色,滤波得到的信号为红色');