• 検索結果がありません。

ラプラス変換と微分方程式の初期値問題

ドキュメント内 Fourier Fourier Gibbs (ページ 86-95)

ただし,F(s) =L(i(t))(s).ここで,F(s)で上の方程式を解くと

F(s) = 1

Ls+R+ 1 Cs

L(v(t))(s) + 1 Ls+R+ 1

Cs Li(0)

ここで,さらに時刻t= 0での電流をi(0) = 0と仮定すると

F(s) = s

Ls2+Rs+C1L(v(t))(s) この伝達関数のラプラス逆変換をL1

( s

Ls2+Rs+C1 )

=w(t)とすると,合成積の(10)より,

L(i(t))(s) =F(s) =L(w∗v)(s) ラプラス逆変換すれば,次を得る.

i(t) = (w∗v)(t) =

t 0

w(τ)v(t−τ)dτ 結論として,まとめれば,次の定理になる.

定理5.31 (Duhamel(デュアメル)の合成定理) (4.26)の電流i(t)i(t) =

t 0

w(τ)v(t−τ)dτ (5.27)

で与えられる.ただしw(t)は伝達関数H(s) = s

Ls2+Rs+C1 のラプラス逆変換である.

注意5.32 定理5.31で述べたことは,(4.1)のRLC回路の微分方程式を,より一般した定係数線形微分方程式に 対しても成立する.

 具体的に,L1

( s

Ls2+Rs+C1 )

=w(t)を求めてみよう.Ls2+Rs+C1= 0が異なる解α, βを持つ場合 と重解を持つ場合に分けて議論する.

(1) 異なる場合Ls2+Rs+C1=L(s−α)(s−β), (α̸=β) 部分分数展開により,

1

L(s−α)(s−β) = 1

α−β = 1 L(α−β)

( α

s−α− β s−β

)

; 逆ラプラス変換をすると,命題4.22の(5),例4.21の(1)より,

w(t) =L1

( 1

L(s−α)(s−β) )

= 1

L(α−β) (

L1 ( α

s−α )

− L1 ( β

s−β ))

= 1

L(α−β)

(αeαt−βeβt)

(i)α, βが実数の場合(CR2>4L),L, R, C0なので,α, β0となる.

w(t) = 1 L(α−β)

(αeαt−βeβt)

(ii)α, βが実数でない場合(CR2<4L),ℜα=−R

2L <0, β=αなので w(t) = 1

L(α−β)

(αeαt−βeβt)

= 1

L(α−α)

(αeαt−αeαt)

= 1

Lℑ(α)(αeαt) ただし,(α), (α)は複素数αのそれぞれ実部,虚部を表す.

(2) 重解の場合Ls2+Rs+C1=L(s−α)2 すなわち,α=−R

2L <0, CR2= 4L 部分分数展開より,

1

L(s−α)2 = 1 L

( 1

s−α+ α (s−α)2

)

; 逆ラプラス変換をすると,命題4.22の(5),例4.21の(1)より,

w(t) =L1

( 1 L(s−α)2

)

= 1 L

( L1

( 1 s−α

) +L1

( α (s−α)2

))

= 1 L

(eαt+teαt)

 いずれの場合も,lim

t→∞w(t) = 0のなり,伝達関数H(s)の解が実数でない場合,振動しながら減衰する.

5.33

x′′+x= 2 cost, x(0) = 3, x(0) = 4 この式をラプラス変換すると(X(s) =L(x(t))(s)で表す),

L(x′′) +L(x) = 2L(cost)

s2X(s)−s·34 +X(s) = 2 s s2+ 1 (s2+ 1)X(s) = 2s

s2+ 1 + 3s+ 4

X(s) = 2s

(s2+ 1)2 +3s+ 4 s2+ 1

= d

ds (

1 s2+ 1

)

+ 3 s

s2+ 1 + 4 1 s2+ 1

x(t) =L1 (d

ds (

1 s2+ 1

))

+ 3L1 ( s

s2+ 1 )

+ 4L1 ( 1

s2+ 1 )

=tsint+ 3 cost+ 4 sint5.34 (積分方程式)

x(t) =t+

t 0

x(τ) sin(t−τ)dτ 式全体をラプラス変換して,畳み込み関数(10)を用いると,

X(s) =L(x(t)) =L(t) +L (∫ t

0

x(τ) sin(t−τ) )

= 1

s2 +L(x(s))L(sint)(s)

= 1

s2 +X(s) 1 s2+ 1; X(s)を解くと,

X(s) =s2+ 1 s4 = 1

s2 + 1 s4

L1(X(s)) =t+t3 6

5.35 (超関数を含む微分方程式 1)

x′′+ 4x+ 5x=δ(t−1), x(0) = 0, x(0) = 3

式全体をラプラス変換すると,

s2X(s)0·s−3 + 4(sX(s)0) + 5X(s) =L(δ(t1)) =es X(s)で解くと,

(s2+ 4s+ 5)X(s) = 3 +es

X(s) = 3

(s+ 2)2+ 1+ es (s+ 2)2+ 1 第一移動(3)を用いると,

x(t) = 3e2tsint+e2tsin(t1)·U(t−1)

5.36 (超関数を含む微分方程式 2)

x′′+ 4x+ 8x=δ(t), x(0) = 1, x(0) = 2 この式をラプラス変換すると,

s2X−sx(0)−x(0) + 4(sX−x(0)) + 8X = 1 s2X−s−2 + 4sX4 + 8X= 1

(s2+ 4s+ 8)X =s+ 7

X = s+ 7

s2+ 4s+ 8 = s+ 2

(s+ 2)2+ 22 + 5 (s+ 2)2+ 22 ラプラス逆変換をすると,

x(t) =L1

( s+ 2 (s+ 2)2+ 22

) +5

2L1

( 2 (s+ 2)2+ 22

)

=e2t (

cos 2t+5 2sin 2t

)

求める解は,次で与えられる.

x(t) =e2t (

cos 2t+5 2sin 2t

)

(5.28)

注: ここで,この解が与えられた微分方程式の条件を満たすかを見てみよう.

x(t) =e2t (

cos 2t+5 2sin 2t

)

; x(t) =2e2t

(

cos 2t+5 2sin 2t

)

+e2t(2 sin 2t+ 5 cos 2t) =e2t(3 cos 2t7 sin 2t) ; x′′(t) =2e2t(3 cos 2t7 sin 2t) +e2t(6 sin 2t14 cos 2t) =e2t(20 cos 2t+ 8 sin 2t) x′′+ 4x+ 8xに代入すると

x′′(t) + 4x(t) + 8x(t) = 0 となり,外力δ(t)に右辺はならない.また,t= 0を解に代入すると,

x(0) = 1, x(0) = 3 初期条件x(0) = 1, x(0) = 2と違ってくる.

 実は解(5.28)は,微分方程式が超関数を含むとき,ラプラス変換して求めた解は超関数としての解であって,普

通の意味での解にはなっていないこのがその理由である.

5.37 連立微分方程式

{x(t) +y(t) = 2 cost, (x(0) = 0, y(0) = 1) x(t) +y(t) = 0

それぞれ式全体をラプラス変換すると,次の連立方程式をえる.

sX(s)−0 +Y(s) = 2 s

s2+ 1, X(s) +sY(s)1 = 0 ただし,X(s) =L(x)(s), Y(s) =L(y)(s) この連立方程式を解くと,

X(s) = 1

s2+ 1, Y(s) = s s2+ 1 ラプラス逆変換すると,

x(t) =L1 ( 1

s2+ 1 )

= sint, y(t) =L1 ( s

s2+ 1 )

= cost

6 補足

6.1 微分と積分の順序交換 , Wierstrass の多項式近似定理

定理6.1 (微分と積分の順序交換) (1) 閉区間[a, b]×[c, d]上のf(x, y)はC1-級の関数とするとき,

d dy

(∫ b a

f(x, y)dx )

=

b a

∂f

∂y(x, y)dx (6.1)

(2) a−∞,あるいは,bがの場合,a < λ, λ < bについて,

λ a

∂f

∂y(x, y) dx+

b λ

∂f

∂y(x, y)

dx λb, λ−→a 0, (変数y∈[c, d]について一様収束) (6.2) ならば,この場合も,(6.1)は成立する.

証明 (1)フビニと積分順序の交換により,

t c

(

b a

∂f

∂y(x, y)dx)dy=

b a

(

t c

∂f

∂y(x, y)dy)dx=

b a

f(x, t)−f(x, c)dx=

b a

f(x, t)dx

b a

f(x, c)dx 右辺の最後の項は定数なので,両端をtで微分すると,求める式が得られる.

b a

∂f

∂y(x, t)dx= d dt(

b a

f(x, t)dx) (2)a=−∞, b=の場合を示す(他も同様).

−∞f(x, y+ ∆y)dx−

−∞f(x, y)dx

∆y

−∞

∂f

∂y(x, y)dx

= ∫

−∞

f(x, y+ ∆y)−f(x, y)

∆y −∂f

∂y(x, y)dx 平均値の定理より,

= ∫

−∞

∂f

∂y(x, y+θ∆y)−∂f

∂y(x, y)dx

(0< θ <1)

−∞

∂f

∂y(x, y+θ∆y)−∂f

∂y(x, y) dx 条件(6.2)より,十分大きいλ >λ0:

|x|

∂f

∂y(x, y)

dx < ε, y∈[c, d]より,

2ε+

λ

λ

∂f

∂y(x, y+θ∆y)−∂f

∂y(x, y) dx

∂f

∂y(x, y)は[−λ, λ]で一様連続なので,δ >0 : |y−y|< δ ∂f

∂y(x, y)−∂f

∂y(x, y)

< ε.従って,

2ε+ 2λε, (|∆y|< δ ) 以上により,

d dy

(∫

−∞

f(x, y)dx )

= lim

∆y0

−∞f(x, y+ ∆y)dx−

−∞f(x, y)dx

∆y =

−∞

∂f

∂y(x, y)dx

(注: (2)の証明方法は,(1)の別証明にもなっている)

定理6.2 (Wierstrassの多項式近似定理) [0,1]上の連続関数f(x)とする.任意のε >0に対して,

max{|f(x))−p(x)|: x∈[0,1]}< ε (6.3) となる多項式p(x)が存在する.

証明 P[0,1]を[0,1]上の多項式全体,C[0,1]を[0,1]上の連続関数全体とする.このとき,一様ノルムでP[0,1]

C[0,1]で稠密であることを,次の方法で示す(この事柄は大事で,証明方法は2項分布を用いたベルンシュタ インの方法が一般的である).

Tn(x) = cos(nCos1x)n次の多項式であることを示す.この多項式をチェビシェフ( ˇCebyˇsev)の多項式と

いう.

cos+isin= (cosθ+isinθ)nの右辺を2項定理で展開して,その実部を見ると,

cos= ( n

k=0

nCkiksinkθcosnkθ )

=

[n2]

m=0

nC2m(1)msin2mθcosn2mθ

=

[n2]

m=0

nC2m(1)m(1cos2θ)mθcosn2mθ ここで,x= cosθとすると,

Tn(x) = cos(nCos1x) =

[n2]

m=0

nC2mxn2m(x21)m 右辺はn次の多項式であるので,Tn(x)はn次の多項式である.

 任意の連続関数f ∈C[0,1]に対して,[−π, π]で偶関数になるようf(x)を拡張しておく.そして,

g(θ) =f(cosθ), (−π < θ≤π)

と定義すると,cosθも偶関数なのでg(θ)は連続な偶関数になる.すると,gの実フーリエ係数bn = 0となり,

フーリエ級数展開はcosnθ, (nN∪ {0})の一次結合で表される.従って,チェサロ平均 σn(θ) =

n k=0

dncos

と表す.ただし,dnC.定理1.36よりチェサロ平均σn(θ)はg(θ)に一様収束する.すなわち,任意のε >0に

対して,

g(θ)−

n k=0

dkcos

< ε (θ∈[−π, π]) となるnが存在する,ここで,θ= Cos1xを代入すると,

f(x)

n k=0

dkcos(

kCos1x)

< ε (x∈[0,1])

f(x)

n k=0

dkTk(x)

< ε (x∈[0,1]) 求めるn次の多項式p(x)p(x) =

n k=0

dkTk(x)である.

索 引

記号

∥ · ∥2. . . 5

A(s), B(s) . . . 27

an, bn. . . 8

cn. . . 8

CW Tf(a, b). . . .38

δ(x) . . . 59, 72 Dn(x) . . . 17

dnT dx . . . 59

f(a0) . . . 10

f(a+ 0) . . . 10

f∗g. . . 28, 66, 77 F(f)(s), fˆ(s) . . . 22

F1(f)(s), fˇ(s) . . . 22

f#(t) . . . 28

Gb. . . 39

Γ(x) . . . 71

H(iω) . . . 64

H(s) . . . 68

h(x) . . . 12

Kn(x) . . . 18

L2(a, b) . . . 3

L(f)(s), F(s) . . . 64, 70 L1(F(s)) . . . 65, 70 S(R). . . 57

pα,β. . . 57

Pn(t) . . . 4

ψj,k(t) . . . 4

S(R) . . . 57

σn(x) . . . 17

Sn(x). . . .11

Tf. . . 58

U(t) . . . 72

<s, t>. . . 38

B Bessel(ベッセル)の不等式. . . 6

Bromwich(ブロムウィッチ)の反転公式. . . 76

Bromwich(ブロムウィッチ)の定理. . . .65

C Cauchy-Schwarzの不等式. . . 5, 34 convolution(対合). . . 28

Cr-級の関数. . . 10

D D’lambert(ダランベール)の公式. . . 47

Duhamel(デュアメル)の合成定理. . . 68, 87 F Fubbini(フビニ)の定理. . . 29

L Lerch(レルヒ)の定理. . . 75

P Parseval(パーセバル)の等式. . . 8

Plancherel(プランシエル)公式. . . 28

R RLC直列回路. . . 57

W Weierstrassの優収束判定法. . . 21

Wierstrassの多項式近似定理. . . 92

Wilbraham-Gibbs定数. . . .16

位数. . . 85

インパルス応答. . . 64

インピーダンス. . . 68

ウェーブレット基底. . . 4

n次元フーリエ級数. . . 37

n次元フーリエ変換とフーリエ逆変換. . . 38

エネルギー保存則. . . 38

応答関数. . . 86

外力. . . 86

重ね合わせの原理. . . 44

緩増加関数. . . 58

(緩増加)超関数(Tempered distribution) . . . 57

Γ-関数. . . 71

ギブスギャップ(Gibbs) . . . 14, 16 ギブス定数. . . 16

逆フーリエ変換. . . 22

逆離散フーリエ変換. . . .37

急減少関数. . . 31, 57 境界条件. . . 40

極. . . 85

区分的に滑らか. . . 10, 22 区分的に連続. . . 10

群上の調和解析. . . 39

計算機断層撮影. . . 37

ゲイン. . . 64

高速フーリエ変換(FFT). . . 37

項別微分. . . 20

固有関数(特性関数). . . 43

固有値. . . 43

混合問題. . . 54

3次元熱方程式. . . 40

3次元波動程式. . . 40

3次元ラプラス方程式. . . 40

自己相関関数. . . 33

指数型. . . 74

実フーリエ係数. . . 8

指標. . . 39

シャノン-染谷のサンプリング定理. . . 34

収束域. . . 65, 74 収束座標. . . 65, 73 周波数伝達関数. . . 64

初期条件. . . 40

スペクトル. . . 43

積分. . . 66, 77 積分方程式. . . 88

絶対可積分. . . 22

絶対収束座標. . . 74

全エネルギー. . . 33

線形システム. . . 57

線形性. . . 66, 77 像(関数). . . 70

相互相関関数. . . 33

相似. . . 66, 77 双対群. . . 39

像の移動. . . 66, 77 像の積分. . . 66, 77 像の第一移動. . . 71

像の微分. . . 66, 77 帯域制限関数. . . 34

第一移動. . . 66, 77 第二移動. . . 66, 77 畳み込み関数. . . 28, 66, 77 妥当性. . . 44

単位インパルス. . . 64

チェサロ平均. . . 17

チェビシェフ( ˇCebyˇsev)の多項式. . . 92

超関数と関数の積. . . 60

超関数の収束. . . 58

超関数と関数の畳み込み. . . 60

超関数の微分. . . 59

超関数のフーリエ変換. . . 60

直交. . . 3

調和関数. . . 54

直交関数系. . . 6

定義セミノルム系. . . 57

ディリクレの積分公式. . . 25

ディリクレ題問. . . 54

テスト関数空間(Test function space) . . . 57

ディリクレ核. . . 17

デルタ関数. . . 59, 72 伝達関数. . . 68, 86 時不変. . . 64

特性関数. . . 68, 86 特性方程式. . . 86

内積. . . 3, 34 内積の性質. . . 3

滑らかな関数. . . 10

軟化子. . . 59

2次元ポアソン方程式. . . 40

2乗可積分関数. . . 33

2乗平均ノルム. . . 5

2乗平均収束. . . 6

熱核(heat kernel). . . 51

ノイマン問題. . . 54

反転公式(Inversion Formula) . . . 27

ピタゴラスの定理. . . 6

左極限. . . 10

左微分可能. . . 10

左微分係数. . . 10

微分. . . 66, 77 微分と積分の順序交換. . . 91

標本化定理. . . 34

フーリエ解析の基本定理. . . 13

フーリエ係数. . . 6

フーリエ積分(実フーリエ変換). . . 27

フーリエ展開. . . 9

フーリエの重積分公式. . . 26

フーリエの単積分公式. . . 26

フーリエ変換. . . 22

フェイエル核. . . 18

ドキュメント内 Fourier Fourier Gibbs (ページ 86-95)

関連したドキュメント