matlab仿真实验报告

示: b) 滤波器指标参考: 通带边缘频率为 4kHz , 阻带边缘频率为 4.5kHz , 阻带衰减大于 50dB ;其实是低通滤波器 c) Matlab 函数 y= awgn(x,snr,'measured') , 首先测量输入信号 x 的功率, 然后对其叠加高斯白噪声; 2 、实验目的: d) 通过本次实验,掌握以下知识: i. FIR 数字滤波器窗口设计法的原理和设计步骤; ii. Gibbs 效应发生的原因和影响; iii. 不同类型的窗函数对滤波效果的影响,以及窗函数和长度 N 的选择。实验分析先确定滤波器的设计指标: 滤波器通带边缘为 4KHZ ,阻带边缘频率为 4.5Khz ,阻带衰减大于 50db 。使用汉明窗。加噪声时使用 awgn 函数。使用作图函数分别画出加噪声后的声音信号和去噪声后的声音信号的时域和频域波形。使用 matlab 库函数 hamming 来进行滤波。再作出滤波后的图形。最后作出滤波器的幅频特性图。代码及注释 fs=32000; [x,fs]=audioread('sound.wav'); snr=20; % 信噪比为 20dB x2=awgn(x,snr,'measured','db'); % 叠加高斯白噪声? audiowrite('sound_noise.wav',x,fs); t=0:1/fs:(size(x2)-1)/fs; wp=8000*pi/32000; ws=9000*pi/32000; wdelta=ws-wp; N=ceil(11*pi/wdelta); % 取整 wn=(ws+wp)/2; b=fir1(N,wn/pi,blackman(N+1)); % 选择窗函数,并归一化截止频率 figure(1) freqz(b,1,512) f2=filter(b,1,x2); title(' 滤波器幅频、相频特性');