数字信号处理实验九剖析.doc

实验报告 课程名称：数字信号处理 实 验 九：信号滤波 班 级：通信1403 学生姓名：强亚倩 学 号：1141210319 指导教师：范杰清 华北电力大学（北京） 实验目的 综合应用信号频谱分析和数字滤波器设计的知识，实现信号的滤波。加深理解信号时域和频域分析的物理概念，理解设计指标的工程概念，认识不同类型滤波器的特性和适用范围。 二、实验原理 首先对待滤波的信号进行频谱分析，观察信号频率分布的规律，从而确定数字滤波器的类型（FIR滤波器、IIR滤波器、自适应滤波器、小波滤波器等）。 在加性噪声的情况下，若信号的频谱与噪声的频谱基本不重叠，可以采用频率选择滤波器（FIR滤波器、IIR滤波器）。 若信号的频谱与噪声的频谱重叠较多，可以采用自适应滤波、小波滤波等。 若为乘性噪声，可以根据同态滤波的原理对信号进行预处理，然后再按照加性噪声的情况处理。 在确定了数字滤波器类型后，还需要根据信号时域特性、频域特性、或时频特性确定滤波器的设计参数，设计出相应的数字滤波器。 最后，利用该数字滤波器对信号进行滤波，在时域和频域观察信号滤波的主观及客观效果。若主观及客观效果满足要求，说明分析过程和滤波方法正确有效，若不满足要求，需要重新分析和设计。 三、 实验内容 1.对周期方波信号进行滤波。 （1）生成一个基频为10Hz的周期方波信号x(t)。 （2）设计一数字滤波器，滤去该周期信号中40Hz以外的频率成分，观察滤波前后信号时域波形及频谱。 （3）若该信号x(t)淹没在噪声中（随机噪声用randn((1,N)生成），试用filter函数滤去噪声。 滤波前的时域波形 Fs=1000; T=1/10; t=0:1/Fs:1; x=square(2*pi*10*t); figure plot(t,x) title(矩形波) xlabel(时间(秒)) ylabel(函数值(福特)) axis([0,0.2,-1.5,1.5]) 滤波前的频谱 N=T*Fs; t1=0:1/Fs:1/Fs*(N-1) x1=square(2*pi*10*t1); f=(-N/2)*Fs/N:Fs/N:((N/2-1))*Fs/N; fx=fftshift(fft(x,N))/N; figure plot(f,abs(fx)); title(方波信号的频谱) xlabel(频率) ylabel(函数值fx) 滤波后的时域波形 [n,Wn]=buttord(40/500,50/500,3,15); [b,a]=butter(n,Wn); y=filter(b,a,x); figure plot(t,y) title(矩形波通过低通滤波器后的时域波形) xlabel(时间(秒)) ylabel(函数值(福特)) axis([0,0.2,-1.5,1.5]) 滤波后的频谱 N=512; f=(-N/2)*Fs/N:Fs/N:(N/2-1)*Fs/N; fy=fftshift(fft(y,N))/N; figure plot(f,abs(fy)) title(滤波后信号的频谱) xlabel(频率) ylabel(函数值fx) figure freqz(b,a); title(低通滤波器的频率特性) 滤除噪声 noise=randn(1,100); HC=x1+noise; f=(-N/2)*Fs/N:Fs/N:(N/2-1)*Fs/N; HC=fftshift(fft(HC,N))/N; figure plot(f,abs(HC)); title(滤除噪声之前) xlabel(频率f) ylabel(幅度) HC1=filter(b,a,HC); figure plot(f,abs(HC1)); title(滤除噪声之后) xlabel(频率f) ylabel(幅度) 2.若原始落信号由5Hz、15Hz、30Hz三个幅度相等的正弦信号构成。分别设计一个FIR和IIR数字滤波器滤除5Hz和30Hz频率成分。 （1）绘制原始信号时域波形和幅度频谱。 （2）分别设计FIR和IIR数字滤波器，滤绘制其幅频特性。 （3）利用设计的滤波器对信号进行滤波，绘制输出信号的时域波形和幅度频谱。 时域波形 （1）k=0:255; x=cos(2*pi*5*k*0.01)+cos(2*pi*15*k*0.01)+cos(2*pi*30*k*0.01); figure stem(x) X=fft(x); figure plot(abs(X)) FI