• ベストアンサー

ラプラス変換です

関数x(t)が  tx''+x'+x=0 x(0)=1、 x'(0)=-1 を満たすものとする。 与えられた微分方程式と初期条件を満足するX(s)を求めたいです。 よろしくお願いします。

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

  • ベストアンサー
  • nubou
  • ベストアンサー率22% (116/506)
回答No.4

X(s)≡L[x(t)](s)≡∫(0<t<∞)dt・x(t)・exp(-s・t)とする L[x’(t)](s)=s・X(s)-x(0)=s・X(s)-1 L[x”(t)](s)=s^2・X(s)-s・x(0)-x’(0)=s^2・X(s)-s+1 X(s)の定義式をsで微分すると X’(s)=-∫(0<t<∞)dt・t・x(t)・exp(-s・t) すなわち L[t・x(t)](s)=-X’(s) とくにx(t)としてx”(t)を選べば L[t・x”(t)](s)=-L[x”(t)]’(s)=-(s^2・X(s)-s+1)’ =-s^2・X’(s)-2・s・X(s)+1 t・x”(t)+x’(t)+x(t)=0 の両辺をラプラス変換すると -s^2・X’(s)-2・s・X(s)+1+s・X(s)-1+X(s)=0 すなわち S^2・X’(s)+(s-1)・X(s)=0 これを解いて X(s)=C・exp(-1/s)/s ラプラス逆変換すればC=1が分かる 従って X(s)=exp(-1/s)/s

noname#4075
質問者

お礼

前回の常微分のときもお世話になりまして どうもありがとうございます。 また何かあったらよろしくお願いします。

全文を見る
すると、全ての回答が全文表示されます。

その他の回答 (4)

  • nubou
  • ベストアンサー率22% (116/506)
回答No.5

わざわざラプラス逆変換するより初期値定理を使った方がスマートですね 1=x(0)=lim(s→∞)・s・X(s)=lim(s→∞)・C・exp(-1/s)=C

全文を見る
すると、全ての回答が全文表示されます。
  • nubou
  • ベストアンサー率22% (116/506)
回答No.3

x(t)=Σ(0≦n<∞)・C[n]・t^nとする x’(t)=Σ(0≦n<∞)・C[n+1]・(n+1)・t^n t・x”(t)=Σ(1≦n<∞)・C[n+1]・(n+1)・n・t^n x(0)=1よりC[0]=1 x’(0)=-1よりC[1]=-1 従って t・x”(t)+x’(t)+x(t)= Σ(1≦n<∞)・(C[n+1]・(n+1)^2+C[n])・t^n 従って1≦nにおいて C[n+1]・(n+1)^2+C[n]=0 従って0≦nにおいて C[n]=(-1)^n/(n!)^2 従って x(t)=Σ(0≦n<∞)・(-1)^n/(n!)^2・t^n x(t)をラプラス変換して X(s)=Σ(0≦n<∞)・(-1)^n/(n!)/s^(n+1) X(s)=exp(-1/s)/s

全文を見る
すると、全ての回答が全文表示されます。
  • nubou
  • ベストアンサー率22% (116/506)
回答No.2

X(s)=C・s・exp(-1/s)-2・s+2 においてx(t)がまともな関数になるにはC=2であるから X(s)=2・s・exp(-1/s)-2・s+2 ラプラス逆変換して x(t)=2・Σ(0≦k<∞)・(-1)^k/(k+1)!/k!・t^k しかし x(0)=1となるがx’(0)=-1/3であるので 結局どっかで間違えたみたいですね

全文を見る
すると、全ての回答が全文表示されます。
  • nubou
  • ベストアンサー率22% (116/506)
回答No.1

X(s)≡L[x(t)](s)≡∫(0<t<∞)dt・x(t)・exp(-s・t)とする X(s)の定義式をsで微分すると X’(s)=-∫(0<t<∞)dt・t・x(t)・exp(-s・t) すなわち L[t・x(t)](s)=-X’(s) L[x’(t)](s)=s・X(s)-x(0)=s・X(s)-1 L[(t・x(t))”] =s^2・L[t・x(t)](s)-s・(0・x(0))-[(t・x(t))’](t=0) =s^2・L[t・x(t)](s)-x(0) =-s^2・X’(s)-1 t・x”(t)+x’(t)+x(t)=0 の両辺をラプラス変換すると -s^2・X’(s)-1+s・X(s)-1+X(s)=0 すなわち S^2・X’(s)=(s+1)・X(s)-2 これを解いて X(s)=C・s・exp(-1/s)-2・s+2 これをラプラス逆変換して初期条件を使えばいいような気がしますが変な式になったのでどっかで間違えたかな?

全文を見る
すると、全ての回答が全文表示されます。

関連するQ&A

  • ラプラス変換

    下の微分方程式を[]内の初期値で、ラプラス変換を用いて解きたいのですが 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) となりました。 この先どのようになりますか? 解答と解説をして頂けると助かります

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

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

  • ラプラス変換についてです。

    ラプラス変換の質問です。 L[x(t)] = X(s) の導関数は, a) L[-tx(t)] = d X(s), ds [ n ] dn b) L (-t) x(t) = dsn X(s) を満たす. の性質をつかって、 t^2×e^λt はどうやって解くのでしょうか?

  • 偏微分方程式の座標変換について

    次の問題がわかりません。 関数u(x,t)は次の偏微分方程式 ∂u/∂t = -2∂u/∂x - u - 2x を満足するものとする。 (1)ξ=x-2t , η=t なる座標変換を考える。関数uが満足するξとηに関する偏微分方程式を求めよ。 (2)v(ξ,η)=u*e^ηとおくとき、(1)の結果を用いて、関数v(ξ,η)の一般解を求めよ。 (3)境界条件u(0,t)=e^(-2t)+4を満足する関数u(x,t)を求めよ。 という問題です。わかる方がいたら回答よろしくお願いします。

  • ラプラス変換

    この問題の解き方を教えて下さい。 微分方程式 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とおいてみましたがこれだと分解できません。 よろしくお願いします。

  • ラプラス変換

    ラプラス変換の初期値問題なのですが、 tx'+x=t のとき方を教えていただけないでしょうか。 初期値は適当に設定してくれてかまいません。 ラプラス変換して (s+1)X'+2X=0 になるとこまでは分かりましたが、そのあとが・・ 勝手な質問かもしれませんが、よろしくお願いします。

  • 大学院入試過去問題(ラプラス変換・微分方程式)

    下記の問題の解法がわかりません。 関数f(x)のラプラス変換F(s)を次のように定義する。 F(s)=∫f(t)*exp(-st)dt   積分区間は0から∞ 以下の問いに答えよ。 (1)関数g(t)が次のように与えられるものとする。 g(t)=(2/π)*∫cos(t*tanz)dz   積分区間は0から(π/2) g(t)のラプラス変換G(s)を求めよ。 (2)関数f(x)が次の方程式 f''+2*f'+3f=g(t)    (f'',f'はtによる微分です) および初期条件 f(0)=1,f'(0)=-2 を満足するものとする。ラプラス変換F(s)を求め、その逆変換f(t)を求めよ。 (1)からさっぱり分かりません。 どなたかよろしくお願いします。

  • ラプラス変換がわかりません

    微分方程式をラプラス変換を用いてとけ 9d^2x/dt^2+x=1 初期値はx(0)=0. x'(0)=0 わかりやすく教えてください

  • ラプラス変換の「s」とは?

    皆さんこんにちは。 解らないことがありましたので質問させていただきます。 今、学校で「ラプラス変換」を使って微分方程式を解いています。 「ラプラス変換は、L[f(t)]=F(s)という操作のことをいって、 微分方程式を解く際には最終的にsの関数を逆変換してtの関数に戻すから、 sとは何か?については何も気にしないで計算を進めて欲しい。」 と、先生から説明されたのですが、どうしても自分の性格上納得がいかないのです。 この「s」って一体何なのでしょう?ネットを使ってかなり調べたのですが、 詳しい説明をされている場所に辿り着くことはできませんでした。 「複素数」だということは解ったのですが・・・・・・ 明確な意味を教えていただきたいな、と思います。 ちなみに高専4年生です。 あまり難しすぎる説明はもしかしたら理解できないかもしれません。申し訳ありません;; それでは、ご教示よろしくお願い致します。

  • ラプラス変換

    微分方程式をラプラス変換使いとけ dx/dt+3x=sin2t 解き方が全然わかりません

このQ&Aのポイント
  • スギ薬局でアルバイトをしている者です。マスクの着用について問題があります。
  • マスクをつけて業務をするのが困難な状況にあります。
  • マスクの着用を任意とすることで、快適な業務環境に改善を求めたいと考えています。
回答を見る