基于matlab的数字信号分析.doc
题目二、已知有一频率为50赫兹幅度为1V的方波信号被频率为1600赫兹幅度为2V的正弦波干扰,请设计方案将方波信号提取出来,要求被衰减的信号要衰减60DB以上,同时对原来的信号和提取后的信号进行谐波分析,并对信号采用MATLAB进行仿真验证。提示用FIR线性相位滤波器FS4000FP100FS500WPFPPI/FS/2WSFSPI/FS/2滤波器边界频率WIDEWSWP过渡带宽NCEIL12PI/WIDE求解滤波器的最小阶数WNWPWS/2截止频率取通带和阻带边界频率的中点BFIR1N,WN/PI,BLACKMANN1FIGURE1滤波器幅频响应曲线H,FFREQZB,1,1024,FS求出滤波器的频率响应,采样频率是20000HZSUBPLOT211PLOTF/FS/2,20LOG10ABSHXLABEL F/HZ YLABEL 振幅/DB TITLE BLACKMAN窗函数设计的滤波器幅频响应曲线 SUBPLOT212PLOTF/FS/2,180/PIUNWRAPANGLEHXLABEL F/HZ YLABEL 相位 TITLE 滤波器相频相应曲线 FIGURE2滤波之前信号的时域和频域图像F150F21600SUBPLOT211T01/FS01XSQUARE2PIF1T,502SIN2PIF2TPLOTT,XXLABEL T/S YLABEL 幅度 TITLE 滤波前的数字信号时域图像 SUBPLOT212N1512FFT采用的数据长度YFFTX,N1进行快速傅里叶变换MAGABSYF0N11FS/N1真实频率PLOTF1N1/2,MAG1N1/22/N1XLABEL F/HZ YLABEL 振幅 TITLE 滤波前数字信号的频率域图像 FIGURE3输出通过滤波器后的时域和频域的图像KFILTERB,1,XSUBPLOT211PLOTT,KXLABEL T/S YLABEL 振幅 TITLE 滤波后的数字信号时域图像 SUBPLOT212N11024FFT采用的数据长度YFFTK,N1进行快速傅里叶变换MAGABSYF0N11FS/N1真实频率PLOTF1N1/2,MAG1N1/22/N1XLABEL F/HZ YLABEL 振幅 TITLE 滤波后数字信号的频率域图像 基于MATLAB设计的数字FIR滤波器图像分析由谐波分析可以知道方波是由多个频率的正弦波叠加而成,而在利用数字滤波器滤除相对于50HZ方波信号的1600HZ正弦波信号时,必然会滤除方波的高次谐波。所以,经过滤波器后的波形并非标准的方波而有波动。这可以从下面的图3中看到,经FFT后得到的幅频响应曲线。且滤波后存在一定的延迟由图3可以看出。这个时间由设计的数字滤波器阶数确定。