- ベストアンサー
微分方程式の問題です。x''=a*x*cos(ωt)
微分方程式の問題です。x''=a*x*cos(ωt) 初めて質問させて頂きます。 a、ωは定数、x=x(t)として、2階の微分方程式x''=a*x*cos(ωt)が解けません。 初期条件をx(0)=x0, x'(0)=vx0としたとき、これには解析的な解はあるのでしょうか? 私なりに考えたのは、 x''=x*exp(iωt) と拡張、s=exp(iωt)と変数変換 (x=x(s))して、 s*x''+x'+(a/ω^2)*x=0 s*x''+x'+k*x=0 (k=(a/w^2)と置く) とし、最後は級数展開という方法...ですが自信がありません。 元の方程式が(見た目は)シンプルですので、もしかしたら解析的な解をご存知の方が いらっしゃるかもしれないと思い質問させて頂いた次第です。 私は物理屋ですので、なるべく純粋数学的な表現は避けて頂けると幸いです。 よろしくお願い致します。
- a0026175
- お礼率62% (5/8)
- 物理学
- 回答数2
- ありがとう数3
- みんなの回答 (2)
- 専門家の回答
質問者が選んだベストアンサー
あまり参考にならないかも知れませんが・・・、 x''=a*x*cos(ωt) Mathieuの微分方程式が連想されます。 Mathieuの微分方程式、Mathieu関数辺りで調べられてみては如何でしょうか?
その他の回答 (1)
なるほど、確かに http://www30.wolframalpha.com/input/ に、 d^2x/dt^2=a*x*cos(ω*t) を入力して実行させると Mathieu関数 とやらが出てきますね。
お礼
ご回答ありがとうございます。 http://www30.wolframalpha.com/input/ という便利なサイトがあるのですね。確かに解はMathieu関数で 表されるようですね。ぜひブックマークしておきます。 ありがとうございました。
関連するQ&A
- 微分方程式の問題です
x(t)''+2ax(t)'+(1+a^2)x=sintを考える。 (1)同次方程式x(t)''+2ax(t)'+(1+a^2)x=0の一般解を求めよ。 (2)同次方程式x(t)''+2ax(t)'+(1+a^2)x=0に対して、初期条件x(0)=0,x(0)'=1を与えた時の特解を求めよ。 (3)a=0の時、非同次方程式x(t)''+2ax(t)'+(1+a^2)x=sin(t)の特解を求めよ。 以上です。 自分で解いたのですが、合っていますでしょうか?また、はずれていれば解説付きで解き方お願いいたします。ただし、(3)については解けませんでした。 (1) x(t)=C1exp{(-a+i)t}+C2exp{-(a+i)t} (2) x(t)=(-i/2)exp{(-a+i)t}+(i/2)exp{-(a+i)t} (3) a=0なので、x(t)''+x=sin(x)となって、未定係数法を使って解いた方がよろしいのでしょうか? その場合、x(t)=x(sin(x)+cox(x))とおいて、微分方程式に代入して解くのですか?
- ベストアンサー
- 数学・算数
- 偏微分方程式のラプラス変換による解法
皆様よろしくお願いいたします。 関数u(x,t)のtに関する偏微分∂u/∂t=u_t、とxに関する2回偏微分∂^2 u/∂x^2=u_xxとおくとき 偏微分方程式 u_t = a*u_xx (aは正の定数) 初期条件:u(x,0) = 0 境界条件:∂u/∂x = u_x = -k (kは正の定数) lim[x→∞]u(x,0) = 0 をラプラス変換して解を求めようとしてますが、ラプラス変換した式が導けません。 偏微分方程式の解は分かっていているので、解をラプラス変換すると答えは次式になるようです。 U(s,x) = k√a・exp( -x*√(s/a) ) / s^(3/2) どのように導けばこうなるのかご教示ください。 ちなみに偏微分方程式の解は次式になります。(上式に入れて成り立つことを確認済み) u(x,t)=2k√(at/π)・exp(-x^2/(4at)) - kx・erfc(x/√(4at)) (※erfcはガウスの余誤差関数です) 【途中までやってみた計算経過】 偏微分方程式を→s、x→yへそれぞれラプラス変換して整理すると U(s,y)=ak/{y(y^2-s/a)} となりました。これをy→xへラプラス逆変換すると U(s,x) = -ka^2/s + ( ka^2/(2s) ) exp(-x√(s/a) ) + ( ka^2/(2s) )exp(x√(s/a) ) となり、答えになりません。 しかもこれだと3項目が境界条件lim[x→∞]u(x,0) = 0に従わず∞に発散してしまいます。
- 締切済み
- 数学・算数
- 減衰振動の微分方程式の解
先生から配られたプリントには減衰振動の微分方程式が「m(dx/dt)^2+2γ・dx/dt+ω^2x=0」の時、解が「x=A・EXP(-γt)cos(ω´t+φ)」って書かれてます。 摩擦:гならг/m=2γ、バネ定数:kならk/m=ω^2、A=√C1^2+C2^2、φ=C1/C2、ω´^2=ω^2-γ^2です。 解の式で、cosじゃなくてsinではないのですか?単振動・強制振動の場合も同様にcosでした。 誰かよろしくお願いいたします。
- ベストアンサー
- 物理学
- 微分方程式の解法。
現在、私は微分方程式が解けなくて困っています。 その微分方程式は次のようになります。 (d^2/dr^2)T+(1/r)(d/dr)T=(1/K)(d/dt)T をラプラス変換した、 T''+(1/r)*T'-(s/K)*T=0 です。 式のsはラプラス演算子で、Kは定数です。 この式の解法を調べたところ、上のような微分方程式はベッセルの変形微分方程式というものであることがわかり、一般解を導出し、計算したのですが、ラプラス逆変換が困難で挫折しました。 なにか他の解法はありませんか? 今、考えているのが解を次のように仮定し、 T=A*exp(-rs)+B*exp(-rs) 上の式に代入し、境界条件によってAとBを決定する方法です。 この方法はまずいですか? 困っているので回答お願いいたします。
- ベストアンサー
- 数学・算数
- 常微分方程式の問題です
二階の常微分方程式の問題です。 P(x)y'' + Q(x)y' +R(x)y =0について P(x) + Q(x) + R(x)=0 のとき 解はy=exp(x)となる事を示せという問題です 必要条件は代入すればよいので簡単ですが 問題は十分条件です。この場合、各yの係数は定数で 無いので、特性方程式はつかえず、 (P,Q,Rを定数とすれば確かにy=exp(x)にはなりますが) またべき級数でとこうとしても、変数が3つとなりとけず、 y=exp(x)という結論を導けません どうかよろしくお願いします。
- ベストアンサー
- 数学・算数
- f(x)が微分方程式を満たすかどうか・・・
f(x)=∫_0^∞(exp(-t^2)cos(xt))dtを定義する y=f(x)が微分方程式 2(dy/dx)+xy=0 の解を満たすことを示せ。 という問題なのですが、これはまずf(x)の式の∫を外すのでしょうか? そうでなければどうやってf(x)の中にあるxだけを取り出すのでしょうか・・・?
- ベストアンサー
- 数学・算数
- 微分方程式と積分
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 このあとどうすればいいのでしょうか。 お願いします。
- ベストアンサー
- 数学・算数
- 微分方程式
次の微分方程式を解けという問題がわかりません。 y''+4y=sin2x 特性方程式s^2+4=0よりs=±2i(虚数解) 補助方程式の一般解はy=Asin2x+Bcos2x 与方程式の右辺を微分して生じる関数は2sin2x,2cos2xであるが、 これらは上の一般解に含まれている。重複度は2なので、 特殊解を求めるために、 y1=ax^2*sin2xとおく y1'=2a(xsin2x+x^2cos2x) y1''=2a(sin2x+4xcos2x-2x^2sin2x) これらを与方程式に代入すると 2asin2x+8axcos2x-4ax^2sin2x+4ax^2sin2x=sin2x となってしまって解けませんでした。どこを直せばいいでしょうか?
- 締切済み
- 数学・算数
お礼
ご回答ありがとうございます。 実は他の質問版でも質問させて頂いており、そちらでもやはり 「Mahieuの微分方程式では?」 との回答を頂きました。私はそれを存じ上げませんでしたので、 昨日図書にて調べて無事、解の挙動(発散するのかしないのか)を 知ることができました。 ありがとうございました。