搜索
32阶数字滤波器设计课程名称：信号检测基础实验名称：32阶数字滤波器设计实验日期：2019年3月18日一、输入:首先进行滤波前信号分析，做出信号1(s1)信号2(s2)原图，并做出信号1与2相加的被污染的信号(s)f=2e6;fm=7e6;fs=40e6;T=15;t=0:1/fs:T/f;s1=sin(2*pi*f*t);s2=sin(2*pi*fm*t);s=s1+s2;plot(s,'b');holdon;plot(s1,'r');holdon;plot(s2,'g');legend('被污染的信号','信号1','信号2');title('滤波前信号分析');二、32阶数字带通滤波器设计：滤波器的阶数（n阶滤波器，specifyorder＝n-1，这里为31），窗函数选用Kaiser，截止频率分别为2MHz和5MHz，采样频率40MHZ（也就是ADC的采样频率）。三、滤波：使用file菜单---generatemfile，命名为bandfilter做出滤波后的图形：figure(2);Hd=bandfilter;output=filter(Hd,s);plot(output);title('滤波后的波形');附录：滤波器源码functionHd=bandfilterFs=40;%SamplingFrequencyN=31;%OrderFc1=2;%FirstCutoffFrequencyFc2=5;%SecondCutoffFrequencyflag='scale';%SamplingFlagBeta=0.5;%WindowParameter%Createthewindowvectorforthedesignalgorithm.win=kaiser(N+1,Beta);%CalculatethecoefficientsusingtheFIR1function.b=fir1(N,[Fc1Fc2]/(Fs/2),'bandpass',win,flag);Hd=dfilt.dffir(b);