• 締切済み

ラプラス変換の微分方程式でy(t)を求める問題

ラプラス変換の微分方程式の問題です。ラプラス変換及びラプラス逆変換を使いY(S)を求めてからy(t)を求めます。 (a)y'(t)-y(t)=(2*t-1)*e^t^2, y(0)=2 (b)t*y''(t)+y'(t)+4*t*y(t)=0 ,y(0)=2,y'(0)=0 (c)y'(t)-3*y(t)+2*∫(0→t) y(u)du=u(t-1) ,y(0)=2,ただしuはヘヴィサイドの階段関数 です。(a)は左側の変換は出来るのですが、右側がよくわかりません。合成積を使う見たいですが。 (b)はtが付くとき微分?をする見たいなことは分かるのですがその後がよくわかりません。 (c)は根本的にわかりません。 よろしくお願いします。

みんなの回答

  • inara1
  • ベストアンサー率78% (652/834)
回答No.1

(a) については別質問(http://sanwa.okwave.jp/qa5564418.html)のANo.3に答えがありますが、どのようにしてそうなるのかの流れだけを紹介します(ラプラス変換の部分はご自分で計算してください)。 微分方程式 y'(t) - y(t) = (2*t-1)*e^(t^2) の右辺はこのままではラプラス変換できません(発散するので)。そこで、両辺を 0 でない e^(t^2) で割ると    y'(t)*e^(-t^2) - y(t)*e^(-t^2) = 2*t-1 --- (1) となります。上の左辺の第2項の y(t)*e^(-t^2) は t の関数なので、それを x(t) とおけば    x(t) = y(t)*e^(-t^2)    → y(t) = x(t)*e^(t^2) --- (2) 式(2)の両辺を t で微分すれば    y'(t) = { x'(t) + 2*t*x(t) }*e^(t^2) --- (3) 式(2), (3) を式(1)に代入すれば     x'(t) + ( 2*t - 1 )*x(t) = 2*t-1 初期条件は式(2)から y(0) = x(0) = 2 となります。これをラプラス変換を使って x(t) について解けば y(t) が求められます。 ちなみにこの問題はここ(http://sanwa.okwave.jp/qa770180.html)に全く同じ問題が出ています。

関連するQ&A

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

    ラプラス変換で微分方程式を解く 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 をラプラス変換により微分方程式で解きたいのですが、部分分数分解がうまくいかないため、解けません。 本日テストですので、途中式、答えをお願いします。

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

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

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

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

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

    微分方程式についてふと思った疑問です。 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)なんだかラプラス変換があれば微分演算子法は いらない子のような気もしなくはないのですが 気のせいでしょうか? 以上です。よろしくお願いいたします。

  • ラプラス変換を用いて微分方程式 - ステップ関数

    下記の問題で、なぜ、いきなりステップ関数 u(t) が出てきたのか理由を教えて下さい。 ラプラス変換を用いて次の微分方程式の解を求めよ。 dy/dt + 3y = f(t) 模範解答 ※ステップ関数をu(t)と記す。 与えられた微分方程式をラプラス変換すると      sY(s) - y(0) + 3Y(s) = F(s) 整理すると      Y(s) = { F(s) + y(0) } / ( s + 3 ) ラプラス逆変換して      y(t) = f(t) * { e^(-3t) * u(t) } + y(0) * e^(-3t) * u(t)         = ∫[0,t] f(τ) * e^{ -3(t-τ) } dτ + y(0) * e^(-3t) * u(t) ・・・と本に書いてあります。(私の回答は u(t) を除けば正解でした。) ただ、ステップ関数はこのラプラス変換の章に入ってすぐにちょっと説明しただけで、 ここ最近の例題の答えにはまったくステップ関数が出てきていませんでした。 例えば、 ラプラス変換を用いて次の微分方程式の解 y(t) を求めよ。      (d^2 y)/(dt^2) - 3 dy/dt + 2y = f(t) という、一つ前の例題の場合、答えは      y = ∫[0,t] f(τ) * [ e^{ 2(t-τ) } - e^(t-τ) ] dτ        + { y'(0) - y(0) } * e^(2t)        + { 2y(0) - y'(0) } * e^(t) でした。似たような問題ですが、こちらにはステップ関数 u(t) がありません。 #今回のメインの問題の答えの左辺はy(0) = …、 #この問題の答えの左辺はy = …ですね。 #(しかも、問題分には「y(t) を求めよ」と書いてあるのに、です)。 #これは誤植でしょうか? ・・・ということで、いつ、どういう場合にステップ関数 u(t) が必要になるのでしょうか? どうか説明をお願いします。

  • 微分方程式を求める問題がわかりません

    ラプラス変換を用いて微分方程式を求める問題がわかりません y''+4y'+5y=δ(t-π) y(0)=y'(0)=0 Y(s)=e^(-sπ)/(s^2+4s+5)までは計算したのですがここからのラプラス変換がわかりません 私の計算は間違っているのでしょうか? 解説をお願いします ちなみに答えはy=U(t-π)e^(2π-2t)sin(t-π)です

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

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