• ベストアンサー

ラプラス変換

この問題の解き方を教えて下さい。 微分方程式 x '' + 4 x = cos 2t x(0) = 0 , x(π/4)=0 <解いたやりかた> L(x(t))=X(s)とする 両辺をラプラス変換 s^2 X(s) - 0 s - 0 + 4X(s) = 2/(s^2 + 4) X(s) = 2/(s^2 +4)^2 ここで右辺の部分分数分解の分解のしかたが分かりません。 2/(s^2 +4)^2 = (As + B)/(s^2 + 4) + (Cs + D)/(s^2 + 4)^2とおいてみましたがこれだと分解できません。 よろしくお願いします。

  • tki-
  • お礼率55% (88/160)

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

  • ベストアンサー
  • info22_
  • ベストアンサー率67% (2650/3922)
回答No.6

#5です。 A#5でミスしましたので訂正します。 最後の方の所で >x(t)=(c/2)sin(2t)+L^-1{-d/ds((1/2)/(s^2+4))} >=(c/2)sin(2t)+tL^-1{(1/2)/(s^2+4)} >=(c/2)sin(2t)+(t/4)L^-1{2/(s^2+4)} >=(c/2)sin(2t)+(t/4)cos(2t) …(▲) 正:=(c/2)sin(2t)+(t/4)sin(2t) …(▲) この影響での訂正が必要になります。 >x(π/4)=0より >0=(c/2)sin(π/2)+(π/8)cos(π/2) ←× >0=c/2 ∴c=x'(0)=0  ←× 正:(c/2)sin(π/2)+(π/8)sin(π/2) 正:0=c/2+π/8 ∴c=x'(0)=-π/4 >(▲)に代入すれば > x(t)=(t/4)cos(2t)  ←× 正:x(t)=(-π/8)sin(2t)+(t/4)sin(2t) ∴x(t)=(1/8)(2t-π)sin(2t)

tki-
質問者

お礼

ありがとうございます!!

その他の回答 (5)

  • info22_
  • ベストアンサー率67% (2650/3922)
回答No.5

>L(x(t))=X(s)とする >両辺をラプラス変換 s^2 X(s) - 0 s - 0 + 4X(s) = 2/(s^2 + 4) 間違い! 正:s^2 X(s) -x(0)s -x'(0) + 4X(s) = s/(s^2 +4) x(0)=0,x'(0)=c(定数)とおくと s^2 X(s) -c + 4X(s) = s/(s^2 +4) X(s)(s^2+4)=c+s/(s^2 +4) X(s)=c/(s^2 +4)+s/(s^2 +4)^2 ...(◆) >X(s) = 2/(s^2 +4)^2 間違い。なので以下は× >ここで右辺の部分分数分解の分解のしかたが分かりません。 >2/(s^2 +4)^2 = (As + B)/(s^2 + 4) + (Cs + D)/(s^2 + 4)^2とおいてみましたが >これだと分解できません。 (◆)ですでに右辺の部分分数分解ができています。 (◆)を逆ラプラス変換すると x(t)=(c/2)sin(2t)+L^-1{-d/ds((1/2)/(s^2+4))} =(c/2)sin(2t)+tL^-1{(1/2)/(s^2+4)} =(c/2)sin(2t)+(t/4)L^-1{2/(s^2+4)} =(c/2)sin(2t)+(t/4)cos(2t) …(▲) x(π/4)=0より 0=(c/2)sin(π/2)+(π/8)cos(π/2) 0=c/2 ∴c=x'(0)=0 (▲)に代入すれば  x(t)=(t/4)cos(2t)

  • 178-tall
  • ベストアンサー率43% (762/1732)
回答No.4

1/(s^2 +4)^2 = 1/{(s + i2)(s - i2)}^2 の部分分数展開? s = -i2 の例。  1/{(s + i2)(s - i2)}^2 = A/(s + i2)^2 + R(s) 両辺に (s + i2)^2 をかけて s → -i2 、  -1/16 = A  R(s) = 1/{(s + i2)(s - i2)}^2 - A/(s + i2)^2 を求めると、1/(s + i2)^2 の項がなくなっているはず。 以下、同じ趣向で続行。   

  • Tacosan
  • ベストアンサー率23% (3656/15482)
回答No.3

あ, 突っ込み忘れた. 本質とはなんにも関係ないけど, そのラプラス変換はたぶん間違ってる.

  • Tacosan
  • ベストアンサー率23% (3656/15482)
回答No.2

第1項の存在意義が....>#1. 複素数, 使っていい?

回答No.1

右辺第2項の分子は、分母が4次だから3次にしないといけません。

関連するQ&A

  • ラプラス変換

    微分方程式 d^2y/dt^2+2c*dy/dt+y=u のt≧0における解を求めるもので、条件としては u(t)=1(t≧0),0(t<0) y(0)=y'(0)=0 0<c<1 をという問題を解こうをしているのですがどうしても途中で止まってしまいます。 自分では・・・ t≧0における解なのでu=1. 両辺をラプラス変換して (左辺)=s^2L[y(t)]-sy(0)-y'(0)+2c(sL[y(t)]-y(0))+L[y(t)] (右辺)=1/s なので y(t)=L^-1[1/(s^3-2cs^2+s)] というところまで行ったのですが、このラプラス逆変換の仕方がよくわかりません。 どなたかご教授お願いします。

  • ラプラス変換

    下の微分方程式を[]内の初期値で、ラプラス変換を用いて解きたいのですが d^2x(t)/dt^2 + 4x(t) = 4 [ x(0+)=0、dx(0+)/dt = 0 ] 両辺をラプラス変換 s^2x(s) - sx(0+) -x'(0+) + 4x(s) = 4/s x(s) = 4/s(s^2 + 4) となりました。 この先どのようになりますか? 解答と解説をして頂けると助かります

  • 『たたみこみ』の逆ラプラス変換

    ”たたみこみ”で逆ラプラス変換の問題を解くものなのですが、 いまいち”たたみこみ”の活用法がわかりません。  S^2/{(S^2+4)^2} という問題で、これを部分分数分解して逆ラプラス変換すると、  (1/2)t・cos2t+(1/4)t・sin2t となる筈(苦)なのですが。 どうも問題を”たたみこみ”で解くことが出来ないのです。  L[cos2t]=S/(S^2+4) という関係式を使うのか、と感じてはいるのですが、そこで止まってしまいます。 ”たたみこみ”について熟知(?)していらっしゃる方々、御回答お願いします。

  • ラプラス変換の問題 教えてください 

    ラプラス変換の問題について分からない問題があります。 次の積分方程式を満たす関数x(t)を求めよ  x(t) = 4∫[0→t] x(τ)×cos2(t-τ)dτ +1 上の問題が分かりません。 最初は右式の∫内を畳み込みを使って展開して、   x(t) = 4{x(t) * cos2t} +1 となりますよね? 次に両辺をラプラス変換して、   X(s) =4{X(s) × (s/s^2 + 4)} + (1/s) 右辺のX(s) を左辺に持ってきて、  X(s) - X(s) × (4s/s^2 +4) = 1/s {1 - (4s/s^2 +4)}X(s) = 1/s {(s^2 -4s +4) / s^2 +4}X(s) = 1/s {(s-2)^2 / (s^2 +4)}X(s) = 1/s X(s) = {(s^2 +4) / (s-2)^2} ×(1/s) というところまでは来ました。 しかしこの後からどうすればいいのか分かりません。 計算しても答えに合わず、もしかしたらどこかで間違っているのかもしれませんが、 どこで間違っているのかもよくわかりません。 お願いいたします。

  • ラプラス変換で微分方程式を解く

    ラプラス変換で微分方程式を解く u"+4u'+4u=e^(-3t) u(0)=0 u'(0)=1 をラプラス変換により微分方程式で解きたいのですが、部分分数分解がうまくいかないため、解けません。 本日テストですので、途中式、答えをお願いします。

  • ラプラス変換で微分方程式を解く

    ラプラス変換で微分方程式を解く u"+2u'+2u=e^(-2t) u(0)=2 u'(0)=2 をラプラス変換により微分方程式で解きたいのですが、部分分数分解がうまくいかないため、解けません。 本日テストですので、途中式、答えをお願いします。

  • ラプラス変換で微分方程式を解く

    ラプラス変換で微分方程式を解く u"+2u'+2u=e^(-2t) u(0)=2 u'(0)=2 をラプラス変換により微分方程式で解きたいのですが、部分分数分解がうまくいかないため、解けません。 本日テストですので、途中式、答えをお願いします。

  • ラプラス変換

    下式を[]内の初期値でラプラス変換すると d^2x(t)/dt^2 - 2dx(t)/dt - 3x(t) = 2 [ x(0+)=1,dx(0+)/dt =0 ] 両辺をラプラス変換 s^2X(s) - sX(0+) - x'(0+) - 2X(s)- 2x(0+) - 3X(s) =2/s となりました初期値を入れてまとめると X(s)はどうなりますか? その後は部分分数展開ですか?

  • ラプラス変換、教えてください

    ラプラス変換を使って、方程式を解く問題です。 (1) ∫cos2(t-τ)y(τ)dτ=e^(2t)sin2t  [0、t] (2) y'+y-5∫e^{-(t-τ)}cos2(t-τ)y(τ)dτ=e^(-t)   [0、t] {y(0)=0} を与えられた条件のなかで解け。 まず(1)は、右辺をラプラス変換をして 2/{(s-2)^2+4} しかし左辺をどのように変換すればよいのか分かりません。畳み込みを使おうと思ったのですが、だめでした。 (2)はラプラス変換をして、   {sY(s)-y(0)}+Y(s)-5~~~=1/(s+1) ∫e^{-(t-τ)}cos2(t-τ)y(τ)dτ の部分のラプラス変換が分かりません。 大まかには分かるのですが、どのような流れで解に持っていけばいいのかも曖昧になっています。 ちなみに答は (1) y(t)=1+e^2t (cos2t+3sin2t) (2) y(t)=(5/2)-(5/2)e^(-2t)-4te^(-t) になります。 どなたか教えてください。 よろしくお願いします。

  • ラプラス変換の初期値問題

    y'+0.2y=0.01t・・・・(1) y(0)=-0.25  解 y=0.05t-0.25 をラプラス変換を使って解くのですが、 L(y')=sL(y)-y(0) Y=L(y)として(1)を変換して sY-y(0)+0.2Y=0.01/s^2・・・(2) Y(s+0.2)=0.01/s^2-0.25・・・・(3) ←Yの係数をまとめる Y=0.01/s^2(s+0.2)-0.25/(s+0.2)・・(4)←両辺に1/(s+0.2)をかける この時点でまず、0.25/(s+0.2)は逆変換した時、0.25exp^(-0.2) になってしまい、この後計算していくと、 解の0.05t-0.25と異なった答えがでてきてしまいます。 どこで間違ってしまったのかわからず困っています。 どなたか教えてください。 それから、文章が長くなって申し訳ないのですが、 初期値問題を解いていると、 部分分数分解や、 部分分数分解を解くための3つ以上の連立方程式がでてくることで、 計算ミスの増加や、 式の整理が不十分な為、遠回りな解答をしてしまったりと一問 にとても時間がかかってしまいます。 部分分数分解のスマートなやり方はないのでしょうか。