北京理工大学信号与系统实验报告5-连续时间系统的复频域分析
实验5连续时间系统的复频域分析
（综合型实验）
一、实验目的
1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。
2）学习和掌握连续时间系统函数的定义及复频域分析方法。
3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。
二、实验原理与方法
1.拉普拉斯变换
连续时间信号x(t)的拉普拉斯变换定义为(s)(t)e st X x dt +∞
--∞=
⎰ （1） 拉普拉斯反变换为1(t)(s)e 2j st j x X ds j σσπ+∞-∞=⎰ （2）
MATLAB 中相应函数如下：
(F)L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。
(F,t)L laplace =用t 替换结果中的变量s 。
()F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。
(,)F ilaplace L x =用x 替换结果中的变量t 。
拉氏变换还可采用部分分式法，当(s)X 为有理分式时，它可以表示为两个多项式之比：
110110
...(s)(s)(s)...M M M M N N N N b s b s b N X D a s a s a ----+++==+++ （3） 上式可以采用部分分式法展成以下形式
1212(s)...N N
r r r X s p s p s p =+++--- （4） 再通过查找常用拉氏变换对易得反变换。
利用residue 函数可将X(s)展成（4）式形式，调用格式为：
[r,p,k]residue(b,a)=其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。
2.连续时间系统的系统函数
连续时间系统的系统函数是指系统单位冲激响应的拉氏变换
(s)(t)e st H h dt +∞
--∞=
⎰
（5）
连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。
(s)(s)/X(s)H Y = （6）
单位冲激响应(t)h 反映了系统的固有性质，而(s)H 从复频域反映了系统的固有性质。
由（6）描述的连续时间系统，其系统函数为s 的有理函数
110110
...(s)...M M M M N N N N b s b s b H a s a s a ----+++=+++ （7） 3.连续时间系统的零极点分析
系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。
通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用⨯表示，这样得到的图形为零极点分布图。
可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下：
r=roots(c)，c 为多项式的系数向量，返回值r 为多项式的根向量。
求取零极点以及绘制系统函数的零极点分布图可以采用pzmap 函数，调用格式如下： pzmap(sys)绘出由系统模型sys 描述的系统的零极点分布图。
[p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。
还有两个专用函数tf2zp 和zp2tf 可实现系统的传递函数模型和零极点增益模型的转换。
调用格式如下：
[z,p,k]=tf2zp(b,a)
[b,a]=tf2zp(z,p,k)
研究系统函数的零极点分布不仅可以了解系统冲激响应的形式，还可以了解系统的频率特性以及判断系统的稳定性。
1） 零极点分布与冲激响应的关系
系统的极点位置决定着系统冲激响应h(t)的波形，冲激响应的幅值是由系统函数的零点和极点共同确定的，系统的零点位置只影响冲激响应的幅度和相位，不影响波形。
2） 零极点分布与系统频率响应的关系
系统函数的零极点分布不仅决定了系统函数H(s)，也决定了系统的频率响应()H ω，根据系统的零极点分布情况，可以由几何矢量法分析系统的频率响应。
3） 零极点分布与系统稳定性的关系
稳定性是系统的固有性质，与激励信号无关，由于系统函数(s)H 包含了系统的所固有的性质，因而可以根据系统函数的零极点分布判断系统的稳定性。
因果系统稳定的充要条件是(s)H 的全部极点位于s 的左半平面。
三．实验内容
（1）已知系统的冲激响应(t)u(t)u(t 2)h =--，输入信号(t)u(t)x =，试采用复频域的方法求解系统的响应，编写MATLAB 程序实现。
代码：
%DFTfifth_2_1.m
syms t。
北京理工大学信号与系统实验实验报告
北京理工大学信号与系统实验实验报告信号与系统实验报告姓名:肖枫学号:1120111431班号:05611102专业:信息对抗技术学院:信息与电子学院12实验1 信号的时域描述与运算一、实验目的1. 掌握信号的MATLAB表示及其可视化方法。
2. 掌握信号基本时域运算的MATLAB实现方法。
3. 利用MATLAB分析常用信号，加深对信号时域特性的理解。
二、实验原理与方法1. 连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。
在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。
表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。
例如一个正弦信号可以表示如下:>> t=0:0.01:10;>> x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。
如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。
例如对于上述正弦信号，可以用符号对象表示如下:>> x=sin(t);>> ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形10.80.60.40.2-0.2-0.4-0.6-0.8-1012345678910Time(seconds)图1 利用向量表示连续时间信号3sin(t)10.5-0.5-1-6-4-20246t图 2 利用符号对象表示连续时间信号常用的信号产生函数函数名功能函数名功能 heaviside 单位阶跃函数 rectpuls 门函数 sin 正弦函数 tripuls 三角脉冲函数 cos 余弦函数 square 周期方波 sinc sinc函数 sawtooth 周期锯齿波或三角波 exp 指数函数2.连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换(尺度伸缩)等。
北京理工大学信号与系统实验讲义电子版
实验1 信号的时域描述与运算 一、实验目的①掌握信号的MATLAB 表示及其可视化方法。
②掌握信号基本时域运算的MA TLAB 实现方法。
③利用MA TLAB 分析常用信号，加深对信号时域特性的理解。
二、实验原理与方法1. 连续时间信号的MATLAB 表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。
在MATLAB 中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
从严格意义上来说，MATLAB 并不能处理连续时间信号，在MATLAB 中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。
表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。
例如一个正弦信号可以表示如下： >> t=0:0.01:10; >> x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。
如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。
例如对于上述正弦信号，可以用符号对象表示如下： >> x=sin(t); >> ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形012345678910-1-0.8-0.6-0.4-0.200.20.40.60.81Time(seconds)图1 利用向量表示连续时间信号常用的信号产生函数 函数名 功能 函数名 功能 heaviside 单位阶跃函数 rectpuls 门函数 sin 正弦函数 tripuls 三角脉冲函数 cos余弦函数square周期方波sinc sinc 函数 sawtooth 周期锯齿波或三角波 exp指数函数-6-4-20246-1-0.50.51t图 2 利用符号对象表示连续时间信号sin(t)2.连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。
北京理工大学信号与系统信号的频域分析
x(t ) a0 ak cos kw0t bk sin kw0t
k 1 k 1
（3）
其中：
a0
1 2 2 T0 x(t )dt , ak T x(t ) cos kw0tdt , bk T x(t ) sin kw0tdt 0 T0 T0 T0 0
（4）
4、离散非周期时间信号的频域分析
非周期序列 x( n) 可以表示成一组复指数序列的连续和
x ( n)
其中
1 2
X (e
2
j
)e jn d
（15）
X (e j )
n
x ( n )e
jn
（16）
式（16）称为 x( n) 的离散时间傅里叶变换，式（15）和式（16）确立了非周期离散时 间信号 x( n) 及其离散时间傅里叶变换 X (e 称为频谱函数，且 X (e
1 .2
1
0 .8
0 .6
0 .4
0 .2
0
-0 .2 -1 .5
-1
-0 .5
0
0 .5
1
1 .5
N=20;T=1;a=0.5;A=1;
1 .2
1
0 .8
0 .6
0 .4
0 .2
0
-0 .2 -1 .5
-1
-0 .5
0
0 .5
1
1 .5
③利用 MATLAB 绘出周期矩形脉冲信号的频谱，观察参数 T 和变 化时对频谱波形的影响。
Answer：频谱包络形状不变，过零点不变，普贤间隔随着 T 变大而缩小。
2、已知 x(t)是如图所示的矩形脉冲信号。 ①求该信号的傅里叶变换； ②利用 MATLAB 绘出矩形脉冲信号的频谱，观察矩形脉冲信号宽 度变化时对频谱波形的影响； ③让矩形脉冲信号的面积始终等于 1，改变矩形脉冲宽度，观察 矩形脉冲信号时域波形和频谱随矩形脉冲宽度的变化趋势。
实验5 连续时间系统的频域和复频域分析
实验5 连续时间系统的频域和复频域分析一．实验目的1.掌握和理解连续时间函数系统频率相应、系统函数的概念和物理意义。
2.学习和掌握连续时间系统频域、复频域的分析方法。
3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。
二．实验原理1.连续时间系统的频率响应系统的频率响应定义为：ττωωτd eh j H j -∞∞-⎰=)()(H （ωj ）反映了LTI 连续时间系统对不同频率信号的相应特性，是系统内在固有的特性，与外部激励无关。
H （ωj ）又可以表示为)()()(ωθωωj ej H j H =其中)(ωj H 称为系统的幅度响应，)(ωθ成为系统的相应响应。
对于由下述微分方程描述的LTI 连续时间系统∑∑===Mm m n Nn n n t xb t ya 0)(0)()()(其频率响应H （ωj ）可以表示为下列式子所示的ωj 的有理多项式1110111...)()(...)()()()()(a j a j a j a b j b j b j b X Y j H N N N N M M M M ++++++++==----ωωωωωωωωωMATLAB 的信号处理工具箱提供了专门的函数freqs ，用来分析连续时间系统的频率响应，该函数有下列几种调用格式：[h,w]=freqs(b,a) 计算默认频率范围内200个频率点上的频率响应的取样值，这200个频率点记录在w 中。
h=freqs （b ，a ，w ） b 、a 分别为表示H （ωj ）的有理多项式中分子和分母多项式的系数向量，w 为频率取样点，返回值h 就是频率响应在频率取样点上的数值向量。
[h ，w]=freqs （b ，a ，n) 计算默认频率范围内n 个频率点上的频率响应的取样值，这n 个频率点记录在w 中。
Freqs （b ，a ，……） 这种调用格式不返回频率响应的取样值，而是以对数坐标的方式绘出来系统的频率响应和相频响应。
连续时间信号与系统的频域分析报告
连续时间信号与系统的频域分析报告1. 引言连续时间信号与系统的频域分析是信号与系统理论中的重要分支，通过将信号和系统转换到频域，可以更好地理解和分析信号的频谱特性。
本报告将对连续时间信号与系统的频域分析进行详细介绍，并通过实例进行说明。
2. 连续时间信号的频域表示连续时间信号可以通过傅里叶变换将其转换到频域。
傅里叶变换将信号分解成一系列不同频率的正弦和余弦波的和。
具体来说，对于连续时间信号x(t)，其傅里叶变换表示为X(ω)，其中ω表示频率。
3. 连续时间系统的频域表示连续时间系统可以通过频域中的频率响应来描述。
频率响应是系统对不同频率输入信号的响应情况。
通过系统函数H(ω)可以计算系统的频率响应。
系统函数是频域中系统输出与输入之比的函数，也可以通过傅里叶变换来表示。
4. 连续时间信号的频域分析频域分析可以帮助我们更好地理解信号的频谱特性。
通过频域分析，我们可以获取信号的频率成分、频谱特性以及信号与系统之间的关系。
常用的频域分析方法包括功率谱密度估计、谱线估计等。
5. 连续时间系统的频域分析频域分析也可以用于系统的性能评估和系统设计。
通过分析系统的频响特性，我们可以了解系统在不同频率下的增益和相位变化情况，进而可以对系统进行优化和设计。
6. 实例分析以音频信号的频域分析为例，我们可以通过对音频信号进行傅里叶变换，将其转换到频域。
通过频域分析，我们可以获取音频信号的频谱图，从而了解音频信号的频率成分和频率能量分布情况。
进一步，我们可以对音频信号进行系统设计和处理，比如对音乐进行均衡、滤波等操作。
7. 结论连续时间信号与系统的频域分析是信号与系统理论中重要的内容，通过对信号和系统进行频域分析，可以更好地理解和分析信号的频谱特性。
频域分析也可以用于系统的性能评估和系统设计，对于音频信号的处理和优化具有重要意义。
总结：通过本报告，我们了解了连续时间信号与系统的频域分析的基本原理和方法。
频域分析可以帮助我们更好地理解信号的频谱特性和系统的频响特性，对系统设计和信号处理具有重要意义。
连续时间信号与系统的频域分析实验报告
《信号与系统》课程实验报告
一．实验原理 1、傅里叶变换 实验原理如下：
傅里叶变换的调用格式
F=fourier(f)：返回关于w 的函数；
F=fourier(f ，v):返回关于符号对象v 的函数，而不是w 的函数。
傅里叶逆变换的调用格式
f=ifourier(F):它是符号函数F 的fourier 逆变换，返回关于x 的函数； f=ifourier(f,u):返回关于u 的函数。
2、连续时间信号的频谱图 实验原理如下：
符号算法求解如下：
ft=sym('4*cos(2*pi*6*t)*(heaviside(t+1/4)-heaviside(t-1/4))'); Fw=simplify(fourier(ft)) subplot(121)
ezplot(ft,[-0.5 0.5]),grid on subplot(122)
ezplot(abs(Fw),[-24*pi 24*pi]),grid on 波形图如下所示：
当信号不能用解析式表达时，无法用MATLAB 符号算法求傅里叶变换，则用MATLAB 的数值计算连续信号的傅里叶变换。
∑⎰
∞
-∞
=-→-∞∞
-==n n j t
j e
n f dt e
t f j F ττωτ
ωτω)(lim
)()(0
若信号是时限的，或当时间大于某个给定值时，信号已衰减的很厉害，可以近似地看成时限信号，设n 的取值为N ，有
1
1()
a jw
++
的分母和分子多项式的系数向量，
1、在调用函数fourier()及ifourier()之前，要用syms命令对所用到的变。
连续时间信号与系统的复频域分析
《信号与系统》课程实验报告一、实验原理的验证 1、 用MATLAB 进行部分分式展开实验原理如下：Residue 函数可以得到复杂有理分式F(s)的部分分式展开式，其调用格式为。
num,den 分别为F(s)的分子和分母多项式的系数向量，r 为部分分式的系数，p 为极点，k 为F(s)中整式部分的系数，若F(s)为有理真分式，则k 为零。
例6-1 用部分分式展开法求F(s)的反变换实验结果如下：理论值分析如下：2、 用MATLAB 分析LTI 系统的特性实验原理如下：系统函数H （s ）通常是一个有理分式，其分子和分母均为多项式。
计算H （s ）的零极点可以用roots 函数，求出分子和分母多项式的根，然后用plot 命令画图。
也可以直接用更简便的pzmap 函数画系统函数H （s ）的零极点分布图，其调用格式为pzmap(sys)，借助tf 函数获得LTI 系统的模型sys ，其调用格式为sys=tf(b,a)，b 和a 分别为系统函数H （s ）的分子和分母多项式的系数向量。
例6-2零极点分布图：[],,(,)r p k residue num den =322()43s F s s s s +=++单位冲激响应频率响应由零极点图分析可知：叉号代表极点，圆圈代表零点。
因为极点均在左半平面，故该系统是稳定系统，且H(s)的收敛域包含虚轴，频率响应函数可直接令s=jw，并且由于极点均在左半平面，故频率响应函数是个衰减的函数。
3、 用MATLAB 进行Laplace 正、反变换实验原理如下：我们采用函数Laplace 和ilaplace 分别计算Laplace 正、反变换,其调用格式分别为 ，上述两式右端的f 和F 分别为时域表示式和s 域表示式的符号表示，可以应用函数sym 实现，其调用格式为S=sym(A)式中，A 为待分析表示式的字符串，S 为符号数字或变量。
例6-3 （1）的拉普拉斯变换；实验结果如下:理论值分析如下：（2）的拉普拉斯反变换。
连续时间信号与系统的频域分析实验报告
《信号与系统》课程实验报告一•实验原理 1傅里叶变换实验原理如下：傅里叶变换的调用格式F=fourier(f):返回关于 W 的函数；F=fourier(f , v):返回关于符号对象V 的函数，而不是W 的函数。
傅里叶逆变换的调用格式f=ifourier(F):它是符号函数F 的fourier 逆变换，返回关于X 的函数; f=ifourier(f,u):返回关于U 的函数。
2、连续时间信号的频谱图实验原理如下： 符号算法求解如下：ft=sym('4*cos(2*pi*6*t)*(heaviside(t+1∕4)-heaviside(t-1∕4))'); FW=SimPlify(fourier(ft))subplot(121)ezplot(ft,[-0.5 0.5]),grid Onsubplot(122) ezplot(abs(Fw),[-24*pi 24*pi]),grid On波形图如下所示：当信号不能用解析式表达时，无法用换，则用MATLAB 的数值计算连续信号的傅里叶变换。
实验步骤或实验方案MATLAB 符号算法求傅里叶变F(j )f(t)ejt dt 叫nf (n )e若信号是时限的，或当时间大于某个给定值时，信号已衰减的很厉 害，可以近似地看成时限信号，设 n 的取值为N ,有4 CO$(12 I )■) (he 如引日环-IMh heaviside(t IeIXW Sin(WM ⅛)yabS(W i -144 >2)3、 用MATLAB 分析LTl 系统的频率特性当系统的频率响应H (jw )是jw 的有理多项式时，有H(S )B(W) b M (jW)Mb Mi (jW)MIL b ι(jw) b oH (jW)NN 1A(W)a N (jw)a ” ι(jw) L α(jw) a °freqs 函数可直接计算系统的频率响应的数值解，其调用格式为H=freqs(b,a,w)其中，a 和b 分别是H(jw)的分母和分子多项式的系数向量，W 定义 了系统频率响应的频率范围，P 为频率取样间隔。
连续时间信号与系统的频域分析实验报告
实验四连续时间信号与系统的频域分析一、实验目的掌握连续时间信号的傅里叶变换及傅里叶逆变换的实现方法，掌握连续时间系统的频域分析方法，熟悉MATLAB 相应函数的调用格式和作用，掌握使用MATLAB 来分析连续时间信号与系统的频域特性及绘制信号频谱图的方法。
二、实验原理（一）连续时间信号与系统的频域分析原理1、连续时间信号的额频域分析 连续时间信号的傅里叶变换为：()()dt e t f j F t j ωω-∞∞-⎰=傅里叶逆变换为：()()ωωπωd e j F t f t j ⎰∞∞-=21()ωj F 称为频谱密度函数，简称频谱。
一般是复函数，可记为：()()()ωϕωωj e j F j F =()ωj F 反映信号各频率分量的幅度随频率ω的变化情况，称为信号幅度频谱。
()ωϕ反映信号各频率分量的相位随频率ω的变化情况，称为信号相位频谱。
2、连续时间系统的频域分析 在n 阶系统情况下，数学模型为：()()()()()()()()t f b dtt df b dt t f d b dt t f d b t y a dtt dy a dt t y d a dt t y d a o m m n m m n o n n n n n n ++++=++++------11111111 令初始条件为零，两端取傅里叶变换，得：()()[]()()()[]()ωωωωωωωωj F b j b j b j b j Y a j a j a j a m n m n n n nn01110111++++=++++----表示为()()()()ωωωωj F j b j Y j a kmk kkn k k∑∑===0则 ()()()()()()()()()∑∑==----=++++++++==nk kk mk kk n n n n m m mm j a j b a j a j a j a b j b j b j b j F j Y j H 0001110111ωωωωωωωωωωω3、系统传递函数 系统传递函数定义为:()()()ωωωj H j Y j H =系统传递函数反映了系统内在的固有的特性，它取决于系统自身的结构及参数，与外部 激励无关，是描述系统特性的一个重要参数。
实验5--连续时间系统的复频域分析
实验5–连续时间系统的复频域分析实验背景在连续时间系统的频域分析中，复频域分析是非常重要的一个方法。
其可以帮助我们更直观地了解系统的频率响应，包括幅频响应和相频响应，对于系统的设计和优化都有非常实际的应用价值。
因此，在本次实验中，我们将通过对一个特定系统的复频域分析来学习这一方法的基本原理和操作流程。
实验目的1.了解连续时间系统的幅频响应和相频响应2.掌握利用MATLAB对系统进行复频域分析的方法3.学会根据复频域图像对系统进行分析和优化实验原理连续时间系统幅频响应和相频响应在连续时间系统的频域分析中，使用的是拉普拉斯变换。
通过对系统的输入信号和输出信号进行拉普拉斯变换，可以得到它们在复平面上的函数，进而求得系统的传递函数H(s):H(s)=Y(s)/X(s)其中，s为复变量。
系统的幅频响应和相频响应分别定义为:H(s)的模和相位:|H(jw)|=sqrt(H(s)H(s)*) (模) arg(H(jw))=tan^-1[Im{H(jw)}]/Re{H(jw)} (相位) 其中，w为实数，j为虚数单位。
利用MATLAB进行系统复频域分析MATLAB提供了众多用于连续时间系统复频域分析的工具。
其中，最基本的是bode命令。
它可以计算和绘制给定系统的幅频响应和相频响应曲线。
常用命令格式如下：[bode(H,w)]其中，H为系统的传递函数，w为频率范围除此之外，MATLAB还提供了很多其他的命令，如nyquist、margin、freqresp 等。
它们可以帮助我们更全面地分析系统的性能和特点。
实验步骤实验环境1.一台已安装MATLAB的计算机实验流程1.根据给定的系统传递函数H(s)，利用MATLAB计算和绘制其幅频响应和相频响应曲线。
%定义系统传递函数H=tf([5+j*10 0.6+0.2*j],[1 2+j 3 4-j 5+j]);%绘制幅频响应和相频响应曲线figure(1)subplot(2,1,1)bode(H);subplot(2,1,2)nyquist(H);2.根据绘制的幅频响应和相频响应曲线，对系统进行分析和优化。