![]() |
中国面包师贴吧-楼主(阅:3137/回:0)数字电路傅里叶变换电路1习题12-7 79.证明三角函数系: 1、cosωx、sinωx、cos2ωr、sin2ωr、...cosωr、sinωx,...,在[-T/2,T/2]上具有正交性,其中T=2π/ω. 将80-87题中周期为2π的周期函数f(x)展开成傅里叶级数,其中f(x)在[-π,π]上的表达式为: 80.f(x)=x 解: n+1 ∞ (-1) x=2 ∑ sin nx (-∞<x<+∞,x≠(2k+1)π,k=0,±1,±2,...) n=1 n π,-π≤x<0 81.f(x)={ x,0≤x<π 解: 3π 2 ∞ cos(2n-1)x ∞ 1 f(x)= - ∑ - ∑ sinnx 4 π n=1 2 n=1 n (2n-1) (-∞<x<+∞,x≠kπ,k=0,±1,±2,...) 2 82.f(x)=3x +1 解: n 2 2 ∞ (-1) cosnx 3x +1=π +1+12 ∑ (-∞<x<+∞) n=1 2 n x 83.f(x)=2sin 3 解: n x 18√3 ∞ (-1) n 2sin = ∑ sinnx 3 π n=1 n (-∞<x<+∞,x≠(2k+1)π,k=0,±1,±2,...) 85.f(x)=│x│ 解: π 4 ∞ cos(2n+1)x │x│= - ∑ (-∞<x<+∞) 2 π n=1 2 (2n-1) ax,-π≤x<0 86.f(x)={ (a,b为不等于零的常数,且a≠b) bx,0≤x<π 解: n+1 π 2(b-a) ∞ cos(2n+1)x (-1) f(x)= (b-a)- ∑ +(b+a) ∑ *sin nx 4 π n=1 2 n (2n-1) (-∞<x<+∞,x≠(2k+1)π,k=0,±1,±2,...) 0,-π≤x<-π/2 87.f(x)={ 1,-π/2≤x<π/2 0,π/2≤x<π 解: 1 2 ∞ 1 nπ f(x)= + ∑ sin cosnx 2 π n=1 n 2 (-∞<x<+∞,x≠(2k+1)π/2,k=0,±1,±2,...) 第八节、周期为T的周期函数的展开 上一节着重研究了将以2π为周期的周期函数展开成傅里叶级数的方法,它有比较普遍的应用价值,下面我们介绍以T(T为任意非零正常数)为周期的周期函数φ(t),在区间[-T/2,T/2)上展开成傅里叶级数问题。为了能按第七节的方法把它展成傅里叶级数,显然首先应当将φ(t)变换成以2π为周期,区间[-T/2,T/2)变换成[-π,π),为此我们作变量代换,令x=2πt/T,即t=Tx/2π,于是φ(t)=φ(Tx/2π)=f(x), 这时,函数f(x)就是以2π为周期的周期函数,假设在区间[-π,π]上满足收敛定理的条件。因此可以将它展成傅里叶级数,并且在连续点上有: a 0 n f(x)= + ∑ (a cosnt+b sinnt) 2 n=1 n n 其中傅里叶系数 1 π a = ∫ f(x)cosnxdx (n=0,1,2,...) n π -π { 1 π b = ∫ f(x)sinnxdx (n=1,2,3,...) n π -π 将变量x再换回成变量t,就得到周期为T的周期函数的傅里叶级数,并且在连续点上有 a 0 n 2nπt 2nπt f(x)= + ∑ (a cos +b sin ) ,(12.8.1) 2 n=1 n T n T 其中傅里叶系数为 2 T/2 2nπ a = ∫ φ(t)cos tdt (n=0,1,2,...) n T -T/2 T { (12.8.2) 1 T/2 2nπ b = ∫ φ(t)sin tdt (n=1,2,3,...) n π -T/2 T 在连续点处收敛于φ(t),在间断点处收敛于 φ(t-0)+φ(t+0) 2 如果以T为周期的周期函数φ(t),在(-T/2,T/2)内是奇函数,那么其傅里叶级数一定是正弦级数,且在连续点处有 n 2nπt φ(t)= ∑ b sin n=1 n T 这时 a =0 (n=0,1,2,...,) n { (12.8.3) 4 T/2 2nπ b = ∫ φ(t)sin tdt (n=1,2,3,...) n T 0 T 同样,以T为周期的周期函数φ(t),在(-T/2,T/2)内是偶函数,那么它的展开式必然是余弦级数,且在连续点处有 a 0 n 2nπt φ(t)= + ∑ a cos 2 n=1 n T 这时 4 T/2 2nπ a = ∫ φ(t)cos tdt (n=1,2,3,...) n T 0 T { (12.8.4) b =0 (n=0,1,2,...,) n 例1,若函数φ(t)以2为周期,在区间[-1,1)上的表达式为 1,-1≤t<0 φ(t)={ 2,0≤t<1 试将其展开成傅里叶级数 解、因为函数φ(t)满足收敛定理条件,且注意到T=2,故可由式(12.8.2)得 2 1 0 1 a = ∫ φ(t)cosnπtdt=∫ cosnπtdt+∫ 2cosnπtdt n 2 -1 -1 0 1 0 2 1 = [sinnπt] + [sinnπt] nπ -1 nπ 0 =0 (n=1,2,3,...,) 2 1 0 1 a = ∫ φ(t) dt=∫ dt+∫ 2dt=3 0 2 -1 -1 0 2 1 0 1 b = ∫ φ(t)sinnπtdt=∫ sinnπtdt+∫ 2sinnπtdt n 2 -1 -1 0 1 0 2 1 =- [cosnπt] - [cosnπt] nπ -1 nπ 0 2 ,n=1,3,5,... 1 n nπ =- [1-(-1) ]={ nπ 0,n=2,4,6,... 因此函数φ(t)的展开式为 3 2 1 1 φ(t)= + (sinπt+ sin3πt+ sin5πt+...) 2 π 3 5 (-∞<t<+∞,t≠k,k=0,±1,±2,...) 例2.试将图12-9中周期性锯齿波展开成傅里叶级数 解:显然首先应当给出这个锯齿波的解析表达式。根据解析几何的知识,我们不难得出周期函数φ(t)在一个周期内的表达式 2At φ(t)= T 由于φ(t)在(-T/2,T/2)内是奇函数,因此只需用式(12.8.3)来计算b n 4 T/2 2nπ 4 T/2 2A 2nπ b = ∫ φ(t)sin tdt= ∫ tsin tdt n T 0 T T 0 T T 4T 2nπ T/2 4A T/2 2nπ = [-tcos t] + ∫ cos tdt nπT T 0 nπT 0 T 2A n+1 = (-1) (n=1,2,3...) nπ 所以,所求的展开式为 2A 2π 1 4π 1 6π φ(t)= (sin t- sin t+ sin t-...) π T 2 T 3 T (-∞<t<+∞,t≠(2k+1)T/2,k=0,±1,±2,...) 当t=(2k+1T/2,(k=0,±1,±2,...)时,级数收敛于0. 在科技问题中,人们常将傅里叶级数用频率ω表示。这只需将ω=2π/T换入,于是上式可改写成 2A 1 1 φ(t)= (sinωt- sin2ωt+ sin3ωt-...) π 2 3 (-∞<t<+∞,t≠(2k+1)π/ω,k=0,±1,±2,...) 例3.若矩形波以T为周期,且在[-T/2,T/2)上表达式为 0,-T/2≤t<-T/4 φ(t)={ A,-T/4≤t<T/4 0,T/4≤t<T/2 试写出前五次谐波 解:这个函数的图形如图12-10所示,φ(t)是偶函数, 所以应利用(12.8.4)式计算傅里叶系数。 4 T/2 2nπ 4 T/2 2nπ a = ∫ φ(t)cos tdt= ∫ Acos tdt n T 0 T T 0 T 2A 2nπ T/4 = [sin ] nπ T 0 2A nπ = sin (n=1,2,3,...) nπ 2 4 T/2 4 T/4 a = ∫ φ(t) dt= ∫ Adt=A 0 T 0 T 0 因此,所求的傅里叶级数为: A 2A 2πt 1 6πt 1 10πt φ(t)= + (cos - cos + cos - ……) 2 π T 3 T 5 T (-∞<t<+∞,t≠(2k-1)T/4,k=0,±1,±2,...) 当t=(2k-t)T/4,(k=0,±1,±2,...)时,该级数收敛于A/2. 令ω=2π/T,利用三角函数公式,则可将φ(t)的傅里叶级数改成谐波的形式。 A 2A π 2A 3π 2A 5π φ(t)= + sin(ωt+ )+ sin(3ωt+ ) + sin(5ωt+ )+… 2 π 2 3π 2 5π 2 如果取五次谐波则有 A 2A π 2A 3π 2A 5π φ(t)≈ + sin(ωt+ )+ sin(3ωt+ ) + sin(5ωt+ ) 2 π 2 3π 2 5π 2 若令ωt=0,即t=0,则 A 2A 2A 2A φ(t)≈ + - + 2 π 3π 5π 这是一个满足莱布尼兹审敛法条件的交错级数的前四项和,显然它的误差 2A │r │< ≈0.09A 4 7π 实际上,由已知的φ(t)表达式知φ(0)=A,因此可以计算出,当取五次谐波时φ(0)≈1.05A,实际误差为0.05A,若想要减小误差,提高精度,就再多取几次谐波。 习题12-8 已知第88-93题中各周期函数在一个周期内的表达式,试将它们展开成傅里叶级数。 2 88.φ(t)=1-t ,(-1/2≤t<1/2) 解: n+1 2 11 1 ∞ (-1) 1-t = + ∑ cos2nπt 12 2 n=1 2 π n (-∞<t<+∞,) t+T/4,-T/2≤t<-T/4 89.φ(t)={ 0,-T/4≤t<T/4 t-T/4,T/4≤t<T/2 解: n+1 nπ nπ(-1) -2sin T ∞ 2 2nπt φ(t)= ∑ sin 2 n=1 2 T 2π n (-∞<t<+∞,t≠(2k+1)T/2,k=0,±1,±2,...) 2,-2≤t<0 90.φ(t)={ 2-t,0≤t<2 解: n 3 4 ∞ 1 (2n-1)πt 2 ∞ (-1) nπt φ(t)= + ∑ cos + ∑ sin 2 2 n=1 2 2 π n=1 n 2 π (2n-1) (-∞<t<+∞,t≠(2k+1)2,k=0,±1,±2,...) 91.φ(t)=t+│t│,(-T/2≤t<T/2) 解: n+1 T 2T ∞ 1 2(2n-1)πt T ∞ (-1) 2πxt t+│t│= - ∑ cos + ∑ +sin 4 2 n=1 2 T π n=1 n T π (2n-1) (-∞<t<+∞,t≠(2k+1)2,k=0,±1,±2,...) A,-T/2≤t<0 92.φ(t)={ -A,0≤t<T/2 解: 4A ∞ 1 2(2n-1)πt φ(t)= ∑ sin π n=1 2n-1 T (-∞<t<+∞,t≠kT/2,k=0,±1,±2,...) T/4,-T/2≤t<-T/4 93.φ(t)={ -t,-T/4≤t<0 t,0≤t<T/4 T/4,T/4≤t<T/2 解: nπ cos -1 3T T ∞ 2 2nπt φ(t)= + ∑ cos 16 2 n=1 2 T π n (-∞<t<+∞) 94将图12-11所示周期性三角波展开为傅里叶级数。 解: A 4A ∞ 1 2(2n-1)πt φ(t)= + ∑ cos 2 2 n=1 2 T π (2n-1) (-∞<t<+∞) 95.若锯齿波在一个周期内的函数表达式为 -2At -A -T/2≤t<0 T 93.φ(t)={ -2At +A 0≤t<T/2 T 试写出其他前五次谐波。 解: 2A 2πt 2A 4πt 2A 6πt 2A 8πt 2A 10πt φ(t)≈ sin + sin + sin + sin + sin π T 2π T 3π T 4π T 5π T 第九节、定义在有限区间上的函数的展开 在有些实际问题中,我们还会遇到定义在有限区间[a,b]上的函数φ(t)展开成傅里叶级数的问题。为了解决这类非周期函数的展开问题,我们根据定义区间[a,b]的不同情况来研究。 一、定义在对称区间[-T/2,T/2]上的情形 假设函数φ(t)仅定义在有限的对称区间[-T/2,T/2]上,而在此区间外无意义。为了能够利用第八节的方法将φ(t)展开成傅里叶级数,我们首先将φ(t)在区间(-∞,+∞)上做周期性延拓,也就是构造一个周期为T的周期函数φ`(t),使得φ`(t)与函数φ(t)在区间[-T/2,T/2]上有φ(t)=φ`(t), 假设φ`(t)满足收敛定理的条件,那么即可按(12.8.1)式将它展开成傅里叶级数,并在连续点上有 a 0 ∞ 2nπt 2nπt φ`(t)= + ∑ (a cos +b sin ) 2 n=1 n T n T 如果函数φ(t)在区间[-T/2,T/2]上连续,则在此[-T/2,T/2]上就有 a 0 ∞ 2nπt 2nπt φ(t)=φ*(t)= + ∑ (a cos +b sin ) (12.9.1) 2 n=1 n T n T 在t=±T/2时,收敛于 1 [φ(T/2-0)+φ(-T/2+0)] 2 在对称区间[-T/2,T/2]上的函数φ(t)的傅里叶展开式,其中的傅里叶系数仍用公式 2 T/2 2nπ a = ∫ φ(t)cos tdt (n=0,1,2,...,) n T -T/2 T (12.9.2) 2 T/2 2nπ b = ∫ φ(t)sin tdt (n=1,2,3,...,) n T -T/2 T 如果函数φ(t)在[-T/2,T/2]内有第一类间断点,在间断点处级数同样收敛于 φ(t-0)+φ(t+0) 2 读者不难看出,定义在对称区间[-T/2,T/2]上的函数φ(t)的展开式的计算,与周期为T的周期函数φ`(t)的展开式的计算完全相同,它们的区别主要在收敛域的确定上,周期函数φ`(t)的收敛域要在整个(-∞,+∞)内考虑,而函数φ(t)仅在有限区间[-T/2,T/2]上来研究。 2 例1.试将定义在[-π,π]上的函数f(x)=x 展开成傅里叶级数。 解.将f(x)在整个数轴上作延拓,如图12-12所示。由于在[-π,π]上f(x)为偶函数,因此 2 π 2 π 2 a = ∫ f(x)cosnxdx= ∫ x cosnxdx= n π 0 π 0 2 2 π 4 π = [x sinx] - ∫ xsinxdx nπ 0 nπ 0 2 2 π 4 π = [x sinx] - ∫ cosxdx 2 0 2 0 n π n π 2 n = (-1) (n=1,2,3,...) 2 n 2 2 π 2 π 2 2π a = ∫ f(x)dx= ∫ x dx= 0 π 0 π 0 3 b =0 (n=1,2,3,...) 0 于是f(x)的展开式(在连续点处)为 2 2 π 1 1 x = -4(cosx- cos2x+ cos3x-...) 3 2 2 2 3 由于f(x)在[-π,π]上连续,经延拓后x=±π为连续点, 2 因此傅里叶级数在收敛域[-π,π]上收敛于x 例2.试求函数φ(t)的傅里叶展开式,设其表达式为 At +A, -T/2≤x<0 T φ(t)={ At -A, 0≤x<T/2 T 其中T与A均为常数。 解.作出φ(t)延拓后的图形,如图12-13所示。当t=(2π+1)T/2时,函数为双值,可以取定其一值使之符合我们的要求,例如取定一值φ(t)为奇函数,但不管取得一值均不影响结果。由于φ(t)是奇函数,因此 a =0 (n=0,1,2,...) n 4 T/2 2nπ b = ∫ φ(t)sin -tdt= n T 0 T 4 T/2 At 2nπ = ∫ ( -A)sin tdt T 0 T T 2 At 2nπ T/2 2A T/2 2nπ = [-( -A)cos ] + ∫ cos tdt nπ T T 0 nπT 0 T A n 2A = (-1) - nπ nπ 3A - ,n=1,3,5,... nπ ={ A - ,n=2,4,6,... nπ 又因为函数φ(t)在[-T/2,T/2]内,当t=0时间断,并且φ(-T/2)≠φ(T/2),即沿拓后t=±T/2为间断点,所以傅里叶展开式为 A 2πt 1 4πt 3 6πt 1 8πt φ(t)= (3sin + sin + sin + sin +…) π T 2 T 3 T 4 T (-T/2<t<0,0<t<T/2) 当t=0和t=±T/2时,级数收敛于0 二、定义在区间[0,T/2]上的情形 如果需要将仅定义在区间[0,T/2]上的函数φ(x)展开成傅里叶级数, 则首先要构造一个定义在对称区间[-T/2,T/2]上的函数Φ(t),使用Φ(t)在区间[0,T/2]上有Φ(t)=φ(t)。然后按照上述的办法,将φ(t)在整个实数轴上周期延拓,再展成傅里叶级数。为了使计算简单,常常是将定义在区间[0,T/2]上的函数φ(t)在区间[-T/2,0]上进行奇延拓或偶延拓。奇延拓是指构造函数Φ(t),使之在区间[-T/2,T/2]上成为奇函数,这只需令 φ(t),0≤t≤T/2 Φ (t)={ s -φ(-t),-T/2≤t≤0 这时函数Φ (t)的展开式为正弦级数,其傅里叶系数为 s a =0,(n=0,1,2,...) n { 4 T/2 2nπ b = ∫ φ(t)sin tdt, (n=1,2,3,...) n T 0 T 在收敛域上确定收敛于什么时,仍用狄利克雷定理来判定 φ(t),0≤t≤T/2 Φ (t)={ s -φ(-t),-T/2≤t≤0 这时展开式为余弦级数,其傅立叶系数为 4 T/2 2nπ 4 2nπ a = ∫ Φ (t)cos tdt= ∫φ(t)cos tdt n T 0 c T T T { (n=0,1,2,...) b =0 (n=1,2,3,...) n 偶延拓的图形12-15,由于是偶函数,所以延拓出去后在端点处都连续,因此如果函数φ(t)在[0,T/2]内连续,展开成余弦级数的收敛域0≤t≤T/2上收敛于φ(t)。 例3.试将函数 2 x πx f(x)= - 4 2 在区间[0,π]上展开成余弦函数。解.按式(12.9.4)计算傅立叶级数,注意T=2π。 2 2 π x π a = ∫ ( - ) cosnxdx n π 0 4 2 2 2 x πx π 2 x π π = [( - ) sinnxdx] + [( - )cosnx] nπ 4 2 0 2 2 2 0 n π 2 π - -[sinnxdx] 2 0 n π 2 = - 2 n 2 2 2 π x πx x a = ∫ ( - )=- (n=1,2,3,...) 0 π 0 4 2 3 由于f(x)在(0,π)上连续,且延拓的函数在x=0,π处连续,因此 2 2 x πx π 1 1 - =- +cosx+ cos2x+ cos3x+… (0≤x≤π) 4 2 6 4 9 例4.试将函数 t,0≤t<T/4 φ(t)={ t-T/4,T/4≤t<T/2 展开成正弦级数 解.按式(12.9.3)计算傅立叶系数如下: 4 T/2 2nπ 4 T/2 2nπ b = ∫ φ(t)cos tdt= ∫ tsin tdt n T 0 T T 0 T 4 T/2 T 2nπ + ∫ (t- )sin tdt T 0 4 T 2 2nπt T/4 2 2nπt T/2 = [-tcos ] + [sin ] nπ T 0 2 2 T 0 n π T T 2nπ T/2 2 2nπt T/2 + [-(t- )cos ] + [sin ] 2nπ 4 T -T/2 2 2 T T/4 n π T nπ T nπ T n+1 T nπ =- cos + sin + (-1) - sin 2nπ 2 2 2 2 2nπ 2 2 2 n π n π T n+1 nπ = [(-1) -cos ] (n=1,2,3,...) 2nπ 2 所以, T 2πt 1 6πt 2 8πt φ(t)= (sin + sin - sin +...) 2nπ T 3 T 4 T 因为当t=T/4时,函数φ(t)间断,且延拓后的函数在x=T/2处间断,故该级数在0≤t<T/4,T/4<t<T/2收敛于φ(t),而当t=T/4时收敛于T/8,当t=T/2时收敛于0. 习题12-9 将第95-99题中所给函数在定义区间内展开成傅立叶级数. -x,-π≤x<0 96.f(x)={ 0,0≤x<π 解: n π 2 ∞ 1 ∞ (-1) f(x)= - ∑ cos(2π-1)x+∑ sinnx x∈(-π,π) 4 π n=1 2 n=1 n (2n-1) -E,-T≤T<0, 97.u(t)={ (E为常量) E,0≤T≤T/2 解: 4E ∞ 1 2(2n-1)πt u(t)= ∑ sin ,t∈(-T/2,0),(0,T/2) π n=1 2n-1 T x,-1≤x<0, 98.f(x)={ 1,0≤x≤1/2 -1,1/2≤x≤1 解: nπ n 1-2cos 1 1 ∞ 1-(-1) 2π nπ 1 ∞ 2 f(x)=- + ∑ { + sin } cosnπx+ ∑ sinnπx 4 2 n=1 2 n 2 π n=1 n π n x∈[-1,0),(0,1/2),(1/2,1] πt 99.φ(t)=cos ,-T/2≤t≤T/2 T 解: n πt 2 4 ∞ (-1) 2nπt cos = + ∑ cos t∈[-T/2,T/2] T π π n=1 2 T 1-4n 将第100、101题中函数展开成正弦级数。 1,0≤x<h 100.f(x)={ 1/2,x=h, (h为常量) 0,h<x≤π 解: 2 ∞ 1-cosnh f(x)= ∑ sinnx, x∈(0,π) π n=1 n 101.φ(t)=(T/4)-t,0≤t≤T/2 解: 4 T ∞ 1 4nπt -t= ∑ sin t∈(0,T/2) T π n=1 2n T 将102、103题中函数展开成余弦级数 -t,0≤t<T/4 102.φ(t)={ -T/4,T/4≤t<T/2 解: nπ 1-cos 3T T ∞ 2 2nπt φ(t)=- + ∑ cos t∈[0,T/2] 16 2 n=1 2 T n x,0≤x<L/2 103.f(t)={ (L为常数) L-x,L/2≤x<L 解: L 2L ∞ 1 nπ n nπx f(x)= + ∑ [2cos -1-(-1) ]cos x∈[0,L] 4 2 n=1 2 2 L π n 2 104.将函数f(x)=1-x (0≤x≤1/2)分别展开成正弦级数和余弦级数 解: 2 ∞ 7 2 1 ∞ 1 1-x = ∑ [ + ] sin2(2π-1)x+ ∑ sin4nx, x∈(0,1/2) n=1 2(2n-1)π 3 2π n=1 2n (2n-1) π n+1 2 11 1 ∞ (-1) 1-x = + ∑ cos2nπx, x∈[0,1/2] 12 2 n=1 2 π n -π/4,-π≤x<0 105.把函数f(x)={ π/4,0≤x≤π 展成傅立叶级数,并由它推出: π 1 1 1 (1) =1- + - +… 4 3 5 7 √3π 1 1 1 1 (2) =1- + - + -… 6 5 7 11 13 解: ∞ 1 f(x)= ∑ sin(2n-1)x n=1 2n-1 令x=π/2及x=π/3 第四部分拉普拉斯变换 下面内容可参见《工程数学》,林益主编,何涛,杨殿生,卢强编,高等教育出版2003年出版。 第三章,拉普拉斯变换 拉普拉斯(Laplace)变换是一种积分变换,它通过无限区间上的广义积分,实现在不同类函数之间的转换,从而达到简化问题,解决问题的目的。变换是数学中常用的方法,拉普拉斯变换在对信号处理,求解微分方程方面是强有力的工具。它被广泛应用到电学、声学、振动力学等学科中去。本章内容涉及到有关复变函数论中的概念与理论,在这里不做深究,对本章出现的运算,只需按一元函数微分、积分的法则进行即可。 3.1拉普拉斯变换 首先给出函数拉普拉斯变换与逆变换的概念。 定义:设函数f(t)在区间[0,+∞]上有定义,如果含复变量s的无穷积分 +∞ -st +∞ -st ∫ e f(t)dt= lim ∫ e f(t)dt, 0 T→+∞ 0 对s的某一取值范围是收敛的,则称 +∞ -st F(t)=£[f(t)]= ∫ e f(t)dt, 0 为函数f(t)的拉普拉斯变换,f(t)称为象原函数,F(t)成为象函数(其中x是复数,s=β+iω)。 -1 记住F(s)=£[f(t)];相应的称f(t)为F(t)的拉普拉斯逆变换,记作f(t)=£ [F(s)]; 注意:为研究方便,本章中讨论的函数f(t)总认为t<0时,f(t)=0。 定理:如果f(t)满足下列条件: (1)在t≥0的任一有限区间上分段连续。 (2)存在实常数a≥0和A>0,使得t充分大时,有 ut │f(t)│≤Ae 则f(t)的拉普拉斯变换在半平面Re(s)>a上存在。这里a称为f(t)的增长指数,当f(t)是有界函数时,可取a=0。 证:设s=β+iω,则Re(s)=β,当β>a时, +∞ -st +∞ -st +∞ -st ∫ │e f(t)│dt=∫ │e │*│f(t)│dt≤∫ e *Ae dt 0 0 0 +∞ -st A =A∫ e dt= 0 β-a 因广义积分 +∞ -st ∫ e f(t)dt 0 绝对收敛,从而收敛,故F(s)存在。 例1,求函数f(t)=1的拉普拉斯变换 解:由拉普拉斯变换的定义,有 +∞ -st F(t)=£[1]= ∫ e f(t)dt 0 1 -st +∞ =- e s 0 1 =- (Re(t)>0) s 例2.求函数f(t)=t的拉普拉斯变换 解:由拉普拉斯变换的定义有 +∞ -st F(t)=£[1]= ∫ e tdt 0 1 +∞ -st =- ∫ tde s 0 1 -st +∞ +∞ -st =- [te - ∫ e dt] s 0 0 1 = (Re(t)>0) 2 s n 例3.求函数f(t)=t 的拉普拉斯变换,其中n是正整数,Re(s)>0 解: n +∞ -st n F(t)=£[t ]= ∫ e t dt 0 1 n -st +∞ n +∞ -st n-1 =- t e - ∫ e t dt s 0 s 0 n n-1 n n-1 n-2 = £[t ]= * £[t ]=... s s s n! n! = £[1]= n n+1 s s at 例4.求函数f(t)=e 的拉普拉斯变换 at +∞ -st at F(t)=£[e ]= ∫ e e dt 0 +∞ -(s-a)t = ∫ e dt 0 t 1 -(s-a) +∞ =- e s-a 0 1 -(β-a)t -iωt = [1- lim (e *e ) ] s-a t→+∞ 1 = (Re(s)>a) s-a 通过上述几个例子的讨论,我们可以看出,对于整个区间[0,+∞]上有定义的不同的函数f(t), n 它们的拉普拉斯变换很可能在s的不同区域上存在,例如f(t)=1,f(t)=t (n是正整数)等函数 2 at t 只对Re(s)>0有定义;而f(t)=e 只对Re(s)>a有定义,此外,对于某些函数,例如f(t)=e 不满足定理的条件,它的拉普拉斯变换对于一切s值都不存在,由此可见,F(s)的定义域是随f(t)而定的。另外,若x为实数,F(s)的定义域可把上述结论中的Re(s)换为s即可,比如,若s为实数,则 +∞ -st F(t)=£[1]= ∫ e dt, 0 的定义域即为s>0。我们通常接触到的函数大部分都能满足定理的条件,其拉普拉斯变换存在,且存在域通常是某个半平面,因此进行拉普拉斯变换时,常常不表出其存在域,只有非常必要时才特别加以注明。一般常用函数的拉普拉斯变换可查拉普拉斯变换简表(见附表5),请读者查表重新计算上述各题。 习题3.1 1.利用拉普拉斯变换的定义,求下列函数的拉普拉斯变换。 2 3 (1)f(t)=t ,解:F(s)=6/s t tln4 (2)f(t)=e-4 =e-e ,解:F(s)=s-1/(s-ln4) -t 2 (3)f(t)=te ,解:F(s)=2/(s+1) 3,当0≤t≤2时 (4)f(t)={ -1,当2≤t≤4时 0,当t≥4时 3,当t≤π/2时 (5)f(t)={ cost,当t>π/2时 2.设f(t)是以2τ为周期的函数,且在一个周期内的表达式为 h,当0≤t<γ时 f(t)={ -h,当γ≤t<2γ时 求f(t)的拉普拉斯变换 3.2拉普拉斯变换的性质 利用定义直接借助于积分的计算来求一个函数的 象函数,在大多数情况下是极其困难的。为了简化各种函数的象函数的求解过程,就有必要介绍拉普拉斯变换的性质,利用这些性质,拉普拉斯变换才能称为解决实际问题的有力工具。 性质1(线性性质)设函数f (t),f (t)满足定理的条件, 1 1 则在它们的象函数的定义域的共同部分上有: £[C f (t)+C f (t)]=C £[f (t)]+C £[f (t)] 1 1 2 2 1 1 2 2 其中,C ,C 是任意常数。 1 2 证:由定义 +∞ -st £[C f (t)+C f (t)]= ∫ e [C f (t)+C f (t)]dt 1 1 2 2 0 1 1 2 2 +∞ -st +∞ -st =∫ e C f (t)dt+∫ e C f (t)dt 0 1 1 2 2 +∞ -st +∞ -st =C ∫ e f (t)dt+C ∫ e f (t)dt 1 0 1 2 0 2 =C £[f (t)]+C £[f (t)] 1 1 2 2 性质1表明:函数的线性组合的拉普拉斯变换,等于各函数拉普拉斯变换的线性组合。 例1:求函数cosωt+isinωt的拉普拉斯变换。 α+iβ α 解:由性质1及欧拉公式[e =e (cosβ+isinβ),得 iωt £[cosωt+isinωt]=£[cosωt]+i£[sinωt]=£[e ] 1 s+iω = = s-iω 2 2 s +ω 再由复数相等的定义,得 s £[cosωt]= 2 2 s +ω ω £[sinωt]= 2 2 s +ω 例2:求f (t)=sht,f (t)=cht的拉普拉斯变换。 1 2 解:由性质1,有 t -t 1 t 1 -t £[sht]=£[(e -e )/2]= £[e ]- £[e ] 2 2 1 1 1 1 1 = - = 2 s-1 2 s-1 2 s +1 t -t 1 t 1 -t £[cht]=£[(e -e )/2]= £[e ]+ £[e ] 2 2 1 1 1 1 1 = + = 2 s-1 2 s-1 2 s -1 n 性质2(原函数的微分性质)如果f`(t),f``(t),...,f (t)均满足定理的条件,则 £[f`(t)]=s£[f(t)]-f(0) 一般地,有 (n) n n-1 n-2 (n-1) £[f (t)]=s £[f(t)]-s f(0)-s f(0)-...-f (0) 证:用数学归纳法证明 当n=1时 +∞ -st +∞ -st -st +∞ +∞ -st £[f`(t)]= ∫ e f`(t)dt=∫ e df(t)=e f(t) +s∫ e f`(t)dt 0 0 0 0 =S£[f`(t)]-f(0) 设n=k时,有 (k) k k-1 k-2 (k-1) £[f (t)]=s £[f(t)]-s f(0)-s f(0)-...-f (0) 成立,现看n=k+1时 (k+1) (k) (k) (k) £[f (t)]=£{[f (t)]`}=s£[f (t)]-f (0) k k-1 (k-1) (k) =s{s £[f(t)]-s f(0)-...-f (0)}-f (0) k+1 k (k-1) (k) =s £[f(t)]-s f(0)-...-sf (0)-f (0) 因此,由数学归纳法,得 (n) n n-1 n-2 (n-1) £[f (t)]=s £[f(t)]-s f(0)-s f`(0)-...-f (0) 性质2表明一个函数导数的拉普拉斯变换等于这个函数的拉普拉斯变换乘以参数,再减去函数的初值。 例3已知f(t)=cost,f(0)=1,f`(0)=0,试求它的二阶导函数的拉普拉斯变换。 解:由性质2知 2 £[f``(t)]=s £[cost]-sf(0)-f`(0) 2 s s =s -s*1-0=- 2 2 s +1 s +1 性质3(象函数的微分性质)如果£[f(t)]=F(s),则 d F(s)=-£[tf(t)] ds 一般地有 n d n n F(s)=(-1) £[t f(t)] n ds 事实上 d d +∞ -st +∞ Ә -st F(s)= ∫ e f(t)dt=∫ e f(t)dt ds ds 0 Әs +∞ -st =∫ e f(t)dt=-£[tf(t)] 0 应用数学归纳法,可得 n d n n F(s)=(-1) £[t f(t)] n ds n n 性质3表明一个函数与t 的乘积的拉普拉斯变换等于其象函数的n阶导数与(-1) 的乘积。 n at 例4,求函数f(t)=t e 的拉普拉斯变换 解:由性质3知 n n n at n d n d at £[t e ]=(-1) F(s)=(-1) £[e ] n n ds ds n d 1 =(-1) ( ) n ds s-a n! = n+1 (s-a) 例5,求£[tcosωt]及£[tsinωt] 解:由性质3知 2 2 d d s s +ω £[tcosωt]=- £[cosωt]=- ( )= 2 2 2 2 2 ds ds s +ω (s +ω ) d ω -2ωs 2ωs £[tsinωt]=- ( )=- = 2 2 2 2 2 2 2 2 ds s +ω (s +ω ) (s +ω ) 性质4(积分性质),如果£[f(t)]=F(s),则 t 1 £[∫ f(t)dt]= F(s) 0 s 证:由于 t +∞ t -st £[∫ f(t)dt]= ∫ [∫ f(τ)dτ]e dt 0 0 0 分部积分,得 -st t e t +∞ 1 t -st £[∫ f(t)dt]= ∫ f(τ)dτ + ∫ f(t)e dt 0 s 0 0 s 0 1 = F(s) s 重复应用性质4,得 t t t 1 £[∫ dt ∫ dt ∫ f(t)dt]= F(s) 0 0 0 n s 性质4表明,一个函数积分后再取拉普拉斯变换,等于这个函数的拉普拉斯变换除以参数s 例6求 t f(t)= ∫ tsin2tdt 0 的拉普拉斯变换 解:由例5知 4s £[tsin2t]= 2 2 (s +4) 再由性质4,得 4 £[f(t)]= 2 2 (s +4) 性质5(象函数的积分性质),若£[f(t)]=F(s),则 +∞ f(t) ∫ F(s)ds=£[ ] s t 一般地,有 +∞ +∞ +∞ f(t) ∫ ds ∫ ds ∫ F(s)ds=£[ ] s s s n t n次 证: +∞ +∞ +∞ +∞ +∞ ∫ F(s)ds=∫ [∫ f(t)e dt]ds=∫ f(t)[ ∫ e ds]dt s s 0 0 s +∞ 1 -st +∞ +∞ f(t) -st =∫ (f(t)[- e ] )dt=∫ e dt t s t f(t) =£[ ] t 反复利用上式可得到 +∞ +∞ +∞ f(t) ∫ ds ∫ ds ∫ F(s)ds=£[ ] s s s n t n次 例7,求 sint f(t)= t 的拉普拉斯变换 解:由 1 £[sint]= 2 1+s 及象函数的积分性质,得 sint +∞ 1 £[ ]= ∫ ds=arccots t s 2 1+s 性质6(位移性质),如果F(s)=£[f(t)],则 at £[e f(t)]=F(s-a) 证:由定义,有 at +∞ -st at £[e f(t)]= ∫ e e f(t)dt s +∞ -(s-a)t = ∫ e f(t)dt=F(s-a) s at 性质6表明,一个函数乘以e 的拉普拉斯变换等于其象函数作位移a, at at 例8,求£[e cosωt]及£[e sinωt] 解:因 s £[cosωt]= 2 2 s +ω ω £[cosωt]= 2 2 s +ω 故由性质6,有 at s-a £[e cosωt]= 2 2 (s-a) +ω at s-a £[e cosωt]= 2 2 (s-a) +ω 性质7(延迟性质),设F(s)=£[f(t)],t<0时,f(t)=0,则对任一非负实数τ,有。 -rs £[f(t-τ)]=e F(s) 证:由定义,有 +∞ -st £[f(t-τ)]= ∫ f(t-τ)e dt s τ -st +∞ -st = ∫ f(t-τ)e dt+∫ f(t-τ)e dt 0 τ =I +I 1 2 因为t<τ时,f(t-τ)=0,故 τ -st I =∫ f(t-τ)e dt 1 s 对I 作变换t-τ=u,则 2 +∞ -s(u+τ) -sτ +∞ -su -sτ I = ∫ f(u)e du=e ∫ f(u)e du=e F(s) 2 0 0 因此 -rs £[f(t-τ)]=e F(s) 从几何上看,f(t-τ)的图形是f(t)的图形沿t轴向右平移τ个单位所得。 性质7表明,时间函数延迟τ的拉普拉斯变换,等于它的象函数乘以指数因子e. 例9,求函数 u(t)-u(t-a) f(t)= a 的拉普拉斯变换,其中u(t)为单位阶跃函数,即 0,当t<0时 u(t)={ 1,当t≥0时 解:由拉普拉斯变换简表得知 £[u(t)]=1/s 根据性质7,有 -as £[u(t-a)]=e /s 故 1 £[f(t)]= [£(u(t))]-£[(u(t-a)] a -as 1 1 -as 1 1-e = ( -e )= a s s as 例10,求 sint,当0≤t≤π时 f(t)={ 0,其它 的拉普拉斯变换。 解,设 sint,当t≥0时 f(t)={ 0,当t<π时 因此, f(t)=f (t)+f (t-π) 1 1 于是 £[f(t)]=£[sint]+£[sin(t-π)] 1 1 -πt = + e 2 2 s +1 s +1 1 -πt = (1+e ) 2 s +1 对于性质7,要特别强调t<0时,f(t)=0的这一约定,因此在利用本性质求逆变换时,应为 -rs £ [e F(s)]=f(t-τ)u(t-τ) 例11,求 t -3t f(t)=t∫ e sin2tdt 0 的拉普拉斯变换。 解:因 2 £[sin2t]= 2 s +4 故 -3t 2 £[e sin2t]= 2 (s+3) +4 t -3t 2 £[∫ e sin2tdt]= 0 2 s[(s+3) +4] 2 t -3t d 2 6s +24s+26 £[t∫ e sin2tdt]= ( )= 0 2 2 ds s[(s+3) +4] s[(s+3) +4] 性质8, £[f(t)*g(t)]=£[f(t)]*£[g(t)] 式中 t t f(t)*g(t)= ∫ f(u)g(t-u)du=∫ f(t-u)g(u)du 0 0 称为函数f(t)和g(t)的褶积(或卷积) 性质9, n n d F(s) £[t f(t)]=(-1) n ds 习题3.2 利用拉普拉斯变换的性质,求下列函数的拉普拉斯变换: 2 (1)f(t)=t +6t-3 (2)f(t)=5sin2t-3cos2t t 3t (3)f(t)=1+te (4)f(t)=e sin4t -4t (5)f(t)=e cos(2t+π/4) (6)f(t)=u(2t-1) t 2 2 (7)f(t)=(sint)te (8)f(t)=sin tcos t 2 (9)f(t)=tcosat (10)f(t)=sin t 2 n at (11)f(t)=cos t (12)f(t)=t e (n为正整数) 3.3拉普拉斯逆变换 前面主要介绍由已知函数f(t)求其象函数£[f(t)]=F(s)。但在很多实际应用中常常遇到与此相反的问题,即已知象函数F(s),要求象原函数f(t)。对此类问题常借助于拉普拉斯变换简表以及拉普拉斯变换的性质来解决. 例1,求 -1 k £ [ ] 2 2 s +k -1 s 和£ [ ](k≠0] 2 2 s +k 解:根据拉普拉斯变换表中公式(5),(6),令k=a,得 -1 k £ [ ]=sinkt 2 2 s +k 同理 -1 s £ [ ]=coskt 2 2 s +k 例2,求 2s-5 F(s)= 的象函数 2 s -5s+6 解:因 2s-5 (s-3)+(s-2) 1 1 F(s)= = = - 2 s -5s+6 (s-3)(s-2) s-3 s-2 故由拉普拉斯变换公式(2)及拉普拉斯变换的线性性质,得 -1 -1 1 -1 1 3t 2t f(t)=£ [F(s)]=£ [ ]+£ [ ] =e +e s-3 s-2 例3,求 2 s F(s)= 2 (s+2)(s +s+2) 的拉普拉斯逆变换 解:因 2 2 s (s +s+2)+(s-2) 1 1 F(s)= = = - 2 2 2 (s+2)(s +s+2) (s+2)(s +s+2) s+2 s +s+2 4 1 1 1 4 7 = - = - 2 7 2 s+2 (s+1/2) +7/4 s+2 (s+1/2) +7/4 故由拉普拉斯变换公式(2)及(15),得 4 -1 -1 1 4 -1 7 f(t)=£ [F(s)]=£ [ ]- £ [ ] 7 2 s+2 (s+1/2) +7/4 -2t 4 -t/2 4 =e - e sin t 7 7 当象函数化成部分分式比较复杂时,可以采用待定系数法求解。 例4,求 s+2 F(s)= 的象函数 2 2 s +6s +9s 解:因 s+2 s+2 F(s)= = 2 2 2 s +6s +9s s(s+3) A B C F(s)= + + 2 s (s+3) S+3 其中A,B,C为待定常数。通分后比较等号两边分子,得 2 A(s+3) +Bs+Cs(s+3)=s+2 令s=0,则有9A=2,于是A=2/9, 令s=-3,则有3B=1,于是B=1/3, 令s=1,则有16A+B+4C=3,于是C=-2/9, s+2 2 1 1 1 2 1 F(s)= = + - 2 2 2 s +6s +9s 9 s 3 (s+3) 9 s+3 于是 -1 s+2 f(t)=£ [ ] 2 2 s +6s +9s 2 -1 1 1 -1 1 2 -1 1 = £ [ ]+ £ [ ]- £ [ ] 2 9 s 3 (s+3) 9 s+3 2 1 -3t 2 -3t = + te e 9 3 9 例5,求 2 2 s +2a F(s)= (τ>0) 2 2 (s +a ) 的拉普拉斯逆变换 解:令 2 2 s +2a F(s)= (τ>0) 2 2 (s +a ) 则 2 2 2 s +2a 1 a F(s)= = + 2 2 2 2 2 2 2 (s +a ) s +a (s +a ) 由拉普拉斯变换公式(5),(29)有 -1 1 -1 a 1 -1 2a £ [F (s)]= £ [ ]+ £ [ ] 1 a 2 2 2a 2 2 2 s +a (s +a ) 3 1 = sinat- tcosat 2a 2 再由拉普拉斯的性质6(延迟性质),得 -1 -1 -τs £ [F(s)]=£ [F (s)e ] 1 3 1 =[ sina(t-τ)- (t-τ)cosa(t-τ)]u(t-τ) 2a 2 3 1 sina(t-τ)- (t-τ)cosa(t-τ),当t≥τ时; 2a 2 ={ 0,当t<τ时 例6,求 2 s -1 F(s)=ln 的拉普拉斯逆变换 2 s 解:因 d 2 1 1 2 F(s)= = + - ds 2 s(s -1) s+1 s-1 s 故由拉普拉斯的性质3(象函数的微分性质),得 -1 d -1 1 -1 1 -1 1 £ [ F(s)]=£ [ ]+£ [ ]-2£ [ ] ds s+1 s-1 s t -t -tf(t)=e +e -2 t -t -1 2-e -e f(t)=£ [F(s)]= t 习题3.3 求下列函数的拉普拉斯逆变换: 3s 1 2 (1)F(s)= = + (s-1)(s+2) s-1 s+2 t -2t 解:f(t)=e +2e 1 (2)F(s)= 4 (s+1) s 2 (3)F(s)= =1- s+2 s+2 -2t 解:f(t)=1/ε-2e 2s+3 (4)F(s)= 2 s +9 4s (5)F(s)= 2 s +16 s+3 (6)F(s)= (s+3)(s-3) 1 (7)F(s)= s(s+1)(s-2) s (8)F(s)= 2 s +4s+20 2s-5 (9)F(s)= 2 s 4s-3 (10)F(s)= 2 s +4 s+2 (11)F(s)= 2 2 (s +10)(s +20) 2 s (12)F(s)= 2 (s+2)(s +2s+2) 2s+3 (13)F(s)= 2 s -2s+5 2 s +1 (14)F(s)= 2 s(s-1) 1 -s (15)F(s)= e 2 s(s-1) s+1 (16)F(s)=ln s-1 3.4拉普拉斯变换的应用 拉普拉斯变换在电路分析和自动控制理论中有着非常广泛的应用,下面介绍这方面的几个例子。 例1:设因果信号 at 0,当t<0时 f (t)=e u(t)={ 1 e ,当t>0(a为实数时) 求其拉普拉斯变换 解:由定义,有 -(s-a)t +∞ at -st e +∞ F(s)= ∫ e e dt= 0 -(s-a) 0 1 ,当Re(s)=β>α时 s-a ={ 不定,当β=α时 无界,当β<α时 因此,对于因果信号,仅当Re(s)=β>α时,其拉普拉斯变换存在。 注:本题可直接查拉普拉斯变换表的公式(2)写出 例2,求矩形脉冲信号 1,当0<t<τ时 f(t)=g (t-τ/2)={ τ 0,其余 的象函数。 解法一:因 -τs +∞ -st τ -st 1-e £[f(t)]= ∫ f(t)e dt=∫ e dt= 0 0 s 故矩形脉冲的象函数为 -τs 1-e £[g (t-τ/2)]= (Re(s)>-∞) t s 解法二:因f(t)=1-u(t-τ), 故, 1 -τs £[f(t)]= -e £[u(t)] s 1 -τs = (1-e ) s 例3,求在t=0时接入的周期性单位冲激序列 ∞ ∑ (t-nT)=δ(t)+δ(t-T)+...+δ(t-nT)+... n=0 由拉普拉斯变换公式及性质6(延迟性质),得 £[δ(t)]=1, -Ts £[δ(t-T)]=e , ……. -nTs £[δ(t-nT)]=e , ………… 因此, ∞ -Ts -nTs 1 £[∑ δ(t-nT)]=1+e +...+e +...= (Re(s)>0) n=0 -Ts 1-e 拉普拉斯变换的另一方面的应用是用于解常系数微分方程,下面来讨论这种解法。 对于给定的常系数n阶线性常微分方程: (n) (n-1) (n-2) a y +a y +a y +...+a y`+a y=f(x) 0 1 2 n-1 n { (n-1) y(0)=y ,y`(0)=y ,...,y (0)=y 0 1 n-1 其中a ≠0,a 为常数(i=1,2,...,n), 0 i f(t)是某些常见函数,如多项式、指数函数、正弦、余弦函数以及这些函数的简单复合。 只要f(t)满足拉普拉斯变换存在定理的条件,并假设y(t)满足拉普拉斯变换微分性质中的条件,我们就可以利用原方程的初始条件,根据拉普拉斯变换的微分性质和线性性质,对原方程两端取拉普拉斯变换,把原函数化为象函数的代数方程,从代数方程中解出象函数,然后取它的拉普拉斯逆变换,从而得到原方程的解。 例4,求微分方程 2 y``+a y=t 满足y(0)=b,y`(0)=c的解,其中a,b,c是常数 解:设 £[f(t)]=F(s) 将所给微分方程两边取拉普拉斯变换,得 2 2 1 s Y(s)-sy(0)-y`(0)+a Y(s)= 2 s 其中,Y(s)=£[y(t)],代入初始条件,得 1 2 bs+c s Y(s)= + 2 2 2 2 s +a s +a bs c 1 1 1 = + + ( - ) 2 2 2 2 2 2 2 2 s +a s +a a s s +a 对上式取拉普拉斯逆变换,便得所求微分方程的解为 c 1 1 y(t)=bcosat+ sinat+ (t- sinat) a 2 a a 例5,求微分方程 t y``-2y`+2y=2e cost 满足y(0)=y`(0)=0的解, 解:对给定的微分方程两边取拉普拉斯变换,得 2 2(s-1) s Y(s)-2sY(s)+2Y(s)= 2 (s-1) +1 解出Y(s),得 2(s-1) d 1 Y(s)= =- ( ) 2 2 ds 2 [(s-1) +1] (s-1) +1 由于 -1 1 t £ [ ]=e sint 2 (s-1) +1 故, -1 d 1 t £ [ ( ) ]=e sint 2 ds (s-1) +1 于是, -1 1 £ [Y(s)]=te sint 因此所求微分方程的解为 t y(t)=te sint 由上面两例可以看出,应用拉普拉斯变换求解微分方程时,已经初始条件用上,因此求出的结果就是满足初始条件的特解了,这样就简化了求微分方程时,先求通解再求特解的步骤。 例6,求微分方程组 t x`+x-y=e { t y`+3x-2y=2e 满足初始条件:x(0)=y(0)=1的解 解:设£[x(t)]=X(s),£[y(t)]=Y(s), 对方程组的每个方程两端取拉普拉斯变换,得 1 (s+1)X(s)-Y(s)-1= s-1 { 2 (s-2)Y(s)+3X(s)-1= s-1 解出象函数X(s)和Y(s): 1 X(s)=Y(s)= s-1 再取拉普拉斯逆变换,求得原方程组的解为 t x(t)=y(t)=e 例7,在如图3-1所示的R,C并联电路内,外加电流为单位脉冲函数δ(t)的电流源, 电容C上初始电压为零,求电路中的u(t) 解:设流经电阻R和电容C的电流分别为i (t)和i (t), 1 2 电学原理知 u(t) du i (t)= , i (t)=C 1 R 2 dt 根据电学上的基尔霍夫定律,有 u(t) u C + =δ(t) R R { u(0)=0 这就是该电路电压应满足的微分方程, 设£[u(t)]=U(s),对微分方程两端取拉普拉斯变换,得 U(s) CSU(s)+ =1 R 解得, 1 1 1 U(s)= = 1 C 1 +C S+ R s RC 取拉普拉斯逆变换,有 -1 1 -t/RC u(t)=£ [U(s)]= e C 解的物理意义是:由于在一瞬间电路受单位脉冲电流的作用,把电容的电压由零跃变到1/C,然后电容C向电阻R按指数衰减规律放电。 习题3.4 1.用拉普拉斯变换解下列微分方程: 3t (1)x``(t)-3x`(t)+2x(t)=2e ,x(0)=0,x`(0)=0 (4) (2) (2)y +2y +y=0,y(0)=y`(0)=y``(0)=0,y```(0)=1 (3)x``(t)-x(t)=4sint+5cost,x(0)=-1,x`(0)=-2 t (4)y``-2y`+2y=2e cost,y(0)=y`(0)=0 2.用拉普拉斯变换解下列微分方程组: y`(t)-2x(t)-y(t)=1, (1){ x(0)=2,y(0)=4, x`(t)-x(t)-2y(t)=1, t x`(t)+x(t)+y(t)=e , (2){ t x(0)=y(0)=1, y`(t)+3x(t)-2y(t)=2e , t y``(t)-x``(t)+x`(t)-y(t)=e -2, (3){ x(0)=x`(0)=0,y(0)=y`(0)=0 2y``(t)-x``(t)-2y`(t)+x(t)=-1, 附表5 拉普拉斯变换简表 序号 f(t) F(t) 1 1 1/s at 2 e 1/(s-a) m Γ(m+1) 3 t ,(m>-1) 注:m+1为正整数时,Γ(m+1)=m! m+1 s m at Γ(m+1) 4 t e ,(m>-1) m+1 (s-a) a 5 sinat 2 2 s +a s 6 cosat 2 2 s +a a 7 shat 2 2 s -a s 8 chat 2 2 s -a 2as 9 tsinat 2 2 2 (s +a ) 2 2 s -a 10 tcosat 2 2 2 (s +a ) 2as 11 tshat 2 2 2 (s -a ) 2 2 s -a 12 tchat 2 2 2 (s -a ) m Γ(m+1) m+1 m+1 13 t sinat,(m>-1) *[(s+ia) -(s-ia) ] 2 2 m+1 2i(s +a ) m Γ(m+1) m+1 m+1 14 t cosat,(m>-1) *[(s+ia) +(s-ia) ] 2 2 m+1 2i(s +a ) -bt a 15 e sinat 2 2 (s+b) +a -bt s+b 16 e cosat 2 2 (s+b) +a -bt (s+b)sinc+acosc 17 e sin(at+c) 2 2 (s+b) +a 2 1 1 s 18 sin t ( - ) 2 s 2 s +4 2 1 1 s 19 cos t ( + ) 2 s 2 s +4 2abs 20 sinatsinbt 2 2 2 2 [s +(a+b) ][s +(a-b) ] at bt a-b 21 e -e (s-a)(s-b) at bt (a-b)s 22 ae -be (s-a)(s-b) 2 2 1 1 b -a 23 sinat- sinbt 2 2 2 2 2 2 (s +a )(s +b ) 2 2 (b -a )s 24 cosat-cosbt 2 2 2 2 [s +(a+b) ][s +(a-b) ] 1 1 25 (1-cosat) 2 2 2 a s(s +a ) 1 1 26 (at-sinat) 3 2 2 2 a s (s +a ) 1 1 2 1 27 (cosat-1)+ t 4 2 2 2 2 a 2a s (s +a ) 1 1 2 1 28 (chat-1)- t 4 2 2 2 2 a 2a s (s -a ) 1 1 29 (sinat-atcosat) 3 2 2 2 2a (s -a ) 2 1 s 30 (sinat+atcosat) 2a 2 2 2 (s -a ) 1 1 1 31 (1-cosat)- tsinat 4 3 2 2 2 a 2a s (s -a ) -at s 32 (1-at)e 2 (s+a) -at s 33 (1-at/2)e 3 (s+a) 1 -at 1 34 (1-e ) a s(s+a) -bt -at 1 1 e e 1 35 + ( - ) ab b-a b a s(s+a)(s+b) 式中a,b,c为不相等的常数 -at -bt -ct e e e 1 36 + + (b-a)(c-a) (a-b)(c-b) (a-c)(b-c) (s+a)(s+b)(s+c) 式中a,b,c为不相等的常数 -at -bt -ct ae be ce s 37 + + (c-a)(a-b) (a-b)(b-c) (b-c)(c-a) (s+a)(s+b)(s+c) 式中a,b,c为不相等的常数 2 -at 2 -bt 2 -ct 2 a e b e c e s 38 + + (c-a)(b-a) (a-b)(c-b) (b-c)(a-c) (s+a)(s+b)(s+c) 式中a,b,c为不相等的常数 -at -bt e -e [1-(a-b)t] 1 39 2 2 (a-b) (s+a)(s+b) 式中a,b,c为不相等的常数 -bt -at [a-b(a-b)t]e -ae s 40 2 2 (a-b) (s+a)(s+b) 式中a,b,c为不相等的常数 2 -at at/2 √3at √3at 3a 41 e -e (cos -√3sin ) 2 2 3 3 s +a 3 4a 42 sinatchat-cosatshat 4 4 s +4a 1 s 43 sinatshat 2 4 4 2a s +4a 1 1 44 (shat-sinat) 3 4 4 2a s -a 1 s 45 (chat-cosat) 2 4 4 2a s -a 1 1 46 √s πt 1 1 47 2 π s√s 1 at 1 48 e (1+2at) (s-a) s-a πt 1 bt at 49 (e -e ) 3 s-a - s-b 2 πt 1 1 -a/s 50 cos2 at e √s πt 1 1 a/s 51 ch2 at e √s πt 1 1 -a/s 52 sin2 at e s √s πt 1 1 a/s 53 sh2 at e s √s πt 1 bt at e-a 54 (e -e ) ln t s-b 2 s+a a 55 shat ln =2Arth =1/ln[(1+a/s)/(1-a/s)] t s-a s 注:arth表示反双曲正切, arthx=1/2ln[(1+x)/(1-x)], 2 2 2 s +a 56 (1-cosat) ln t 2 s 2 2 2 s -a 57 (1-chat) ln t 2 s 1 a 58 sinat arctg t s 2 2 1 s +b 59 (chat-cosbt) ln t 2 1 s -a 1 a 60 sin(2a√t) erf( ) πt s 2 2 x -t 注:erf(x)= ∫ e dt,称为误差函数 √π 0 2 1 -2a√t 1 a /s a 61 e e erf( ) πt √s s a 1 -a√s 62 erfc( ) e 2√t s 2 2 t 1 a s 63 erfc( ) e erfc(as) 2a s 1 -2√at 1 a/s a 64 e e erfc( ) √πt √s s 1 1 as 65 e erfc( as ) π(t+a) √s 1 1 66 erfc( at ) a s s+a 1 at 1 67 e erfc( at ) a √s(s-a) 1 68 u(t) s 1 69 tu(t) 2 s m 1 70 t u(t),(m>-1) Γ(m+1) m+1 s 注:m+1为正整数时,Γ(m+1)=m! 71 δ(t) 1 72 δ`(t) s 73 sgnt 2/s 1 74 J (at) 0 2 2 s +a -n 注:I (x)=i J (jx),J 称为第一类n阶贝塞尔(Bessel)函数, n n n I 称为第一类n阶变形的贝塞尔函数,或称为虚宗量的贝塞尔函数。 n 1 75 I (at) 0 2 2 s -a 1 -a/s 76 J (2 at ) e s -bt 1 77 e I (at) 0 2 2 (s+b) -a s 78 tJ (at) 0 2 2 s/2 (s +a ) s 79 tI (at) 0 2 2 s/2 (s -a ) 2 2 1 b(s- s +a ) 80 J (a t(t+2b) e 0 2 2 s +a 第五部分拉普拉斯运算法 下面的电路可参见《数字电路》上册,上海师范大学物理系,编,上海人民出版社1975年出版。 附-2-2拉普拉斯运算法 用频谱变换法求线性电路输出信号,在实际工作中并不是简单可行的,有时往往无法进行计算。如果把福立叶变换推广到拉普拉斯变换,许多复杂问题就变得比较简单可行了。 这种方法通常又称为运算法。下面我们先来讨论福立叶变换推广到拉普拉斯变换的问题。 如从严格的数学分析出发,对一非周期脉冲函数f(t)进行福立叶变换,除了这个函数必须满足狄里赫利条件外,还必须满足绝对可积条件: +∞ ∫ │f(t)│dt<∞ (附-2-14) 0 bt 这个条件对诸如阶跃函数H(t)、斜波函数at(a>0)及指数函数e (b>0)等一类非周期信号就不 -Ct -Ct 能满足。如果这些函数乘上一个指数衰减函数e (C>0),从而构成另一个新函数f`(t)=f(t)e , 使这个新函数f`(t)能满足绝对可积条件,这个指数衰减函数常称为收敛因子。对不同的函数f(t),收敛因子中的实系数C应有不同的数值。 bt (b-C)t 例如:当f(t)=e 时,应取C>b,才能使新函数f`(t)=e 满足绝对可积条件。 因此,收敛因子中的实系数C应是一个参变量。这样,新函数f`(t)的福立叶变换应为: +∞ -Ct -jωt F(C,ω)= ∫ f(t)e e dt (附-2-15) -∞ 或 +∞ -(C+jω)t F(C+jω)= ∫ f(t)e dt -∞ 考虑到非周期脉冲函数f(t)总起始时刻t=0的实际情况,令t<0时,f(t)=0,于是上式的积分范围可改写为0→∞。另外,令复变量p=C+jω,(附-2-15)式即可改写为: +∞ -pt F(p)= ∫ f(t)e dt (附-2-16) -∞ 上式即为所谓拉普拉斯变换式。其中,f(t)常称为原函数;F(p)称为象函数,它是复变量p的函数。如果把复变量p当作复频率看待,F(p)也可称为复频谱密度函数。 (附-2-16)式常用符号式F(p)=£[f(t)]表示。由福立叶反变换式,不难推出拉普拉斯反变换 -Ct 式。新函数f(t)e 的福立叶反变换应为: -Ct 1 +∞ jωt f(t)e = ∫ F(C+jω)e dω 2π -∞ Ct 上式两边乘以因子e ,并且 1 dω= d(C+jω) j 于是得: 1 C+jω (C+jω)t f(t)= ∫ F(C+jω)e d(C+jω) 2πj C-jω 或 1 C+jω pt f(t)= ∫ F(p)e d(p) (附-2-17) 2πj C-jω 这个积分是沿着复平面的右半平面上一垂直线进行的线积分,这条作为积分路径的垂直线应选取在被积函数的所有极点的右边。可以证明,作为拉普拉斯反变换的(附-2-17)式可改写成为: 1 pt f(t)= ∮F(p)e d(p) (附-2-18) 2πj 上式积分号上的小圆圈表示积分式围线积分,这个围线必须包含被积函数的所有极点在内。 -1 (附-2-18)式常用符号式f(t)=£ [F(p)]表示。必须指出,在实际应用中,拉普拉斯变换与反变换一般都不进行具体的验算。而可以通过f(t)与F(p)的变换对照表(附录3)求得。 大家知道,线性电路的输入和输出信号的时间函数v (t)和u (t)之间的关系, i 0 一般可用n阶线性微分方程表示: n n-1 d v d v dv 0 0 0 a +a +...+a +a v =v (t) (附-2-19) n n-1 n-1 1 0 0 i dt dt dt 式中,系数a ,a ,...,a ,a 与网络的元件与结构有关。 n n-1 1 0 图附-2-1(a)所示的简单阻容电路,微分方程为: n dv 0 RC +v =v (t) (附-2-20) 2πj 0 i 设£[v (t)]=V (p);£[v (t)]=V (p) i i 0 0 dv (1) 0 一阶导数v (t)= 的变换可求得为: 0 dt dv (t) ∞ -pt 0 £[v (t)]= ∫ e dtpV (p)-v (0) (附-2-21) 0 0 dt 0 0 式中,v (0)决定于线性电路输出的初始状态。如果在输入信号作用前, 0 线性电路中的储能元件C或L未充磁,则v (0)=0, 0 2 d v (2) 0 二阶导数v (t)= 的变换为: 0 2 dt (2) (1) (1) 2 (1) £[v (t)]=p£[v (t)]-v (0) =p V (p)-pv (0)-v (0) (附-2-22) 0 0 0 n阶导数的变换即为: (n) n n-1 n-2 (1) (n-1) £[v (t)]=p V (p)-p v (0) -p v (0)-...-v (0) (附-2-23) 0 0 0 0 0 式中, n d v (n) 0 v (0)= lim t→0 n dt 由于拉氏变换是一种线性运算,(附2-19)式左边的变换等于各项分别进行变换后之和,因而得 V (p)=K(p)[V (p)+Q(p)] (附-2-24) 0 i 式中, 1 K(p)= (附-2-25) n n-1 a p +a p +...+a p+a n n-1 1 0 n-1 n-2 n-2 n-3 (1) Q(p)=(a p +a p +...+a p+a )v (0) +(a p +a p +...+a p +a )v (0)+... n n-1 2 1 0 n n-2 2 1 0 (n-2) (n-1) +(a p+a )v (0)+a v (0) (附-2-26) n n-1 0 n 0 K(p)是线性网络的传输函数,它完全取决于线性网络的元件和结构形式;Q(p)是与网络 (i) 的初始状态密切有关的函数,如网络的初始状态为零,v (0)和有关的v (0)均等于零, 0 0 从而Q(p)=0,这样,(附-3-24式)则为: V (p)=K(p)V (p) (附-2-27) 0 i 此式说明,在线性网络的初始状态为零的情况下,网络的输出信号v (t)的象函数等于输入信号的象函数与网络的传输函数的乘积。这里的传输函数K(p)与频谱法中的K(ω)具有同一形式,只要将K(ω)中的jω用p代替。如果图附-2-1(a)所示的阻容电路的初始状态为: C充有负极性的电压——θ,即v (0)=-θ, 0 如图附-2-4所示,电路的输入电压为阶跃电压E,即 v (t)=E (t>0), i 这样 dv 0 RC +v =E dt 0 上式的拉普拉斯变换为: 1 E V (p)= ( -θCR) (附-2-28) 0 dt P 查表(附录3)求得原函数v (t)为: 0 -t/RC -t/RC v (t)=E(1-e )-θe =v` (t)+v`` (t) (附-2-29) 0 0 0 式中,v` (t)是由输入电压v (t)=E所产生的,常称为强迫相应; 0 I v`` (t)是由电路的初始值v (0)=-θ所产生的,常称为自由响应。 0 0 波形如图附-2-5所示,如θ=0,则v`` (t)=0,得 0 -t/RC v (t)=E(1-e ) 0 此式与用富立叶变换法所得结果完全相同。由(附-2-21)式可以推导出两个很有用的定理。当p→∞时,如果p的实部为正数,(附-2-21)式左边的积分等于零,因此,得 v (0)= lim v (t)= lim pV (p) (附-2-30) 0 t→∞ 0 t→∞ 0 这是所谓始值定理。当p→0时,(附-2-21)式左边的积分为: dv ∞ dv lim £[ ]= ∫ dt=v(∞)-v(0) p→0 dt 0 dt 因而得: v (0)= lim v (t)= lim pV (p) 0 t→∞ 0 t→∞ 0 这是所谓终值定理。 始值定理和终值定理可以帮助我们直接从网络的输出象函数V (p)求的输出原函数v (t)的 0 0 瞬态值v (0)与稳态值v (∞),而不必等到v (t)求得后再进行推导,例如,如图附-2-1(a) 0 0 0 所示的RC电路,由(附-2-28)式得: E θCR v (0)= lim pV (p)= lim ( - )=-θ 0 p→∞ 0 p→∞ 1+pCR 1 +CR p E θCR v (∞)= lim pV (p)= lim ( - )=E 0 p→0 0 p→0 1+pCR 1 +CR p 可见,由始值、终值定理求得的结果与前面结果完全一致,但要方便很多。在实际工作中,往往只需要知道始值和终值。例如:当矩形脉冲v (t)输入RC微分电路, i 其输出是正负相间的尖脉冲,如图附-2-6所示,这点大家都是很熟悉,就没有必要去取输出函数v (t),但是尖脉冲的起始值和稳定值是大家需要知道的。设电路的初始状态为零, 0 由(附-2-7)式可得: V (p)=K(p)V (p) 0 i 式中, p K(p)= 1 p+ p V (p)可查附录3得: i -pτ E u V (p)= (1-e ) p 因此,得 -pτ 1 1 u V (p)=E( - e )=V (p)-V (p) 1 1 01 02 p+ p+ CR CR 由始值、终值定理很快就求得: -pτ E u v (0)= lim pV (p)= lim [ (1-e )]=E 0 p→0 0 p→0 1+ 1 CR v (∞)= lim pV (p)=0 0 p→0 0 附-2-3卷积积分和重迭积分法 1.卷积积分和重迭积分原理 如图-2-7所示,图(a)为跳变沿是ε的阶跃函数;(b)为跳变沿ε→0的理想阶跃函数;(c)为阶跃函数发生于t=τ,表示式分别为: 0,t<0 H (t)={ t/ε,0≤t≤ε ε 1,t>ε 0,t<0 H(t)= lim H (t)={ 附-2-32 ε→0 ε 1,t>0 0,t<τ H(t-τ)={ 1,t>τ 一个矩形脉冲可以看作两个阶跃函数,H(t)和-RH(t-τ)迭加而成,如图附-2-8(a)所示;而图附2-8(b)所示说明一个阶梯脉冲可看作四个阶跃函数迭加。同理,任何一个脉冲信号都可看成由若干个发生在不同时刻的阶跃函数迭加而成。图附-2-9所示为单位脉冲波形,图(a)的脉宽为ε何幅度为1/ε;图(b)为脉宽ε→0的理想单位脉冲函数;图(c)为发生于时间t=τ的单位脉冲函数,表示式分别为: 0,t>ε,t<0 δ (t)={ ε t/ε,0≤t≤ε 0,t≠0 δ(t)= lim δ (t)={ 附-2-33 ε→0 ε ∞,t=0 0,t≠τ δ(t-τ)={ ∞,t=τ 单位脉冲函数还具有下列特性: ∞ ε ∫ δ(t)dt=∫ δ(t)dt=1 -∞ 0 ∞ τ+ε ∫ δ(t-τ)dt=∫ δ(t-τ)dt=1 附-2-34 -∞ 0 δ(t)= lim H` (t)=H`(t) 附-2-35 ε→0 ε ∞ ∫ f(t)δ(t-τ)dt=f(τ) 附-2-36 -∞ (附-2-34)式表明单位脉冲函数曲线与时间轴t围城的面积等于1,因而成为单位脉冲。(附-2-35)式表明阶跃函数的导数是单位脉冲函数。(附-2-36)式的物理意义说明单位脉冲起着扫描取样脉冲的作用,也就是说,如果在τ时刻出现单位脉冲δ(t-τ),函数f(t)的值f(τ)便被取出。从这一特性出发可以把任何时间函数f(t)都可看作若干个发生在不同时刻的δ (t-τ)和一定常数的乘积组合而成。 ε 设线性网络的传输函数为N(p),网络对阶跃函数H(t)的输出响应为h(t);对脉冲函数δ(t)的输出响应为k(t),如图附-2-10所示。再令: H(p)=£[H(t)];h(p)=£[h(t)];δ(p)=£[δ(t)];k(p)=δ[k(t)] 根据(附2-27)式可得: h(p)=£[h(t)]=N(p)H(p); (附-2-37) k(p)=£[k(t)]=N(p)δ(p); (附-2-38) H(p)和δ(p)可查表得到,或者直接用(附-2-16)式计算出: ∞ -pt H (p)=£[H(t-τ)]= ∫ H(t-τ)e dt τ 0 τ -pt 1 -pt = ∫ e dt= e (附-2-39) 0 p 当τ=0时,H(p)=1/p,同样,由(附-2-16)式,并考虑到(附-2-30)式的关系,可得 ∞ -pt -pt δ(p)=£[δ(t-τ)]= ∫ δ(t-τ)e dt=e (附-2-40) 0 -pt 当τ=0时,δ(p)=1,上列两式中的e 为位移因子,它反应原函数在时间上向后位移了τ值。将(附-2-39),(附-2-40)式分别代入(附-2-37),(附-2-38)式,可得: 1 h(p)=£[h(t)]= N(p) (附-2-41) p k(p)=£[k(t)]=N(p) (附-2-42) 可见,阶跃响应函数h(t)和脉冲响应函数k(t)完全决定于线性网络的传输函数N(p). 这样,输入给网络的信号为一任意函数x(t), 其输出函数y(t)可由若干个对应的阶跃响应函数或脉冲响应函数线性迭加而成,这就是下面要导出重迭积分和卷积积分的物理基础。 2.重迭积分 设有某一线性网络N(p),输入函数为x(t),具有任意波形,但当t<0时,x(t)=0,对应的输出函数为y(t). 我们将输入函数x(t)用阶跃函数序列之和来近似,如图附2-11所示,因此得 x(t)=x(0)H(t)+△x H(t-△τ)+△x H(t-2△τ)+..+△x H(t-n△τ)+... 1 2 n 式中:△x =x(i△τ)-x[(i-1) △τ] i 根据线性迭加原理,线性网络的输出函数y(t)可近似为: y(t)=x(0)h(t)+△x h(t-△τ)+...+△x h(t-n△τ)+.. 1 n △x n i =x(0)h(t)+ ∑ h(t-i△τ)△τ i=1 △τ 令:△x→0;n→∞;n△τ→τ,则得 △x n i y(t)=x(0)h(t)+ lim ∑ h(t-i△τ)△τ i=1 △τ 即, t y(t)=x(0)h(t)+ ∫ x`(t)h(t-τ)dτ 0 或 (附-2-43) t y(t)=x(0)h(t)+ ∫ x`(t-τ)h(t)dτ 0 上两式就是所求的重迭积分式。重迭积分式有多种形式,除(附-2-45)式外,下列两种也是常用的重迭积分式: t y(t)=h(0)x(t)+∫ x(τ)h`(t-τ)dτ 0 (附-2-44) t y(t)=h(0)x(t)+∫ x(t-τ)h`(t)dτ 0 下面仍举矩形脉冲通过RC微分电路之例,说明重迭积分的应用。 用x(t)=v (t),y(t)=v (t)代入(附-2-43)式,则为 i 0 t v (t)=v (0)h(t)+∫ v` (τ)h(t-τ)dτ (附-2-45) 0 0 0 i 由图附-2-8(a)可知: v (t)=EH(t)-EH(t-τ ) I 0 再根据(附-2-32)式,得: v (0)=0 i 由(附-2-35)式可得: v` (t)=E*H`(t)-E*H`(t-τ )=E[δ(t)-δ(t-τ )] i u u RC电路的传输函数的象函数应为: R p N(p)= = 1 1 R+ p+ pC CR 因此 1 £[h(t)]= 1 p+ CR 查表可得: -1/RC e ,t>0 h(t)={ 0,t<0 或 -1/RC h(t)=e *H(t) 将v (0),v` (t),h(t)代入(附-2-45)式得: I I T -(t-τ)/RC v (t)= ∫ E[δ(τ)-δ(τ-τ )]H(t-τ)e dτ 0 0 u t -(t-τ)/RC t -(t-τ)/RC =E[∫ δ(τ)H(t-τ)e dτ- ∫ δ(τ-τ )H(t-τ)e dτ] 0 0 u -1/RC -(t-τ)/RC =E[e H(t)-e H(t-τ )] (附-2-46) u 上式中的第一项为发生在t=0处的正尖脉冲;第二项则为发生在t=τ 处的负尖脉冲,如图附-2-6所示。 3.卷积积分 与重迭积分的导出方法相似,y(t)可近似地表示为: y(t)=x △τk(t)+x △τk(t-△τ)+...+x △τk(t-n△τ)+... 0 1 n n =∑ x △τk(t-i△τ) i=0 i 令:△τ→0,n→∞;n△τ=τ,得: n y(t)= lim ∑ x △τk(t-i△τ) △x→0 i=0 i n→∞ 即, t y(t)= ∫ x(τ)k(t-τ)dτ (附-2-47) 0 或 t y(t)= ∫ x(t-τ)k(τ)dτ (附-2-48) 0 上面两式就是卷积积分式,它是求网络的输出函数的一种方法。其实,卷积积分式可以很方便的从拉普拉斯变换法求得。设网络的初始状态为零,根据(附-2-27)式可得: Y(p)=N(p)X(p) (附-2-49) 因此: -1 -1 y(t)=£ [Y(p)]=£ [N(p)X(p)] 1 C+j∞ pt = ∫ N(p)X(p)e dp 2πj C-j∞ 由(附-2-12)式可得: ∞ -pτ N(p)=£[k(τ)]= ∫ k(τ)e dτ 0 将N(p)代入(附-2-50)式,则得: 1 C+j∞ pt ∞ -pτ y(t)= ∫ X(p)e [∫ k(τ)e dτ]dp 2πj C-j∞ 0 ∞ 1 C+j∞ p(t-τ) = ∫ k(τ)[ ∫ X(p)e dp]dτ 0 2πj C-j∞ ∞ = ∫ k(τ)x(t-τ)dτ 0 t = ∫ k(τ)x(t-τ)dτ (附-2-48) 0 可见,卷积积分式和运算法的基本公式(附-2-19)式是相对应的。前者是在时域中表出了网络与输入、输出函数之间关系;而后者则在复频域中表现了这三者的关系。 再应用卷积积分来求解图附-2-6所示的例题,由(附-2-47)式得: t v (t)=∫ v (τ)k(t-τ)dτ (附-2-51) 0 0 i 式中, v (τ)=E*H(τ)-E*H(τ-τ ) i u 由(附-2-42)式得 1 p RC £[k(t)]=N(p)= =1- 1 1 p+ p+ RC CR 查表得: 1 -t/RC k(t)=δ(t)- e H(t) RC 将k(t)和v (τ)的结果代入(附-2-51)式: i t 1 -(t-τ)/RC v (t)= ∫ E[H(τ)-H(τ-τ )][δ(t-τ)- e H(t-τ)dτ 0 u RC t t 1 -(t-τ)/RC = ∫ E[H(τ)-H(τ-τ )]δ(t-τ)dτ-∫E[H(τ)-H(τ-τ )] e H(t-τ)dτ 0 u 0 u RC -(t-τ )/RC -t/RC u =v (t)-{v (t)-E[e H(t)-e H(t-τ )]} i i u -(t-τ )/RC -t/RC u =E[e H(t)-e H(t-τ )] (附-2-52) u 从这个简单的例子可直接清楚的看到,应用卷积积分式求解网络远没有运算法简便有效。因此它不是求解线性网络的一种常用的方法,但它常常作为运算法中求解拉普拉斯反变换等的一种辅助方法,例如,求解下列象函数 1 F(p)= (p+α)(p+β) 的拉普拉斯反变换,可将F(p)写成: F(p)=F (p)*F (p) 1 2 式中: 1 F (p)= 1 p+α 1 F (p)= 2 p+β 因此 -1 -αt f (t)=£ [F (p)]=e 1 1 -α(t-τ) f (t-τ)=e 1 -1 -βt f (t)=£ [F (p)]=e 2 2 应用卷积积分式则得: -1 t f(t)=£ [F(p)]= ∫ f (t-τ)f (τ)dτ 0 1 2 t -α(t-τ) -βt -αt t -(β-α)τ =∫ e e dτ=e ∫ e dτ 0 0 1 -βt -αt = (e -e ) α-β 3.脉冲信号的拉普拉斯变换表 序号 F(p) f(t)(脉冲波形) 1 1/p -ap 2 e /p 3 1 -ap 4 e -ap 5 pe 2 6 1/p -ap 3 7 e /p -ap 1-e 8 p -ap -bp e -e 9 p -ap 2 (1-e ) 10 p -ap 2 (1-e ) 11 2 p -ap -bp 2 (e -e ) 12 2 p -ap be 13 p(p+b) 1 14 p+b -ap e 15 p+b 1 16 (p+a)(p+b) -ap -bp e -e 17 2 p -ap 1-e 18 2 p -ap e +ap-1 19 2 ap -2cπp/a a(1-e ) 20 2 2 p +a -2cπp/a p(1-e ) 21 2 2 p +a -πp/a a 1+e 22 2 2 -πp/a p +a 1-e a 1 23 2 2 -πp/a p +a 1-e 1 24 -ap p(1+e ) -ap e +ap-1 25 -ap ap (1+e ) 1 26 ap p(1+e ) -ap/v 1-e 27 -ap p (1+e ) -ap 1-e 28 -ap p (1+e ) -ap e -1 29 -ap p (1+e ) -ap 4-e 30 -ap p(4+2e ) -ap 1-e 31 ap -ap p(e +e ) -ap 1-e 32 2 -ap ap (1+e ) -ap 2 (1-e ) 33 2 -4ap ap (1+e ) -ap/2v 2 2v(1-e ) 34 2 -ap ap (1+e ) -ap 2 ap/v v(v-1)+ve -v e 35 2 -ap (v-1)ap (1-e ) a ap+1-e p 36 2 -ap ap (1-e ) -ap 2-ap-(2+ap)e 37 2 -ap ap (1-e ) -ap 2(1-e ) 1 38 - 2 -ap ap (1-e ) p 1 39 p p(e -1) 1 40 ap p(e -1) 1 41 -p p(1-e ) 1 42 -ap p(1-e ) 1 43 2 2 p +a p 44 2 2 p +a 1 45 2 2 p +a p 46 2 2 p +a 1 47 2 2 p +a 1 48 2 2 2 2 p +a (p+ p +a ) 1 49 2 2 (p+a) +b -1 a 50 tg ( ) p 第六部分拉普拉斯变换AD采样电路 根据拉普拉斯变换 +∞ -st F(s)=£[f(t)]= ∫ e f(t)dt, 0 |
| 发帖须知: 1,发帖请遵守《计算机信息网络国际联网安全保护管理办法》、《互联网信息服务管理办法》、 《互联网电子公告服务管理规定》、《维护互联网安全的决定》等法律法规。 2,请对您的言论负责,我们将保留您的上网记录和发帖信息。 3,在此发帖表示认同我们的条款,我们有权利对您的言论进行审核、删除或者采取其他在法律、地方法规等条款规定之内的管理操作。 |