電気学生なら容易に解ける三角関数を含む一階微分方程式

単純な1階微分方程式であるが三角関数が入ると面倒そうである。ラプラス変換をすれば解けそうだが，それでも面倒そうな気がして解くのを躊躇してしまう。でも，電気工学科の学生ならば，比較的簡単に例題の解を導くことができるのだ。さあ，解いてみよう。

例題の式を電気回路で置き換えると，実は$LR$回路に交流電源を印加した場合の電流特性を見ていることと等価になる。つまり，
$$L\frac{di}{dt}+Ri=V\sin \omega t但しt=0のときi=0(1)$$
を解くことと同じである。電気回路に落とし込めそうならば，電気工学科ならお手のものだ。

(1)式を解くにあたり，まず，一般解(過渡解)を考える。つまり，
$$L\frac{di}{dt}+Ri=0(2)$$
を解く。これは，下記のように比較的容易に解くことができる。
$$L\frac{di}{dt}=-Ri$$
$$L{\log }_{e}i=-Rt+K_1$$
$${\log }_{e}i=-\frac{R}{L}t+K_2$$
$$i=K{e}^{-\frac{R}{L}t}ただし，Kは積分定数(3)$$

次に，特殊解(定常解)を求める。さて，どう求めるか？この特殊解の導出が面倒なことが例題の式を解くことを困難にしているのである。でも大丈夫，なぜなら2年生の時にフェイザーを習得しているから。つまり，(1)式は，
$$j\omega Li+Ri=V\underset{―}{I_m\left[e^{j\omega t}\right]}(4)$$
と表せる。(アンダーラインの部分は普通は表記しないが，ここでは，分かりやすさのためにあえて残してみた。) (4)式からは$\sin$の項が消えており，見通しが良さそうだ。しかし，そもそも(4)式は正しいのか？答えはもちろん正しい。なぜならば，$\sin \omega t$はオイラーの公式より$e^{j\omega t}$の虚数項であり，実軸と虚軸は直交しているため，互いに作用しないと考えれば(4)式のように表しても何ら問題はない。もっと言えば，$L$や$R$の線形素子に正弦波の電圧を与えているのだから，電流$i$の周波数は変わるはずもなく，$i$も$e^{j\omega t}$を用いて表されるはずである。(だから，$L\frac{di}{dt}$を$j\omega L$($j\omega$は$e^{j\omega t}$を微分した時に現れる係数成分)と表されるとも言える。)したがって，電気工学科の学生ならば(1)式は通常は，
$$j\omega Li+Ri=V(5)$$
と表すのである。そして，$i$を解いた後に $\sin \omega t$を戻してあげれば答えにたどり着くはずだ。では，早速，(4)式を解いてみよう。
$$\begin{array}{rcl}i& =& \frac{V}{R+j\omega L}{e}^{j\omega t}\\ & =& \frac{V}{R^2+{\omega }^2{L}^2}(R-j\omega L)e^{j\omega t}(6)\end{array}$$
ここで，$R-j\omega L$は大きさが$\sqrt{R^2+{\omega }^2{L}^2}$であり，位相角は$\varphi =-{\mathrm{T}\mathrm{a}\mathrm{n}}^{-1}\frac{\omega L}{R}$であるので，
$$\begin{array}{rcl}i& =& \frac{V}{\sqrt{R^2+{\omega }^2{L}^2}}\sin (\omega t+\varphi )ただし，\varphi =-{\mathrm{T}\mathrm{a}\mathrm{n}}^{-1}\frac{\omega L}{R}(7)\end{array}$$
となり，あっさりと特殊解を求められる。ここまでくれば，あとは簡単。一般解と特殊解を纏めると，
$$\begin{array}{rcl}i(t)& =& K{e}^{-\frac{R}{L}t}+\frac{V}{\sqrt{R^2+{\omega }^2{L}^2}}\sin (\omega t+\varphi )(8)\end{array}$$
$t=0$で$i=0$なので，
$$i(0)=K+\frac{V}{\sqrt{R^2+{\omega }^2{L}^2}}\sin \varphi ＝0(9)$$
となり，
$$K=-\frac{V}{\sqrt{R^2+{\omega }^2{L}^2}}\sin \varphi (10)$$
したがって，
$$i(t)=\frac{V}{\sqrt{R^2+{\omega }^2{L}^2}}[\sin (\omega t+\varphi )-\sin \varphi e^{-\frac{R}{L}t}](11)$$
と解くことができる。

念のため,例題の解は，$L=1$，$R=A$，$V=B$，$\omega =C$と考えればよく，
$$y=\frac{B}{\sqrt{A^2+C^2}}[\sin (Cx+\varphi )-\sin \varphi e^{-Ax}]ただし，\varphi =-{\mathrm{T}\mathrm{a}\mathrm{n}}^{-1}\frac{C}{A}(12)$$
となる。めでたし。めでたし。