1120108218王维韬.docVIP

《DSP原理与应用》课程 大作业 姓名 王维韬 学号 1120108218 班级 电科1112 成绩 2014年5月29日 FIR滤波器的设计 一、实验目的 熟悉数字滤波的基本原理和实现方法； 熟悉线性相位FIR数字滤波器特性； 二、实验设备 PC兼容机一台，操作系统为 Windows2000( 或 Windows98 ， WindowsXP ，以下默认为 Windows2000)，安装 Code Composer Studio 2.0 软件。 实验内容 1、利用汉宁窗设计线性相位FIR低通滤波器，要求通带截止频率=rad，N=21. 加窗得到FTR滤波器单位脉冲响应h（n）N=21;a=(N-1)/2;Wc=pi/4; 2、用MATLAB仿真得： n=[0:1:(N-1)]; m=n-a+eps; hd=sin(Wc*m)./(pi*m) W_han=(hanning(N)); h2=hd.*W_han; [H2,W]=freqz(h2,1); figure(2); subplot(211); stem(n,h2);title(汉宁窗实际脉冲响应); xlabel(n);ylabel(h(n)); subplot(212); plot(W/pi,20*log10(abs(H2)/max(H2))); xlabel(频率);ylabel(对数幅度响应); title(加汉宁窗时的频谱图); 3、ccs仿真： （1)代码 #include myapp.h #include csedu.h #include scancode.h #include math.h #define FIRNUMBER 21 #define SAMPLEF 10000 #define PI 3.1415926 float InputWave(); float FIR(); float fHn[FIRNUMBER]={0,0,0,0,-0.022,-0.026,0,0.062,0.147,0.221,0.25,0.221,0.147,0.062,0,-0.026,-0.022,0,0,0,0 }; float fXn[FIRNUMBER]={ 0.0 }; float fInput,fOutput; float fSignal1,fSignal2; float fStepSignal1,fStepSignal2; float f2PI; int i; float fIn[256],fOut[256]; int nIn,nOut; main() { nIn=0; nOut=0; f2PI=2*PI; fSignal1=0.0; fSignal2=PI*0.1; fStepSignal1=2*PI/30; fStepSignal2=2*PI*1.4; while ( 1 ) { fInput=InputWave(); fIn[nIn]=fInput; nIn++; nIn%=256; fOutput=FIR(); fOut[nOut]=fOutput; nOut++; /* break point */ if ( nOut=256 ) { nOut=0; } } } float InputWave() { for ( i=FIRNUMBER-1;i0;i-- ) fXn[i]=fXn[i-1]; fXn[0]=sin((double)fSignal1)+cos((double)fSignal2)/8.0; fSignal1+=fStepSignal1; if ( fSignal1=f2PI ) fSignal1-=f2PI; fSignal2+=fStepSignal2; if ( fSignal2=f2PI ) fSignal2-=f2PI; return(fXn[0]); } float FIR() { float fSum; fSum=0; for ( i=0;iFIRNUMBER;i++ ) { fSum+=(fXn[i]*fHn[i]); } return(fSum); } (2)仿真图 四、实验步骤 1．实验准备： -设置软件仿真模式， -启动 CCS，。 2．打开工程，工程目录光盘program里面FIR\Fir.pjt。 3．编译并下载程序。 4．打开观察窗口。 *选择菜