![]() ![]() *** 大区分 *** 数Ⅰ・A数Ⅱ・B数Ⅲ高卒・大学初年度 *** 中区分 *** ベクトル・行列連立方程式複素数関数・数列微分積分微分方程式統計maxima ※高校から大学初年度レベルの「微分方程式」について,このサイトには次の教材があります.
この頁へGoogleやYAHOO ! などの検索から直接来てしまったので「前提となっている内容が分からない」という場合や「この頁は分かったがもっと応用問題を見たい」という場合は,他の頁を見てください. が現在地です. |
← PC用は別頁
次の形の常微分方程式を1階線形常微分方程式といいます.
(解説)y'+P(x)y=Q(x)…(1) 方程式(1)の右辺:Q(x)を0とおいてできる同次方程式 (この同次方程式は,変数分離形になり比較的容易に解けます) y'+P(x)y=0…(2) の1つの解をu(x)とすると,方程式(1)の一般解は y=u(x)(∫ ![]() で求められます.
参考書には
上記のu(x)の代わりに,e−∫P(x)dxのまま書いて
y=e−∫P(x)dx(∫ Q(x)e∫P(x)dxdx+C)…(3')
と書かれているのが普通です.この方が覚えやすい人は,これで覚えるとよい.ただし,赤と青で示した部分は,定数項まで同じ1つの関数の符号だけ逆のものを使います.
筆者は,この複雑な式を見ると頭がクラクラ(目がチカチカ)して,どこで息を継いだらよいか困ってしまうので,上記の(3)のように同次方程式の解をu(x)として,2段階で表すようにしています. 同次方程式(2)は,次のように変形できるので,変数分離形です. y'+P(x)y=0 ![]() ![]() 両辺を積分すると ∫ ![]() log|y|=−∫ P(x)dx |y|=e−∫P(x)dx+A=eAe−∫P(x)dx=Be−∫P(x)dxとおく y=±Be−∫P(x)dx=Ce−∫P(x)dx…(4) (4)式は,Cを任意定数とするときに(2)を満たすが,そのままでは(1)を満たさない.
このような場合に,
同次方程式y'+P(x)y=0 の一般解の定数Cを関数に置き換えて, 非同次方程式y'+P(x)y=Q(x) の解を求める方法を定数変化法という.
なぜ,そんな方法を思いつくのか?自分にはなぜ思いつかないのか?などと考えても前向きの考え方にはなりません.思いついた人が偉いと考えるとよい.
(4)式において,定数Cを関数z(x)に置き換えて定数変化法は,数学史上に残るラグランジェの功績ですが,後からついていく我々は,ラグランジェが発見した方法のおいしいところをいただいて,節約できた時間を今の自分に必要なことに当てたらよいと割り切るとよい. ただし,この定数変化法は2階以上の微分方程式において,同次方程式の解から非同次方程式の解を求める場合にも利用できるなど適用範囲の広いものなので,「今度出てきたら,真似してみよう」と覚えておく値打ちがあります. u(x)=e−∫P(x)dxは(2)の1つの解 y=z(x)u(x)…(5) とおいて,関数z(x)を求めることにする. 積の微分法により:y'=(zu)'=z'u+zu'だから,(1)式は次の形に書ける. z'u+zu'+P(x)y=Q(x)…(1') ここでu(x)は(2)の1つの解だから u'+P(x)u=0 zu'+P(x)zu=0 zu'+P(x)y=0 そこで,(1')において赤で示した項が消えるから,関数z(x)は,またしても次の変数分離形の微分方程式で求められる. z'u=Q(x) ![]() dz= ![]() したがって z=∫ ![]() (5)に代入すれば,目的の解が得られる. y=u(x)(∫ ![]()
理論の上では上記のように解けますが,実際の積分計算
が難しいかどうかは
u(x)=e−∫P(x)dxや∫ ![]() になるかによります. すなわち,P(x)や ![]() 筆算では手に負えない問題になることがあります. |
||||
【例題1】
微分方程式y'−y=2xの一般解を求めてください.
この方程式は,(1)において,P(x)=−1, Q(x)=2xという場合になっています.
(解答)♪==定数変化法の練習も兼ねて,じっくりやる場合==♪ はじめに,同次方程式y'−y=0の解を求める.
【指数法則】…よく使う
ex+C1=exeC1 ![]() ![]() ∫ ![]() log|y|=x+C1 |y|=ex+C1=eC1ex=C2ex(eC1=C2とおく) y=±C2ex=C3ex(2=C3とおく) ±C 次に,定数変化法を用いて,3=z(x)とおいて Cy=zex(zはxの関数)の形で元の非同次方程式の解を求める. y=zexのとき y'=z'ex+zexとなるから 元の方程式は次の形に書ける. z'ex+zex−zex=2x z'ex=2x ![]() dz= ![]() ∫ dz=2∫ xe−xdx z=2∫ xe−xdx
∫ fg' dx=fg−∫ f 'g dxにより ∫ xe−xdx=−xe−x+∫ e−xdx=−xe−x−e−x+C4 z=2(−xe−x−e−x+C4) yに戻すと y=2(−xe−x−e−x+C4)ex y=−2x−2+2C4ex=−2x−2+Cex…(答) ♪==(3)または(3')は公式と割り切って直接代入する場合==♪ P(x)=−1だから,u(x)=e−∫P(x)dx=ex Q(x)=2xだから,∫ ![]() ![]() =2(−xe−x−e−x)+C したがって y=ex{ 2(−xe−x−e−x)+C}=−2x−2+Cex…(答) |
||||
【例題2】
微分方程式y'+2y=3e4xの一般解を求めてください.
この方程式は,(1)において,P(x)=2, Q(x)=3e4xという場合になっています.
(解答)♪==定数変化法の練習も兼ねて,じっくりやる場合==♪ はじめに,同次方程式y'+2y=0の解を求める. ![]() ![]() ∫ ![]() log|y|=−2x+C1 |y|=e−2x+C1=eC1e−2x=C2e−2x(eC1=C2とおく) y=±C2e−2x=C3e−2x(2=C3とおく) ±C 次に,定数変化法を用いて,C3=z(x)とおいてy=ze−2x(zはxの関数)の形で元の非同次方程式の解を求める. y=ze−2xのとき y'=z'e−2x−2ze−2xとなるから 元の方程式は次の形に書ける. z'e−2x−2ze−2x+2ze−2x=3e4x z'e−2x=3e4x ![]() dz=3e4xe2xdx=3e6xdx ∫ dz=3∫ e6xdx z=3∫ e6xdx = ![]() yに戻すと y=( ![]() y= ![]() ♪==(3)または(3')は公式と割り切って直接代入する場合==♪ P(x)=2だから,u(x)=e−∫2dx=e−2x Q(x)=3e4xだから,∫ ![]() = ![]() したがって y=e−2x{ ![]() ![]() |
※正しい番号をクリックしてください.
それぞれの問題は暗算では解けませんので,計算用紙が必要です.
※ブラウザによっては,番号枠の少し上の方が反応することがあります.
【問題1】
微分方程式y'−2y=e5xの一般解を求めてください. 1y= ![]() ![]() 3y= ![]() ![]()
≪同次方程式の解を求めて定数変化法を使う場合≫
同次方程式を解く: ![]() ![]() ∫ ![]() log|y|=2x+C1 |y|=e2x+C1=eC1e2x=C2e2x y=±C2e2x=C3e2x そこで,元の非同次方程式の解をy=z(x)e2xの形で求める 積の微分法によりy'=z'e2x+2e2xzとなるから z'e2x+2e2xz−2ze2x=e5x z'e2x=e5x 両辺をe2xで割ると z'=e3x z= ![]()
≪(3)または(3')の結果を使う場合≫
P(x)=−2だから,u(x)=e−∫(−2)dx=e2x Q(x)=e5xだから,∫ ![]() ![]() = ![]()
y=e2x(
![]() ![]() |
||||||||
【問題2】
微分方程式y'cosx+ysinx=1の一般解を求めてください. 1y=sinx+Ccosx 2y=cosx+Csinx 3y=sinx+Ctanx 4y=tanx+Csinx ヒント1ヒント2解答
≪同次方程式の解を求めて定数変化法を使う場合≫
元の方程式は y'+ytanx= ![]() そこで,同次方程式を解くと: ![]()
tanx=
![]() ![]() だから ∫tanx dx=−∫ ![]() =−log|cosx|+C ![]() ∫ ![]() log|y|=log|cosx|+C1 =log|eC1cosx| |y|=|eC1cosx| y=±eC1cosx y=C2cosx そこで,元の非同次方程式の解をy=z(x)cosxの形で求める 積の微分法によりy'=z'cosx−zsinxとなるから z'cosx−zsinx+zcosx tanx= ![]()
(tanx)'=(
z'cosx=![]() ![]() だから ∫ ![]() ![]() z'= ![]() ![]() ![]() dz= ![]() ∫ dz=∫ ![]() z=tanx+C
≪(3)または(3')の結果を使う場合≫
【元に戻る】…よく使う
P(x)=tanxだから,elog A=A log eA=A u(x)=e−∫tanxdx=elog|cos x|=|cos x| その1つはu(x)=cos x Q(x)= ![]() ![]() ![]() =tanx+C
y=(tanx+C)cosx=sinx+Ccosxになります.→1
|
||||||||
【問題3】
微分方程式xy'−y=2x2+xの一般解を求めてください. 1y=x(x+log|x|+C) 2y=x(2x+log|x|+C) 3y=x(x+2log|x|+C) 4y=x(x2+log|x|+C) ヒント1ヒント2解答
≪同次方程式の解を求めて定数変化法を使う場合≫
元の方程式は y'− ![]() 同次方程式を解く: ![]() ![]() ![]() ![]() ∫ ![]() ![]() log|y|=log|x|+C1=log|x|+log eC1=log|eC1x| |y|=|eC1x| y=±eC1x=C2x そこで,元の非同次方程式の解をy=z(x)xの形で求める 積の微分法によりy'=z'x+zとなるから z'x+z− ![]() z'x=2x+1 両辺をxで割ると z'=2+ ![]() z=2x+log|x|+C
≪(3)または(3')の結果を使う場合≫
元の方程式は y'− ![]() P(x)=− ![]() その1つはu(x)=x Q(x)=2x+1だから,∫ ![]() ![]() ![]() =2x+log|x|+C
y=(2x+log|x|+C)xになります.→2
|
||||||||
【問題4】
微分方程式y'+y=cosxの一般解を求めてください. 1y=( ![]() 2y=( ![]() 3y= ![]() 4y= ![]()
≪同次方程式の解を求めて定数変化法を使う場合≫
I=∫ excos x dxは,次のよう
に部分積分を(同じ向きに)2回行うことによりIをIで表すことができ,これを「方程式風に」解くことによって求めることができます.
同次方程式を解く:
−{−excos x+∫ excos x dx} =exsin x+excos x−I 2I=exsin x+excos x I= ![]() ![]() ![]() ∫ ![]() log|y|=−x+C1 =log e−x+C1=log(eC1e−x) |y|=eC1e−x y=±eC1e−x=C2e−x そこで,元の非同次方程式の解をy=z(x)e−xの形で求める 積の微分法により y'=z'e−x−ze−xとなるから z'e−x−ze−x+ze−x=cos x z'e−x=cos x z'=excos x z=∫ excos x dx 右の解説により z= ![]()
≪(3)または(3')の結果を使う場合≫
P(x)=1だから,u(x)=e−∫P(x)dx=e−x Q(x)=cos xだから,∫ ![]() = ![]()
y=
![]() |
○ 微分方程式の解は,y=f(x)の形のyについて解かれた形(陽関数)になるものばかりでなく,x2+y2=Cのような陰関数で表されるものもあります.もちろん,x=f(y)の形でxがyで表される場合もありえます. そうすると,場合によってはxをyの関数として解くことも考えられます.
【例題3】
微分方程式(y−x)y'=1の一般解を求めてください. ![]()
この方程式は,y'=
![]() しかし, ![]() ![]() ![]() ![]() ![]() ![]() 同次方程式: ![]() ![]() ∫ ![]() log|x|=−y+C1 |x|=e−y+C1=eC1e−y x=±eC1e−y=C2e−y 非同次方程式の解をx=z(y)e−yの形で求める 積の微分法によりx'=z'e−y−ze−yとなるから,元の微分方程式は z'e−y−ze−y+ze−y=y z'e−y=y
I=∫ yeydxは,次のよう
に部分積分で求めることができます.
両辺にeyを掛けると
z'=yey z=∫ yeydy =yey−ey+C したがって,解は x=(yey−ey+C)e−y =y−1+Ce−y |
||||
≪同次方程式の解を求めて定数変化法を使う場合≫
(y2+x) ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() と変形すると,変数yの関数xが線形方程式で表される. 同次方程式を解く: ![]() ![]() ![]() ![]() ∫ ![]() ![]() log|x|=log|y|+C1=log|y|+log eC1=log|eC1y| |x|=|eC1y| x=±eC1y=C2y そこで,元の非同次方程式(1)の解をx=z(y)yの形で求める 積の微分法により x'=z'y+zとなるから z'y+z−z=y z'y=y z'=1 z=∫ dy=y+C
≪(3)または(3')の結果を使う場合≫
P(y)=− ![]() Q(y)=yだから,∫ ![]() (u(y)=y (y>0)の場合でもu(y)=−y (y<0)の場合でも,結果は同じになります.)
x=(y+C)y=y2+Cyになります.→2
|
||||
≪同次方程式の解を求めて定数変化法を使う場合≫
(ey−x) ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() と変形すると,変数yの関数xが線形方程式で表される. 同次方程式を解く: ![]() ![]() ![]() ![]() ∫ ![]() ![]() log|x|=−log|y|+C1 log|x|+log|y|=C1 log|xy|=C1 |xy|=eC1 xy=±eC1=C2 そこで,元の非同次方程式(1)の解をx= ![]() 商の微分法により x'= ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() z'=ey z=∫ eydy=ey+C
≪(3)または(3')の結果を使う場合≫
P(y)= ![]() ![]() 1つの解はu(y)= ![]() Q(y)= ![]() ![]()
x=
![]() |
||||
【問題7】
微分方程式(x+2ylogy)y'=y (y>0)の一般解を求めてください. 1x= ![]() ![]() 3x=y(logy+C) 4x=y((logy)2+C) ヒント1ヒント2解答
≪同次方程式の解を求めて定数変化法を使う場合≫
(x+2ylogy) ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() と変形すると,変数yの関数xが線形方程式で表される. 同次方程式を解く: ![]() ![]() ![]() ![]() ∫ ![]() ![]() log|x|=log|y|+C1 log|x|=log|y|+eC1 log|x|=log|eC1y| |x|=|eC1y| x=±eC1y=C2y
∫
そこで,元の非同次方程式(1)![]() おく置換積分で計算できます. t=logy ![]() ![]() dy=y dt ∫ ![]() ![]() =∫ t dt= ![]() = ![]() の解をx=z(y)yの形で求める 積の微分法により x'=z'y+zとなるから z'y+z−z=2logy z'y=2logy z'= ![]() z=2∫ ![]() =2( ![]() =(logy)2+C
≪(3)または(3')の結果を使う場合≫
P(y)=− ![]() 1つの解はu(y)= ![]() Q(y)=2log yだから,∫ ![]() ![]() =2( ![]()
x=y(logy)2+C)になります.→4
|
1階線形微分方程式[例と解]
において
1.(2)
同次方程式の1つの解は 元の方程式の解は (検算) とおくと だから,微分方程式を満たす. |
1.(3)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
1.(4)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
1.(5)
同次方程式の1つの解は 元の方程式の解は (検算) のとき だから,微分方程式を満たす. |
1.(6)
同次方程式の1つの解は 元の方程式の解は
この積分は,部分積分2回で求められる.
(検算) のとき だから,微分方程式を満たす. |
1.(7)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
2.(1)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
2.(2)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
2.(3)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
2.(4)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
2.(5)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
2.(6)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
3.(2)
同次方程式の1つの解は 元の方程式の解は (検算) のとき だから,微分方程式を満たす. |
3.(3)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
3.(4)
同次方程式の1つの解は 元の方程式の解は (検算) のとき だから,微分方程式を満たす. |
3.(5)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
3.(6)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
4.(2)
同次方程式の1つの解は 元の方程式の解は 部分積分 のとき だから,微分方程式を満たす. |
4.(3)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
4.(5)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
5.(2)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
6.(2)
同次方程式の1つの解は 元の方程式の解は
この形の積分
(検算)とおくと だから,微分方程式を満たす. |
![]() ![]() |
■[個別の頁からの質問に対する回答][1階線形 微分方程式について/17.8.16]
特に気にすることではないと思いましたが、微分方程式の1回線形微分方程式のページにおいて、例題3の解答の積の微分によりx&apos=…… となっています。
もしかすると私のスマートフォンが原因かもしれませんが、一応アンケートに書かせていただきました。
■[個別の頁からの質問に対する回答][1階線形 微分方程式について/17.8.13]
=>[作者]:連絡ありがとう.訂正しました.
' と書くべきところを &apos と書いていたようです.(セミコロンが1つ足りない)
以前参考書で読んだのは、
exp(∫P(x)dx) ※不定積分でCを出さない
を両辺にかける事で
左辺=exp(∫P(x)dx) *y'+exp(∫P(x)dx) *y=(exp(∫P(x)dx) *y)'
右辺=exp(∫P(x)dx) *Q(x)
となるので、両辺をxで積分して最後に両辺にexp(-∫P(x)dx)を掛ければyが出てきます。
これでやっていることは茶色字で書かれている公式と中身は同じですが、
暗記事項はexp(∫P(x)dx)を両辺に掛けることだけで良く、
0からの導出は厳しいですが結果的に正しいことを示すのは容易で
変形の検算も比較的容易なのでオススメです。ご参考までに。
まぁ結局の所慣れた方法で計算するのが一番な気はしますけども…
=>[作者]:連絡ありがとう.なるほど,そのやり方でできます・・・ただし,あなたの式では左辺第2項に P(x) が抜けているようです. |
隨�ソス邵コ阮呻ソス郢ァ�オ郢ァ�、郢昜コ・�ス邵コ�ョGoogle隶諛�スエ�「隨�ソス |