数字信号处理课程设计以及MATLAB程序

离散信号与系统的时域分析

摘要：通过对数字信号处理的深刻理解，在有MATLAB对信号进行仿真前提下，使得离散信号输出更准确更形象。利用MATLAB工具箱中filter函数对差分方程求解，防真出系统响应和系统单位脉冲响应，用conv函数完成线形卷积。在判断系统稳定性时，巧妙的将输入信号改为单位阶跃序列，观察稳态输出是否趋于常数（包括零）来判断，如果输出趋于数，则系统稳定，否则系统不稳定。然而，完整的模拟也会存在细小的瑕疵，所以本文还需要改进。尤其是在卷积编程时，冗余度较大。

报告要求：

（1）简述离散系统时域分析方法；

答：离散系统在时域中，是通过差分方程，单位脉冲响应或系统函数对系统特性进行描写。所以在输入信号已知的情况下，知道其中之一就可以利用MATLAB就可以对其进行求解。倘若一直差分方程就适合用工具箱函数filter求差分方程的解。也可以用conv函数计算输入信号核系统的单位脉冲响应的线形卷积，求出系统响应。

（2）简述通过实验判断系统稳定性的方法；

答：判断系统稳定性，从时域看可以通过检查系统单位脉冲响应是否绝对可和来判断。但实际中，只要用单位阶跃序列作为输入信号，如果稳态输出趋于常数或零，则系统稳定，否则系统不稳定。

（3）完成设计实验，并对结果进行分析和解释；

（4）打印程序清单和要求画出的信号波形；

（5）写出本次课程设计的收获和体会。

程序清单和波形图

定义基础函数：

1  uDT为单位阶跃序列

function y=uDT(n)   

y=n>=0;            %当参数大于等于0时阶跃为1，否则为0

1. impDT为单位采样序列

function y=impDT(n)

y=(n==0);              %当参数为0时冲激为1,否则为0

第一题：

1 给定一个低通滤波器的差分方程为



输入信号分别为

时输出波形如下：

程序一：

%系统响应%

N=8;

a=[1 -0.9];

b=[0.05 0.05];

n=0:30;

subplot(211)

x1=uDT(n)-uDT(n-N);         

y1=filter(b,a,x1);

stem(n,y1,'fill','g'),grid on

xlabel('n'),title('系统响应y1(n)')

subplot(212)

x2=uDT(n);

y2=filter(b,a,x2);

stem(n,y2,'fill','r'),grid on

xlabel('n'),title('系统响应y2(n)')

2 系统的单位脉冲响应波形图如下：

程序二：

%系统单位脉冲响应%

a=[1 -0.9];

b=[0.05 0.05];

n=0:30;

impz(b,a,30),grid on

title('系统单位脉冲响应h(n)')

第二题：

给定系统的单位脉冲响应为

用线性卷积法求x1(n)=R8(n)分别对系统h1(n)和h2(n)的输出响应，画出波形如下：

程序三：

%线性卷积%

nx=-1:5;             %x(n)向量显示范围(添加了附加的零值)

nh=-2:10;            %h(n)向量显示范围(添加了附加的零值)

x=uDT(nx)-uDT(nx-8);

h1=uDT(nh)-uDT(nh-10);

h2=impDT(nh)+2.5.*impDT(nh-1)+2.5.*impDT(nh-2)+impDT(nh-3);

y1=conv(x,h1);

y2=conv(x,h2);

ny1=nx(1)+nh(1);  %卷积结果起始点

%卷积结果长度为两序列长度之和减1,即0到(length(nx)+length(nh)-2);

%因此卷积结果的时间范围是将上述长度加上起始点的偏移值

ny=ny1+(0:(length(nx)+length(nh)-2));

subplot(321)

stem(nx,x,'fill'),grid on

xlabel('n'),title('x(n)')

axis([-4 16 0 3])

subplot(322)

stem(nh,h1','fill'),grid on

xlabel('n'),title('h1(n)')

axis([-4 16 0 3])

subplot(323)

stem(nh,h2','fill'),grid on

xlabel('n'),title('h2(n)')

axis([-4 16 0 3])

subplot(324)

stem(ny,y1,'fill'),grid on

xlabel('n'),title('y(n)=x(n)*h1(n)')

axis([-4 16 0 3])

subplot(325)

stem(ny,y2,'fill'),grid on

xlabel('n'),title('y(n)=x(n)*h2(n)')

axis([-4 16 0 3])

第三题：

给定一谐振器的差分方程为



令b0=1/100.49，谐振器的谐振频率为0.4 rad。

① 用实验方法检查系统是否稳定。输入信号为u(n)时，画出系统输出波形如下：

从频域看：

程序四：

B=[1,1.8237,-0.9801];

A=[1/100.49,0,-1/100.49];

zplane(B,A),grid on

legend('零点','极点')

title('零极点分布图')

（由于系统所有极点在单位圆上，则系统稳定）

从时域看：

程序五：

%稳定性判断

a=[1 -1.8237 0.9801];

b=[1/100.49 0 -1/100.49];

n=0:300;

x1=uDT(n);

y1=filter(b,a,x1);

stem(n,y1,'fill''b'),grid on

xlabel('n'),title('系统响应y1(n)')

（由于系统输出趋于零，则系统稳定）

（2） 给定输入信号为



求出系统的输出响应，画出其波形图如下：

程序六：

%差分方程对应的系统函数

a=[1 -1.8237 0.9801];

b=[1/100.49 0 1/100.49];

n=0:30;

subplot(211)

x1=uDT(n);

y1=filter(b,a,x1);

stem(n,y1,'fill','g'),grid on

xlabel('n'),title('系统响应y1(n)')

subplot(212)

x2=sin(0.014.*n)+sin(0.4.*n);

y2=filter(b,a,x2);

stem(n,y2,'fill','r'),grid on

xlabel('n'),title('系统响应y2(n)')

 

参考文献：

[1] 吴礼斌，李伯年，数学实验与建模 北京，国防工业出版社2007

[2] 张威 MATLAB基础与编程入门（第二版）西安， 西安电子科技大学  2008

[3] 高西全，丁玉美，数字信号处理（第三版）西安，西安电子科技大学  2008

[4] Vinay K Lngle,John G Proakis.数字信号处理及其MATLAB实现 陈怀琛 王朝英,高西全, 译.北京:电子工业出版社,1988

课程设计论文撰写的内容和要求

课程设计论文要求每个人写一份，字数要求3000～5000字。

涉及到计算机软件：学生要独立完成一个软件或较大软件中的一个模块，要有足够的工作量；要写出软件说明书，能够进行计算机演示和给出运算结果。

涉及到计算机硬件或电控装置：学生要独立完成一个完整的实验，要有完整的测试结果和实验数据，实验要有探索性，要写出使用说明书。

1、标题

标题应该简短、明确、有概括性。标题字数要适当，不宜超过20个字。如果有些细节必须放进标题，可以分成主标题和副标题。

2、论文摘要

论文摘要简要陈述研究课题的内容，创新见解和主要论点，中文摘要在200字左右。摘要的最后应列出关键词。

[3、目录

目录按三级标题编写（即：[一、…，1. …，（1）…]），要求标题层次清晰。目录中的标题应于正文中的标题一致。每一章节后应标明起始页码。课程设计的全部标题层次应有条不紊，整齐清晰。]

3、正文

正文包括绪论、正文主体与结论，其内容分别如下：

• 绪论：说明本课题的意义、目的、研究范围及要达到的技术要求；简述本课题在国内外的发展概况及存在的问题；说明本课题的指导思想；阐述本课题应解决的主要问题。
• 方案：说明设计原理并进行方案选择，阐明为什么要选择这个设计方案以及所采用方案的特点。
• 过程论述：指作者对自己的研究工作的详细表述。要求理论正确、论据确凿、逻辑性强、层次分明、表达确切，
• 结果分析：对研究过程中所获得的主要的数据、现象进行定性或定量分析，得出结论或推论。
• 结论：对整个研究工作进行归纳和综合，阐述本课题研究中尚存在的问题及进一步开展研究的见解和建议。结论要写的概括、简短。

4、致谢

谢辞应以简短的文字对在课程设计过程中曾给予帮助的人（如指导教师、答疑5、5、参考文献

参考文献是课程设计论文中不可缺少的组成部分，所引用的文献必须是本人真正阅读过的，近期发表的与设计或论文工作直接相关的文献。它反映课程设计论文的取材来源、材料的广博程度和材料的可靠程度，列入的主要的文献在3篇以上。

[

1. 附录

附录是对于一些不宜放在正文中，但又直接反映完成工作的成果内容。如图纸、实验数据、计算机程序等材料附于课程设计论文之后，附录所包括的材料是课程设计论文的重要组成部分。]