• ベストアンサー

y=(exp(-at)-exp(-bt))/(b-a)をt=…の式にす

y=(exp(-at)-exp(-bt))/(b-a)をt=…の式にするには? yがある値をとる時のtの値を算出したいのですが、式が変換できずに困っています。 どなたか解る方がいらっしゃいましたら、解答お願いいたします。

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

  • ベストアンサー
  • info22_
  • ベストアンサー率67% (2650/3922)
回答No.2

ab<0の時 yは単調増加関数または単調減少関数でyとtが1:1に対応しますので、逆関数が存在します。しかし数値係数でなく、a,bのような文字定数に対しては tをyを使って一般的に表すことは無理でしょう。a,bが具体的な数値係数の場合はtの範囲を場合分けすればニュートン法(数値計算法)を使って yに対するtを求めることが出来ます。 ab>0の時 単調増加関数または単調減少関数にならず、単一の極大値または極小値をもちますので 全ての実数 tに対する逆関数は存在しません。なぜなら tは yの2価関数になるからです。 yの変域を極大値または極小値の上または下の範囲に限定すれば1価関数に出来ますので 逆関数が存在します。しかし a,bのような文字定数に対しては tをyを使って一般的に解析的には求めることは無理でしょう。 a,bが具体的な数値係数の場合はtの範囲により場合分けすればニュートン法(数値計算法)を使って yに対するtを1個または2個(2価関数の範囲)、求めることが出来ます。

vallisneria
質問者

お礼

各条件における算出の仕方まで回答していただき、ありがとうございます。 具体的な条件を記載していないために、お手を煩わせてしまい申し訳ありませんでした。 再度質問になってしまうのですが、もしよろしければ以下のことについて教えていただけますでしょうか。 実際の条件はa>0、b>0、t>0で、tを算出する際には、a,bに数値を当てはめて計算する予定です。 上記の式で求めたかったのは、yがある値に到達するまでの時間(t)です。 (極大値がそれよりも小さい場合はtは解なしと判断し、2つ解がある場合は値が小さい方を選択します。) a,bの値を変動させた場合のt(時間)を求めたかったため、t=・・・の式にしたかったのですが無理なようですね。 この場合、 (1) a,bに数値を当てはめる。 (2) y=(exp(-at)-exp(-bt))/(b-a)の極大値およびその時のt(Tmとする)を求める。 (3) 極大値が設定した値よりも小さければ解なしとする。 (4) 大きい場合、0<t<Tmの範囲で数値計算法によりtを求める。 という手順で良いでしょうか? また、数値計算法かそれを詳しく解説しているサイトを教えていただけますでしょうか。 拙い文章のため分かり辛いかもしれませんが、よろしくお願いいたします。

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

その他の回答 (2)

  • info22_
  • ベストアンサー率67% (2650/3922)
回答No.3

#2です。 A#2の補足質問の回答 >この場合、 (1) a,bに数値を当てはめる。 (2) y=(exp(-at)-exp(-bt))/(b-a)の極大値およびその時のt(Tmとする)を求める。 (3) 極大値が設定した値よりも小さければ解なしとする。 (4) 大きい場合、0<t<Tmの範囲で数値計算法によりtを求める。 という手順で良いでしょうか? ↑の手順で良いですよ。 >また、数値計算法かそれを詳しく解説しているサイトを教えていただけますでしょうか。 ニュートン法(ニュートン=ラプソン法ともいう)は高校の教科書に載っている方法です。 参考URL http://computation.cside.com/math/math008.html エクセルを使ったニュートン法 http://homepage1.nifty.com/gfk/excel_newton.htm

参考URL:
http://ja.wikipedia.org/wiki/ニュートン法
vallisneria
質問者

お礼

ありがとうございました。 おかげ様で、tの値を算出することができました。 ちょうどtを算出するエクセルファイルを作成していたところだったので、 紹介いただいたサイトがとても参考になりました。

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

もし x = exp(t) 式は y = (x^(-a) - x^(-b))/(b-a) a,bは具体な値があれば、簡単に解けるが、総括的に解けるかどうか、まだ分かりません。

vallisneria
質問者

お礼

ご回答ありがとうございます。 やはり、a,bが文字係数のままで式を変換するのは難しいのですね。

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

関連するQ&A

  • y=-a(a+b)t+b*exp(-(a+b)t) (a>0,b>0)

    y=-a(a+b)t+b*exp(-(a+b)t) (a>0,b>0) yがある値をとる時のt(t>0)を算出したいのですが、 上記の式をt=・・・の式にすることはできるでしょうか? 実際にtを算出する時にはa,bに数値を当てはめて計算を行うのですが、 a,bの値を変えた場合のtも求めたいので、文字係数のままで式を変換したいのです。 どなたか解る方がいらっしゃいましたら、解答お願いいたします。 補足 上記の式は以下の式と単純化したものです。 もしできるならば、こちらの式でt=・・・にしていただけると助かります。 d-(c*b^2)/(a+b)^2-abct/(a+b)+((c*b^2)/(a+b)^2)*exp(-(a+b)t)=0 (a>0,b>0,c>0,d>0) ↓ -a(a+b)t+b*exp(-(a+b)t)=b-(d*(a+b)^2)/bc b-(d*(a+b)^2)/bc=y と置いて単純化しています。 計算ミス等ありましたらご指摘下さい。

  • 1 y=at^2 x=bt+c のとき dy/dxを求めよ。

    1 y=at^2 x=bt+c のとき dy/dxを求めよ。 2 V=exp(j2π/3) のとき V^2をA+jBの形で求めよ。 3 関数y=tcost+1 の導関数を求めよ。 またこの式の点(1.0)における接線の方程式を求めよ。 4 微分方程式 dy/dt+2y=1 の関数yを求めよ。ただしt=0のときy=0とする 5 t=1における関数 y=t*eの-st乗 の接線の方程式を求めよ。 基本的な問題かもしれませんがお願いします><

  • a,bは定数とし、a>0とする。関数y=a(x^2+2x+3)^2-2a(x^2+2x+3)+bの-2≦x≦2における最大値は14、最小値は3である。

    a,bは定数とし、a>0とする。関数y=a(x^2+2x+3)^2-2a(x^2+2x+3)+bの-2≦x≦2における最大値は14、最小値は3であるとする。このとき、a,bの値を求めよ。 という問題で解法、解釈で分からない部分があります。 以下は解答に沿って自分なりに解釈した結果を書いています。 y=a(x^2+2x+3)^2-2a(x^2+2x+3)+bの-2≦x≦2における最大値は14、最小値は3である・・・(1) ⇒x^2+2x+3=tとおく・・・(2) ⇒y=at^2-2at+bの-2≦x≦2における最大値は14、最小値は3である・・・(3) ⇒-2≦x≦2におけるtの変域を求める・・・(4) ⇒y=at^2-2at+bの2≦t≦11における最大値は14、最小値は3である・・・(5) こう考えると(3)まではいいのですが、(4)からうまく納得できません。 -2≦x≦2だったらtの変域は確かに2≦t≦11なのですが、 そしたらy=at^2-2at+b(2≦t≦11)ですが、元((1))のグラフがどういうものかは分かりませんが、これは元のグラフとは異なっているので、2≦t≦11と求めたは良いものの、y=at^2-2at+bでこの範囲における最大・最小が何故、元の-2≦x≦2における最大・最小の14,3と同じなのかイメージが沸きません。

  • 三角関数と指数関数のフーリヘ変換について

    f(t)=exp(-at)cos(bt) (t>=0) t<0は、0。をフーリヘ変換したいのですが、どうすればいいでしょうか? どういうやり方をすればいいのでしょうか? 自分でやってみたのですが、自信がありません。 自分でやってみた考え方はこうです。 cos(bt)=exp(ibt)として、後でReal partをとってやればいいのでは?と思いました。 よって、f(t)=exp(-at)exp(ibt)としました。 ここで、F(ω)=∫[0 to ∞]exp(-at)exp(-iωt)を計算しました。 計算結果は、1/(a+iω)となりました。 フーリヘ変換の性質として、exp(ibt)f(t)=F(ω-b)という性質があるので、 先ほど、得られた、1/(a+iω)を、-bだけ推移させるという考え方で 1/{a+i(ω-b)}となりました。 今回求めているのは、exp(-at)cos(bt)のフーリヘ変換であり、 つまり、1/{a+i(ω-b)}のReal partを取ることにして、 Re[1/{a+i(ω-b)}]=a/{a^2+(ω-b)^2}となりました。 F(ω)=a/{a^2+(ω-b)^2} となったのですが、このReal partを取る考え方はあっていますでしょうか? フーリヘ逆変換して答えを確かめようと思ったのですが、 逆変換が分からなく、答えを確かめられないです。 よければ、逆変換もしていただけるとありがたいです。

  • y=A・x・exp(B・√x)の解法

    y=A・x・exp(B・√x) の B・√x≠0のときのxについての解法を教えてください。

  • 下記(a)(b)の式からt及びσ0の値を求めたいのですが、

    下記(a)(b)の式からt及びσ0の値を求めたいのですが、 どのように計算すればよいのでしょうか? 3.807=lnσ0-600/t    ...(a) 3.689=lnσ0-10,800/t  ...(b) 以上です。よろしくお願い致します。

  • 両辺にyを含んだexpの式をグラフ化したいのですが

    y=A*(exp((x-y*B)/C)-1)+(x-y*B)/D 上記の式で、定数A,B,C,Dを任意に設定して エクセル上(絶対にエクセルである必要はないんですが)で x-yのグラフを描きたいです。 expの中にyが入っているため、 y=f(x) の形に変形することができません。 よってグラフ化できず、困っています。 A~Dの任意定数と変数xの大まかな範囲は以下の通りです。 A:10^(-10) ~ 10^(-7) B:10 ~ 200 C:0.02 ~ 0.2 D:10^5 ~ 10^7 x:0 ~ 2 A~Dの任意定数を変化させたときの、 グラフの変化の傾向が知りたいです。 よってyは厳密な解である必要はなく、近似値で構いません。 質問は以上です。 何卒力を貸して下さい。よろしくお願いします。 以下、参考になるかわかりませんが、 問題の式がどこから出てきたかの説明です。 現在、ショットキー型ダイオードの特性評価をしており、 サンプルとなるダイオードを以下に示す回路図でモデル化しました。                ダイオード                  ↓               ┌──>l──┐               │        │               │        │ ──/\/\/──┤        ├────               │        │               │        │               └/\/\/-┘ 直列抵抗R1↑   並列抵抗R2↑ (フォントが崩れてないといいんですが・・・見づらくてすいません) 測定対象のダイオードに対して、 直列および並列に抵抗が入っていると考えます。 直列抵抗:ダイオードを形成する各材料による抵抗成分(接触抵抗など) 並列抵抗:ショットキー接合における欠陥によるリーク成分      (低電圧時はR2を通って電流が流れてしまうと考えます) この回路に電圧Vを印加したときに流れる電流Iは I=I0*(exp(q*(V-I*R1)/nkT)-1)+(V-I*R1)/R2 I0:逆方向飽和電流 q:電気素量 n:ダイオードの電気特性を示すnパラメータ k:ボルツマン定数 T:ケルビン温度(ここでは室温:T=300としています) と定義されます。 最初に示した問題の式は I→y V→x I0→A R1→B nkT/q→C R2→D と置き換えを行ったものです。 実際にサンプルを測定して得られたVとIのデータがあり、 そのデータに対してフィッティングをかけて、 A~Dのどのパラメータが支配的なのかを探すことが目的です。

  • y=a*exp(b*x)+cによるフィッティングについて

    測定時間(分)  0  5  10  15  20  25  30 平均血圧(mmHg) 92  86  82  82  80  80  77 というデータが手元にあります。  横軸を測定時間t(分)、縦軸を平均血圧BP(mmHg)にして、 BP=a*exp(-b*t)+cという指数曲線で近似したいと考えていますが、 できないでいます。  ExcelやSPSSを用いると、BP=a*exp(-b*t)でしか近似できません。 お手数なのですが、どなたかご存知の方がいらっしゃったら、 ご教授お願いいたします。  私が使えるソフトはExcelとSPSS(統計用ソフト)です。

  • y(t)=exp(-t)*sin(t)を相平面に描

    質問内容  減衰振動 y(t) = exp(-t) * sin(t) を相平面上に描いて、らせん状になることを確認したい。  50歳の数学が趣味の会社員です。  いま勉強している「技術者のための高等数学-1 常微分方程式」の第5版 (ISBN4-563-00561-4)の 131頁に載っている下記の問題が難しくてどうしても 解けません。  かれこれ 1週間ほど悩んでいます。 12.減衰振動 y(t) = exp(-t) * sin(t) を考える。対応する曲線を相平面上に 描き、この曲線がらせん状をしていることを見よ。 y(t) = exp(-t) * sin(t) から特性方程式 の解はλ = -1±i。(i = √(-1)) ゆえに特性方程式は、 (λ - (-1 + i)) * ((λ - (-1 - i)) = ((λ + 1) - i) * ((λ + 1) + i) = (λ + 1) ^ 2 - i ^ 2 = λ ^ 2 + 2 * λ + 1 + 1 = λ ^ 2 + 2 * λ + 2 特性方程式に対応する微分方程式は、 y'' + 2 * y' + 2 * y = 0(' は d/dt) … (1) v = y' = dy/dt とおくと、 y'' = v' = dv/dt = dv/dy * dy/dt = dv/dy * v ゆえに (1) は、 y'' + 2 * y' + 2 * y = v * dv/dy + 2 * v + 2 * y = 0 dv/dy + 2 + 2 * y / v = 0 … (2) u = v / y, v = u * y とおくと、dv/dy = du/dy * y + u ゆえに (2) は、 dv/dy + 2 + 2 * y / v = du/dy * y + u + 2 + 2 * y / (u * y) = du/dy * y + u + 2 + 2 / u = 0 du/dy * y = -(u + 2 + 2 / u) = -(u ^ 2 + 2 * u + 2) / u u / (u ^ 2 + 2 * u + 2) * du = -dy / y ∫(u / (u ^ 2 + 2 * u + 2))du = -∫(dy / y) + B (B は積分定数) … (3) (3) の右辺について -∫(dy / y) + B = -log|y| + B (3) の左辺について ∫(u / (u ^ 2 + 2 * u + 2))du = (1 / 2) * ∫((2 * u + 2 - 2) / (u ^ 2 + 2 * u + 2))du = (1 / 2) * ∫((2 * u + 2) / (u ^ 2 + 2 * u + 2))du - ∫(1 / (u ^ 2 + 2 * u + 2))du = (1 / 2) * log|u ^ 2 + 2 * u + 2| - ∫(1 / ((u + 1) ^ 2 + 1))du = (1 / 2) * log(u ^ 2 + 2 * u + 2) - atan(u + 1)(atan は tan の逆関数) ゆえに (3) は、 (1 / 2) * log(u ^ 2 + 2 * u + 2) - atan(u + 1) = -log|y| + B log(u ^ 2 + 2 * u + 2) - 2 * atan(u + 1) = -2 * log|y| + C(C = 2 * B) = -log(y ^ 2) +C u = v / y なので、 log(v ^ 2 / y ^ 2 + 2 * v / y + 2) - 2 * atan(v / y + 1) = -log(y ^ 2) + C log((v ^ 2 + 2 * v * y + 2 * y ^ 2) / (y ^ 2)) - 2 * atan(v / y + 1) = -log(y ^ 2) + C log(v ^ 2 + 2 * v * y + 2 * y ^ 2) - log(y ^ 2) - 2 * atan(v / y + 1) = -log(y ^ 2) + C log(v ^ 2 + 2 * v * y + 2 * y ^ 2) - 2 * atan(v / y + 1) = C … (4) ここで手詰まり。これ以上、式を簡単にできません。  ここで yv 平面を極座標に変換してみると、(4) は、 y = r * cos(θ), v = r * sin(θ) とおく log(r ^ 2 * sin(θ) ^ 2 + 2 * r ^ 2 * sin(θ) * cos(θ) + 2 * r ^ 2 * cos(θ) ^ 2) - 2 * atan((r * sin(θ)) / (r * cos(θ)) + 1 ) = C log(r ^ 2 * (sin(θ) ^ 2 + 2 * sin(θ) * cos(θ) + 2 * cos(θ))) - 2 * atan(tan(θ) + 1 ) = C log(r ^ 2 * (1 + 2 * sin(2 * θ) + cos(θ) ^ 2)) - 2 * atan(tan(θ) + 1 ) = C log(r ^ 2 * (1 + 2 * sin(2 * θ) + cos(θ) ^ 2)) = C + 2 * atan(tan(θ) + 1) r ^ 2 * (1 + 2 * sin(2 * θ) + cos(θ) ^ 2) = exp(C + 2 * atan(tan(θ) + 1)) = C * exp(atan(tan(θ) + 1))(exp(C) を Cに再定義) r ^ 2 = C / (1 + 2 * sin(2 * θ) + cos(θ) ^ 2) * exp(atan(tan(θ) + 1))  らせん状ということは、r = α * exp(β * θ) の形にならないといけないと 思うのですが、どうしてもその形に持って行けません。  減衰しない振動 y(t) = sin(t) の場合は、対応する曲線を相平面上に描くと 円になるので、減衰振動がらせん状になることは予想できます。  私がどこかで計算を誤っているのでしょうか? それとも、そもそも解き方が 根本的に誤っているのでしょうか ?  最終学歴は工業高校卒業で大学には行っていません。  45歳から数学を独学で勉強していて、周りには質問に答えてくれる人はいません。勉強している本にも答えは載っていません。  ご教示いただければ幸いです。

  • exp{-a(x^2+y^2)}のフーリエ変換

    exp{-a(x^2+y^2)}のフーリエ変換を行うとどのような関数になるのか教えてください。いろいろ調べたのですが、変数が一つ(xだけ)の場合はいろいろなとこで載っていたのですが、これは見つかりませんでした。よろしくお願いします。