• ベストアンサー

ラプラス変換法について(微分方程式)

nubouの回答

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

重大なミスが有りましたので報告します (h(t)・g(t))’=δ(t)・g(t)+h(t)・g’(t) はいいのですが (h(t)・g(t))”=(δ(t)・g(t)+h(t)・g’(t))’ =δ’(t)・g(t)+2・δ(t)・g’(t)+h(t)・g”(t) は間違いです 正しくは (h(t)・g(t))”=(δ(t)・g(0)+h(t)・g’(t))’ =δ’(t)・g(0)+δ(t)・g’(t)+h(t)・g”(t) =δ’(t)・g(0)+δ(t)・g’(0)+h(t)・g”(t) です つまり δ(t)・g(t)=δ(t)・g(0) ですから (δ(t)・g(t))’=δ’(t)・g(t)+δ(t)・g’(t) ではなく (δ(t)・g(t))’=(δ(t)・g(0))’=δ’(t)・g(0) なのです 以上を元にNo.1を修正します 微分方程式:x”(t)+a・x’(t)+b・x(t)=g(x)・・・(1) 初期条件:x(0)=α,x’(0)=β を解いてみましょう h(t)・x(t)の両側ラプラス変換をL(h(t)・x(t))(s)≡X(s)とし h(t)・g(t)の両側ラプラス変換をL(h(t)・g(t))(s)≡G(s)とする (h(t)・x(t))’=δ(t)・x(0)+h(t)・x’(t) であるから s・L(h(t)・x(t))=x(0)+L(h(t)・x’(t)) ∴L(h(t)・x’(t))=s・X(s)-α (h(t)・x(t))”=(δ(t)・x(0)+h(t)・x’(t))’ =δ’(t)・x(0)+δ(t)・x’(0)+h(t)・x”(t) であるから s^2・L(h(t)・x(t))(s)=s・x(0)+x’(0)+L(h(t)・x”(t))(s) ∴L(h(t)・x”(t))(s)=s^2・X(s)-s・α-β (1)式の両辺にh(t)をかけて両側ラプラス変換をすると L(h(t)・x”(t))(s)+a・L(h(t)・x’(t))(s)+b・L(h(t)・x(t))(s) =L(h(t)・g(t))(s) すなわち (s^2+a・s+b)・X(s)=s・α+β+a・α+G(s) ∴X(s)≡ (s・α+β+a・α)/(s^2+a・s+b)+G(s)/(s^2+a・s+b) このあと部分分数展開と両側ラプラス変換表からx(t)を求める

関連するQ&A

  • ラプラス変換による微分方程式について

    ラプラス変換による微分方程式について u'+3u=cost u(0)=-1 について微分方程式をラプラス変換で解きたいのですが、答えがわかりません。部分分数分解でつまずいています。 明日テストの為、途中式、答えをお願いします。

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

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

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

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

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

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

  • ラプラス変換・微分方程式

    微分方程式についてふと思った疑問です。 y'' + (w^2)y' = f(t), y(0)=a, y'(0)=v という微分方程式で、両辺ラプラス変換してやって、 (s^2)Y(s) - sa - v + (w^2)Y(s) = F(s) ・・・ と計算を進めていきますが、 y(0) = a, y'(0) = v という初期条件が、仮に y(1) = a, y'(1) = v だとしたらどのように計算すればよいのでしょうか。 今のところはラプラス変換を形式的にしか理解していないので、もしこうなったらどうするんだろうと疑問に思ってしまいました。 (手持ちの参考書では、全てy(0)=..., y'(0)=...,という条件になっています)

  • ラプラス変換で連立微分方程式を解くとき

    お願いします。 連立微分方程式をラプラス変換で解くとき、 たとえばx'をラプラス変換すると sL(x) - x(0) のようにx(0)が出てきますよね。 ラプラス変換の問題集の場合たいてい初期条件が付いているのですが、 初期条件がない場合はこのままx(0)を答えに使用してもよいのでしょうか。 たとえば演算子法で解く問題の場合、 x' = x - 4y y' = x + 5y となっていて、問題集の回答の通り微分演算子で解けば 答えは x = {(C2 - 2C1) - 2C2t}exp(3t) y = (C1 + C2t)exp(3t) (C1,C2は任意定数) となります。一方ラプラス変換で解くと x = (x0 - (2x0 + 4y0)t)exp(3t) y = (y0 + (x0 + 2y0)t)exp(3t) (x0 = x(0),y0 = y(0)) となります。 これは実は C1 = y0, C2 = x0 + 2y0 と置き直すと同じになります。ここで質問です。 (1)このような問題でふつうは任意定数を使うべきでしょうが、 x(0),y(0)を使ったら不正解なのでしょうか。 (2)そもそもx(0),y(0)は任意定数になるのでしょうか。 (3)なんだかラプラス変換があれば微分演算子法は いらない子のような気もしなくはないのですが 気のせいでしょうか? 以上です。よろしくお願いいたします。

  • 微分方程式 ラプラス変換

    微分方程式を解く際のラプラス変換の問題で分からない問題があり、教えていただきたいです。 (1)y' + y = 2y (y(0)=1) (2)y" + 2y' + y = 0 (y(0)=y'(0)=1) どちらも部分分数展開の係数を決定するところがよくわかりません。 申し訳ありませんがよろしくお願いします。

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

    皆様よろしくお願いいたします。 関数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(t) のうち、 初期条件 y(0)=1 を満足するものを求めよう。      dy/dt + 3y = 0 y(t) のラプラス変換を Y(s) とすると、dy/dt のラプラス変換は      sY(s) - y(0) = sY(s) - 1     ←sY(s)はどうやって出てきたの? であるから、微分方程式の両辺のラプラス変換を作ると次の式を得る。      sY(s) - 1 + 3Y(s) = 0 したがって、      Y(s) = 1/(s+3) これをラプラス変換すれば、      y(t) = e^(-3t) ・・・と書いてあるんですが、sY(s)のところが分かりません。 y(t) のラプラス変換を Y(s) とすると、dy/dt のラプラス変換は      sY(s) - y(0) = sY(s) - 1 となる、の sY(s) はどうやって出てきたんですか? 最初の s の出所が知りたいです。 ちなみに、ラプラス変換の表では 「基本的な関数は f(t) で表し、そのラプラス変換をF(s)と表す」そうで、 f(t)     F(s) ------------------- e^(at) f(t) F(s-a) のように書かれています。 sY(s)のようなのは書かれていないと思うんですけど…。 どうか sY(s) を得るまで解説お願いします。