==ラプラス変換の基本==
【1. はじめに】
ラプラス変換とは何かということを,数式を使って正確に定義してみても,初めてそれを学ぶ人にはそもそもどのようなものなのかというイメージがないので,伝わらない.ここでは,具体例でイメージをつかんでもらって,簡単な作業ができるところまでを目指します.(あまり高度な内容は扱いません).このような変換は, 積分変換の例として,フーリエ変換やラプラス変換があるが,ここではラプラス変換を扱う. ラプラス変換では,核 この教材では,ラプラス変換の具体例とそれを利用して微分方程式を解くところまでを扱う. おそらく,この教材を眺めているだけでは身に着かないので,鉛筆を持って簡単な問題を幾つか解くことが習得のための秘訣ですのでよろしく. |
【2. ラプラス変換とは】
元の関数が
【3. ラプラス変換で使われる記号】
関数のとき と書く.( このとき,元の関数は独立変数,従属変数とも実数であるが,変換後の関数は独立変数,従属変数とも複素数で考える.(ただし,このページで取り上げる問題のレベルでは,複素数であることを意識しなくてもできる) と書く. |
【4. ラプラス変換を用いた微分方程式の解き方】
[簡単な図解でイメージ作り]![]()
元の微分方程式
変換された方程式
解の方程式
変換された解
(注1) yやy',その他よく使う関数のラプラス変換は,あらかじめ一覧表にしておく. (注2) ラプラス逆変換も,個別に計算をすることはめったになく,あらかじめ作っておいた表の右欄から左欄を読むという使い方を考える. [最低限覚えるべき手順] (1) (2)
(証明)
(3)のとき 次の部分積分を行う ラプラス変換は となる関数に対して定義される. もしそうでなければ は収束しないから また,次の式が成り立つ したがって (4) (*) |
なぜ,
ラプラス変換 が収束するためには,この被積分関数の極限値が0でなければならない. ![]() もしそうでなければ,たとえば正の値 結局,ラプラス変換が収束するためには,被積分関数の極限値が0になることが必要条件である. 微積分でとてもよい働きをするのが指数関数であるが,たとえば は発散してしまう.指数関数の内で, また,この関数は 結局,
【用語】 収束座標
正確に言えば, |
【5. よく使うラプラス変換とその証明】
(1)次の部分積分を行う だから,このラプラス変換が収束するには 次々に,この部分積分を繰り返すと この式は, |
(4) 次の部分積分を行う ここで, さらに,次の部分積分を行う ここで, 証明は(3)とほぼ同様[略] |
(6)推移則[シフト公式]
入力関数に
入力関数に
(証明)のとき だから,成り立つ. |
推移則[シフト公式]を使うと,各々の公式が次のように関連していることが分かる. (1)→(3) だから (2)→(*3) (*3)は後出 だから (4)→(*1) (*1)は後出 だから (5)→(*2) (*2)は後出 だから |
【6. ラプラス変換表】
(入力変数を*** 多項式 *** |
*** 双曲線関数 ***
|
【7. ラプラス変換の基本性質】
(1) 線形性
(証明)
(2) 自己相似性 (
(証明)
に対して |
【例1】
(解答)次の関数のラプラス変換を求めてください.ただし,出力変数は (1) (2) (3) (1)ラプラス変換の線形性により,次のように分けられる. 次に,ラプラス変換表から (2)ラプラス変換の線形性により,次のように分けられる. 次に,ラプラス変換表から (3)ラプラス変換の線形性により,次のように分けられる. 次に,ラプラス変換表から |
【例2】
(解答)次の関数のラプラス変換を求めてください.ただし,出力変数は (1) (2) (3) (1) (*1)の公式にa=−2, b=3を代入する (2) (*3)の公式にa=−2, b=3を代入する (3) 積や累乗の形になっているとラプラス変換の線形性が利用できないので,半角公式を使って三角関数の1次式にしておく |
【例3】
(解答)次の関数のラプラス逆変換を求めてください.ただし,変数は (1) (2) (3) (1) 部分分数に分解する (2) 部分分数に分解する (3) 部分分数に分解する とおいて定数 より この連立方程式を解くと |
線形な微分方程式(
【1階線形微分方程式】
(3.1)微分方程式
≪A:ラプラス変換による解き方≫初期条件 を満たす 両辺をラプラス変換する
右辺を部分分数に分解する.
とおいて,両辺の係数を比較すると ≪B:1階線形微分方程式の筆算による解き方≫
同次方程式
この公式を使う |
同次方程式 を解く 非同次方程式の1つの特別解は, 次の部分積分を行う ≪C:フリーソフトwxMaximaによる解き方≫ メニューをたどって,「方程式」→「微分方程式を解く」
siki1:ode2('diff(y,x)-y=2*x, y, x);
のように方程式に名前を付けておく y=(2*(-x-1)*%e^(-x)+%c)*%e^x という形で一般解が得られる ic1(siki1,x=0,y=0); という形で初期条件を入れる y=2*%e^x-2*x-2 が返される |
(3.2)微分方程式
≪A:ラプラス変換による解き方≫初期条件 を満たす 両辺をラプラス変換する 右辺を部分分数に分解する ≪B:1階線形微分方程式の筆算による解き方≫ 同次方程式 を解く |
非同次方程式の1つの特別解を求める.右辺の形から したがって,非同次方程式の1つの特別解は 元の微分方程式の一般解は 初期条件を入れると 結局,求める方程式は ≪C:フリーソフトwxMaximaによる解き方≫ メニューをたどって,「方程式」→「微分方程式を解く」
siki1:ode2('diff(y,x)-y=2*exp(3*x), y, x);
のように方程式に名前を付けておく という形で一般解が得られる ic1(siki1,x=0,y=1); という形で初期条件を入れる が返される |
【2階線形微分方程式】
(3.3)微分方程式
≪A:ラプラス変換による解き方≫初期条件 を満たす 両辺をラプラス変換する 初期条件を入れると
右辺を部分分数に分解する.
とおいて,両辺の係数を比較すると 5次: 4次: 3次: 2次: 定数項: 求まるものから順に求めると ≪B:2階線形微分方程式の筆算による解き方≫
非同次方程式
同次方程式 非同次方程式 |
非同次方程式 ( 同次方程式 の特性方程式 元の微分方程式の一般解は ≪C:フリーソフトwxMaximaによる解き方≫ メニューをたどって,「方程式」→「微分方程式を解く」
siki1:ode2('diff(y,x,2)-y=sin(x)+x, y, x);
のように方程式に名前を付けておく ※2階導関数を表現するには,'diff(y,x,2)とするとよい という形で一般解が得られる ic2(siki1,x=0,y=0,'diff(y,x)=0); という形で初期条件を入れる (2階微分方程式で初期条件に対応する特殊解を求めるには…initial value problems for 2nd order differential equations・・・ic2( )関数が使える.y'=0という条件は,'diff(y,x)=0で書く・・・先頭のアポストロフィーに注意) |
(3.4)微分方程式
≪A:ラプラス変換による解き方≫初期条件 を満たす 両辺をラプラス変換する 初期条件を入れると
右辺を部分分数に分解する.
とおいて,両辺の係数を比較すると 3次: 2次: 1次: 定数項: この連立方程式を解くと ≪B:2階線形微分方程式の筆算による解き方≫
非同次方程式
同次方程式 非同次方程式 |
非同次方程式 ( 同次方程式 の特性方程式 元の微分方程式の一般解は ≪C:フリーソフトwxMaximaによる解き方≫ メニューをたどって,「方程式」→「微分方程式を解く」
siki1:ode2('diff(y,x,2)-3*'diff(y,x,1)+2*y
=2*cos(x), y, x); のように方程式に名前を付けておく ※2階導関数を表現するには,'diff(y,x,2)とするとよい という形で一般解が得られる ic2(siki1,x=0,y=1,'diff(y,x)=1); という形で初期条件を入れる (2階微分方程式で初期条件に対応する特殊解を求めるには…initial value problems for 2nd order differential equations・・・ic2( )関数が使える.y'=1という条件は,'diff(y,x)=1で書く・・・先頭のアポストロフィーに注意) |
(3.5)微分方程式
≪A:ラプラス変換による解き方≫初期条件 を満たす 両辺をラプラス変換する 初期条件を入れると
右辺を部分分数に分解する.
とおいて,両辺の係数を比較し,連立方程式を解くと
ここで,次のラプラス変換公式を使う
≪B:2階線形微分方程式の筆算による解き方≫ 非同次方程式 ( |
同次方程式 の特性方程式 元の微分方程式の一般解は ≪C:フリーソフトwxMaximaによる解き方≫ メニューをたどって,「方程式」→「微分方程式を解く」
siki1:ode2('diff(y,x,2)+2*'diff(y,x,1)+4*y
=2*x^2, y, x); のように方程式に名前を付けておく ※2階導関数を表現するには,'diff(y,x,2)とするとよい という形で一般解が得られる ic2(siki1,x=0,y=0,'diff(y,x)=0); という形で初期条件を入れる (2階微分方程式で初期条件に対応する特殊解を求めるには…initial value problems for 2nd order differential equations・・・ic2( )関数が使える.y'=0という条件は,'diff(y,x)=0で書く・・・先頭のアポストロフィーに注意) |
推移則[シフト公式]の補足問題 次の図で示された関数のラプラス変換を求めよ.
(1) 矩形波
(解答)… ラプラス変換の定義に当てはめる![]() { }内は,公比 (別解)… 推移則[シフト公式](6.2)を使う に対して,2つの目の矩形はtの正の向きに2aだけ平行移動したものだから,積分は第1項に 以下,次々に |
(2) 鋸波
(別解)… 推移則[シフト公式](6.2)を使う![]() 2つの目の鋸形はtの正の向きにaだけ平行移動したものだから,積分は第1項に 以下同様にして,公比 |
(3) 矩形波
(別解)… 推移則[シフト公式](6.2)を使う![]() 次々に周期 ※この形でもよいが,双曲線関数 |
(4) 三角波
(別解)… 推移則[シフト公式](6.2)を使う![]() (途中経過略) 次々に周期 前問(3)と同様に,この答も双曲線関数を使って表せる. |
笆�縺薙�繧オ繧、繝亥�縺ョGoogle讀懃エ「笆� |