• ベストアンサー

ラプラス変換・微分方程式の一般解

x'' + x = 2sin2t の一般解が求められません><;解いていくうちに、4/{(s^2+4)(s^2+1)}という項が出来てしまい、そこでストップしてしまいます。ご教授ください。

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

  • ベストアンサー
  • 151A48
  • ベストアンサー率48% (144/295)
回答No.2

♯1です x(0)=a , x'(0)=b, L(x(t))=X(s)とする。与式をラプラス変換して s^2 X(s)-as-b+X(s)=4/(s^2 +4) 整理して X(s)=4/(s^2+1)(s^2+4) +as/(s^2+1) +b/(s^2 +1) =(4/3){1/(s^2 +1)}-(2/3){2/(s^2+4)}+as/(s^2+1) +b/(s^2+1) よって x(t)=(4/3)sint -(2/3)sin2t +acost+bsint

evlm4a1
質問者

お礼

できました!分かりやすい解答をどうも有難うございました!

その他の回答 (1)

  • 151A48
  • ベストアンサー率48% (144/295)
回答No.1

4/(s^2 +4)(s^2+1) =(4/3){1/(s^2+1) -1/(s^2+4)} =(4/3){1/(s^2+1) -(1/2)・2/(s^2+4)} =(4/3){1/(s^2+1)} -(2/3){2/(s^2+4)} ではどうでしょう。 初期条件x(0), x'(0)の値はないのですか?

evlm4a1
質問者

補足

ご回答どうもありがとうございます。初期条件は一般解なので無いようなのです… 何卒よろしくお願いします…

関連するQ&A

  • 偏微分方程式の一般解

    偏微分方程式:∂u(x,t)/∂x=A (Aは定数) の解き方と一般解を教えてください。

  • ラプラス変換を使った微分方程式です

    dx/dt + 2y = cos(t) , x(0)=1 x - dy/dt = sin(t) , y(0)=-√2 を解きなさい。なお、x(t),y(t)にはiは含まれない形にすること。 という問題なのですが、両辺をラプラス変換して、X(s),Y(s)を出せてもそこからがうまく逆変換できません。どなたか教えてください。

  • 1階線形偏微分方程式の一般解

    数学のことでちょっと皆様のお知恵を拝借いたしたく質問します。 次の偏微分方程式の一般解の求め方を教えてください。 ∂T(x,t)/∂t + (q(t)/S)(∂T(x,t)/∂x) = c(T_w(x,t) - T(x,t)) c,S:定数 僕の所有する参考書によるとこの種の方程式は ラグランジュの偏微分方程式と呼ばれていて、 ちょっとだけ一般解の求め方が書いてありました。 しかし、どうしても一般解にたどりつけません。 その方法とは、偏微分方程式 P(x,y,z)(∂z/∂x) + Q(x,y,z)(∂z/∂y) = R(x,y,z) に対して連立補助方程式 dx/P = dy/Q = dz/R を解いた解を f(x,y,z) = a, g(x,y,z) = b (a,bは積分定数) とする。φを任意の関数として、一般解は  φ(f,g) = 0 である。 という解法です。しかし、T_wが邪魔でうまくいかないです。 詳しい参考書を手に入れようにも近くに本屋がないのでお手上げです。 どなたかご教授お願いしますm( _ _ )m

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

    d^2x(t)/dt^2 + 4dx(t)/dt + 3x(t) = 12 この微分方程式の一般解は何になるのでしょうか? また3x(t)が4x(t)だとどうなりますか? よろしければ教えてください。

  • 微分方程式の一般解

    微分方程式の一般解を求める問題なのですが、どうしてもよく分かりません。 y'6+4y'2=40x^3  (ここで'○は微分の回数を示すとします。また、以下ではD=d/dxのことです) 同時方程式(D^6+4D^2)y=D^2(D^4+4)y=D^2{(D^2+2)^2-4D^2}y=D^2(D^2+2+2D)(D^2+2-2D)y=D^2{(D+1)^2+1}{(D-1)^2+1}y=0 の基本解は{1、x、e^(-x)cosx、e^(-x)sinx、e^xcosx、e^xsinx} 次に特殊解Y(x)を求める。 非同次項R(x)=40x^3は同時方程式40D^4y=0の解だから、 特殊解Y(x)の式Yは同時方程式 40D^4・D^2{(D+1)^2+1}{(D-1)^2+1}y=0の解である。 この基本解{1、x、x^2、x^3、x^4、x^5、e^(-x)cosx、e^(-x)sinx、e^xcosx、e^xsinx}から与式の基本解を除いたx^2、x^3、x^4、x^5の一次結合として Y(x)=Ax^2+Bx^+Cx^4+Dx^5とおく。 与式の左辺に代入… と続いていくのですが、どうにもしっくりきません。 答えも y=(20/3)x^3+c1+c2x+e^(x/√2){c3cos(x/√2)+c4sin(x/√2)}+e^(-x/√2){c5cos(x/√2)+c6sin(x/√2)} となり、私の解からでは到底結びつくとは思えないです。 気になるのが 「非同次項R(x)=40x^3は同時方程式40D^4y=0の解だから」 としていますが、本当にこれで良いのか自信もありません。 もし間違えていたら解説をお願いします。 また、他に違うというようなところがあったら指摘してください。 回答、よろしくお願いします。

  • 偏微分方程式のラプラス変換による解法

    皆様よろしくお願いいたします。 関数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に従わず∞に発散してしまいます。

  • 微分方程式

    次の微分方程式を解けという問題がわかりません。 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 となってしまって解けませんでした。どこを直せばいいでしょうか?

  • ラプラス変換で微分方程式の一般解を求めるには限界がある?

    ラプラス変換を覚えて、微分方程式を簡単に解いてしまおうと思い勉強していたのですが、 y' = (1+x)y という問題において、 y(0) = a , L[y(t)] = Y(s) , L[y'(t)] = sY(s) - y(0) とし、与式の両辺のラプラス変換を取って sY(s) - a = Y(s) - Y'(s) <-像関数の微分法則より となると思います。このY'(s)の処理の仕方が分かりません。 答えは y = Cexp(x+x^2/2) (Cは定数) らしいのですが、これはラプラス変換では難しいのでしょうか。 (s-1)Y(s)が出てくるのでexp(x)は納得できるのですが、何故xを積分したと思われる値がexp()内に出るのか分かりません。

  • 偏微分方程式の解について。

    現在、私は3変数(x、y、z)2階の偏微分方程式を解いています。 その同次解を導いています。 まず、変数の一般解をΣX(r)*(cosmθ)、ΣY(r)*(cosmθ)、ΣZ(r)*(cosmθ)と仮定し元の式に代入したのち、r=exp(s)と変数変換します。 そして同次解の形をX=X'exp(λs),Y=Y'exp(λs),Z=Z'exp(λs)のように仮定し代入することによって、自明でない解をもつ次の特性方程式を得ました。 p^3+d*p+f=0 このときp=(λ^2-A)とします。 またAとdとfは定数です。 ここから解を導くのですが λ^2=p+A>0のときは、 X=F*exp(λs)+S*exp(λs)  =F*r^λ+S*r^(-λ) このときのF,Sは勝手においた未知数です。 とまずおきました。 次にXを既知だと仮定し、YとZの関係を求めるのですが、 関数型はXと同様のために、F=1として 同次解を仮定して代入した式で計算してYとZの関係を導きました。 (簡単な2次方程式を解く作業です) 同様にS=1としても行いました。 そこで以下の解を得ました。 Y=G(λ)*F*r^λ+G(-λ)*S*r^(-λ) Z=H(λ)*F*r^λ+H(-λ)*S*r^(-λ) G(λ)とH(λ)は2次方程式を解いて出した関係式です。 次がわからないところです。 λ^2=p+A<0の場合、つまりλの根が複素数の場合です。 上と同様に係数を比較して求めるのですが、 X=F*cos(λs)+S*sin(λs) と仮定するところまではわかりますが、 その仮定によって Y={Re[G(j*λ)]cos(λs)-Im[G(j*λ)]sin(λs)}*F +{Im[G(j*λ)]cos(λs)+Re[G(j*λ)]sin(λs)}*S となるのがわかりません。Zについても式の形は同様です。 本当に困っています。 意味がわからない文章かもしれませんが、汲み取っていただけると幸いです。 ヒントでもいいのでください。 ちなみに 実部については G(j*λ)=G(j*-λ)が成り立ち      虚数部については G(j*λ)=-G(j*-λ)が成り立っております。

  • 連立一次方程式の一般解

     現在、大学で線形代数を学んでいる者です。     問: 特殊解を見出すことによって、次の連立一次方程式の一般解を求めよ。       ⅹ1+4ⅹ2-ⅹ3+ⅹ4=1     2ⅹ1+5ⅹ2-2ⅹ3+ⅹ4=1              -ⅹ1+2ⅹ2+ⅹ3+ⅹ4=1      3ⅹ1+9ⅹ2-3ⅹ3+2ⅹ4=2      *以下の()内は列ベクトルを示す     (解答)ⅹ4の係数を見れば、特殊解ⅴ0=(0000)にとれることが見える。     以下の計算は省略      :      :   求める一般解は   (ⅹ1ⅹ2ⅹ3ⅹ4)=s(1010)+t(011-3)+ⅴ0           =s(1010)+t(011-3)+(0000)     (s、tは任意のKの元) というように教科書には書かれているのですが、なぜ、特殊解ⅴ0=(0000)と決められるのかがわかりません。わかる方、ぜひ、理由を教えてください。