「制御と振動の数学/第一類/Laplace 変換による解の吟味/t0≠0で初期値が与えられている場合」の版間の差分

提供: testwiki
ナビゲーションに移動 検索に移動
imported>MathXplore
 
(相違点なし)

2022年11月23日 (水) 15:26時点における最新版

Laplace 変換を用いて微分方程式を解くとき,初期値は t=0 で与えられている必要があった. しかし,前節で述べた定常性の原理を用いると,t=t00 で初期値が与えられている場合,すなわち, テンプレート:制御と振動の数学/equation テンプレート:制御と振動の数学/equation の解法を与えることができる.

[定理 3.3]

式 (3.12) の解は, テンプレート:制御と振動の数学/equation テンプレート:制御と振動の数学/equation の解を y(t) とするとき,y(tt0) で与えられる.


証明

式 (3.13) の解を y(t) とする. そうすれば,定常性の原理 Ⅱにより, テンプレート:制御と振動の数学/equation式 (3.12) を満足し,しかも, テンプレート:制御と振動の数学/equation となる[1]


例68 テンプレート:制御と振動の数学/equation テンプレート:制御と振動の数学/equation を解け.

テンプレート:制御と振動の数学/equation を Laplace 変換すると テンプレート:制御と振動の数学/equation となる.この原像は, テンプレート:制御と振動の数学/equation テンプレート:制御と振動の数学/equation ここで,τ+t0 を改めて τ とおくと[2]テンプレート:制御と振動の数学/equation となる.これが求める結果である.

一般的な解放も同様である.式 (3.13) を Laplace 変換すれば, テンプレート:制御と振動の数学/equation となる.ここに q(s) は高々 n1 次の s の多項式である. テンプレート:制御と振動の数学/equation の原像を求めるため, テンプレート:制御と振動の数学/equation とおけば, テンプレート:制御と振動の数学/equation テンプレート:制御と振動の数学/equation ここで τ+t0τ とおくと, テンプレート:制御と振動の数学/equation を得る.

原像を求める際,(3.13b) で, テンプレート:制御と振動の数学/equation とおいても誤りではない.このとき,y(tt0) は, テンプレート:制御と振動の数学/equation となる.


例69

テンプレート:制御と振動の数学/equation を直接示せ.

解答例

I=t0tg(tτ)f(τ)dτ
にて,θ=tτ とおいて,τθ の変数変換を行う.
τ:t0t のとき,θ:tt00,また τ=tθ,dθ=dτ
I=tt00g(θ)f(tθ)dθ
=0tt0g(θ)f(tθ)dθ=0tt0g(τ)f(tτ)dτ


例70

x+ax=f(t),x(t0)=x0 の解は,

テンプレート:制御と振動の数学/equation または テンプレート:制御と振動の数学/equation となることを示せ.

解答例

例 69 と同じ推論で解く.

y+ay=f(t+t0),y(0)=x0 にて
Y(s)y(t) とおき両辺をラプラス変換すると,
sYx0+aY=[f(t+t0)]
(s+a)Y=x0+[f(t+t0)]
Y=x0s+a+1s+a[f(t+t0)]
y(t)=x0eat+eat*f(t+t0)…①


①からの展開は二通りある,まず一つは,
y=x0eat+0tea(tτ)f(τ+t0)dτ
x(t)=y(tt0)=x0ea(tt0)+0tt0ea(tt0τ)f(τ+t0)dτ
t0+τ=θ とおいて,τθ の積分変数変換を行う.
このとき dτ=dθτ:0tt0 のとき θ:t0t
すなわち,x(t)=x0ea(tt0)+t0tea(tθ)f(θ)dθ
これが解の表現の一つ.


もう一つは①から, y=x0eat+0teaτf(t+t0τ)dτ
よって x=y(tt0)=x0ea(tt0)+0tt0eaτf(tt0+t0τ)dτ
=x0ea(tt0)+0tt0eaτf(tτ)dτ


例71

テンプレート:制御と振動の数学/equation テンプレート:制御と振動の数学/equation を解け.

解答例

直前の x(t)=x0(tt0)+0tt0g(τ)f(tτ)dτ を使う. t=t0 における初期値はすべて 0 より過渡解 x00
p(s)=s42s3+2s22s+1
=(s1)2(s2+1)
1p(s)=12s1+12(s1)2+12ss2+1
1p(s)g とおくと,g(t)=12{et+tet+cost}
x=0tt0g(τ)f(tτ)dτ
=120tt0(eτ+τeτ+cosτ)f(tτ)dτ



  1. yp(D)x=f(t+t0) の解,すなわち p(D)y=f(t+t0) であるならば, y(tt0)p(D)x=f{(t+t0)t0}=f(t) の解である.
  2. したがって置き換える前の τ について,τ:0tt0 のとき,τ+t0:t0t