1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
| | %indent
大学で習う部分分数分解は、高校で習う剰余定理である。
;,物理・工学の諸分野で線形微分方程式を解く必要がある。
;,その中、ラプラス変換で線形微分方程式を手軽に解く分野がある。
;,例えば、電気工学や制御工学、振動工学などが特にラプラス変換を好む。
;,ラプラス変換自体は以下の積分式により定義されるが、
;,ラプラス変換で微分方程式を解く際は予め計算されるラプラス変換表を使う。
#ceq(e)
正変換： $$ \mathcal{L}[f(t)] $ = $ F(\:s) $ = $ \int_0^\infty $ f(t) $ e^{-\:st} $ dt $$(($$ \:s $$は複素数。))
#ceq(e)
逆変換：$$ \mathcal{L}^{-1}[F(\:s)] $ = $ f(t) $ = $ \ffd{1}{2\pi \:i} $ \int_{\sigma - \:i\infty}^{\sigma + \:i \infty} $ F(s) $ e^{st} $ ds $$((\:i は複素数単位。$$ \:i^2 $ = $ -1 $$。))
#ceq(d)
|*ラプラス変換表|<|<|<|<|<|<|h
|*原関数$$ f(t)$$|$$ 1 $$|$$ t $$|$$ t^n $$|$$ e^{at} $$|$$ \sin(at) $$|$$ \cos(at) $$|
|*像関数$$ F(\:s)$$|$$ \ffd{1}{\:s} $$|$$ \ffd{1}{\:s^2} $$|$$ \ffd{n!}{\:s^{n+1}} $$|$$ \ffd{1}{\:s-a} $$|$$ \ffd{a}{\:s^2+a^2} $$|$$ \ffd{\:s}{\:s^2+a^2} $$|
;,ここで、良く使う関数のラプラス変換結果は$$ s $$の有理式((多項式を分母と分子とする分数式))になっているのが分かる。
;,このため、ラプラス変換で線形微分方程式を解く際には、複雑な有理式を表内の有理式に分解する計算が大量発生する。
;,例えば、$$ F(\:s) $ = $ \ffd{1}{\:s^2 + 3\:s + 2} $ = $ \ffd{1}{\:s+1} $ - $ \ffd{1}{\:s+2} $$。
;,表内の分数の線形結合にさえ分解できれば、直ちに微分方程式の解が分かる。
;,例の場合は、$$ e^{at} $ \Longrightarrow $ \ffd{1}{\:s-a} $$より、$$ e^{-t} $ \Longrightarrow $ \ffd{1}{\:s+1} $$、$$ e^{-2t} $ \Longrightarrow $ \ffd{1}{\:s+2} $$が分かり、
;,方程式の解はラプラス変換が$$ F(\:s) $$になる関数$$ f(t) $ = $ e^{-1} $ - $ e^{-2t} $$となる。
;,分母の因数分解で$$ \ffd{1}{\:s^2 + 3\:s + 2} $ = $ \ffd{1}{(\:s+1)(\:s+2)} $$までは良いが、
;,その先は恒等式$$ \ffd{1}{(\:s+1)(\:s+2)} $ = $ \ffd{A}{\:s+1} $ + $ \ffd{B}{\:s+2} $$を解くことになる。
;,この分母が高次多項式の有理式を、分母が低次多項式の有理式の線形和に分解する式変形を部分分数分解と言う。
;,部分分数分解は分母を払ってから展開して係数比較しても地道に解けるが、次のテクニックが良く知られている。
#ceq(e)
$$ A $ = $ (\:s+1) $ F(\:s) $ |_{s=-1} = $ \cancel{(\:s+1)} $ \left.\ffd{1}{\cancel{(\:s+1)}(s+2)}\right|_{s=-1} $ = \ffd{1}{-1+2} $ = $ +1 $$
#ceq(e)
$$ B $ = $ (\:s+2) $ F(\:s) $ |_{s=-2} = $ \cancel{(\:s+2)} $ \left.\ffd{1}{(s+1)\cancel{(\:s+2)}}\right|_{s=-2} $ = \ffd{1}{-2+1} $ = $ -1 $$
#ceq(d)
%bodynote
;,一般に、多項式$$ P(x) $$を１次式$$ x - a $$で割った商を$$ Q(x) $$、余りを$$ R $$と置くと、
#ceq(e)
$$ P(x) $ = $ (x - a) $ Q(x) $ + $ R $$
#ceq(d)
;,と表せて、
#ceq(e)
$$ P(a) $ = $ \cancelto{0}{(a - a)}\;\;\; $ Q(a) $ + $ R $ = $ R $$
#ceq(d)
;,が成り立つ。
;,これが高校で学ぶ剰余定理である。
;,そこで、$$ F(\:s) $ = $ \ffd{A}{\:s+1} $ + $ \ffd{B}{\:s+2} $$を眺めてみると、
;,$$ B $$の分母を払って、$$ P_{2}(\:s) $ = $ (\:s+2) $ F(\:s) $ = $ (\:s+2) $ \ffd{A}{s+1} $ + $ B $$と置くと、剰余定理とそっくりな式を作り出せる。
;,そのため、$$ P_{2}(-2) $ = $ \cancelto{0}{(-2+2)}\;\;\; $ \ffd{A}{-2+1} $ + $ B $ = $ B $$が導ける(($$ Q(\:s) $ = $ \ffd{A}{\:s+1} $$が多項式ではないが、係数$$ a+2 $$が$$ 0 $$になるのが本質で、$$ Q(-a) $$が$$ \ffd{1}{0} $$にならない限り、問題なく成り立つ。))。
;,同様に$$ A $$の分母を払えば、$$ P_{1}(\:s) $ = $ (\:s+1) $ F(\:s) $ = $ A $ + $ (\:s+1) $ \ffd{B}{\:s+2} $$と置くと、
;,$$ P_{1}(-1) $ = $ A $ + $ \cancelto{0}{(-1+1)}\;\;\; $ \ffd{B}{-1+2} $ = $ A $$が得られる。
;,そこで、$$ P_{1}(\:s) $$や$$ P_{2}(\:s) $$の計算に使う$$ F(\:s) $$に$$ \ffd{1}{(\:s+1)(\:s+2)} $$を使えば、
;,前節で紹介した部分分数分解のテクニックとなる。
%bodynote
|
![[PukiWiki]](image/NekoPunch.png "[PukiWiki]")
