如何用matlab进行滤波器设计和仿真.docVIP

华 中 科 技 大 学 本科课程考试答题本 考生姓名 考生学号 班 级 联系电话 Email： 考试科目 考试日期 一、报告摘要： 本报告首先论述如何用matlab进行滤波器设计和仿真，并且将参数转换成DSP可用的滤波器参数。然后，将这组参数在CCS上进行纯软件仿真。随之是本次实验过程中代码拼接的问题与解决方法。软硬件流程图及核心代码的陈述。将实际滤波的结果进行展示。与理想滤波器的比较之间的差异与改良。最后是对于本课程的感想 二、滤波器基础知识以及基于matlab的参数设计： 设计一个能够将1000Hz一下保留，1500Hz以上的频率滤出的滤波器，利用理想的fir滤波器截止频率可以设定为（1000+1500）/2=1250Hz，利用哈明窗设计本滤波器，阶数为21阶， function[a]=ccc() clear all; clf; N=1024;%设置采样点数 fs=6000;%设置采样频率 dt=1/fs; for k=1:N; f1=1500; f2=1000; y(k)=sin(2*pi*f1*k*dt)+sin(2*pi*f2*k*dt)+sin(2*pi*(f1+200)*k*dt);%+sin(2*pi*1250*k*dt); end wp=pi*1000*2/6000;%归一化的通带频率 ws=pi*1500*2/6000;%归一化的阻带频率 wc=(wp+ws)/2; %理想的滤波器的截止频率 M=21; %设置阶数 b=fir1(M-1,wc/pi,hamming(M))%;利用fir1函数得到滤波系数b [H,F]=freqz(b,1,N); % figure(1); plot(F/pi,20*log10(abs(H)));grid;xlabel(f(Hz));ylabel(幅度);%画出幅频特性 yy1=conv(b,y); %做卷积即为滤波 y=fft(y,N); %做fft运算得到频域特性 pyy=y.*conj(y); %乘以共轭，得到的是y的幅值的平方 f=(0:(N/2-1)); figure(2); plot(f*fs/N,pyy(1:N/2)); y=fft(yy1,N); pyy=y.*conj(y); f=(0:(N/2-1)); figure(3); plot(f*fs/N,pyy(1:N/2)); 执行以上函数后得到以下图像： 可以看到在（2*1000）/6000=0.33处的衰减基本不变，而在（1000+1500）/6000=0.42处衰减值已经大于3dB了，而在（2*1500）/6000=0.5衰减值已经接近20dB了。显然本滤波器与理想的滤波器有些差异，当然随着阶数的提高会更加接近理想的滤波器。 滤波前后的变化，可以清楚地看到1000Hz处的幅值并没有什么衰减，而在1500Hz以上的分量基本全部被滤除（在1500Hz只衰减不到20dB）。 滤波系数： A0 A1 A2 A3 A4 A5 A6 A7 A8 A9 A10 0.0013 -0.0026 -0.0058 0.0032 0.0211 0.0089 -0.0470 -0.0608 0.0726 0.3007 0.4169 A11 A12 A13 A14 A15 A16 A17 A18 A19 A20 0.3007 0.0726 0.0608 -0.0470 0.0089 0.0211 0.0032 -0.0058 -0.0026 0.0013 转换后得到的滤波系数如下： a0 002AH a1 -0054H a2 -00BEH a3 0068H a4 02B4H a5 0124H a6 -0605H a7 -07C9H a8 094CH a9 267DH a10 355CH a11 267DH ; a12 094CH a13 -07C9H a14 -0605H a15 0124H a16 02B4H a17 0068H a18 -00BEH a19 -0054H a20 002AH 三、算法验证： 将matlab得到的的滤波器系数转换成十六进制数，并设计一个21阶的fir滤波器在ccs2.0