基于MATLAB的FIR滤波器的设计及应用(图像去噪).docVIP

数字信号处理课程设计报告书 课题名称 MATLAB的FIR滤波器的设计及应用（图像去噪） 姓 名 学 号 院、系、部 专 业 指导教师 基于MATLAB的FIR滤波器的设计及应用（图像去噪） 一、实验目的 1．学会MATLAB的使用，掌握MATLAB的程序设计方法； 2．掌握数字信号处理的基本概念、基本理论和基本方法； 3. 在MATLAB环境下产生噪声信号 4．掌握MATLAB设计FIR数字滤波器的窗函数设计方法； 5．学会用MATLAB对进行分析和处理。FIR滤波器的原理 如果所希望的滤波器的理想的频率响应函数为，则其对应的单位脉冲响应为 （.1） 窗函数设计法的基本原理是用有限长单位脉冲响应序列逼近。由于往往是无限长序列，而且是非因果的，所以用窗函数将截断，并进行加权处理，得到： （.2） 就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为 （.3） 式2.3)中，N为所选窗函数的长度。 我们知道，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见表。 各种窗函数的基本参数 窗函数 旁瓣峰值幅度/dB 过渡带宽 阻带最小衰减/dB 矩形窗 -13 4π/N -12 三角形窗 -25 8π/N -25 汉宁窗 -31 8π/N -44 哈明窗 -41 8π/N -53 莱克曼窗 -57 12π/N -74 凯塞窗(α=7.865) -57 10π/N -80 这样选定窗函数类型和长度N之后，求出单位脉冲响应，并按照式（.3）求出。是否满足要求，要进行演算。一般在尾部加零使长度满足2的整数次幂，以便用FFT计算。如果要观察细节补零点数增多即可 2．FIR滤波器的原理 用表示希望逼近的幅度特性函数，要求设计线性相位FIR数字滤波器时，必须满足线性相位约束条件。用表示实际设计的滤波器幅度特性函数。定义加权误差函数为式中，成为误差加权函数，用来控制不同频段（一般指通带和祖带）的逼近精度。等波纹最佳逼近基于切比雪夫逼近，在通带和阻带以的最大值和最小值为准则，采用Remez多重交代迭代算法求解滤波器系数。在Remez多重交换迭代过程中是确知函数。 利用等波纹最佳逼近准则设计线性相位FIR数字滤波器数学模型的建立及其求解算法的推导复杂，所以求解必须借助MATLAB信号处理工具箱函数remezord和remez，只要简单的调用这两个函数就可以完成线性相位FIR数字滤波器的等波纹最佳逼近设计。 如果要求线性相位特性，则还必须满足 根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类，例如，要设计线性相位低通特性，可以选择这一类，而不能选择这一类。三、实验仪器及材料 计算机、MATLAB71教学版。 四、实验内容 1噪声信号的频谱分析2.设计数字滤波器和画出频率响应 利用窗函数的方法设计FIR滤波器；最后，利用MATLAB中的函数freqz画出各滤波器的频率响应。 3用滤波器对噪声进行滤波 4.比较几种滤波器去滤波前后噪声的波形及频谱 五、程序设计 1.先原始图像，然后加噪由各项参数求出单位脉冲响应和其频率响应用布莱克曼窗函数法等波纹最佳逼近法I = imread(eight.tif); J = imnoise(I,gaussian,0.1); figure(1) ;imshow(I);title(原始图像); figure(2);imshow(uint8(J));title(加噪图像); %窗函数法 fp=1000;fs=2000;Fs=10000; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; B=ws-wp;N=ceil(12*pi/B);wc=(wp+ws)/2/pi; hn=fir1(N-1,wc); n=0:N-1; figure(3) subplot(1,2,1);stem(n,hn,.);title(h(n)波形); grid on; [hw,w]=freqz(hn,1); subplot(1,2,2); plot(w/pi,20*log10(abs(hw)));title(频率响应hw);grid on; axis([0 1 -100 4]); yt=filter2(hn,J); figure(4);imshow(yt,[]);title(滤波后图像1); %用等波纹最佳逼近方法绘制滤波器损耗函数以