matlab在信号与系统中的应用,Matlab在信号与系统中的应用

广西轻工业

ＧＵＡＮＧＸＩＪＯＵＲＮＡＬＯＦＬＩＧＨＴＩＮＤＵＳＴＲＹ

计算机与信息技术

２００８年３月第３期(总第１１２期)

【作者简介】官宗琪，３０岁，女，山东平度人，讲师，研究方向：移动通信。

１引言

信号与系统是通信电子类专业的一门重要的专业基础课，

其特点是理论性较强，概念难于理解，计算繁琐，系统分析时出现的时域图和频谱图较难绘制。Ｍａｔｌａｂ软件的出现给信号与系统的计算机辅助教学带来了便捷，它不仅能把我们从以往复杂的计算中解脱出来，更重要的是，通过几句简单的程序，能准确的绘制出各种图形，这有助于同学们对理论知识的理解，使这门理论性很强的课程更易于接受。下面简述一些Ｍａｔｌａｂ在信号与系统中的应用。

２用Ｍａｔｌａｂ实现周期信号的傅利叶级数分解与综合

系统的时域分析中，常以正弦函数为基本信号，分析工程

上常用的周期信号的一些基本特性以及信号在系统中的传输问题。周期信号到正弦信号的分解是利用傅里叶级数，根据这种正弦分解思想，系统的响应可看作各不同频率正弦信号产生的响应的叠加。信号的分解是很繁琐的数学计算，运用Ｍａｔｌａｂ的符号运算功能，能很容易的实现信号的分解。

三角形式的傅里叶级数：ｆ

(ｔ)＝ａ

０２

＋

∞

ｎ＝１

!

ａｎｃｏｓ(ｎｗ０ｔ)＋ｂｎｓｉｎ(ｎｗ０ｔ"

#)其中，傅里叶系数为：ａｎ＝２

Ｔ

Ｔ

０$

ｆ(ｔ)ｃｏｓ(ｎｗ０ｔ)ｄｔ

ｂｎ＝２

Ｔ

Ｔ０

$ｆ(ｔ)ｓｉｎ(ｎｗ０

ｔ)ｄｔ

公式表明，任何满足狄利赫利条件的周期信号都可以分解成不同频率正弦分量和余弦分量的叠加。ｎ＝０时，为直流分量，

ｎ＝１时为一次谐波分量，以此类推。可见，周期信号可分解成各

次谐波分量的叠加，反过来，各次谐波分量的叠加可以合成周期信号。

下面，我们用Ｍａｔｌａｂ实现上面的信号分解与合成。给出信号为周期矩形脉冲信号，幅度为１，脉冲宽度为，周期为，用

Ｍａｔｌａｂ语言编程求出该信号三角形式的傅里叶系数，并绘出各

次谐波叠加的傅里叶综合波形图。

实现流程为：(１)编写子函数ｘ＝ｔｉｍｅ＿ｆｕｎ＿ｘ(ｔ)，用符号表达式表示出周期信号在第一个周期内的符号表达式。(２)编写子函数ｙ＝ｔｉｍｅ＿ｆｕｎ＿ｅ(ｔ)，求出该信号在绘图区间内的信号样值。(３)编写求解信号傅里叶系数及绘制合成波形图的普通函数ＣＴＦＳｈｃｈｓｙｍ，该函数流程为：①调用函数ｘ＝ｔｉｍｅ＿ｆｕｎ＿ｘ(ｔ)，获取周期信号的符号表达式；②求信号的傅里叶系数；③求各次谐波；④绘制各次谐波叠加波形图；⑤调用信号ｔｉｍｅ＿ｆｕｎ＿ｅ

(ｔ)，绘制原信号波形图。

程序运行结果如图１所示：

图１周期信号的分解与合成

３用Ｍａｔｌａｂ实现连续时间系统Ｓ域零极点分析

根据系统的系统函数，利用Ｍａｔｌａｂ绘制出系统的零极点

分布图，以及系统的冲激响应和频率响应，加深对系统频域分析的理解。假设系统的系统函数为Ｈ(ｓ)＝ｓ２

＋ｓ＋２３ｓ３＋５ｓ２

＋４ｓ－６

。Ｍａｔｌａｂ中利用系统函数的分子分母多项式系数就可以描

述出其对应的连续系统。因此设参数ｎｕｍ＝［１３２］，ｄｅｎ＝［３５４

－６］，ｓｙｓ＝ｔｆ(ｎｕｍ，ｄｅｎ)。绘制零极点图的函数为ｐｚｍａｐ(ｓｙｓ)，绘制

冲激响应的函数为ｉｍｐｕｌｓｅ(ｎｕｍ，ｄｅｎ，ｔ)，绘制频率响应的函数为ｆｒｅｑｓ(ｎｕｍ，ｄｅｎ)。根据以上重要的函数，在选择合适的坐标范围，可以很方便的绘制出该连续时间系统的零极点分布图、冲激响应和频率响应。在学生绘制零极点分布图的过程中，有的同学根据例题的坐标范围，绘制出来系统的极点只有一个，显然左边范围取的偏小，经过更正后，得到了完整的零极点分布图。此外还有系统的冲激响应和频率响应如图２所示：

图２系统的冲激响应与频率响应

Ｍａｔｌａｂ在信号与系统中的应用

官宗琪

(浙江万里学院电子信息学院，浙江宁波３１５１００)

【摘

要】针对电子信息类专业以及信号与系统课程的特点，提出了将ＭＡＴＬＡＢ引入信号与系统课程的教与学中，通过

简单的程序及运行结果，形象的描述了课程中难于理解的概念，这样既能加强学生对理论知识的掌握及提高解决实际问题的能力，又能为课堂教学及教学方法和手段的改革增添活力。

【关键词】傅里叶级数；分解；合成；谐波【中图分类号】ＴＰ３９１．６

【文献标识码】Ａ

【文章编号】１００３－２６７３(２００８)０３－７７－０１

(下转第３１页)