• ベストアンサー

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

inara1の回答

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

Excelのマクロで y を計算する方法を以下に紹介します。 セルに = func1(A,B,C,D,x) と書けば y(電流)を計算してくれます。 x(電圧)がマイナスでも計算できます(デフォルトは-10Vまで)。 y の相対精度は10^(-10)です(変更可)。動作確認済みです。  (1) Excelシートのメニューの[ツール] → [マクロ] → [Visual Basic Editor]  (2) Visual Vasic Editorのメニューの[挿入] → [標準モジュール]  (3) 出てきた空白画面に以下の文を貼り付ける     ↓ここから Function fnc1(A As Double, B As Double, C As Double, D As Double, x As Double) As Double ' If x = 0 Then fnc1 = 0 Exit Function End If ' Dim xmin As Double, ymax As Double Dim y As Double, y1 As Double, y2 As Double Dim s As Double, s1 As Double, s2 As Double, eps As Double ' xmin = -10 '    xの最小値 [V] ymax = 10 '    yの最大値 [A] eps = 10 ^ (-10) ' 収束判定 ' y1 = (xmin / D - A) / (1 + B / D) y2 = ymax s1 = A * (Exp((x - B * y1) / C) - 1) + (x - B * y1) / D - y1 s2 = A * (Exp((x - B * y2) / C) - 1) + (x - B * y2) / D - y2 ' While Abs((y1 - y2) / y2) > eps y = (y1 + y2) / 2 s = A * (Exp((x - B * y) / C) - 1) + (x - B * y) / D - y If Sgn(s1 * s) < 0 Then y2 = y Else y1 = y End If Wend fnc1 = y End Function    ↑ここまで  (4) その画面を閉じるか、Visual Vasic Editorのメニューの[ファイル] → [終了して Microsoft Excelへ戻る]  (5) ユーザー定義関数 func1(A,B,C,D,x) として使用する 最初、逐次近似法で計算してみましたが、x が大きくなると解が振動してしまいました。これは2分法で y の値を探す方法です。この電流-電圧特性は単調増加なので解は1つしかありません。

関連するQ&A

  • 式からxとyを求めたい.

    以下の二つの式から、xとyを求めたいのですがどうしても解けません。 (x^2)*(1-a)^2+[1+(1/2)*(r*a*x)^2]y^2 = 2*[1+(r*a*x)^2] 3*(x^2)+y^2 = 6*(1-a) aとrは定数です. 解き方を教えてください。

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

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

  • ダイオードの特性式、順方向降下電圧について

     LEDの特性について理解を深めようとしている者です。 ダイオードの特性式と順方向降下電圧について質問させてください。  ダイオードの順方向降下電圧はダイオードの種類によって様々で、 シリコンダイオードなら0.6Vから0.7V,LEDになると2.0Vから3.5V程度になる、ということを学びました。 また、ダイオードの電流-電圧特性は I=I0{exp(qV/nkT)-1} で与えられる、とさまざまなところで見たのですが、 この式でq,kはそれぞれ電気素量、ボルツマン定数で一定、 nは1から2までの値を取る、 Tは絶対温度、およそ300K、となるので、 結局大きく変わる変数としては逆方向飽和電流I0のみとなると思うのですが、 この式で、例えば順方向降下電圧3.5Vのダイオードの電流-電圧特性のグラフにフィッティングさせようとすると、 I0がとんでもなく小さな値(1.0×10^(-40)A程)になり、いくら逆方向に電流を流さないとはいっても違和感を覚えてしまいます。 pn接合部の抵抗を考え、 V'=V+RI (V'はダイオードと抵抗にかかる電圧の合計、Rは抵抗) の式を先ほどのダイオード特性式に当てはめ、 I=I0{exp(q(V'-RI)/nkT)-1} という式に変形させフィッティングしてみても、順方向降下電圧にはほぼ影響がなく、やはりI0が低い値になってしまいました。 LEDにおいて、順方向降下電圧が上がってしまうのはやはりI0が極端に低いからなのでしょうか、 それとも他に要因があるのでしょうか? LEDだけでなく、順方向降下電圧に同一温度下でも幅が出るのは、すべてI0によるものなのでしょうか? 物性面で違いが出るのはもちろんわかるのですが、それが特性式のどこに関わってくるのかが理解できていません…。  他に要因があるのでしたら、ダイオードの特性式に絡めて説明していただけると非常に助かります。 よろしくお願いします。

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

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

  • 複素積分

    I1=∫[-∞,-∞]cos(a*x)/(x^2+b^2)dxを複素積分を使って求めます。 まず ∫[-∞,-∞]cos(a*x)/(x^2+b^2)dxの被積分関数の分子にi*sin(a*x)を (iは虚数単位)を加えても加えた部分が奇関数でかわらないので加え ると ∫[-∞,-∞]{cos(a*x)+i*sin(a*x)}/(x^2+b^2)dxとなります するとI=∫[-∞,-∞]exp(i*a*x)/(x^2+b^2)dxです。 ここで複素積分 I=∫exp(i*a*z)/(z^2+b^2)dz (積分路は実軸と虚軸の正の部分を通る 反時計回りの半径Rの半円) またI2=∫exp(i*a*z)/(z^2+b^2)dz (積分路は虚軸の正の部分のみを通 る反時計回りの半径Rの半円)を考えるとRが十分大きいとき I=I1+I2・・・(1)になります。 Iは留数定理よりI=2*π*i*Res[f,i*b]=π*exp(-a*b)/b・・・(2) I2はz=R*exp(i*θ)とおき I2=∫[0,π]exp(i*a*R*exp(i*θ))/(R*exp(i*θ)^2+b^2)dθ =∫[0,π]exp(-a*R*sinθ+)*exp(i*a*R*cosθ)*i*R*exp(i*θ)/(R^2*exp (2*i*θ)+b^2)dθ 三角不等式より 0<|I2|<∫[0,π]|exp(-a*R*sinθ+)*exp(i*a*R*cosθ)*i*R*exp(i*θ)|/|(R^2*exp(2*i*θ)+b^2)|dθ<π*R*exp(-a*R*sinθ)/|-R^2+b^2|・・・(3) ここでsinθ >0よりでexp(-a*R*sinθ)<1なので π*R*exp(-a*R*sinθ)/|-R^2+b^2|<π*R/|-R^2+b^2|となり π*R/|-R^2+b^2|はR-->∞で0なので結局 |I2|-->0   なので(1)より I1=π*exp(-a*b)/bが答えです。  これはわかるのですが、スタートで i*sin(a*x)ではなく-i*sin(a*x)を加えても変わらないですよね? そこで-i*sin(a*x)を加えて実際にやってみると (2)の部分はπ*exp(a*b)/bに変わってしまい、また (3)の部分はπ*R*exp(a*R*sinθ)/|-R^2+b^2|となってしまいこれでは R-->∞で発散するように思えます。 どこがまちがっているのでしょうか

  • PCでのグラフの作図法

     まず、状況について書かせて頂きます。  現在、実験でシリコンを用いたショットキーダイオードを作成し、その電流-電圧特性を調べているのですが、回路中の抵抗がどれだけ特性に影響を与えるか見積もるため、理論式に基づきPC上でグラフ化しシミュレーションを行おうと考えています。  テキストによると、回路の直列抵抗をRs、並列抵抗をRpとし、端子間の電圧と電流をV、I、またダイオードにかかる電圧、流れる電流をVd、Idとすると、  オームの法則から  Vd=V-RsI I=Id+Vd/Rp ショットキー接触の特性の式から  Id=Io{exp(qVd/kT)-1}  (Ioは逆方向飽和電流、qは電子の電荷、kはボルツマン定数、Tは絶対温度) が導かれます。これらの式から  I=Io{exp(q(V-RsI)/kT)-1}+(V-RsI)/Rp とIとVの関係が得られます。  本題の質問なのですが、上の式を実際にPCで図示するにはどのようにすればよろしいでしょうか。上式のままではうまく扱うことができず困っています。さらに何らかの変形を行う必要があるのか、PC上で何か別の手法を用いれば描くことができるのか…。  単に数学の基礎的な問題であるのかもしれませんが、念のため背景についても書かせて頂きました。長々と申し訳ありません。  勉強不足な素人でお恥ずかしいのですが、何かヒントを頂けたらと思います。

  • グラフの形について

    グラフの形について ある関数 y=ax^(2)exp(bx^(2)) について横軸x、縦軸yのグラフを書くとするとどのようなグラフの外形が描かれるのでしょうか。 ちなみに a、b:定数 x^(2)というのはxの二乗 という条件でお願いします!! 分かる方がいましたら教えてください。よろしくお願いします

  • sumとexpについて?

    =sum(exp(x))とは、どのように計算すればいいのでしょうか? 例題)xに入る数字が、 A: 0, B: 2.4, C:4.9, D: 12 というデータセットに対し、 exp(x)/sum(exp(x))は、どのように解けばいいのでしょうか? Aの場合は、 exp(0)/sum(exp(19.3))かなと思いましたが、解が合いません。ご存知の方、教えてください。

  • ホイートストンブリッジにテブナンの定理を

    図のホイートストンブリッジで抵抗R5に流れる電流I5をテブナンの定理を用いて求めよという問題なのですが、端子bc間の電位差は{R2/(R1+R2)-R4/(R3+R4)}Eとなるのはわかります。 抵抗R5を取り去ったとして端子bcから見た回路の等価抵抗R0は、起電力Eを短絡したときの合成抵抗であるから、わかりやすくみれば図(右)のような回路の端子bc間の抵抗である。 よって,R1とR2の並列抵抗とR3とR4の並列抵抗が直列に接続された回路が求める等価回路である。 したがって、R0=R1R2/(R1+R2)+R3R4/(R3+R4) この起電力Eを短絡したときの合成抵抗というのはわかりますが、だからといってなぜR1とR2の並列抵抗とR3とR4の並列抵抗が直列に接続されたものになるのかがわかりません。 そもそもわかりやすくいえばa→dに導線を引いた図になるといっていますが、なぜこのような図になるのか理解できません。ただ単に起電力Eがないものと考えればいいんではないのですか? 私の考えではa→dに導線をひいていない図(左)についてc→b方向に電流を流したときに流れる電流の道から、R1とR3の直列抵抗とR2とR4の直列抵抗が並列に接続されたもののように思えるのですが、なぜこのようにならないんですか? このことについてわかる方がいらっしゃいましたら是非教えていただけないでしょうか。 よろしくお願いします。

  • 行列とexp

    行列Yが A 0 B 0 で与えられた時、 exp(Y)が行列 exp(A) 0 0 exp(B) で与えられるのはなぜでしょうか? 終えてください。宜しくお願い致します。