第二章-卷积积分与积分变换1PPT课件

《第二章-卷积积分与积分变换1PPT课件》由会员分享，可在线阅读，更多相关《第二章-卷积积分与积分变换1PPT课件（84页珍藏版）》请在装配图网上搜索。

1、工程应用数学工程应用数学主主 讲讲北京理工大学机械与车辆学院北京理工大学机械与车辆学院 李晓雷李晓雷第二章第二章 卷积积分与积分变换卷积积分与积分变换 本章主要复习傅立叶级数，介本章主要复习傅立叶级数，介绍傅立叶级数的复数表示，脉冲响绍傅立叶级数的复数表示，脉冲响应与卷积积分、傅立叶变换和拉普应与卷积积分、傅立叶变换和拉普拉斯变换及其性质。拉斯变换及其性质。2.1 2.1 傅立叶级数傅立叶级数 2.2 2.2 脉冲响应与卷积积分脉冲响应与卷积积分 2.3 2.3 傅立叶变换傅立叶变换 2.4 2.4 拉普拉斯变换拉普拉斯变换 2.1 2.1 傅立叶级数傅立叶级数 1 1 周期函数的傅立叶级数周

2、期函数的傅立叶级数 2 2 复数形式的傅立叶级数复数形式的傅立叶级数 1 1 周期函数的傅立叶级数周期函数的傅立叶级数 如果如果g(t)是以是以T 为周期的周期函数，即为周期的周期函数，即 g(t+T) = g(t) 并且并且g(t)还满足下列条件还满足下列条件(Dirichlet条件条件) 1. 在在- -T/2，T/2上连续上连续,或者只有有限个一类或者只有有限个一类间断点间断点(即存在极限的间断点即存在极限的间断点); 2. 只有有限个极值；只有有限个极值； 3. 在在- -T/2，T/2上绝对可积，即上绝对可积，即 /2/2| ( )|dTTg tt 则在则在-T/2，T/2上上g(t

3、)可以展成傅立叶级数，可以展成傅立叶级数，即有即有 这里：这里：000)sincos()(ppptpbtpatg00cos()pppcpt称为周期函数称为周期函数g(t)的基频的基频T20 p0(p=2，3，)称为基频称为基频的的p次倍频。次倍频。 a0=c0， b0=0 0=0 p=1,2, 2021( )dTTcg ttT2022( )cosdTpTag tpt tT2022( )sindTpTbg tpt tT222pppcabarctanpppbacospppacsinpppbcap、bp、cp和和p(p=0,1,2,)称为周期函数称为周期函数g(t)的傅立叶系数或谐波系数，的傅立叶系

4、数或谐波系数，cp是谐波的振幅是谐波的振幅,p是谐波的初相位。是谐波的初相位。 c1cos(t1)称为周期函数称为周期函数g(t)的基波；的基波； cpcos(ptp) （p1）称为周期函数）称为周期函数g(t)的的基波的基波的p次谐波。次谐波。2 2 复数形式的傅立叶级数复数形式的傅立叶级数单边傅立叶级数单边傅立叶级数 根据欧拉公式根据欧拉公式 有有 代入代入g(t)的傅立叶级数表达式的傅立叶级数表达式,有有 iecosisin0i0cosReeptpt0i()000( )cos()Reepptpppppg tcptc0i0ReeptppA这里这里 0 = 0= ap ibp = cp(co

5、sp isinp) 20021( )dTTAcg ttTiepppAc20022( )(cosisin)dTTg tptpttT02i22( )edTptTg ttT p=1,2, 00.511.520123456-505双边傅立叶级数双边傅立叶级数 由于对任意复数由于对任意复数z有有 因此，当因此，当p0时时 Re 22zzz 000iii1Reeee2ptptptpppAAA002ii21 2( )ed e2TptptTg ttT002ii21 2( )ed e2TptptTg ttT002ii21( )ed eTptptTg ttT002ii21( )ed eTptptTg ttT则则

6、而而令令02i211( )e(i)2TptpppTdg tdtabT02i21( )d TptpTdg t etT02i211( )e(i)2TptppTg tdtabT02i()21( )eTptpTdg tdtT 02i21( )eTptpTg tdtdT因而可以得到因而可以得到 即即000Reeee0iptiptiptpppAddp00ii0( )Reeeptptppppg tAd在频率轴的正半轴上，双边谱的系数在频率轴的正半轴上，双边谱的系数dp与单边谱与单边谱的系数的系数Ap之间有如下关系之间有如下关系 Ap= 2dpReAp= 2RedpImAp= 2Imdp通常用频谱图来直观显示

7、周期函数所包含的频通常用频谱图来直观显示周期函数所包含的频率成分及其大小。率成分及其大小。以频率以频率f (很少用很少用)为横轴，分别以为横轴，分别以cp(或或Ap)和和p为纵轴作图，并称为纵轴作图，并称f cp图为图为g(t)的幅频图，的幅频图，f p图为相频图。图为相频图。一般来说，频谱图多为单边的，只画出一般来说，频谱图多为单边的，只画出f 0的的部分。部分。周期函数频谱图的特点是只在离散点周期函数频谱图的特点是只在离散点0，f ，2f ，上有值，被称为离散谱，有时也形象，上有值，被称为离散谱，有时也形象地称为谱线图。地称为谱线图。 例：求下图所示周期方波信号例：求下图所示周期方波信号g

8、(t)的傅立叶级数的傅立叶级数102( )102Ttg tTt a0=c0=0 2022( )cosd0TpTag tpt tT2022( )sindTpTbg tpt tT0200202( 1)sindsindTTpt tpt tT02020000211cos( cos)TTptptTpp000211 cos() 1 cos()22TTppT p 利用利用 0T = 2 ，有，有21 cospbpp41,3,5,02,4,6,ppp即即000411( )sinsin3sin535g tttt014sin(21)nntn另外另外0002ii202( 1)eded TptptpTAttT00ii

9、022000211()iiptpt TTeeTpp 00ii22021(1)(1)iTTppeeT pii4(1)i22ppeep因此因此 2(1 cos)ipp41,3,5,i02,4,6,ppp000ii3i5444( )ReeReeReeii3i5tttg t0i(21)14Reei(21)ntnn同样同样 1(i)2pppdab21,3,5,i02,4,6,ppp0i(21)2( )ei(21)ntng tn-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.500.51-1.5-1-0.500.51

10、1.5-1-0.500.51-1-0.500.51-1-0.500.51-1.5-1-0.500.511.5-1-0.500.51-1-0.500.51-1-0.500.51-1.5-1-0.500.511.5-1-0.500.51-1-0.500.51-1-0.500.51-1.5-1-0.500.511.5-1-0.500.51-1-0.500.51-1-0.500.51-1.5-1-0.500.511.5-1-0.500.51-1-0.500.51-1-0.500.51-1.5-1-0.500.511.5-1-0.500.51-1-0.500.51周期函数的均方值周期函数的均方值 如果如

11、果g(t)是简谐函数，即是简谐函数，即 g(t)=Acos0t则则g(t)的均方值为的均方值为 22222021( )cosd2TTAgtAt tT显然，简谐函数的均方值只与它的振幅有关，显然，简谐函数的均方值只与它的振幅有关，与它的频率无关。与它的频率无关。如果如果g(t)是周期函数，则可展为傅立叶级数是周期函数，则可展为傅立叶级数g(t)的均方值为的均方值为 000( )(cossin)pppg taptbpt22221( )( )dTTgtgttT2200201(cossin) dTppTpaptbpttT200200001(cossin)cos()sin()dNTppTNpNpNpap

12、tbptTaNptbNptt根据三角函数的性质根据三角函数的性质2002cossind0TTmtnt t2002coscosd0TTmtnt t(mn) 2002sinsind0TTmtnt t(mn) 2202cosd2TTTpt t2202sind2TTTpt tp=1,2, p=1,2, 有有 2222220020( )(cossin)dTppTpgtaptbptt222011()2pppaab22011|2ppAA2|ppd上式为周期函数的巴塞伐上式为周期函数的巴塞伐(Parserval)公式。公式。它的右端为一无穷级数，级数的每一项是周期它的右端为一无穷级数，级数的每一项是周期函数的

13、谐波的均方值。函数的谐波的均方值。巴塞伐公式表明，周期函数的均方值是它的基巴塞伐公式表明，周期函数的均方值是它的基波和各次谐波的均方值之和。波和各次谐波的均方值之和。可以解释为，周期函数在时域内的总能量等于可以解释为，周期函数在时域内的总能量等于在频域内的总能量。在频域内的总能量。巴塞伐公式表明，在频域中，运动的总能量是巴塞伐公式表明，在频域中，运动的总能量是各谐波能量之和。换句话说，在频域中能量可各谐波能量之和。换句话说，在频域中能量可按谐波直接相加。按谐波直接相加。 2.2 2.2 脉冲响应与卷积积分脉冲响应与卷积积分 1 1 脉冲函数脉冲函数 2 2 脉冲响应脉冲响应 3 3 卷积积分卷

14、积积分 1 1 脉冲函数脉冲函数 脉冲函数又称为脉冲函数又称为函数函数,定义如下：定义如下： 它表示在它表示在t = 0时刻作用的一个幅值无穷大，但时刻作用的一个幅值无穷大，但冲量为冲量为1的脉冲力。力学上称为单位脉冲力。的脉冲力。力学上称为单位脉冲力。 1d)(00)(tttt-2.5-2-1.5-1-0.500.511.522.5-1-0.500.511.521t在在时刻作用的单位脉冲力可表示为：时刻作用的单位脉冲力可表示为： 1d)(0)(tttt在任意时刻在任意时刻作用的幅值为作用的幅值为 的脉冲力可以表示的脉冲力可以表示为为 意义：在意义：在时刻的一个力值无限大，但作用时时刻的一个力

15、值无限大，但作用时间为零的脉冲力。间为零的脉冲力。 F( )()F tFt其冲量为：其冲量为： ()dFttF函数性质：函数性质：(a) 如果如果F(t)是一个连续函数，则是一个连续函数，则 ( ) ()d( )F tttF(b)(t)为偶函数为偶函数，即，即 (t) = ( t)2 2 脉冲响应脉冲响应 例：设如图所示单自由度系统在例：设如图所示单自由度系统在t = 0以前静止，以前静止，在在t = 0时受到脉冲力时受到脉冲力(t)的激励。的激励。 这里：这里：0表示小于零但无限接近于零的时刻，表示小于零但无限接近于零的时刻， 0+表示大于零但无限接近于零的时刻。表示大于零但无限接近于零的时

16、刻。 运动微分方程运动微分方程 ( )( )(0 )0, (0 )0mxcxkxf tFtxx根据动量定理，在根据动量定理，在0到到0+这段时间内系统动量的这段时间内系统动量的改变为改变为mv(0+)mv(0) = 即在即在t=0时的脉冲力作用下，质量的速度时的脉冲力作用下，质量的速度v由由v(0)=0变为变为v(0+)= /m，而位移没有变化。，而位移没有变化。 FF当当t0后，系统不受外力，是自由振动。后，系统不受外力，是自由振动。 系统的运动微分方程为系统的运动微分方程为 =1时，即系统受到单位脉冲力作用时，系统响时，即系统受到单位脉冲力作用时，系统响应称为系统脉冲响应，用应称为系统脉冲

17、响应，用h(t)表示：表示： 0(0 )0,(0 )/mxcxkxxxF m( )esin(0)ntddFx tttmF1( )esin(0)ntddh tttm 在在t时刻以前静止的系统，在时刻以前静止的系统，在t时受一个单时受一个单位脉冲力激励后的响应为：位脉冲力激励后的响应为： 012345678910-8-6-4-202468x 10-3()1()esin()()ntddh tttm3 3 卷积积分卷积积分 以上例说明卷积积分。以上例说明卷积积分。 在振动系统受任意持续激励时，可把激励看为一在振动系统受任意持续激励时，可把激励看为一系列脉冲力的迭加。设在系列脉冲力的迭加。设在t = 时

18、刻附近的冲量。时刻附近的冲量。 ( )dFF它对它对t 上一定存在上一定存在,并且积分并且积分|( )|etg tM0( )( )edstG sg tt绝对可积且一致收敛，绝对可积且一致收敛，G(s)为解析函数。为解析函数。2 拉普拉斯变换的性质拉普拉斯变换的性质 1.线性线性 2.微分微分 LLa1g1(t)+a2g2(t)=a1LLg1(t)+a2LLg2(t)LL1a1G1()+a2G2()=a1LL1G1()+a2LL1G2() LLg(t) = sLLg(t) g(0)LLg(t) = s2LLg(t) sg(0) g(0) LLg(n)(t) = snLLg(t) = snLLg(

19、t) sn 1g(0) sn 2g(0) sg(n 2)(0) g(n 1)(0)(1)1(0)nnjjjsg3.积分积分 LL =s 1LLg(t)0( )dtg tt4.延时延时 LLg(t )= LLg(t) es5.位移位移 LL eatg(t) = G(s a)6.尺度尺度 LL g(at)=1( )sGaa7.初值与终值定理初值与终值定理 LLg(t) = G(s) 且且 存在存在 则则 (a)初值定理初值定理 如果如果 lim( )sG ss(0)lim( )lim( )0gg tsG sst(b)终值定理终值定理 如果如果 LLg(t) = G(s) 且且 存存在在 lim(

20、)g tt则则 ( )lim( )lim( )0gg tsG sts 8.卷积卷积 LLg1(t)*g2(t)=G1(s)G2(s)LL1G1(s)G2(s)=g1(t)*g2(t)LLg1(t)g2(t)=G1(s)*G2(s)g1(t)g2(t)=LL1G1(s)*G2(s)它的拉普拉斯变换为它的拉普拉斯变换为 解：如果解：如果g(t)是周期函数是周期函数,即有即有 例例 周期函数的拉普拉斯变换周期函数的拉普拉斯变换 g(t +T) = g(t) (t0)LLg(t) = 0( )edstg tt(1)0( )edjTstjTjg tt令令u=tjT，则，则 t=u+jT，dt=du。t

21、=jT时，时，u = 0；t = (j+1)T时，时，u = T，因此，因此 由于由于g(t)是周期函数，是周期函数，g(u+jT)g(u)，因此，因此 LLg(t) = ()00()edTs ujTjg ujTu00()eedTsujsTjg ujTu00e()edTjsTsujg ujTuLLg(t) = 00e( )edTjsTsujg uu01( )ed1 eTsusTg uu3 拉普拉斯逆变换的求法拉普拉斯逆变换的求法 由于由于g(t)的拉普拉斯变换就是的傅立叶变换的拉普拉斯变换就是的傅立叶变换， 根据傅立叶逆变换公式根据傅立叶逆变换公式,有有ii1( )e( )=( )e( )ed

22、 ed2ttg tu tgu (i)i01=( )ed ed2tg i1=(i )ed2tG由此可得由此可得 i1( ) ( )=(i )ede2ttg t u tG(+i)1=(i )ed2tG令令s=+i，则有，则有ds = id。当当时时，si，所以得到，所以得到 i-i1( )=( )e d(0)2 istg tG sst 例例 求求g(t) = sint的拉普拉斯变换的拉普拉斯变换 解：由于解：由于g(t)满足满足 LLg(t)=s2LLg(t) = 2LLg(t) 用同样的方法可以得到用同样的方法可以得到 g(0) = 0, g(0) =, g(t) = 2g(t) 因此有因此有 LLg(t)=LLsint =s22LLcost= ss22部分资料从网络收集整理而来，供大家参考，感谢您的关注！