课 程 设 计
课程设计名称：_数字信号处理课程设计 专 业 班 级 ： 学 生 姓 名 ： 学 号 ： 指 导 教 师 ： 课程设计时间：2014-6-16至2014-6-20
说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页
一、需求分析： 滤波器从广义上来说对特定的频点或频点以外的频率进行有效滤波的电路，这种电路保留输入信号中的有用信息，滤除不需要的信息，从而达到信号的检测、提取、识别等不同的目的。如果处理的信号是时域离散信号，那么相应的处理系统就称为数字滤波器，由于在实际工作中被处理的信号都是幅度量化的数字信号，因此，数字滤波器实际上是用有限精度的算法实现一个线性时不变的时域离散系统。目前，数字滤波器的应用越来越广泛，它已深入到很多领域，如图象处理、医学生物信息处理、地质信号处理和模式识别处理等。
数字滤波器的种类很多，分类方法也不同，可以从功能上分类，也可以从实现方法上及设计方法上来分类等等。滤波器在功能上总的可分为四类，即低通(LP)、高通(HP)、带通(BP)、带阻(BS)滤波器等，从实现方法上，由有限长冲激响应的数字滤波器被称为FIR滤波器，具有无限长冲激响应的数字滤波器称作IIR滤波器。
切比雪夫滤波器的幅度特性具有等纹波特性。他有两种类型，一种是通带内为等纹波的，在阻带内是单调的成为切比雪夫Ι型滤波器；一种是通带内单调，阻带内等纹波的，称为切比雪夫ΙΙ型。本实验采用切比雪夫Ι型滤波器。
利用双线性变换法将模拟传输信号Ha（s
）变换为数字传输函数G（z），从而是z域的数字传输函数保留s域的模拟传输函数的基本性质。设计成的IIR数字低通滤波器能够去掉信号中不必要的高频成分，降低采样频率，避免频率混淆，去掉高频干扰。
二、设计原理及设计思路：
滤波器的设计流程图如图1所示
图1 滤波器 设计流程方框图
包括在此说明每个部分的算法设计说明（可以是描述算法的流程图）；每个程序中使用的存储结构设计说明（如果题目已经指定了数据存储的，按照指定的设计，并且写出该存储结构的定义）。也就是简单介绍一下你程序中自定义的函数（函数名，参数，以及功能）等。
具体实现方法如下：
（1）确定数字低通滤波器的技术指标：通带截止频率ωp、通带衰减?p、阻带截止频率ωs、阻带衰减?s
切比雪夫滤波器的振幅平方特性如图2所示：
图2 切比雪夫滤波器的振幅平方特性
（2）将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。如果采用脉冲响不变法，边界频率的转换关系为: ??? T
如果采用双线性变换法，边界频率的转换关系为
21 ??tg(?)T2
(3) 按照模拟低通滤波器的技术指标设计模拟低通滤波器。
(4) 利用双线性变换法将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。
（5）数字低通技术指标为：
ωp=0.4πrad, ? p=1dB; ωs=0.5πrad, ? s=40Db
（6）模拟低通的技术指标为：
21tan?p,T?1T2?P?2tan0.1??0.65rad/s,?p?1dB?p??s?2tan0.15??1.019rad/s,?s?15dB
归一化截止角频率wp=2pi*Fs/Ft; ws=2pi*Fs/Ft
(7)利用模拟切比雪夫滤波器设计数字滤波器。通带截止频率为：wp=0.4*pi; 阻带截止频率为：ws=0.5*pi;通带最大衰减为：Rp=1;阻带最大衰减为：As=15;设定周期为1s；模拟低通滤波器的生成：[b,a]=cheby1(n,1,Wn,'low','s'); 满足设计指标的最小阶数n和截止频率为：
Wn[n,Wn]=cheb1ord(OmegaP,OmegaS,1,40,'s')。最后实现输入输出、幅频特性、相频特性的图形，理论结合实验对比分析，得出心得体会。
三、运行环境：
PC机，windows2000及其以上，matlab软件。
Matlab一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序（M文件）后再一起运行。使之更利于非计算机专业的科技人员使用。而且这种语言可移植性好、可拓展性极强。
MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用
的算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。MATLAB的这些函数集包括从最简单最基本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数。函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等。
四、开发工具和编程语言：
MATLAB 软件，编程语言为DSP设计语言。
五、程序源代码：
%实验设计程序如下：
clf;
ft=2000;fp=400;fs=500;
wp=2*pi*fp/ft;ws=2*pi*fs/ft;
omegp=tan(wp/2);omegs=tan(ws/2);
omega=1;omegb=omegp/omegs;
[n,wn]=cheb1ord(omega, omegb, 1, 40,'s');
[b,a] = cheby1(n,1, wn, 's');
[num,den]=bilinear(b,a,0.5);
[z,p,k]=tf2zp(num,den);
[h,w]=freqz(num,den,256);
g = 20*log10(abs(h));
figure(1)
subplot(2,1,1);
plot(w/pi,g);grid on;
axis([0 1 -60 1]);
xlabel('频率，Hz'); ylabel('增益， dB');
title('切比雪夫1型低通滤波器');
subplot(2,1,2);
zplane(z,p);
xlabel('z'); ylabel('p');
title('传输零极点'); f1=300,f2=600;
n=0:600;
t=n/10000;
x1=sin(2*pi*f1*t); x2=sin(2*pi*f2*t); x=x1+x2;
figure(2)
subplot(2,2,1);
plot(x1);grid on;
axis([0,50*pi,-3,3]); xlabel('t');ylabel('x1'); title('x1的波形'); subplot(2,2,2);
plot(x2);grid on;
axis([0,50*pi,-3,3]); xlabel('t');ylabel('x2'); title('x2的波形'); subplot(2,2,3);
plot(x);grid on;
axis([0,50*pi,-3,3]); xlabel('t');ylabel('x'); title('输入x的波形'); y=filter(num,den,x); subplot(2,2,4)
plot(y);grid on;
axis([0,50*pi,-3,3]); xlabel('t');ylabel('y'); title('滤波器输出y的波形');
六、调试分析：
编写程序有一定难度，调试是不断出错，更多的错误是：函数未定义，函数功能调用错误等，不过经过耐心的学习，最终还是调试除了比较理想的结果，如有不妥，还望指正。调试好后，发现没有输入的x1、与x2波形，才知道我们编程的时候，忘记了将其实现。最后我们又编程实现了此功能，见后文。
七、测试结果：
图三 实验程序截图
运行结果如下：
图四 切比雪夫1型低通滤波器及其传输零极
图五 滤波器输入x、x1、x2与输出y的波形
八、参考资料：
参考文献格式如下：([1]书籍 [2] 文章例)
[1] 边肇祺，模式识别（第二版），北京：清华大学出版社,1988，25～35
[2] 李永忠，几种小波变换的图像处理技术，西北民族学院学报（自然科学版），2001.6,22（3），15～18
[3]程佩青著，《数字信号处理教程》，清华大学出版社，2001
[4]Sanjit K. Mitra著，孙洪，余翔宇译，《数字信号处理实验指导书（MATLAB版）》，电子工业出版社，20xx年1月
[5]郭仕剑等，《MATLAB 7.x数字信号处理》，人民邮电出版社，20xx年
[6]胡广书，《数字信号处理 理论算法与实现》，清华大学出版社，20xx年
[7]刘泉，阙大顺编 数字信号处理原理与实现，电子工业出版社，20xx年6月
[8]苏金明， 王永利编 Matlab7.0使用指南， 电子工业出版社，20xx年11月
[9]薛定宇 陈阳泉编 基于matlab/simulink的系统仿真技术与应用，清华大学出版社 20xx年
[10]余成波等编，数字信号处理及MATLAB实现（第二版），北京：清华大学出版社，2008
[11]陈怀琛编著，数字信号处理教程——MATLAB释义与实现，北京：电子工业出版社2004
九、心得体会:
本次课程设计中，我的任务是用双线性变换法设计切比雪夫Ι型低通数字滤波器，并应用MATLAB软件进行仿真分析。这是数字信号处理的知识与软件应用相结合的一次训练，我受益匪浅。
切比雪夫滤波器包括两种，即切比雪夫1型滤波器和2型滤波器，通带内为等纹波的，在阻带内是单调的即为切比雪夫Ι型滤波器。通过本次课程设计，我对平时所学的知识有了更深入的了解，对MATLAB的应用也较以前熟练，知识与应用相结合提高了我的学习热情，在实验中，我确实遇到了一些问题，但通过我们的共同努力，最终还是解决了问题，收获了宝贵的经验。我相信，在下一阶段的学习中，我会取得更优异的成绩。
第二篇：数字信号处理课程设计报告(样本)
数字信号处理
课程设计报告
课题名称：语音信号的处理与滤波
系 别：
学 号：
姓 名：
班 级：
指导教师：
20##年10月25日
目录
摘录…………………………………………………………….3
正文…………………………………………………………….4
一、设计内容…………………………………………………..4
二、设计原理…………………………………………………..4
三、设计过程………………………………………………..…7
1、语音信号的录制…………………………………………....7
2、取点频谱分析……………………………………………...8
3、低通滤波器设计…………………………………………..10
4、滤滤前后信号对比………………………………………..12
结论…………………………………………………………...14
参考文献……………………………………………………...15
摘要
随着计算机和信息科学的飞速发展，信号处理逐渐发展成为一门独立的学科，成为信息科学的重要组成部分，在语音处理、雷达、图像处理、通信、生物医学工程等众多领域中得到广泛应用。 Matlab语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，Matlab功能强大、简单易学、编程效率高，深受广大科技工作者的喜爱。特别是Matlab还具有信号分析工具箱，不需具备很强的编程能力，就可以很方便地进行语音信号分析、处理和设计。 数字信号处理课程在现代科学中具有很大重要性及自身特点，理解与掌握课程中的基本概念、基本原理、基本分析方法，对用Matlab进行数字信号处理课程设计的思路，具有很大帮助。 语音信号的处理与滤波是数字信号处理课程中常出现的课题，也是现代科学中值得深入究研的一个课题。 语音信号的处理与滤波的设计主要是用Matlab作为工具平台，设计中涉及到声音的录制、播放、存储和读取，语音信号的抽样、频谱分析，滤波器的设计及语音信号的滤波，通过数字信号处理课程的理论知识的综合运用。从实践上初步实现对数字信号的处理。
关键词：抽样频率；滤波器；抽样点；频谱分析
正文
一、设计内容：
1、熟悉并掌握MATLAB中有关声音（wave）录制、播放、存储和读取的函数。
2．在MATLAB环境中，使用声音相关函数录制2秒左右自己的声音，抽样率是8000Hz/s。（考虑如何解决一个实际问题：录制刚开始时，常会出现实际发出声音落后录制动作半拍的现象，如何排除对这些无效点的采样？）
3．分别取8000个和16000个数据进行频谱分析，得到幅度和相位谱，比较二者异同并分析原因。
4．针对电话信道（最高3500Hz），设计一个FIR或IIR滤波器进行滤波，把抽样率转变为7000Hz/s，并进行频谱分析，得到幅度和相位谱。
5．把处理后的所有数据储存为声音文件，与原始声音进行比较。
二、设计原理
1、 语音信号的采样
l 语音信号的录制
语音信号的录制可通过MATLAB工具也可以通过windows自带的录音机录制一段自己的声音。在MATLAB里面，用函数wavrecord录制一段时间自定的声音，用函数wavwrite存储原始语音信号。其调用格式如下：
X=wavrecord(n ,fs ,ch ,dtype) ；
式中n 为抽样点数，决定了录音长度；
fs为抽样频率，默认值为11025Hz，还可根据要求自己选择合适的抽样率；
ch为声道数，默认值为1，表示单声道，如果指定为2，则抽样为双声道立体声数据；
dtype为抽样数据的存储格式，用字符串指定，可以是‘double’、‘single’、‘int16’、‘int8’，指定存储格式的同时也就规定了每个抽样值量化的精度，int8 对应8 位精度抽样，其它都是16 位抽样精度。
语音信号的存储可用函数wavwrite实现,其调用格式如下:
wavwrite(X,Fs,'f:\课程设计\sound.wav');
式中X为所录制的语音信号，Fs为抽样频率，引号内表示存储路径
l 语音信号的频谱分析
对语音信号进行频谱分析，是认识语音信号和处理语音信号的重要方法。下面将采用频谱来分析语音信号。
对一个语音时域离散信号{x(n)}，其频谱函数
是x(n)的傅里叶变换。傅里叶变换定义为：
其中j为虚数单位，
称为数字频率。
由于处理的信号是有限长的，即
，故实际采用的是离散傅里叶变换DFT 。 长度为N的序列x(n)，其DFT定义为：
，其中
如果x(n)为一个周期序列，得到的X(k)为x(n)的周期频谱。若x(n)不是周期序列，先对x(n)进行周期性扩展，把它看成某个周期性信号的一个周期，得到的X(k)是x(n) 频谱在一个周期的采样。
DFT有一种快速算法FFT，称为快速傅里叶变换。MATLAB中采用的就是FFT算法。
2、滤波器设计
滤波器的设计一般可分为无限冲激响应数字滤波器（IIR）和有限冲激响应数字滤波器（FIR）。本课题设计主要用FIR滤波器来设计。
l IIR数字滤波器
IIR数字滤波器的系统函数是Z的有理函数，可表示为系统函数的设计就是要确定系数或者零、极点，以使滤波器满足给定的性能要求。这种设计方法如下：
（1）用模拟滤波器理论来设计数字滤波器。在IIR数字滤波器的设计中较多采用这种方法。
（2）用计算机辅助设计，优化技术设计。即要求所设计的数字滤波器的输入x(n)等于x(t)的抽样信号x(Tn)时，输出y(n)也恰好等于y(t)的抽样信号。
l FIR数字滤波器
（1）FIR数字滤波器的算法描述
FIR数字滤波器又称卷积滤波器，因为它在时域上的输入x(n)和输出信号y(n)的数学关系是卷积运算：y(n)= hd(n)*x(n)。其中hd(n)是FIR滤波器的单位冲激响应函数，其表达式是由理想滤波器频响的幅度函数H（w）进行傅里叶逆变换得到的，式中wc为滤波器的截止频率。
（2）FIR滤波器的窗函数设计法
窗函数法的基本思想是，由于滤波器的时间响应函数hd(n)是无限长的，不能用FIR滤波的方法实现，所以要截取hd(n)幅度较大的部分，舍弃hd(n)幅度较小的部分来近似表示hd(n)。截取方法是用某一种窗函数和截取hd(n)的一段进行卷积以得到实际滤波器的相应hd(n)。窗函数不同，窗口宽度不同，实际频响会有较大区别。如加海明窗的频响曲线的通带和阻带特性比矩形窗好；窗口宽度N的阶数高，也能提高滤波器频响性能，但实时信号处理的时间会相应增加。因此，应根据滤波参数要求确定选用窗函数的类型和阶数。
（3）逼近
确定了技术指标后，就可以建立一个目标数字滤波器模型。通常采用理想的数字滤波器模型。之后，利用数字滤波器的设计方法，设计一个实际滤波器模型来逼近给定的目标。
（4）性能分析和计算机仿真
上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器。根据这个描述就可以分析其频率特性和相位特性，以验证设计结果是否满足指标要求，或利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。
（5）线性相位FIR滤波器具体的窗函数法设计方法
线性相位FIR滤波器通常采用窗函数法设计。窗函数法设计FIR滤波器的基本思想是：根据给定的滤波器技术指标，选择滤波器长度N和窗函数w(n)，使其具有最窗宽度的主瓣和最小的旁瓣。其核心是从给定的频率特性，通过加窗确定有限长单位脉冲响应序列h(n)。工程中常用的窗函数共有6种，即矩形窗、巴特勒特(bartlett)窗、汉宁(hanning)窗、海明(hamming)窗、布莱克曼(blackman)窗和凯泽(kaiser)窗。
三、设计过程（程序实现）
l 语音信号的录制
l 语音信号的录制
Fs=8000; %抽样频率
time=3; %录音时间
fprintf('按Enter键录音%ds',time); %文字提示
pause; %暂停命令
fprintf('录音中......');
x=wavrecord(time*Fs,Fs,'double'); %录制语音信号
fprintf('录音结束'); %文字提示
fprintf('按Enter键回放录音');
pause; %暂停命令
wavplay(x,Fs); %播放语音信号
wavwrite(x,Fs,'f:\课程设计\sound.wav'); %存储语音信号
在录音过程中，出现了这样一个现象：录制刚开始时，常会出现实际发出声音落后录制动作半拍的现象。出现这种现象的原因主要是录音开始时，人的反应慢了半拍，导致出现了一些无效点。要去掉这些无效点，可用函数wavread()来滤除。调用格式如下：
y= wavread('f:\课程设计\sound.wav',[3423 20050]);
subplot(2,1,1);
plot(y);
title('滤除无效点后的语音信号')
l 取8000点进行频谱分析
y=wavread('f:\课程设计\sound.wav'); %读取原始语音信号
u=fft(y,8000); %对语音信号做8000点的FFT变换
subplot(2, 2, 1); %分割窗口
plot(abs(u)); %绘出语音信号的幅度谱
axis([0, 4000 ,0,80]); %显示X轴、Y轴取值范围
title('8000点幅度谱'); %加图形说明
grid; %添加网格
subplot(2, 2, 2); %分割窗口
plot(angle(u)); %绘出语音信号的相位谱
axis([0, 4000 ,-4,4]); %显示X轴、Y轴取值范围
title('8000点相位谱'); %加图形说明
grid; %添加网格
l 取16000点进行频谱分析
y=wavread('f:\课程设计\sound.wav'); %读取原始语音信号
u=fft(y,16000); %对语音信号做16000点的FFT变换
subplot(2, 2, 3); %分割窗口
plot(abs(u)); %绘出语音信号的幅度谱
axis([0, 8000 ,0,150]); %显示X轴、Y轴取值范围
title('16000点幅度谱'); %加图形说明
grid; %添加网格
subplot(2, 2, 4); %分割窗口
plot(angle(u)); %绘出语音信号的相位谱
axis([0, 8000 ,-4,4]); %显示X轴、Y轴取值范围
title('16000点相位谱'); %加图形说明
grid %添加网格
8000点和16000点数的频谱分析图
取8000点和16000点进行频谱分析，可以发现两幅度谱和相位谱图形状大致相同，但在取数点的始段，16000点比8000点的振幅大，因为取16000点比取8000点数多，16000点显示得更详细，更接近真实值，而8000点相对于16000点来说只是大致趋向真实值而已。
l 滤波器的设计
function hd = ideallp(Wc,N); %自定义函数
% hd = ideallp(Wc,N)
% ----------------------------------------------
% hd = 点0 到 N-1之间的理想脉冲响应
% wc = 截止频率(弧度)
% N = 理想滤波器的长度
%
tao = (N-1)/2; % 理想脉冲响应的对称中心位置
n = [0: (N-1)]; % 设定脉冲响应长度
m = n - tao + eps; % 加一个小数以避免零作除数
hd = sin(Wc*m) ./ (pi*m); % 理想脉冲响应
function [db,mag,pha,W]=freqz_m(b,a); %自定义函数
%db为相对振幅(dB)；mag为绝对振幅；pha为相位响应；W频率样本点向量。
[H,W]=freqz(b,a,1000, 'whole');
%b和a分别为H(z)分子和分母多项式系数;
H=(H(1:501))';
W=(W(1:501))';
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);
%。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
Wp=2* 3000/7000; Ap=1; %滤波器性能指标
Wr=2* 2400/7000; Ar=50;
Width= Wp-Wr; %过渡带宽度
Wc= (Wr+Wp)/2; %理想滤波器的截止频率
N1=ceil(6.6*pi/Width); %滤波器长度
N=N1+mod(N1+1,2); %确保N为奇数
n=0:1:N-1; %步长
hd=ideallp(Wc,N); %理想单位脉冲响应
window=hamming(N)'; %海明窗
h=window.*hd; %截取得到实际单位脉冲响应
[db,mag,pha,W]=freqz_m(h,1); %计算实际滤波器的幅频响应
subplot(2,1,1); %分割窗口
stem(n,h); %绘出实际单位脉冲响应
title('实际脉冲响应h(n)'); %加标题说明
grid; %添加网格
subplot(2,1,2); %分割窗口
plot(W/pi,db); %绘出实际滤波器的幅频响应
title('幅度响应(db)'); %加标题说明
grid %添加网格
%滤波后信号的频谱分析
fs=8000; %抽样频率
y=wavread('f:\课程设计\sound.wav'); %读取原始语音信号
y1=fftfilt(h,y); %将语音信号送入滤波器滤波，y1为输出
y2=fft(y1); %将滤波后的信号做FFT变换
subplot(2,1,1); %分割窗口
plot(abs(y2)); %绘出幅度图
axis([0, 12000 ,0,200]); %显示X轴、Y轴取值范围
title('滤波后信号的幅度谱图'); %加标题说明
grid; %添加网格
subplot(2,1,2); %分割窗口
plot(angle(y2)); %绘出相位图
axis([0, 12000 ,-4,4]); %显示X轴、Y轴取值范围
title('滤波后信号的相位谱图'); %加标题说明
grid; %添加网格
wavwrite(y1,fs,'f:\课程设计\sound1.wav'); %储存滤波后的语音信号
实际滤波器的幅度响应
滤波器的长度为121，最小阻带衰减53dB，满足滤波器设计要求。
滤波后信号的幅度、相位谱图
l 原始信号与滤波后的信号相比较
ys=wavread('f:\课程设计\sound.wav'); %读取原始语音信号
ys1= wavread('f:\课程设计\sound1.wav'); %读取滤波后语音信号
fprintf('按Enter键播放原始语音信号'); %文字提示
pause; %暂停命令
wavplay(ys,fs); %播放原始语音信号
fprintf('按Enter键播放滤波后语音信号'); %文字提示
pause; %暂停命令
wavplay(ys1,fs); %播放滤波后语音信号
yf=fft(ys); %对原始声音进行FFT变换
subplot(2,1,1); %分割窗口
plot(abs(yf)); %绘出原始信号的振幅图
axis([0, 12000 ,0,200]); %显示X轴、Y轴取值范围
title('滤波前的振幅图'); %加图形标题说明
yf1=fft(ys1); %对滤波后的声音进行FFT变换
subplot(2,1,2); %分割窗口
plot(abs(yf1)); %绘出滤波后信号的振幅图
axis([0, 12000 ,0,200]); %显示X轴、Y轴取值范围
title('滤波后的振幅图'); %加图形标题说明
通过滤波前后两信号的比较，通过上面的振幅图可以很清晰地观察到经过滤波后与滤波前存在明显差别，滤波后的图比较平稳，原始语音信号比滤波后的信号清晰，滤波后的声音显得比较低沉，因为原始信号经过低通滤波器后，低通滤波器已经把信号高频部分滤掉了，只剩下低频部分。
结论
经过一个星期的设计，本人基本按质按量完成了课程设计的各项工作。通过本次课程设计，使我加深了对数字信号处理这门课和MATLAB这个软件工具的认识。语音信号处理是以MATLAB为工具平台、数字信号处理为基础而形成的一门综合性学科，处理的目的是要得到一些语音参数以便高效的传输或存储，或者通过处理的某种运算以达到某种用途的要求。而MATLAB是一款多功能多特点的的软件，对于现代科学工程的计算和分析具有重要的意义。语音信号处理与滤波的设计对于所学的数字信号处理课程的理论得到了综合的应用，并从实践上初步实现对数字信号的处理。
在设计过程中遇到了很多从未遇到过的问题，小到一些标点符号，大到程序的设计，但我还是有耐心，经过多方面查资料和请教老师、同学，逐渐把一些不懂的问题解决，最终完成本次课程设计。在本次设计过程中，使我对数字信号处理的设计和MATLAB软件的认识更进一层，更重要的是使我学会了如何面对困难。在设计中，总结了以下几点经验，这也是自己在实践中一步一步摸索出来。
1、开始设计时，不要太过了想象完美，有思路就马上记录下来，在查阅资料也要注意记录重点。
2、遇到问题，不要着急，保住沉着的心态，通过查阅资料、请教老师或同学找出解决途径。
3、做设计，单凭自己的力量是不够的，要学会学习别人好的方面，多查阅资料，多与老师、同学交流，把别人的知识转变为自己的知识。
4、养成良好的工作态度，对以后工作具有具有重要意义。
在完成设计的过程中，由于所学知识的狭窄，本次设计不免有很多不理想的地方，这些都是有待提高的。
参考文献
(1)程佩青.数字信号处理教程(第三版).清华大学出版社.20##年2月
(2)赵红怡,张常年.数字信号处理及其MATLAB实现.化学工业出版社.20##年1月
(3)钱同惠.数字信号处理.机械工业出版社.20##年2月
(4)李海涛,邓樱.MATLAB6.1基础及应用技巧国防工业出版社20##年3月