• ベストアンサー

ラプラス変換で・・・

未知関数y=y(t)は0<=tでは値が0で、 y(t)-3y(t-1)=t+1 という関係があるy(t)をラプラス変換で求めたいのですが、どうもよくわかりません。 指数倍の平行移動の性質を使うのかなぁとは思うのですが・・・

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

  • ベストアンサー
  • mmky
  • ベストアンサー率28% (681/2420)
回答No.1

参考程度まで 未知関数y=y(t)は0≧t(tがゼロより小さい)では値が0で、 y(t)-3y(t-1)=t+1 のラプラス変換ですね。 y(t-1)のラプラス変換が出来れば変換できますね。。 ∫[0→∞]y(t)e^-st dt=Y(s) ∫[0→∞]te^-st dt=1/s^2 ∫[0→∞]y(t-1)e^-st dt= t-1=x と置くとy(t-1)=y(x), dt=dx, t=x+1 =∫[-1→∞]y(x)e^-s(x+1)dx= {e^-s}*{∫[0→∞]+∫[-1→0]}y(x)e^-sx dx ={e^-s}Y(s)+ ∫[0→1]}y(-x)e^sx dx ={e^-s}Y(s) :0≧t, y(t)=0 故: y(t)-3y(t-1)=t+1 →Y(s)-3(e^-s)Y(s)=1/s^2 +1/s Y(s){1-3(e^-s)}=1/s^2 +1/s y(s)={1/{1-3(e^-s)}}{1/s^2 +1/s} 式の解釈があっているかどうかわかりませんが、 y(t-1)であればちょっと逆変換は反転公式を使わないといけないかも。

その他の回答 (1)

  • stomachman
  • ベストアンサー率57% (1014/1775)
回答No.2

y(t-1) = y(t)*δ(t-1) (*は畳み込み積分) ですから、 y(t)-3(y(t)*δ(t-1))=t+1 ラプラス変換すると Y(s)-3Y(s)exp(-s)=1/(s^2)+1/s すなわち、 Y(s)=(1+s)/(s^2)/(1-3exp(-s)) と、まあここまでは良いでしょうけど、このあとが面倒ですねえ。分母にある (1-3exp(-s)) というのが、「exponentialに増大するということを除けば、周期性を持つような関数」とでも言うしかないようなヨクワカラン代物を暗示しております。(これが(1-exp(-s))なら、単に周期性を言っているに過ぎない訳ですが。) それで、ラプラス変換のことは忘れてみたらどうでしょうか。つまり方程式 y(t)-3y(t-1)=t+1 を満たすy(t)をとにかく見つけてみようという訳です。 ところがこの方程式はどうも曖昧です。どう解釈すれば良いのでしょう? 方程式がt<0の時でも成り立つ、ということを要求すると、 t≦0のときy(t)=0 であることと両立しません。なぜなら、 y(-10)=0, y(-9)=0 なら y(-9)-3y(-10)=-9+1 となって 0=-8 ですもんね。 だから、 y(t)-3y(t-1)=t+1 はt≧0のときだけ成り立つと考えるべきでしょう。従って、 1>t≧0のとき、y(t)=t+1 でなくてはなりません。(∵y(t-1)=0ですから。) 一方、 y(t)=f(t)+bt+c f(t)=3f(t-1) としてみると、方程式は f(t)+bt+c-3f(t-1)-3b(t-1)-3c=t+1 ゆえに -2bt+3b-2c=t+1 でなくてはならず、よって、 b=-1/2, c=-5/4 と決まります。 これだけだと、f(t)は f(t)=3f(t-1) を満たしさえすれば何でも良い。たとえば f(t)=a(3^t)  (aは任意の定数) なんてのを思いつきます。a=0でも構わない。或いは[t]をtを越えない最大の整数とするとき f(t)=a(3^[t]) でも良い。これは階段状になっています。はたまた f(t)=a(3^[t])sin(2πt) だって構わない訳で、いやいやそれどころか tが有理数のとき、f(t)=(3^[t])sin(2πt) tが無理数のとき、f(t)=-100(3^[t])((sin(100πt))^99) なんてハチャメチャの関数でもアリです。 ところが 1>t≧0のとき、y(t)=t+1 という条件を考慮しますと、f(t)を好き勝手に選ぶ訳には行かなくなります。 1>t≧0のとき、y(t)=f(t)-(1/2)t-5/4=t+1 より 1>t≧0のとき、f(t)=(3/2)t+9/4 でなくてはなりません。そして f(t)=3f(t-1) なのだから t=1のときf(t)は不連続になります。なぜなら t→1のとき(3/2)t+(9/4)→15/4 f(1)=3f(0)=3(9/4)=27/4 ですからね。同様にtが自然数のときf(t)は不連続点を持つことになります。 つまり、 f(t)=(3^[t])((3/2)(t-[t])+9/4) と表されるでしょう。まとめると y(t)=(3^[t])((3/2)(t-[t])+9/4)-(1/2)t-5/4 ということになります。ひえ~ 問題の解釈が合っているかどうかよく分からないから「自信なし」です。

関連するQ&A

  • ラプラス変換の問題について

    ラプラス変換を用いて問題で次の問題が分かりません。 次の関係を同時に満たす関数x(t) y(t)を求めよ。  x' = 3s - y y' = x + y x(0) = 0, y(0) = 1 やったことのない型なのでどうすればいいのかわかりません。 お願いします。

  • ラプラス変換、教えてください

    ラプラス変換を使って、方程式を解く問題です。 (1) ∫cos2(t-τ)y(τ)dτ=e^(2t)sin2t  [0、t] (2) y'+y-5∫e^{-(t-τ)}cos2(t-τ)y(τ)dτ=e^(-t)   [0、t] {y(0)=0} を与えられた条件のなかで解け。 まず(1)は、右辺をラプラス変換をして 2/{(s-2)^2+4} しかし左辺をどのように変換すればよいのか分かりません。畳み込みを使おうと思ったのですが、だめでした。 (2)はラプラス変換をして、   {sY(s)-y(0)}+Y(s)-5~~~=1/(s+1) ∫e^{-(t-τ)}cos2(t-τ)y(τ)dτ の部分のラプラス変換が分かりません。 大まかには分かるのですが、どのような流れで解に持っていけばいいのかも曖昧になっています。 ちなみに答は (1) y(t)=1+e^2t (cos2t+3sin2t) (2) y(t)=(5/2)-(5/2)e^(-2t)-4te^(-t) になります。 どなたか教えてください。 よろしくお願いします。

  • ラプラス変換について

    1.ラプラス変換Y(s)の求め方を教えて下さ い。 (1) y(t) = 10 + 5t (2) y(t) = 5e^-t ー e^-5t (3) y(t) = Asinat + Bcosbt 2.ラプラス変換して、伝達関数G(S) (X(S)/F (s))の求め方を教えて下さい。初期値 は0です。 (1) 2x(t) +10dx(t)/dt=f(t) (2) x(t) +1/5∫x(t)dt=f(t) (3) 2d^2 x(t)/dt^2 +10dx(t)/dt +5x(t)=f (t) 3.逆ラプラス変換のやり方を教えて下さ い。 (1) F(s)=17s/(2sー1)(s^2 +4) (2) F(s)=2s+5/s^2+4s+13 (3) F(s)=1/(sー1)^3 以上教えていただけないでしょうか?

  • ラプラス変換

    f(t)=cost (0<t<π) ,0 (π<t) という関数を階段関数を用いてラプラス変換を求めてください。なるべくわかりやすくお願いします。

  • ラプラス変換について

    こんにちは。 今現在制御工学を勉強しているんですがラプラス変換で一箇所わからないところがあるんで、回答お願いいただけますでしょうか m*d^2y/dt^2=A*g^2(t)のラプラス変換です。そこでわからないのはg(t)の2乗のラプラス変換はどうなるかということです。 よろしくお願いいたします。

  • デルタ関数をラプラス変換すると何故1になるか?

    デルタ関数をラプラス変換すると何故1になるか? わかり易く説明お願いします。 デルタ関数はt=0の時、∞になり       t≠0の時、0である のは理解しているのですが、どうもラプラス変換して何故1になるかが分かりません。 そういうものとして丸暗記するべきことなのでしょうか?

  • ラプラス変換による解法

    ラプラス変換の問題で教えてください。 次の方程式を解け。ただし初期値は0とする。 y'-3y+2∫[τ=0からtまで]y(τ)dτ=e^t という問題です。 ラプラス変換による解法で教えてください。

  • ラプラス変換

    y'+(1/r)y=(γ/r)Θ(t) という微分方程式をラプラス変換を用いて解いてください

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

    ラプラス変換の質問です。 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 はどうやって解くのでしょうか?

  • tが短い短形波のラプラス変換について

    ある力Yが働いたときに,短形波が発生します.この短形波のY軸の値をY0とします.短形波のX軸の値を⊿tとします. このYをラプラス変換するとY(s)となります. けれど,⊿tが十分小さいときY(s)はY0⊿tで近似的に表せると言う記述があるのですが よく分かりません. これについて分かる方がいれば教えてください. よろしくお願いします.