Introduction to Computor (by R. Saito) Electronics Eng.
UEC Education System

Maple 初級 の 微分方程式 (2)

初期値、境界条件付の微分方程式も、簡単に解けます。
以下では RC 回路に V sin(wt) の電圧を加えた時の問題を解きましょう。
> deq:= R*diff(Q(t),t) + Q/C = V*sin(w*t); <= 微分方程式 deq を定義。

           /  d      \
         R |---- Q(t)| + Q/C = V sin(w t)
           \ dt      /

> init_cond:= Q(0)=0, D(Q)(0)=0; <== 初期条件です。D(Q)はQの微分。

         init_cond := Q(0) = 0, D(Q)(0) = 0

> dsolve({deq,init_cond},Q(t));  <== 初期条件付きで、微分方程式を解きます。
bytes used=1000032, alloc=720764, time=2.13 <== 考えています。

                                                    t          2
                 2                           exp(- ---) R V w C
          R V w C  cos(w t)   V C sin(w t)         R C
 Q(t) = - ----------------- + ------------ + -------------------
                  2  2  2          2  2  2           2  2  2 2
             1 + w  R  C      1 + w  R  C      (1 + w  R  C )

               t    3    3  4
        exp(- ---) R  V w  C
              R C
      + ---------------------
                 2  2  2 2
           (1 + w  R  C )


境界条件も同じ様に与えることができます。

コメントまたはアドバイスなどがあれば以下のアドレスへどうぞ。


tutor1@edu.cc.uec.ac.jp