無限領域での波動方程式の計算に出てくる偏微分方程式

このQ&Aのポイント
  • 偏微分方程式の解の計算方法が分からない
  • フーリエ変換を使った解法を用いればよい
  • フーリエ変換を利用して解を求める
回答を見る
  • ベストアンサー

無限領域での波動方程式の計算に出てくる偏微分方程式

波動方程式の計算に出てくる、偏微分方程式の解の計算方法が分かりません。 本から引用します: ここで、弦を伝わる波の問題などで使われる波動方程式 { (∂^2) u(x,t) } / (∂t^2) - c^2 * { (∂^2) u(x,t) } / (∂x^2) = 0 (式7.33) を考えてみよう。ここで、u(x,t)は座標xの位置での時刻tにおける弦の変位を表わし、cは正の定数とする。そして、∞に長い弦を考え(すなわち、-∞<x<∞の範囲で考え)、境界条件は、すべての t>=0 に対して u(x,t)→0 (式7.34) (x→±∞) を満たすとする。つまり、無限遠では波が存在しないとする。更に初期条件は u(x,0) = f(x) { ∂u(x,t) } / ∂t |t=0 = 0 (式7.35) とし、ここでf(x)は x→±∞ で0に近付く絶対可積分な関数であるとする。また、上式の縦棒(|)の後のt=0は、「t=0での偏微分の値」という意味である。(式7.35)のように初期条件として2つの式を与えるのは、(式7.33)がtについて2階の微分方程式だからである。今の場合、xの無限領域での関数u(x,t)を取り扱うので、フーリエ変換を使った解法を用いればよい。 例題 初期条件(式7.35)と境界条件(式7.34)を満たす(式7.33)の解を求めよ。 [解] u(x,t)のxについてのフーリエ変換を F(k,t) = ∫[-∞,∞] u(x,t) e^(-ikx) dx (式7.36) と表す。(式7.33)にe^(-ikx)を掛け、xについて-∞から∞まで積分すると、熱伝導方程式(式7.20)を導いたときと同様な考え方から、 { (∂^2)F(k,t) } / (∂t^2) + (c^2) * (k^2) * F(k,t) = 0 (式7.37) ←質問箇所 を得る。この微分方程式の解は、 F(k,t) = C[1](k) e^(ickt) + C[2](k) e^(-ickt) (式7.38) ←これをどう導いたのかが不明 であることが、代入すれば確かめられる。ここで、C[1](k)、C[2](k)は任意のkの関数で ある。 ・・・以上、引用終わり。 私は偏微分方程式自体、変数分離とかいう方法でサラッとやっただけで、上記の方法は見たことがありません。ネットで検索しましたが、同様の式を見つけることが出来ませんでした。そんな私が敢えて解こうとすると: { (∂^2)F(k,t) } / (∂t^2) + (c^2) * (k^2) * F(k,t) = 0 第2項を右辺に移項する { (∂^2)F(k,t) } / (∂t^2) = - (c^2) * (k^2) * F(k,t) 左辺の(∂t^2)と右辺のF(k,t)を交換する { (∂^2)F(k,t) } / F(k,t) = - (c^2) * (k^2) * (∂t^2) 両辺をtで積分する(もう既に未知の領域…きっと2乗が減って1乗になるのでしょう…) ln{F(k,t)} * {∂F(k,t)} / F(k,t) = - (c^2) * (k^2) * ∫(1)(∂t^2) ln{F(k,t)} * {∂F(k,t)} / F(k,t) = - (c^2) * (k^2) * t (∂t) + C[1](k) もう一度両辺をtで積分するだろう雰囲気を漂わせたところでやめておきます。 もしかしたらln{F(k,t)}を積分しなければならないのでは、と思ったら思考が停止しました。多分、既に間違っているのでしょう。 …ということで、この偏微分方程式の解き方を教えて下さい。お願いします。

質問者が選んだベストアンサー

  • ベストアンサー
  • jcpmutura
  • ベストアンサー率84% (311/366)
回答No.1

{(∂^2)F(k,t)}/(∂t^2)+(c^2)*(k^2)*F(k,t)= 0 で D=∂/∂t とすると (D^2)F(k,t)+(c^2)*(k^2)*F(k,t)= 0 {D^2+(c^2)(k^2)}F(k,t)=0 (D-ick)(D+ick)F(k,t)=0 D(D+ick)F(k,t)-ick(D+ick)F(k,t)=0 D[e^{ickt}e^{-ickt}(D+ick)F(k,t)]-ick(D+ick)F(k,t)=0 ick(D+ick)F(k,t)+e^{ickt}D[e^{-ickt}(D+ick)F(k,t)]-ick(D+ick)F(k,t)=0 e^{ickt}D[e^{-ickt}(D+ick)F(k,t)]=0 ↓両辺にe^{-ickt}をかけると D[e^{-ickt}(D+ick)F(k,t)]=0 ↓両辺をtで積分すると e^{-ickt}(D+ick)F(k,t)=C ↓両辺にe^{ickt}をかけると (D+ick)F(k,t)=Ce^{ickt} DF(k,t)+ickF(k,t)=Ce^{ickt} De^{-ickt}e^{ickt}F(k,t)+ickF(k,t)=Ce^{ickt} -ickF(k,t)+e^{-ickt}De^{ickt}F(k,t)+ickF(k,t)=Ce^{ickt} e^{-ickt}De^{ickt}F(k,t)=Ce^{ickt} ↓両辺にe^{ickt}をかけると De^{ickt}F(k,t)=Ce^{2ickt} ↓両辺をtで積分し,C/(2ick)=C[1]とすると e^{ickt}F(k,t)=C[1]e^{2ickt}+C[2] ↓両辺にe^{-ickt}をかけると F(k,t)=C[1]e^{ickt}+C[2]e^{-ickt}

futureworld
質問者

お礼

ありがとうございます。 すみません、 D[e^{ickt}e^{-ickt}(D+ick)F(k,t)]-ick(D+ick)F(k,t)=0 ↓ ick(D+ick)F(k,t)+e^{ickt}D[e^{-ickt}(D+ick)F(k,t)]-ick(D+ick)F(k,t)=0 の計算過程を教えて下さい。 e^{ickt}だけDの外に出ているので、e^{ickt}の方だけtで偏微分したのかなとは思いますが、それで+ick(D+ick)F(k,t)が副次的に出てくる計算過程が分かりません。そういう一般式があれば教えて下さい。よろしくお願いします。

futureworld
質問者

補足

ベストアンサーを差し上げます。 「積の偏微分」で理解できました。なんか部分積分を思い起こさせますね。 1 → e^{ickt}e^{-ickt}が出てくるところなどは多少天下り(下から逆算したのでは)と思えるところがありますが、計算過程が垣間見れたので私的には大満足です。 ありがとうございました。

その他の回答 (1)

  • f272
  • ベストアンサー率46% (7998/17100)
回答No.2

Fはkとtの2変数関数ですが,式7.37はtのみに関する偏微分方程式なので,tの関数のように扱えばよい。つまりF’’=-(ck)^2 Fです。 最も簡単な微分方程式F'=a*Fの解としてF=exp(at)があるということは常識的に知っているでしょう。またF''=b^2*Fの解としてF=exp(bt)とF=exp(-bt)があることも常識と言ってよいでしょう。さらに2階微分方程式の解は,1次独立な解の線形結合になることはどこかで習ったことと思います。 http://www.geisya.or.jp/~mwm48961/kou3/differ_eq2.htm 従ってF’’=-(ck)^2 Fの解は,F=exp(ickt)とF=exp(-ickt)が1次独立であることから F=C1*exp(ickt)+C2*exp(-ickt) です。この場合の積分定数であるC1とC2はtに関しては定数ですが,kに関しては関数になるので F(k,t)=C1(k)*exp(ickt)+C2(k)*exp(-ickt) と書くことができます。

futureworld
質問者

お礼

ご回答、ありがとうございます。 1次従属と1次独立は線形代数の時にやった記憶があります。なるほど、1次独立と判った時点であの一般解の式が使えるのですね。そしてC1とC2はkに関する関数にする、と。 実は、一般的にはこちらの回答の方が正しいような気がしてきましたが、No.1さんには補足もしていただいた労力も考えて、今回は見送らさせていただきます、すみません。でも大変勉強になりました。またのご回答をよろしくお願いします。 ありがとうございました。

関連するQ&A

  • 微分方程式と積分

    1.次の微分方程式を解け。 (1)y''+2y'+y=3sin2x 同次微分方程式の一般解はu(x)=(C₁+C₂x)exp(-x) と求められるのですが、非同次微分方程式の特殊解u₀(x)が求められません。 どうやって求めればいいのでしょうか。 (2)y''-5y'+6y=x(exp(x)) 非同次微分方程式の特殊解u₀(x)はどうやって求めたらいいのでしょうか。 2.置換積分によって、次の定積分を求めよ。 1.∫[0→π/2] 1/(1+cosx)dx tanx/2=tと置いた後、どうすればいいのでしょうか。 2.∫[0→a] x^2(√a^2-x^2)dx(a>0) x=asintとおくと、dx=acost dt .∫[0→a] x^2(√a^2-x^2)dx=∫[0→π/2] a^2sin^2t*acos^2t dt このあとどうすればいいのでしょうか。 お願いします。

  • 偏微分方程式(その2)

    http://oshiete1.goo.ne.jp/kotaeru.php3?q=104706 の続きです。 『(2) (**)式で次の形の解を求めよ。 φ(x, t) = Φ(ξ), ξ = x - ct ここで、cは正の定数で、Φ(ξ)のξ = 0における条件が次のように与えられている。 Φ(ξ = 0) = 1, (dΦ)/(dξ)|_ξ=0 = -c/(2ν)』 一応解けたつもりなので添削をお願いします。 ∂φ/∂t = Φ'(ξ)(∂ξ)/(∂t) = -cΦ'(ξ) ∂φ/∂x = Φ'(ξ)(∂ξ)/(∂x) = Φ'(ξ) ∂^2φ/∂x^2 = Φ''(ξ)(∂ξ)/(∂x) = Φ''(ξ) よって(**)式より Φ''(ξ) = -(c/ν) Φ'(ξ) 両辺を積分して Φ'(ξ) = -(c/ν) Φ(ξ) + d (dは定数) 初期条件Φ(0) = 1, Φ'(0) = -c/(2ν) より d = c/(2ν) よって Φ'(ξ) = -(c/ν) Φ(ξ) + c/(2ν) この非斉次方程式に対応する斉次方程式 Φ'(ξ) = -(c/ν) Φ(ξ) の解はΦ(ξ) = e^{-(c/ν)ξ}なので上非斉次方程式の解を Φ(ξ) = f e^{-(c/ν)ξ} + g (f, gは定数) とおくと f(c/ν) e^{-(c/ν)ξ} = -(c/ν){f e^{-(c/ν)ξ} + g} + c/(2ν) ∴g = 1/2 Φ(ξ) = f e^{-(c/ν)ξ} + 1/2 初期条件Φ(0) = 1よりf = 1/2、よって求める解は Φ(ξ) = (1/2)[e^{-(c/ν)ξ}] 合ってますでしょうか?

  • 微分方程式の解き方が分からず、困っています。

     現在、試験に向けて微分方程式の勉強をしているのですが、下記の問題の解き方が分かりません。  教科書を参考に(1)は変数分離系、(2)は同次形、(3)は線形で解こうとしましたが、どの問題も積分するところで複雑な式になってしまい、解けれません。  分かる問題だけでも良いのでアドバイス、解き方を教えてください。よろしくお願いします。     (1)次の微分方程式の一般解を求めよ dy/dx=y^2+1 (2)次の微分方程式の一般解を求めよ y'=(y/x)(log(y/x)+1) (3)次の微分方程式の解でt=0のときx=1の条件を満たすものを求めよ x'cost+xsint=1

  • 微分方程式について。

    微分方程式の一般解をもとめます。 (1)dy/dx=(y^2)+y これは、線形微分方程式を使ってとくのでしょうか?? (2)(x-y)y'=2y 同次形で解きましたが 途中の式、 ∫du(1-u)/(u+u^2)=∫1/xでの右辺の積分がわかりません。 両者の解答の導き方を教えてください。お願いします。

  • 偏微分方程式と常微分方程式

    物質濃度をC、時間をt、座標をx、物質の分子拡散係数をνとすると分子拡散による物質濃度の時空間変化は以下の偏微分方程式によって記述される。これについて以下の問いに答えよ。 ∂C/∂t=ν((∂^2)C/∂x^2) (1)C(x,t)=X(x)T(t)と仮定することにより、X(x)およびT(t)に関する常微分方程式をそれぞれ導出せよ。 (2)(1)での2つの常微分方程式の一般解をそれぞれ求めよ。 (3)上記拡散方程式は一般に放物型と言われる偏微分方程式に分類される。これとは別の楕円型と言われる偏微分方程式を1つ、数式で記述せよ。 困っているのは(2)の問題です。 以下のようなwebサイトを見つけました。 http://www12.plala.or.jp/ksp/mathInPhys/partial/ これに沿って問題を解いていったとき、一般解をどのようにするべきか迷いが生じました。今回の問題では初期条件や境界条件はないため、一般解はλが正、ゼロ、負のとき全ての場合の一般解を求めなければならないということですか? 後もう1点、もしよければ、楕円型の微分方程式として有名な物理現象、あるいは式を教えていただけないでしょうか? ヨロシクお願いしますm(_ _)m 特に(2)の問題に関する質問、ヨロシクお願いします。。。

  • 1・2次元の波動方程式

    ∂^2u/∂^t2=c^2∂^2u/∂x^2 を以下の境界条件の下で解け。 (1)x=0でu=0、x=Lでu=0 (2)x=0でu=0、x=Lで∂u/∂x=0 という問題をやっているのですが、この微分方程式の解き方がわかりません。1、2階の線形、非線形微分方程式は習ったのですが、この微分方程式は、左辺はtで微分していて、右辺はxで微分していて、どういうことなのかわかりません。また、これが2次元になった場合はどのようにすればいいのでしょうか?

  • 微分方程式の解き方

    自分の趣味で、{f(x)}^2-f'(x)=0 という微分方程式が解けるかどうかやってみました。 解答 (1) f(x)=0は、与えられた微分方程式を満たす。 (2) f(x)=a (aは0以外の任意の実数の定数)は与えられた微分方程式を満たさないのでf(x)≠0、f'(x)≠0とする。 {1/f(x)}^2=1/f'(x)…(A) {1/f(x)}'=-f'(x)/{f(x)}^2 より {-1/f(x)}'=1とすると、{-1/f(x)}'=f'(x)/{f(x)}^2 f'(x)/{f(x)}^2=1 1/{f(x)}^2=1/f'(x) よって(A)と同じ式になる。 なので{-1/f(x)}'=1の両辺を積分して -1/f(x)=x+C (Cは任意定数) f(x)=-1/(x+C) となる。 (1),(2)より、一般解はf(x)=-1/(x+C)、特殊解はf(x)=0である。 これでOKでしょうか? この解き方が正しいか教えていただきたいですm(__)m

  • 微分方程式について

    微分方程式について。 yやdy/dxの形ならば解けるのですが ちょっと変わった形になると解けずに困っております。 回答お願いします。 1 未知関数x(t),y(t)に関する微分方程式 x´(t)=y(t), y´(t)=-x(t)を 初期条件x(0)=a, y(0)=bの下で解け。 2 x=x(t)を変数tのC^∞級関数とする。 このとき、 d^2x/dt^2 +(dx/dt)^2 -4=0 を解け。 3 tの関数x(t)が次の微分方程式を満たすとする x´+x^2+a(t)x+b(t)=0 ただしx´=dx/dtである。 ・x(t)=u´(t)/u(t)のとき、関数u(t)の満たす微分方程式を求めよ。 ・微分方程式 x´=x(1-x)の一般解を求めよ。 長いですが回答お願いします

  • 微分方程式の解法を教えてください!

    常微分方程式の解法はどんなものがあり、どのような場合に適用すれば解けるでしょうか。 解法を覚えても、それが適用される場合についての判断ができません。教えてください! 以下の場合だとどのように解けばよいでしょうか。 (1)d^2x/dt^2+ω^2x=0の一般解の求め方。(ωは定数) (2)dx/dt=-c^2y、 dy/dt=c^2x の一般解の求め方。(cは定数) (3)dx/dt=u、    du/dt=-kx-cu+f(t) (k,cは定数)  のとき  (1)f(t)=0のとき、t=0でx=x0のもとでの解を求め   る。  (2)f(t)=cosωtのときの解。

  • 微分方程式について

    解析で習う、微分方程式についてですが、 一般解を求めますよね? そのときに、特殊解を求めますよね? その特殊解の置き方(例えば、f(x)= Ax^2+Bx+cなど)は式によって、どう置けばいいか わかりますか? 教えて下さい。