跳到主要内容

東京工業大学 工学院 情報通信系 2017年8月実施 S1 Laplace変換と微分方程式

Author

思齐塾, 祭音Myyura (co-authored with GPT 5.6 SOL)

Description

S1.

以下の問に答えよ. ただし, tt の関数 f(t)f(t) のラプラス変換 F(s)F(s)

F(s)=0f(t)estdtF(s) = \int_{0}^{\infty} f(t)e^{-st} dt

と定義され, ss は複素変数である. また, 以下のすべての tt の関数は t<0t < 0 においては零であるとする. なお, ラプラス変換に関する諸性質は導出過程を示さずに用いてよい.


  1. 次の 3 種類の関数のラプラス変換を求めよ. ただし, aa は正の定数である.

    a) f1(t)=eatf_1(t) = e^{-at} b) f2(t)=costf_2(t) = \cos t c) f3(t)=eatsintf_3(t) = e^{-at} \sin t

  2. x(t)=sintx(t) = \sin t のとき, 以下の式を満足する関数 y(t)y(t) をラプラス変換を用いて求めよ. ただし, y(0)=0y(0) = 0 とする.

dy(t)dt+2y(t)=x(t)\frac{dy(t)}{dt} + 2y(t) = x(t)
    1. の関数 y(t)y(t) は関数 x(t)(=sint)x(t) (= \sin t) と関数 h(t)h(t) によって以下のように表すことができる. 関数 h(t)h(t) を求めよ.
y(t)=h(tτ)x(τ)dτy(t) = \int_{-\infty}^{\infty} h(t - \tau)x(\tau)d\tau
  1. x(t)=sintx(t) = \sin t のとき, 以下の式を満足する関数 z(t)z(t)tt を十分大きくすると, z(0)z(0)dz(t)dtt=0\left.\frac{dz(t)}{dt}\right|_{t=0} に無関係にある関数に近づく. tt を十分大きくしたときの z(t)z(t) を求めよ.
d2z(t)dt2+2dz(t)dt+z(t)=x(t)\frac{d^2z(t)}{dt^2} + \sqrt{2}\frac{dz(t)}{dt} + z(t) = x(t)

题目描述

函数 f(t)f(t) 的 Laplace 变换定义为

F(s)=0f(t)estdt,F(s)=\int_0^\infty f(t)e^{-st}\,dt,

其中 ss 是复变量。以下出现的所有时间函数在 t<0t<0 时均为零;可以直接使用 Laplace 变换的各项性质,无须推导。

  1. a>0a>0,分别求下列三个函数的 Laplace 变换:

    1. f1(t)=eatf_1(t)=e^{-at}
    2. f2(t)=costf_2(t)=\cos t
    3. f3(t)=eatsintf_3(t)=e^{-at}\sin t
  2. x(t)=sintx(t)=\sin t 时,用 Laplace 变换求满足

dy(t)dt+2y(t)=x(t),y(0)=0\frac{dy(t)}{dt}+2y(t)=x(t), \qquad y(0)=0

的函数 y(t)y(t)。 3. 第 2 问的 y(t)y(t) 可写成 x(t)x(t) 与某个函数 h(t)h(t) 的卷积:

y(t)=h(tτ)x(τ)dτ.y(t)=\int_{-\infty}^{\infty} h(t-\tau)x(\tau)\,d\tau.

h(t)h(t)。 4. 仍令 x(t)=sintx(t)=\sin t。满足

d2z(t)dt2+2dz(t)dt+z(t)=x(t)\frac{d^2z(t)}{dt^2} +\sqrt2\frac{dz(t)}{dt} +z(t)=x(t)

的解在 tt 充分大时会趋近于一个与 z(0)z(0)dz(t)/dtt=0\left.dz(t)/dt\right|_{t=0} 无关的函数。求这一大时间下的 z(t)z(t)

Kai

1)

F1(s)=1s+a,F2(s)=ss2+1,F3(s)=1(s+a)2+1.F_1(s)=\frac1{s+a},\quad F_2(s)=\frac{s}{s^2+1},\quad F_3(s)=\frac1{(s+a)^2+1}.

それぞれの収束域は s>a\Re s>-a, s>0\Re s>0, s>a\Re s>-a

2)

y(0)=0y(0)=0 より (s+2)Y(s)=1/(s2+1)(s+2)Y(s)=1/(s^2+1)。部分分数分解すると

Y(s)=15(1s+2ss2+1+2s2+1).Y(s)=\frac15\left(\frac1{s+2}-\frac{s}{s^2+1}+\frac2{s^2+1}\right).

よって

y(t)=15(e2tcost+2sint)(t0).\boxed{y(t)=\frac15(e^{-2t}-\cos t+2\sin t)\quad(t\ge0)}.

3)

Y(s)=H(s)X(s)Y(s)=H(s)X(s) より H(s)=1/(s+2)H(s)=1/(s+2)。したがって

h(t)=e2t(t0),h(t)=0(t<0).\boxed{h(t)=e^{-2t}\quad(t\ge0)},\qquad h(t)=0\quad(t<0).

4)

同次解の特性根は 1/2±i/2-1/\sqrt2\pm i/\sqrt2 であり,同次解は減衰する。定常解を zp=Asint+Bcostz_p=A\sin t+B\cos t とおくと

zp+2zp+zp=2(AcostBsint)=sint.z_p''+\sqrt2z_p'+z_p=\sqrt2(A\cos t-B\sin t)=\sin t.

したがって A=0A=0, B=1/2B=-1/\sqrt2 であり,

z(t)12cost(t).\boxed{z(t)\sim-\frac1{\sqrt2}\cos t\quad(t\to\infty)}.