心得体会：

课程设计时间虽短，但是这次我也基本熟悉了一种新的集成开发环境CCS，学习新的知识的过程也是自己学习能力培养与提升的过程。

ＴＤＳ５１０ＵＳＢ 仿真器驱动程序的安装以及相应的配置流程也有了一定的了解。仿真环境的配置，到工程的建立，文件的加载，到程序的仿真，与目标板的链接与调试，整个过程在摸索中逐渐熟悉。对已有程序进行修改，重复相应的过程也能实现预定的功能，在短短的时间里能掌握这些基本就差不多了，由于有的需要配置的文件的缺失，无法完成对相应工程的配置设置，所以采用的参考例程里的程序，完成整个过程，这也是一个学习的过程。做项目不是一个人的事，每个成员都应积极的参与，为整个项目的完成提供保障，团队的协作，尽可能的去发挥每个成员的专长，在整个项目的完成都能有所收获，这才应该是做项目的真正目的，加强同学之间的交流，用心付出，共同享受带给大家的成功的喜悦。相关课程的学习只是个基础，在此基础之上对相应的硬件与软件结合，切实去体验一个芯片所能实现的各种功能，去发现所学的知识会在哪些方面用到，是如何应用，有怎样可以改进的方法，更深层次去掌握跟其他相关课程的交叉点，提升学习能力，从近期来看，可以为我们将要开始的毕业设计做准备，当做是一次练手，争取出色完成毕业设计，为四年的大学交出一份完美的答卷。从长远看，为自己以后的工作也在一定程度的奠定基础，学习能力强了，自己就能比较快的接受新知识，更能适应社会对人的要求

。

相信团队的力量，同时也要提高个人解决问题的能力，让自己在团队中发挥的作用，将个人融入团队中，才能让自己有更大的收获。好好珍惜每次锻炼学习的机会，不断提升自己，不断超越自己，成就人生美好的梦想！

专业班级：微电子学

姓 名

学 号：

第二篇：DSP课程设计2

一、DSP设计题目

    设计合适的巴特沃斯、切比雪夫IIR滤波器、窗函数FIR滤波器分别对同一段语音进行信号处理。

二、 DSP设计要求

1.运用双线性变换法设计高通巴特沃斯滤波器：

性能指标：fst=4000 Hz, fc=5000 Hz, Rst=30 dB, Rp＝5 dB

2.运用冲击响应不变法设计带通切比雪夫滤波器：

     性能指标：fc1=2000 Hz, fc2=3000 Hz, fst1=1000 Hz, fst2=4000 Hz, Rst=30 dB, Rp＝5 dB。  

3.选择合适的窗函数设计带阻滤波器：

     性能指标：fst1=2000 Hz, fst2=3000 Hz, fc1=1000 Hz, fc2=4000 Hz, Rst=30

dB, Rp=5 dB

4.了解和掌握用MATLAB实现IIR和FIR滤波器的设计方法、过程，为以后的设计打下良好基础。

三、DSP设计原理

（1）双线性变换法高通巴特沃斯数字滤波

     通过所给的技术指标我们可得数字域技术指标：通带截止频率wc和阻带截止频率wst，通过阻带和通带的最小衰减，运用buttord函数,得到返回阶数N及3dB的频率Wn。运用butter就可以得到高通巴特沃斯滤波器的系统函数系数AB, 用freqz对系统函数进行抽样，就可以得到高通巴特沃斯滤波器的抽样频谱，

（2）冲击响应不变法带通切比雪夫数字滤波器

     通过所给的技术指标可以得到模拟滤波器的技术指标，从而用过cheb1ord函数可以得到N和Wn，cheby1函数可以得到得到模拟域的系统函数，通过冲击响应不变法impinvar就可以得到带通切比雪夫数字滤波器的系统函数的系数

（3）选择合适的窗函数设计FIR数字滤波器

    通过阻带最小衰减30dB及过渡带宽要求选择汉宁窗及N的的大小。通过fir1函数当中选择’stop’得到窗函数的系统函数的系数，从而就可以得到FIR带阻数字滤波器

四、DSP源程序清单

（1）语音信号的采集

[y,fs]=wavread('d:\333.wav',[1000 60000]); %利用函数wavread对语音信号进行采样，记住采样频率和采样点数。

Sound(y,fs);                          %通过sound函数播放该文件的声音：

Y=fft(y,1000);                       %对语音号进行快速傅里叶变换，得到信号的频谱

%特性Y

Subplot(432);plot(y);title('滤波前的信号波形');

Subplot(433);plot(abs(Y));title('滤波前的信号频谱');

（2）双线性变换法高通巴特沃斯数字滤波

wp=2*pi*5000/fs;                  %数字域的通带截止频率

wst=2*pi*4000/fs;                  %数字域的阻带截止频率

Rp=5;

Rst=30;                           %阻带最小衰减

[N,Wn]=buttord(wp/pi,wst/pi,Rp,Rst);  %返回阶数N及3dB的频率Wn.

[B,A]=butter(N,Wn,'high');           %根据N,阻带最小衰减Rp及3dB的频率%Wn,求高通数字巴特沃斯滤波器的系统%函数.得到系统函数的系数A.B

[h,w]=freqz(B,A,256);               %对系统函数进行抽样

Subplot(434);

plot(w/pi,20*log10(abs(h)));     

title('高通巴特沃斯滤波器的抽样频谱');%标题

xlabel('f in pi unit');                  %横坐标的名字

ylabel('gain in dB');                  %纵坐标的名字

axis([0,1,-50,10]);                   %对图形X轴 (0,1),Y轴(-50,10)部分进行%放大

x=filter(B,A,y);                     %利用函数filter对信号进行滤波

X=fft(x,1000);                    %对滤波后的语音号再次进行快速傅里叶%变换，得到信号的频谱特性X

Subplot(435);plot(x);title('高通巴特沃斯滤波后的信号波形');

Subplot(436);plot(abs(X));title('高通巴特沃斯滤波后的信号频谱');

 (3)冲击响应不变法带通切比雪夫数字滤波器

 Wc=[2*pi*2000 2*pi*3000];            %模拟的通带截止频率

Wst=[2*pi*1000 2*pi*4000];            %模拟的阻带截止频率

Rst=30;                             %最小衰减

Rp=5;

 [N,Wn]=cheb1ord(Wc,Wst,Rp,Rst,'s');    %切比雪夫模拟滤波器的技术指标

 [B,A]=cheby1(N,Rp,Wn,'s');

 [b,a]=impinvar(B,A,fs);                %冲击响应不变法求系统函数的系数

 [h,w]=freqz(b,a,256);

Subplot(437);

plot(w/pi,20*log10(abs(h)),'-'); title('切比雪夫滤波器的抽样频谱');

xlabel('f in pi');ylabel('gain in db');axis([0,1,-50,10])

x=filter(b,a,y);

X=fft(x,1000);

Subplot(438);plot(x);title('切比雪夫滤波后的信号波形');

%利用函数filter对信号进行滤波

Subplot(439);plot(abs(X));title('切比雪夫滤波后的信号频谱');

% 对滤波后的语音号再次进行快速傅%里叶变换，得到信号的频谱特性X

(4) 选择合适的窗函数设计FIR数字滤波器

wst1=2*pi*2000/fs;                         %数字域的阻带截止

wst2=2*pi*3000/fs;       

wc1=2*pi*1000/fs;                          %数字域的通带截止

wc2=2*pi*4000/fs;

N=max(ceil(3.1*pi/(wst1-wc1))*2+1,ceil(3.1*pi/(wc2-wst2))*2+1);

wst=[(wst1+wc1)/2,(wst2+wc2)/2];

b=fir1(N-1,wst/pi,'stop',Hanning(N));          %根据阻带最小衰减30dB,选择

%汉宁窗。且’stop’表示为设计为

%带阻，返回fir滤波器的系统函

%数

[h,f]=freqz(b,1,500);              

title('汉宁窗带阻FIR');

Subplot(4310);

plot(f/pi,20*log10(abs(h)),'-');ylabel('20log|H|');xlabel('w/pi');title('汉宁窗带阻FIR');                                     

>> axis([0 1 -100 10]);

>> grid

x=filter(b,1,y);

X=fft(x,1000);

Subplot(4311);plot(x);title('经汉宁窗FIR滤波后的信号波形');

Subplot(4312);plot(abs(X));title('经汉宁窗FIR滤波后的信号频谱');

五、DSP图形

                                                                                                                                                                                                                                                                                                                

（6）收获和体会

     在这次DSP设计主要是设计滤波器对音频信号进行滤波处理。这次的任务是通过所给的技术指标设计三个不同的滤波器：高通巴特沃斯滤波器、带通切比雪夫滤波器、FIR数字滤波器。

将所要处理的音频信号分别通过三个滤波器，得到滤波后的信号图形和信号频谱。进行二次比较：(1)滤波前和滤波后的信号频谱和图形。(2)信号通过三种滤波器后的信号频谱和图形。

我们通过图形可以看到通过高通巴特沃斯数字滤波器后滤波效果特别明显，信号图形基本滤完，基本只剩高频部分。通过切比雪夫数字滤波器的效果一般，基本能把低频和高频滤去。FIR滤波器效果就不是特别明显。

在设计FIR滤波器时，我们用到fir1(n,wn,'type',window)实现加窗线性相位FIR滤波器。在设计高通或带阻滤波器时若N-1为奇数时，fir1函数会自动加一使其成为偶数，这样就保证了h(n)有奇数个，从而满足线性相位高通和带阻FIR滤波器的设计要求。

在这次操作当中，开始我在快速变换的函数(Y=fft(y,50000))我取了50000个点，看图时发现语音信号通过各个设计好的滤波器后的效果不是很容易分辨，后来听老师说可以试着把快速傅立叶变换的点数取小一点，这样效果比较明显，我就把点数取成1000，确实减少抽样点数后效果明显了。

我个人觉得这次的设计让我们学的课程有了实际性的作用，比如录歌曲的时候如果有了杂音，就可以自己试着设计一个滤波器，提高歌曲的音质。我想在破案电视上的那个滤去杂音获取破案人员所需的声音并且放大，其中肯定也有一部分用到了滤波器。

  

（7）参考文献

(1)MATLAB信号处理  刘波  电子工业出版社

(2)Matlab信号处理与应用  董长虹  国防工业出版社

(3)数字信号处理原理及其MATLAB实现 从玉良  电子工业出版社

(4)数字信号处理基础及MATLAB实现 周辉 中国林业出版社

《DSP课程设计》

实验报告

姓    名:         

学    号:                

专    业:  电子信息科学类    

班    级:            

指导老师：