• ベストアンサー

ラプラス変換と微分方程式について

下記の式、 da(t)/dt = kA - f・a(t) + f/p・C(t) - lamda・a(t) kA, f, f/p, lamdaは定数。 をラプラス変換で解くと、 a(t) = -kA/(f+lamda)・(1 - exp((f+lamda)・t) - f/p・exp((f+lamda)・t) * C(t) *はコンボリューション。 となったのですが、合ってますでしょうか。 当方数学が苦手なので、得意な方に是非助けて頂きたいです。 よろしくお願いいたします。

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

  • ベストアンサー
  • bran111
  • ベストアンサー率49% (512/1037)
回答No.1

コンボリューションがこの場合何を意味するのか分かりませんが C(t)が具体的に与えられていないのでC(t)の部分については積分のまま残るはず。 元の微分方程式は定数係数の線形1次方程式で完全解が求められます。それは a(t)=kA/(f+lamda)+c'exp((f+lamda)・t) +(f/p)exp(-(f+lamda)・t)∫C(t)exp((f+lamda)・t)dt c'は積分定数 です。解を元の微分方程式に代入してOKかどうか確認してください。

makomakoman2003
質問者

お礼

符号の違い等、認識いたしました。 ご指摘いただき大変ありがとうございました。 積分定数も初期値が分かるので計算できそうです。 とても助かりました!!

makomakoman2003
質問者

補足

コンボリューションは畳み込み積分です。

関連するQ&A

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

    お願いします。 連立微分方程式をラプラス変換で解くとき、 たとえば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)なんだかラプラス変換があれば微分演算子法は いらない子のような気もしなくはないのですが 気のせいでしょうか? 以上です。よろしくお願いいたします。

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

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

  • 以下の微分方程式をラプラス変換で解いてください

    以下の微分方程式をラプラス変換で解いてください。 答えがわからなくて困っています。よろしくお願いします。 dS(t) / dt = -(S(t) / A) + Σ δ(t - k) dS(t) / dt = -(S(t) / B) + A(1 - S(t)) * Σ δ(t -k) (k,A,Bは定数) よろしくお願いします。

  • ラプラス変換を常微分方程式に応用

    ラプラス変換を用いて、次の微分方程式の解 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) を得るまで解説お願いします。

  • 微分方程式でラプラス変換を用いるやりかた(急いでます)

    急いでます!!次の計算などのやり方を答えを含めて教えてください。 1.t^(-a)のラプラス変換をガンマ関数を用いて表せ。 2.f(τ)=(e^(-2τ)-e^(-τ))/τの逆ラプラス変換を求めよ。 3.x'(t)-x(t)=(2t-1)e^t^2,x(0)=2をラプラス変換を用いて解け。 4.Hebisideの階段関数H(t)={1(t>0),0.5(t=0),0(t<0)}としたとき、H(t-a)のラプラス変換を求めよ。 以上ですが、1は自分は、-aτ^(a-1)Γ(-a)となりました。あっているでしょうか? 2はお手上げです。変換表を見てもそれらしいのは乗っていませんし。 3.はe^t^2のラプラス変換ができませんでした。 4.はt>aのとき1/τ,t=aのとき、t<aのときと場合分けするだけでいいのでしょうか? なにとぞよろしくお願いします。

  • 指数関数の入れ子構造のラプラス変換

    f(t) = a/exp(b*exp(-ct)) (a,b,cは定数、tは時間[s]) という関数を、F(s) = ○ の形にラプラス変換したいのですが、誰かわかる方いらっしゃいますか?

  • ラプラス変換による回路解析

    RC直列回路を2パラでつないだ回路をラプラス変換で解析しているのですが、 次のような形の式が出てきました。(余計な係数は省略してます。)   1/(s+A) + 1/{s(s+B)}   (A,Bは定数) この式を逆変換したいのですが、第二項は分母がsでくくられているので そのまま公式に入れられませんでした。 そこで、ラプラス逆変換の定義式   f(x) = (1/2πj)∫F(s)exp(st)ds  ( s:[c+jω,c-jω]) をtで不定積分 (積分定数は0) してみると、  ∫f(x)dt = (1/2πj)∫F(s)exp(st)ds × (1/s) といった感じでうまいこと分母にsが飛び出してきてくれたので、 これを参考に、先ほどの逆変換も次のように計算してみました。  L-1[ 1/{s(s+B)} ]  =  ∫{ L-1 [1/(s+B)] }dt               =  ∫{ exp(-Bt) }dt               =  -exp(-Bt)/B 個人的にはこれで間違っていないような気はするのですが、いきなり積分定数の無い不定積分を使うのは何か気持ち悪いですし、それにこの解法で導かれる答えもどこかしっくりきません。(直感ですが。) しかし、自分では正解かどうか確かめる術が無いので、もしラプラス変換に詳しい方おりましたらご助言願いたいと思います。宜しくお願い致します。

  • ラプラス変換で解けるのでしょうか

     次の連立微分方程式は、ラプラス変換を使って解けますか。あるいは他の方法があれば教えて下さい。 dx/dt=(a+bx+cy)x dy/dt=(p+qx+ry)y x,yは時間tの関数、a,b,c,p,q,rは定数。

  • 微分方程式の求め方について

    はじめまして、微分方程式の解き方がわからず困っております。 問題は以下となります。 dΔf/dt=-5Δf-0.5 から Δf(t=0)=0として上記の微分方程式を解くと以下の式になる Δf(t)=-0.1{1-exp(-t/0.2)} このΔf(t)=-0.1{1-exp(-t/0.2)}がどうやって導いたのかがわからないです。 (ラプラス変換でも解けるみたいですが・・・) ご教授のほどお願いいたします。

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

    ラプラス変換を覚えて、微分方程式を簡単に解いてしまおうと思い勉強していたのですが、 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()内に出るのか分かりません。