为什么u(t)与1的拉普拉兹变换都是1/s
如何解释u(t)与1的拉普拉兹变换都是1/s一、问题发现通过matlab求解拉普拉兹变换我们可以发现,无论是单位阶跃函数1还是常数1它的拉普拉兹变换都是1/s,而1/s的拉普拉兹逆变换却是常数1,并不是单位阶跃函数。那么问题来了时域的不同函数,为什么映射到s域是同一个函数。通过拉普拉兹正变换的公式,我们很容易发现,常数1的s域变换是不存在的,因为它在实数域上不收敛。于是我们推测,matlab将常数
如何解释u(t)与1的拉普拉兹变换都是1/s
一、问题发现
通过matlab求解拉普拉兹变换我们可以发现,无论是单位阶跃函数1还是常数1它的拉普拉兹变换都是1/s,而1/s的拉普拉兹逆变换却是常数1,并不是单位阶跃函数。那么问题来了时域的不同函数,为什么映射到s域是同一个函数。
通过拉普拉兹正变换的公式,我们很容易发现,常数1的s域变换是不存在的,因为它在实数域上不收敛。于是我们推测,matlab将常数1与u(t)看作了同一个函数,也就是得出的结果1就是u(t)为了证实这一猜想,我们使用拉普拉兹逆变换来寻求答案。
二、推导过程
- 知识准备
-
通过查阅资料,我们可以获得拉普拉兹逆变换的公式如下:
L T − 1 ( 1 s ) = 1 2 π i ∫ σ − R j σ + R j e s t s d s 其 中 R → + ∞ L{T^{ - 1}}(\frac{1}{s}) = \frac{1}{{2\pi i}}\int_{\sigma - Rj}^{\sigma + Rj} {\frac{{{e^{st}}}}{s}ds} 其中R→+∞ LT−1(s1)=2πi1∫σ−Rjσ+Rjsestds其中R→+∞ -
为了求出该积分的结果,我们先令:
∫ σ − R j σ + R j e s t s d s = g ( t ) \int_{\sigma - Rj}^{\sigma + Rj} {\frac{{{e^{st}}}}{s}ds = g(t)} ∫σ−Rjσ+Rjsestds=g(t) -
我们发现被积函数具有奇点,我们通过复变函数的知识构建一个闭合回路,包裹奇点,复数域积分回路如下图,我们令r→∞且R→∞即可求得积分值g(t):
-
这时候沿着A→B→C→D→A的积分,就等效于被积函数零点处的“留数”,故存在等式:
g ( t ) ∣ A B + ∫ σ + R j − r + R j e s t s d s ∣ B C + ∫ − r − R j σ − R j e s t s d s ∣ D A + ∫ − r + R j − r − R j e s t s d s ∣ C D = 2 π i R e s ( f ( t ) , 0 ) g(t)\left| {AB} \right. + \int_{\sigma + Rj}^{ - r + Rj} {\frac{{{e^{st}}}}{s}ds} \left| {BC} \right. + \int_{ - r - Rj}^{\sigma - Rj} {\frac{{{e^{st}}}}{s}ds} \left| {DA} \right. + \int_{ - r + Rj}^{ - r - Rj} {\frac{{{e^{st}}}}{s}ds\left| {CD = 2\pi i{\rm{Res}}(f(t),0)} \right.} g(t)∣AB+∫σ+Rj−r+Rjsestds∣BC+∫−r−Rjσ−Rjsestds∣DA+∫−r+Rj−r−Rjsestds∣CD=2πiRes(f(t),0) -
对中间两项BC与DA积分进行变量替换可以得到:
∫ − r σ − e ( x + R j ) t x + R j d x ∣ x = s − R j + ∫ − r σ e ( x − R j ) t x − R j d x ∣ x = s + R j {\left. {\int_{ - r}^\sigma { - \frac{{{e^{(x + Rj)t}}}}{{x + Rj}}dx} } \right|_{x = s - Rj}} + {\left. {\int_{ - r}^\sigma {\frac{{{e^{(x - Rj)t}}}}{{x - Rj}}dx} } \right|_{x = s + Rj}} ∫−rσ−x+Rje(x+Rj)tdx∣∣∣∣x=s−Rj+∫−rσx−Rje(x−Rj)tdx∣∣∣∣x=s+Rj -
合并后我们可以发现:
∫ − r σ ( x + R j ) e ( x − R j ) t − ( x − R j ) e ( x + R j ) t x 2 + R 2 d x \int_{ - r}^\sigma {\frac{{(x + Rj){e^{(x - Rj)t}} - (x - Rj){e^{(x + Rj)t}}}}{{{x^2} + {R^2}}}dx} ∫−rσx2+R2(x+Rj)e(x−Rj)t−(x−Rj)e(x+Rj)tdx -
当R→∞时,被积函数趋向于0:
lim R → ∞ ( x + R j ) e ( x − R j ) t − ( x − R j ) e ( x + R j ) t x 2 + R 2 = 0 \mathop {\lim }\limits_{R \to \infty } \frac{{(x + Rj){e^{(x - Rj)t}} - (x - Rj){e^{(x + Rj)t}}}}{{{x^2} + {R^2}}} = 0 R→∞limx2+R2(x+Rj)e(x−Rj)t−(x−Rj)e(x+Rj)t=0 -
所以原积分值为0,步骤4中公式变为:
g ( t ) − ∫ − r + R j − r − R j e s t s d s = 2 π i R e s ( f ( t ) , 0 ) g(t) - \int_{ - r + Rj}^{ - r - Rj} {\frac{{{e^{st}}}}{s}ds = 2\pi i{\rm{Res}}(f(t),0)} g(t)−∫−r+Rj−r−Rjsestds=2πiRes(f(t),0) -
由于s=σ+jw由积分上下限可以知道,σ为常数-r,因此CD积分可化简为:
∫ − r + R j − r − R j e s t s d s = − ∫ − R R e − r t ( cos ( w t ) + j sin ( w t ) ) j − r + w t d w \int_{ - r + Rj}^{ - r - Rj} {\frac{{{e^{st}}}}{s}ds = - \int_{ - R}^R {\frac{{{e^{ - rt}}(\cos (wt) + j\sin (wt))j}}{{ - r + wt}}dw} } ∫−r+Rj−r−Rjsestds=−∫−RR−r+wte−rt(cos(wt)+jsin(wt))jdw -
当r→∞时,被积函数也为0,所以CD积分为0:
lim r → ∞ e − r t ( cos ( w t ) + j sin ( w t ) ) j − r + w t = 0 \mathop {\lim }\limits_{r \to \infty } \frac{{{e^{ - rt}}(\cos (wt) + j\sin (wt))j}}{{ - r + wt}} = 0 r→∞lim−r+wte−rt(cos(wt)+jsin(wt))j=0 -
这样我们就得到了,积分g(t)的值:
g ( t ) = 2 π i lim s → 0 s f ( s ) = 2 π i g(t) = 2\pi i\mathop {\lim }\limits_{s \to 0} sf(s) = 2\pi i g(t)=2πis→0limsf(s)=2πi -
带入第一个式子就得到了:
L T − 1 ( 1 s ) = 1 L{T^{ - 1}}(\frac{1}{s}){\rm{ = 1}} LT−1(s1)=1 -
确实得到了常数1,但是常数1的拉普拉兹双边正变换确实是不存在的,这是为什么呢。
三、问题解释
-
这个问题主要来自于,当我们进行上述推导时,已经默认了t≥0了,何以见得?
-
在第10步推导中,"当r→∞时,被积函数也为0"这句话成立的前提,就是t不能为负数,一旦t<0,通过洛必达法则我们就可以推知,被积函数是趋于负无穷的,g(t)的积分也就趋于正无穷,积分是不存在的,因此g(t)积分存在的前提就是t≥0这是在课本中没有提及的内容。
-
所以结论就是,这里的常数1指代的就是u(t)
开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!
更多推荐
所有评论(0)