典型信号的拉普拉斯变换和拉普拉斯逆变换
- 格式:doc
- 大小:217.00 KB
- 文档页数:15
419附录A 拉普拉斯变换及反变换1.表A-1 拉氏变换的基本性质 齐次性)()]([s aF t af L =1线性定理叠加性)()()]()([2121s F s F t f t f L ±=±一般形式=−=][′−ٛ−=−=−−−−=−∑11)1()1(1222)()()0()()(0)0()(])([)0()(])([k k k k nk k n n n n dt t f d t f f s s F s dt t f d L f sf s F s dt t f d L f s sF dt t df L M )( 2微分定理初始条件为0时)(])([s F s dtt f d L n nn =一般形式 }}∑∫∫∫∫∫∫∫∫∫∫∫==+−===+=++=+=nk t n n k n n n n t t t dt t f ss s F dt t f L sdt t f s dt t f s s F dt t f L sdt t f s s F dt t f L 101022022]))(([1)(])()([]))(([])([)(]))(([])([)(])([个共个共L L M3积分定理初始条件为0时}nnn s s F dt t f L )(]))(([=∫∫个共L4 延迟定理(或称t 域平移定理))()](1)([s F e T t T t f L Ts−=−− 5 衰减定理(或称s 域平移定理))(])([a s F e t f L at +=−6 终值定理 )(lim )(lim 0s sF t f st →∞→=7 初值定理 )(lim )(lim 0s sF t f s t ∞→→=8 卷积定理)()(])()([])()([21021021s F s F d t f t f L d f t f L tt =−=−∫∫τττττ4202.表A-2 常用函数的拉氏变换和z 变换表 序号 拉氏变换E(s)时间函数e(t) Z 变换E(z)1 1δ(t)12 Ts e −−11∑∞=−=0)()(n T nT tt δδ 1−z z3 s1 )(1t 1−z z 421s t2)1(−z Tz5 31s 22t32)1(2)1(−+z z z T6 11+n s!n t n )(!)1(lim 0aTn n n a e z z a n −→−∂∂− 7 as +1ate −aTe z z−− 8 2)(1a s + atte− 2)(aT aT e z Tze −−−9 )(a s s a +ate−−1 ))(1()1(aT aT e z z z e −−−−− 10 ))((b s a s ab ++− btate e−−−bTaT e z z e z z −−−−− 11 22ωω+st ωsin 1cos 2sin 2+−T z z Tz ωω12 22ω+s st ωcos1cos 2)cos (2+−−T z z T z z ωω 13 22)(ωω++a s t eatωsin −aTaT aT eT ze z T ze 22cos 2sin −−−+−ωω 14 22)(ω+++a s a st e at ωcos −aTaT aTe T ze z T ze z 222cos 2cos −−−+−−ωω15aT s ln )/1(1−T t a/az z −4213. 用查表法进行拉氏反变换用查表法进行拉氏反变换的关键在于将变换式进行部分分式展开,然后逐项查表进行反变换。
常用拉普拉斯变换及反变换在数学和工程领域中,拉普拉斯变换是一种非常有用的工具,它可以将时域中的函数转换为复频域中的函数,从而使一些复杂的微分方程和积分方程的求解变得更加简单。
接下来,让我们一起深入了解一下常用的拉普拉斯变换及反变换。
拉普拉斯变换的定义是对于一个实值函数\(f(t)\),其拉普拉斯变换\(F(s)\)定义为:\F(s) =\int_{0}^{\infty} f(t) e^{st} dt\其中\(s =\sigma +j\omega\)是一个复变量,\(\sigma\)是实部,\(\omega\)是虚部,\(j\)是虚数单位。
下面我们来看一些常见函数的拉普拉斯变换:单位阶跃函数\(u(t)\),当\(t < 0\)时,\(u(t) = 0\);当\(t \geq 0\)时,\(u(t) = 1\)。
它的拉普拉斯变换为:\\mathcal{L}u(t) =\frac{1}{s}\指数函数\(e^{at}\),其拉普拉斯变换为:\\mathcal{L}e^{at} =\frac{1}{s a}\正弦函数\(sin(\omega t)\)的拉普拉斯变换为:\\mathcal{L}sin(\omega t) =\frac{\omega}{s^2 +\omega^2}\余弦函数\(cos(\omega t)\)的拉普拉斯变换为:\\mathcal{L}cos(\omega t) =\frac{s}{s^2 +\omega^2}\这些常见函数的拉普拉斯变换在解决实际问题中经常会用到。
那么,拉普拉斯反变换又是什么呢?拉普拉斯反变换就是将复频域中的函数\(F(s)\)转换回时域中的函数\(f(t)\)。
拉普拉斯反变换的计算通常比较复杂,但是对于一些常见的形式,我们可以通过一些方法来求解。
例如,对于形如\(F(s) =\frac{A}{s a}\)的函数,其反变换为\(f(t) = Ae^{at}\)。
成绩评定表课程设计任务书目录1.Matlab介绍............... 错误!未定义书签。
2.利用Matlab实现信号的复频域分析—拉普拉斯变化和拉普拉斯逆变换的设计 (5)2.1.拉普拉斯变换曲面图的绘制 (5)2.2.拉普拉斯变化编程设计及实现 (7)2.3.拉普拉斯逆变化编程设计及实现 (8)3.总结 (14)4.参考文献 (15)1.Matlab介绍MATLAB语言是当今国际上在科学界和教育界中最具影响力、也最具活力的软件;它起源于矩阵运算,现已发展成一种高度集成的计算机语言;它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、丰富的交互式仿真集成环境,以及与其他程序和语言便捷接口的功能。
经过多年的开发运用和改进,MATLAB已成为国内外高校在科学计算、自动控制及其他领域的高级研究工具。
典型的用途包括以下几个方面:1)数学计算;2)新算法研究开发;3)建模、仿真及样机开发;4)数据分析、探索及可视化;5)科技与工程的图形功能;6)友好图形界面的应用程序开发。
1.1Matlab入门Matlab7.0介绍Matlab7.0比Matlab的老版本提供了更多更强的新功能和更全面、更方便的联机帮助信息。
当然也比以前的版本对于软件、硬件提出了更高的要求。
在国内外Matlab已经经受了多年的考验。
Matlab7.0功能强大,适用范围很广。
其可以用来线性代数里的向量、数组、矩阵运算,复数运算,高次方程求根,插值与数值微商运算,数值积分运算,常微分方程的数值积分运算、数值逼近、最优化方法等,即差不多所有科学研究与工程技术应用需要的各方面的计算,均可用Matlab来解决。
MATLAB7.0提供了丰富的库函数(称为M文件),既有常用的基本库函数,又有种类齐全、功能丰富多样的的专用工具箱Toolbox函数。
函数即是预先编制好的子程序。
在编制程序时,这些库函数都可以被直接调用。
无疑,这会大大提高编程效率。
拉普拉斯变换、连续时间系统的S 域分析基本要求通过本章的学习,学生应深刻理解拉普拉斯变换的定义、 收敛域的概念:熟练掌握拉普拉斯变换的性质、卷积定理的意义及它们的运用。
能根据时域电路模型画出S 域等效电路模型,并求其冲激响应、零输入响应、零状态响应和全响应。
能根据系统函数的零、极点分 布情况分析、判断系统的时域与频域特性。
理解全通网络、最小相移网络的概念以及拉普拉 斯变换与傅里叶变换的关系。
会判定系统的稳定性。
知识要点1. 拉普拉斯变换的定义及定义域(1) 定义 单边拉普拉斯变换:st正变换 [f(t)] F(s) 0 f(t)e dt双边拉普拉斯变换:的收敛域。
0与函数f(t)的性质有关。
2. 拉普拉斯变换的性质逆变换[F(s)] f(t)stF(s)e正变换F B(S )f(t)edt1 jst逆变换 f(t)2 jjF B(s)eds(2)定义域若0 时,lim f (t)et0则St 「 ” ”t ”f(t)e 在0的全部范围内收敛,积分0就是f(t)的单边拉普拉斯变换st[f2(t)] F2(S) , 1 , 2 为常数(2 ) 原函数微分若[f (t)] F(s)则[響]sF(s) f(0 ) dt[d df>] s n F(s) n1s nr1f(r)(0 ) dt r 0r式中f⑴(0 )是r阶导数在0时刻的取值。
dt r(3)原函数积分(4)延时性F (s),则[f(t t°)u(t t。
)] e st0F(s)(5)s域平移at若[f (t)] F (s),则[f(t)e ] F(s a)(6)尺度变换1 s若[f (t)] F (s),则[f (at)] F( )( a 0)a a(7)初值定理lim f (t) f(0 ) limsF(s)to s(8)终值定理lim f (t) lim sF(s)t s(9)卷积定理若[f1(t)] F1(s),[f2(t)] F2(S),则有[f1(t) f2(t)] F1(S)F2(S) (1) 线性性[仏“⑴]1 1肓[h(s) F2(s)] = ^-j.h(p)F2(s p)dpj[i f l(t) 2f2(t)] 1F1(S) 2F2(S)t若[f (t)] F (s),则[f(t)dt] F(s)s3 式中f(D(0)s f(t)dt若[f (t)]3.拉普拉斯逆变换(1 ) 部分分式展开法首先应用海维赛展开定理将F (s)展开成部分分式,然后将各部分分式逐项进行逆变换,最后叠加起来即得到原函数 f (t)。
第十二章 拉普拉斯变换及逆变换拉普拉斯(Laplace)变换就是分析与求解常系数线性微分方程得一种简便得方法,而且在自动控制系统得分析与综合中也起着重要得作用。
我们经常应用拉普拉斯变换进行电路得复频域分析。
本章将扼要地介绍拉普拉斯变换(以下简称拉氏变换)得基本概念、主要性质、逆变换以及它在解常系数线性微分方程中得应用。
第一节 拉普拉斯变换在代数中,直接计算328.957812028.6⨯⨯=N 53)164.1(⨯就是很复杂得,而引用对数后,可先把上式变换为164.1lg 53)20lg 28.9lg 5781(lg 3128.6lg lg ++-+=N然后通过查常用对数表与反对数表,就可算得原来要求得数N 。
这就是一种把复杂运算转化为简单运算得做法,而拉氏变换则就是另一种化繁为简得做法。
一、拉氏变换得基本概念定义12、1 设函数()f t 当0t ≥时有定义,若广义积分()pt f t e dt +∞-⎰在P 得某一区域内收敛,则此积分就确定了一个参量为P 得函数,记作()F P ,即dte tf P F pt ⎰∞+-=)()( (12、1)称(12、1)式为函数()f t 得拉氏变换式,用记号[()]()L f t F P =表示。
函数()F P 称为()f t 得拉氏变换(Laplace) (或称为()f t 得象函数)。
函数()f t 称为()F P 得拉氏逆变换(或称为()F P 象原函数),记作)()]([1t f P F L =-,即)]([)(1P F L t f -=。
关于拉氏变换得定义,在这里做两点说明:(1)在定义中,只要求()f t 在0t ≥时有定义。
为了研究拉氏变换性质得方便,以后总假定在0t <时,()0f t =。
(2)在较为深入得讨论中,拉氏变换式中得参数P 就是在复数范围内取值。
为了方便起见,本章我们把P 作为实数来讨论,这并不影响对拉氏变换性质得研究与应用。
常用拉普拉斯变换及反变换拉普拉斯变换在工程和数学中是个非常实用的工具。
它不仅能帮助我们解决微分方程,还能简化许多复杂的问题。
今天我们就来聊聊常用的拉普拉斯变换和反变换,看看它们是如何发挥作用的。
一、拉普拉斯变换的基本概念1.1 定义拉普拉斯变换是一个积分变换,它将时间域的函数转换为复频域的函数。
简单来说,它把一个函数从“时间的世界”带到了“频率的世界”。
公式上,拉普拉斯变换可以表示为:\[ \mathcal{L}\{f(t)\} = F(s) = \int_0^{\infty} e^{-st} f(t) dt \]这里的 \( s \) 是复数变量,\( f(t) \) 是我们要变换的时间域函数,\( F(s) \) 则是变换后的结果。
1.2 性质拉普拉斯变换有几个重要的性质,比如线性性、时间延迟和微分等。
这些性质使得在实际应用中,可以灵活地对待不同类型的函数。
例如,线性性让我们可以把两个函数的变换简单相加,这对于解决复杂问题很有帮助。
二、常用的拉普拉斯变换2.1 单位阶跃函数单位阶跃函数 \( u(t) \) 是拉普拉斯变换中最常用的函数之一。
它的变换结果是:\[ \mathcal{L}\{u(t)\} = \frac{1}{s} \]这个简单的公式为很多工程应用奠定了基础,因为很多信号和系统可以用阶跃函数来描述。
2.2 指数函数另一个常见的函数是指数函数 \( e^{at} \)。
它的拉普拉斯变换结果为:\[ \mathcal{L}\{e^{at}\} = \frac{1}{s - a} \]这在处理自然衰减或增长的过程时特别有用,比如在电子电路中,我们经常会遇到这种情况。
2.3 正弦和余弦函数正弦和余弦函数的拉普拉斯变换也很重要。
它们分别为:\[ \mathcal{L}\{\sin(\omega t)\} = \frac{\omega}{s^2 + \omega^2} \] \[ \mathcal{L}\{\cos(\omega t)\} = \frac{s}{s^2 + \omega^2} \]这些变换结果在振动分析和控制系统中应用广泛,帮助我们理解系统的频率响应。
完整版拉普拉斯变换表拉普拉斯变换是一种重要的数学工具,它能将时间域上的函数转换为频率域上的函数,为信号处理、电路分析等领域的数学建模和分析提供了极大的便利。
下面是完整版的拉普拉斯变换表,列出了常用的函数及其对应的拉普拉斯变换公式。
1. 常数函数:f(t) = 1,其拉普拉斯变换为:F(s) = 1/s2. 单位阶跃函数:f(t) = u(t),其拉普拉斯变换为:F(s) = 1/s3. 指数函数:f(t) = e^-at,其拉普拉斯变换为:F(s) = 1/(s + a)4. 正弦函数:f(t) = sin(ωt),其拉普拉斯变换为:F(s) = ω/(s^2 + ω^2)5. 余弦函数:f(t) = cos(ωt),其拉普拉斯变换为:F(s) = s/(s^2 + ω^2)6. 指数衰减正弦函数:f(t) = e^-at sin(ωt),其拉普拉斯变换为:F(s) = ω/( (s+a)^2 + ω^2 )7. 指数衰减余弦函数:f(t) = e^-at cos(ωt),其拉普拉斯变换为:F(s) = (s+a)/( (s+a)^2 + ω^2 )8. 阻尼正弦函数:f(t) = e^-αt sin(ωt),其拉普拉斯变换为:F(s) = ω/( (s+α)^2 + ω^2 )9. 阻尼余弦函数:f(t) = e^-αt cos(ωt),其拉普拉斯变换为:F(s) = (s+α)/( (s+α)^2 + ω^2 )10. 给定函数f(t)的导数Laplace变换:f'(t) 的Laplace 变换 F(s) 为:F(s)= s*F(s) - f(0)11. 给定函数f(t)的不定积分Laplace变换:∫f(t)dt 的 Laplace 变换 F(s) 为:F(s)= 1/s*F(s)12. Laplace变换与乘法定理:L{f(t) g(t)} = F(s)G(s)13. Laplace变换与移位定理:L{f(t-a) u(t-a)} = e^-as F(s)14. Laplace变换与初值定理:f(0+) = lims→∞ sF(s)f'(0+) = lims→∞ s^2F(s) - sf(0+)f''(0+) = lims→∞ s^3F(s) - s^2f(0+) - sf'(0+)15. Laplace变换与终值定理:limt→∞ f(t) = lims→0 sF(s)limt→∞ f'(t) = lims→0 s^2F(s) - sf(0+)limt→∞ f''(t) = lims→0 s^3F(s) - s^2f(0+) -sf'(0+)这是完整版的拉普拉斯变换表,其中列出了常用的函数及其对应的拉普拉斯变换公式,以及常见的拉普拉斯变换定理和公式。
第四章 拉普拉斯变换—连续信号s 域分析一、考试内容(知识点)1.拉普拉斯变换的定义及其性质、拉普拉斯逆变换; 2.系统的复频域分析法; 3.系统函数)(s H ;4.系统的零极点分布决定系统的时域、频域特性; 5.线性系统的稳定性;6.拉普拉斯变换与傅里叶变换之间的关系。
二、内容(知识点)详解1.拉普拉斯变换的定义、收敛域(1)变换式与反变换式dt e t f t f s F st -∞⎰-==0)()]([)(L ds e s F js F t f stj j ⎰∞+∞--==σσπ)(21)]([)(1L )(s F 称为)(t f 的象函数,)(t f 称为)(s F 的原函数。
下限值取-0,主要是考虑信号)(t f 在t =0时刻可能含有冲激函数及其导数项也能包含在积分区间之内。
(2)收敛域在s 平面上,能使式0)(lim =-→∞t t e t f σ满足和成立的σ的取值范围(区域),称为)(t f 或)(s F 的收敛域。
2.常用时间函数的拉普拉斯变换(1)冲激函数 )()(t t f δ= 1)(=s F)()()(t t f n δ= n s s F =)((2)阶跃函数 )()(t u t f = ss F 1)(= (3)n t (n 是正整数) t t f =)( 21)(s s F =2)(t t f = 32)(s s F =n t t f =)( 1!)(+=n s n s F(4)指数信号 t e t f α-=)( α+=s s F 1)(t te t f α-=)( ()21)(α+=s s F t n e t t f α-=)( ()1!)(++=n s n s F αt j e t f ω-=)( ωj s s F +=1)( (5)正弦信号、余弦信号系列)sin()(t t f ω= 22)(ωω+=s s F)cos()(t t f ω= 22)(ω+=s ss F)sin()(t e t f t ωα-= 22)()(ωαω++=s s F)cos()(t e t f t ωα-= 22)()(ωαα+++=s s s F )sin()(t t t f ω= 222)(2)(ωω+=s ss F )cos()(t t t f ω= 22222)()(ωω+-=s s s F )()(t sh t f ω= 22)(ωω-=s s F )()(t ch t f ω= 22)(ω-=s ss F (6) ∑∞=-=0)()(n nT t t f δ sT e s F --=11)(∑∞=-=00)()(n nT t f t f sTes F s F --=1)()(0 3.拉普拉斯变换的基本性质象函数)(s F 与原函数)(t f 之间的关系为:)]([)(t f s F L = (1)线性(叠加性)∑∑===⎥⎦⎤⎢⎣⎡ni i i n i i i s F a t f a 11)()(L ,其中i a 为常数,n 为正整数。
成绩评定表课程设计任务书目录1.Matlab介绍............... 错误!未定义书签。
2.利用Matlab实现信号的复频域分析—拉普拉斯变化和拉普拉斯逆变换的设计 (5)2.1.拉普拉斯变换曲面图的绘制 (5)2.2.拉普拉斯变化编程设计及实现 (7)2.3.拉普拉斯逆变化编程设计及实现 (8)3.总结 (14)4.参考文献 (15)1.Matlab介绍MATLAB语言是当今国际上在科学界和教育界中最具影响力、也最具活力的软件;它起源于矩阵运算,现已发展成一种高度集成的计算机语言;它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、丰富的交互式仿真集成环境,以及与其他程序和语言便捷接口的功能。
经过多年的开发运用和改进,MATLAB已成为国内外高校在科学计算、自动控制及其他领域的高级研究工具。
典型的用途包括以下几个方面:1)数学计算;2)新算法研究开发;3)建模、仿真及样机开发;4)数据分析、探索及可视化;5)科技与工程的图形功能;6)友好图形界面的应用程序开发。
1.1Matlab入门Matlab7.0介绍Matlab7.0比Matlab的老版本提供了更多更强的新功能和更全面、更方便的联机帮助信息。
当然也比以前的版本对于软件、硬件提出了更高的要求。
在国内外Matlab已经经受了多年的考验。
Matlab7.0功能强大,适用范围很广。
其可以用来线性代数里的向量、数组、矩阵运算,复数运算,高次方程求根,插值与数值微商运算,数值积分运算,常微分方程的数值积分运算、数值逼近、最优化方法等,即差不多所有科学研究与工程技术应用需要的各方面的计算,均可用Matlab来解决。
MATLAB7.0提供了丰富的库函数(称为M文件),既有常用的基本库函数,又有种类齐全、功能丰富多样的的专用工具箱Toolbox函数。
函数即是预先编制好的子程序。
在编制程序时,这些库函数都可以被直接调用。
无疑,这会大大提高编程效率。
第十二章 拉普拉斯变换及逆变换拉普拉斯(Laplace)变换是分析和求解常系数线性微分方程的一种简便的方法,而且在自动控制系统的分析和综合中也起着重要的作用。
我们经常应用拉普拉斯变换进行电路的复频域分析。
本章将扼要地介绍拉普拉斯变换(以下简称拉氏变换)的基本概念、主要性质、逆变换以及它在解常系数线性微分方程中的应用。
第一节 拉普拉斯变换在代数中,直接计算328.957812028.6⨯⨯=N 53)164.1(⨯是很复杂的,而引用对数后,可先把上式变换为164.1lg 53)20lg 28.9lg 5781(lg 3128.6lg lg ++-+=N然后通过查常用对数表和反对数表,就可算得原来要求的数N 。
这是一种把复杂运算转化为简单运算的做法,而拉氏变换则是另一种化繁为简的做法。
一、拉氏变换的基本概念定义12.1 设函数()f t 当0t ≥时有定义,若广义积分()pt f t e dt +∞-⎰在P 的某一区域内收敛,则此积分就确定了一个参量为P 的函数,记作()F P ,即dte tf P F pt ⎰∞+-=0)()( (12.1)称(12.1)式为函数()f t 的拉氏变换式,用记号[()]()L f t F P =表示。
函数()F P 称为()f t 的拉氏变换(Laplace) (或称为()f t 的象函数)。
函数()f t 称为()F P 的拉氏逆变换(或称为()F P 象原函数),记作)()]([1t f P F L =-,即)]([)(1P F L t f -=。
关于拉氏变换的定义,在这里做两点说明:(1)在定义中,只要求()f t 在0t ≥时有定义。
为了研究拉氏变换性质的方便,以后总假定在0t <时,()0f t =。
(2)在较为深入的讨论中,拉氏变换式中的参数P 是在复数范围内取值。
为了方便起见,本章我们把P 作为实数来讨论,这并不影响对拉氏变换性质的研究和应用。
成绩评定表课程设计任务书目录1.Matlab介绍.............. 错误!未定义书签。
2.利用Matlab实现信号的复频域分析—拉普拉斯变化和拉普拉斯逆变换的设计 (5)2.1.拉普拉斯变换曲面图的绘制 (5)2.2.拉普拉斯变化编程设计及实现 (7)2.3.拉普拉斯逆变化编程设计及实现 (8)3.总结 (14)4.参考文献 (15)1.Matlab介绍MATLAB语言是当今国际上在科学界和教育界中最具影响力、也最具活力的软件;它起源于矩阵运算,现已发展成一种高度集成的计算机语言;它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、丰富的交互式仿真集成环境,以及与其他程序和语言便捷接口的功能。
经过多年的开发运用和改进,MATLAB已成为国内外高校在科学计算、自动控制及其他领域的高级研究工具。
典型的用途包括以下几个方面:1)数学计算;2)新算法研究开发;3)建模、仿真及样机开发;4)数据分析、探索及可视化;5)科技与工程的图形功能;6)友好图形界面的应用程序开发。
1.1Matlab入门Matlab7.0介绍Matlab7.0比Matlab的老版本提供了更多更强的新功能和更全面、更方便的联机帮助信息。
当然也比以前的版本对于软件、硬件提出了更高的要求。
在国内外Matlab已经经受了多年的考验。
Matlab7.0功能强大,适用范围很广。
其可以用来线性代数里的向量、数组、矩阵运算,复数运算,高次方程求根,插值与数值微商运算,数值积分运算,常微分方程的数值积分运算、数值逼近、最优化方法等,即差不多所有科学研究与工程技术应用需要的各方面的计算,均可用Matlab来解决。
MATLAB7.0提供了丰富的库函数(称为M文件),既有常用的基本库函数,又有种类齐全、功能丰富多样的的专用工具箱Toolbox函数。
函数即是预先编制好的子程序。
在编制程序时,这些库函数都可以被直接调用。
无疑,这会大大提高编程效率。
MATLAB7.0的基本数据编程单元是不需要指定维数的复数矩阵,所以在MATLAB环境下,数组的操作都如数的操作一样简单方便。
而且,MATLAB7.0界面友好,用户使用方便。
首先,MATLAB具有友好的用户界面与易学易用的帮助系统。
用户在命令窗里通过help 命令可以查询某个函数的功能及用法,命令的格式极为简单。
其次,MATLAB 程序设计语言把编辑、编译、连接、执行、调试等多个步骤融为一体,操作极为简单。
除此之外,MATLAB7.0还具有强大的图形功能,可以用来绘制多姿多彩的图形,直观而形象。
综上,在进行信号的分析与仿真时,MATLAB7.0无疑是一个强大而实用的工具。
尤其对于信号的分析起到了直观而形象的作用,非常适合与相关课题的研究与分析2 利用Matlab 实现信号的复频域分析—拉普拉斯变化和拉普拉斯逆变换的设计2.1 拉普拉斯变换曲面图的绘制连续时间信号)(t f 的拉普拉斯变换定义为:⎰+∞-=0)()(dt e t f s F st(6-1)其中ωσj s +=,若以σ为横坐标(实轴),ωj 为纵坐标(虚轴),复变量s 就构成了一个复平面,称为s 平面。
显然,)(s F 是复变量s 的复函数,为了便于理解和分析)(s F 随s 的变化规律,可以将)(s F 写成:)()()(s j es F s F ϕ= (6-2)其中,)(s F 称为复信号)(s F 的模,而)(s ϕ则为)(s F 的幅角。
从三维几何空间的角度来看,)(s F 和)(s ϕ对应着复平面上的两个平面,如果能绘出它们的三维曲面图,就可以直观地分析连续信号的拉普拉斯变换)(s F 随复变量s 的变化规律。
上述过程可以利用MATLAB 的三维绘图功能实现。
现在考虑如何利用MATLAB 来绘制s 平面的有限区域上连续信号)(t f 的拉普拉斯变换)(s F 的曲面图,现以简单的阶跃信号)(t u 为例说明实现过程。
我们知道,对于阶跃信号)()(t u t f =,其拉普拉斯变换为ss F 1)(=。
首先,利用两个向量来确定绘制曲面图的s 平面的横、纵坐标的范围。
例如可定义绘制曲面图的横坐标范围向量x1和纵坐标范围向量y1分别为:x1=-0.2:0.03:0.2; y1=-0.2:0.03:0.2;然后再调用meshgrid()函数产生矩阵s ,并用该矩阵来表示绘制曲面图的复平面区域,对应的MATLAB 命令如下:[x,y]=meshgrid(x1,y1); s=x+i*y;上述命令产生的矩阵s 包含了复平面2.02.0<<-σ, 2.02.0<<-ωj 范围内以时间间隔0.03取样的所有样点。
最后再计算出信号拉普拉斯变换在复平面的这些样点上的值,即可用函数mesh()绘出其曲面图,对应命令为:fs=abs(1./s); mesh(x,y,fs); surf(x,y,fs);title('单位阶跃信号拉氏变换曲面图'); colormap(hsv);axis([-0.2,0.2,-0.2,0.2,0.2,60]); rotate3d;执行上述命令后,绘制的单位阶跃信号拉普拉斯变换曲面图如图1所示。
2.2 拉普拉斯变化编程设计及实现已知连续时间信号)()sin()(t u t t f =,求出该信号的拉普拉斯变换,并利用MATLAB 绘制拉普拉斯变换的曲面图。
解:该信号的拉普拉斯变换为:11)(2+=s s F利用上面介绍的方法来绘制单边正弦信号拉普拉斯变换的曲面图,实现过程如下:绘制单边正弦信号拉普拉斯变换曲面图程序图2 单边正弦信号拉氏变换曲面图clf;a=-0.5:0.08:0.5; b=-1.99:0.08:1.99; [a,b]=meshgrid(a,b);d=ones(size(a)); c=a+i*b;%确定绘制曲面图的复平面区域c=c.*c; c=c+d; c=1./c; c=abs(c);%计算拉普拉斯变换的样值 mesh(a,b,c);%绘制曲面图surf(a,b,c);axis([-0.5,0.5,-2,2,0,15]);title('单边正弦信号拉氏变换曲面图'); colormap(hsv);上述程序运行结果如图2所示。
2.3 拉普拉斯逆变化编程设计及实现连续信号)(t f 的拉普拉斯变换具有如下一般形式:∑∑====Li iiKj jjsd s c s D s C s F 11)()()(若L K ≥,则)(s F 可以分解为有理多项式与真分式之和,即∑∑==+=+=+=Ni i iMj jjs a sb s P s A s B s P s R s P s F 11)()()()()()()(其中,)(s P 是关于s 的多项式,其逆变换可直接求得(冲激信号及其各阶导数),)(s R 为关于s 的有理真分式,即满足N M <。
以下进讨论N M <的情况。
设连续信号)(t f 的拉普拉斯变换为)(s F ,则∏=-==Ni ip s s B s A s B s F 1)()()()()(在满足N M <情况下,有以下几种情况(1)极点均为单重情况下,可对其直接进行部分分式展开得:NNp s r p s r p s r s F -++-+-= 2211)(其中,),,2,1()()(N i s F p s r ip s i i =-==称为有理函数)(s F 的留数。
则)(s F 的拉普拉斯逆变换为:)()(1t u e r t f Ni ti p i ∑==(2)有k 重极点,设为1p ,则部分分式展开为)()()()()()(111112111s D s E p s K p s K p s K s F k k k +-++-+-=-i K 1可用下式求得[]11111)()()!1(1p s ki i i s F p s ds d i K =----= 则)(s F 的拉普拉斯逆变换为:)()()!()(211t u e r t u e t j k K t f N i ti p i kj ti p j k j ∑∑==-+-= (3)有共轭极点N Nt f p s r p s r p s r p s r s F -++-+-+-=32)(22211)(设)(s F 有一对共轭极点βαj p ±-=2,1,则θj p s er s F p s r 1111)()(=-==*12r r =由共轭极点所决定的两项复指数信号可以合并成一项,故有)()cos(2)(12t u t e r t f tθβα+=-从以上分析可以看出,只要求出)(s F 部分分式展开的系数(留数)i r ,就可直接求出)(s F 的逆变换)(t f 。
上述求解过程,可以利用MATLAB 的residue()函数来实现。
令A 和B 分别为)(s F 的分子和分母多项式构成的系数向量,则函数:[r,p,k]=residue(B,A)将产生三个向量r 、p 和k ,其中p 为包含)(s F 所有极点的列向量,r 为包含)(s F 部分分式展开系数i r 的列向量,k 为包含)(s F 部分分式展开的多项式的系数行向量,若N M <,则k 为空。
例:已知连续信号的拉普拉斯变换为:ss s s F 442)(3++=试用MATLAB 求其拉普拉斯逆变换)(t f 。
解:MATLAB 命令如下:a=[1 0 4 0];b=[2 4];[r,p,k]=residue(b,a)运行结果:r =-0.5000 - 0.5000i-0.5000 + 0.5000i1.0000p =0 + 2.0000i0 - 2.0000ik =[] 由上述结果可以看出,)(s F 有三个极点22,1j p ±=,03=p ,为了求得共轭极点对应的信号分量,可用abs()和angle()分别求出部分分式展开系数的模和幅角,命令如下:abs(r)ans =0.70710.70711.0000angle(r)/pians =-0.75000.75000 由上述结果可得)()]432cos(21[)(t u t t f π-+=。
例:求下式函数的逆变换3)1(2)(+-=s s s s F解:MATLAB 程序如下:a=[1 3 3 1 0];b=[1 -2];[r,p,k]=residue(b,a)运行结果:r =2.00002.00003.0000-2.0000p =-1.0000-1.0000-1.0000k =[] 则ss s s s F 2)1(3)1(2)1(2)(32-+++++=,对应的逆变换为)(]2)2223[()(2t u e t t t f t -++=-3. 总结通过本次综合实践让我们在学习“信号与系统”课程的同时,掌握MATLAB 的应用,对MATLAB 语言在中的推广应用起到促进作用。