ラプラス変換をやろう
[size=150][b]このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][br][/b]前回まで、常微分方程式のタイプ分けをタイプに応じた解法を見てきた。[br]今回からは、視点を変えて解法の種類を広げてみよう。まずはラプラス変換。[b][br][size=150][b]<ラプラス変換>[/b][/size][br][/b][/size]ラプラス変換は、微分方程式を代数式に変えて、代数的に単純化する、[br]それを逆変換するだけで、微分方程式が解けてしまう。[br][br]つまり、微分積分計算を代数計算に置き換えて、手計算を減らすことができるというおいしい手順の1つだ。ラプラスさんの編み出した方法は、電気技師のヘビサイドという人物が演算子を使って電気回路の設計に使われて広がったといわれている。制御システムを作るのにも使われる。[br][br]もちろんどんな微分方程式でも解けるというわけではないにしろ幅広いタイプに使える。[br]・もとの関数f(t)に対する[b]ラプラス変換 L(f)=integral(exp(-[color=#0000ff]s[/color]t) f(t), dt, [0, ∞])= ∫e[sup]-st[/sup] f(t) dt[br][br][/b]不定積分ではなく定積分の形になったいるが、[0, ∞]が固定なので、最後の数値代入以外は省略する。[br]たて棒の上下に積分範囲をかくスタイルもあるが、行の節約から[∞,0]閉区間のように記すことにする。[br]積分範囲に∞があるので特異な積分だけど、数Tを入れてからTを∞にすると考えればよいね。[br]積分変数はtだが、指数関数の[b]変数がtの-s倍[/b]になっていることに着目しよう。[br]指数関数exp(-st)をf(t)にかけることで、f(t)を[b]急速に0に近づける[/b]ことで、片側無限区間の積分値が[br]収束するようにしているね。[br][br]・ラプラス変換は積分の1種なので[b][color=#0000ff]線形性がある[/color][/b]。だから、基本となる関数の変換公式の組み合わせで[br]変換ができてしまう。逆変換も同様。[br][color=#0000ff](例) [/color][br]L([b]1[/b])=∫e[sup]-st[/sup] 1dt= -1/s e[sup]-st[/sup][∞,0]=0-(-1/s)=[b]1/s=F(s)[/b][br]L([b]e[sup]at[/sup][/b])=∫e[sup](a-s)t[/sup]dt= 1/(a-s) e[sup](a-s)t[/sup][∞,0]=0-(1/(a-s))=[b]1/(s-a)=F(s)[br][/b]ステップ関数θ(非負で1、負で0)L([b]θ[/b])=∫e[sup]-st [/sup]θdt=∫e[sup]-st[/sup] 1dt=1/s[br]デルタ関数(δ(t)=0 , ∫g(t)δ(t)[∞,-∞] dt = g(0) ) L([b]δ[/b])=∫e[sup]-st [/sup]δ(t)dt=e[sup]-s 0[/sup] =1[br]・第一推移定理L(f(t))=F(s)ならば、L([b]e[sup]at[/sup][/b]f(t))=F(s[b]-a[/b]) ,L([b]e[sup]-at[/sup][/b]f(t))=F(s+[b]a[/b]) [br]・相似定理 L(f(t))=F(s)ならば、aが正のときL(f([b]a[/b]t))=[b]1/a[/b]F(s[b]/a[/b])[br]・[b]基本変換集( f(t)→L(f)) [br] べき関数1→1/s , t → 1/s[sup]2[/sup], t[sup]2[/sup] →2!/s[sup]3[/sup], nが整数のとき [/b][math]t^n→\frac{n!}{s^{n+1}}[/math][b][sup][/sup]、aが正数のとき[/b][math]t^a→\frac{Γ(a+1)}{t^{a+1}}[/math][b] [br] 指数関数e[sup]at[/sup]→[/b][math]\frac{1}{s-a}[/math][b] (aが複素数でも)[br] 三角関数cosωt→[/b][math]\frac{s}{s^2+\omega^2}[/math][b] , sinωt→[/b][math]\frac{\omega}{s^2+\omega^2}[/math][b] ,cosh at=(e[sup]at[/sup]+e[sup]-at[/sup])/2→[/b][math]\frac{s}{s^2-a^2}[/math][b], sinh at=(e[sup]at[/sup]-e[sup]-at[/sup])/2→[/b][math]\frac{a}{s^2-a^2}[/math][b] [br] 振幅変化三角関数 e[sup]at[/sup] cosωt→[/b][math]\frac{s-a}{\left(s-a\right)^2+\omega^2}[/math][b] , e[sup]at[/sup] sinωt→ [/b][math]\frac{\omega}{\left(s-a\right)^2+\omega^2}[/math] [br]・[b]微分[/b]の変換集 (sはあるkより大きい、そこではfは連続関数)初期値[b]f(0)=z[/b]とする。[b]f→Y[/b]とする。[br] 微分を1階するごとにYに[b][color=#0000ff]sを1回かける[/color][/b]。[br] L(f')=s[b]Y[/b]-z (sの1次式)[br] L(f'')=s[sup]2[/sup][b]Y[/b]-sz-z' (sの2次式)[br] ....................................... (n解のLはn次式だけど、最高次以外はzのn-1階以下の微分の初期値が係数)[br] L(f[sup](n)[/sup])=s[sup]n[/sup][b]Y[/b]-(s[sup]n-1[/sup]z+s[sup]n-2[/sup]z'+......+ z[sup](n-2)[/sup]) sの指数とzの微分階数の和がn-1のsのn-1次多項式をひく。[br]・[b]積分[/b]の変換(微分のs倍の[b][color=#0000ff]逆だからsで割る[/color][/b])[br] L(∫[b]f(r)[/b] dr, [0,t] ) = [b]1/s Y[/b][br]・[b]部分分数展開[/b]の例[br] [math]\frac{1}{\left(s-a\right)\left(s-b\right)}=\frac{\left(s-b\right)-\left(s-a\right)}{\left(s-a\right)\left(s-b\right)}\cdot\frac{1}{a-b}=\frac{1}{a-b}\cdot\left\{\frac{1}{s-a}-\frac{1}{s-b}\right\}[/math][br]・ラプラス変換の手順[br](ステップL)[br] 微分方程式y''+ay'+by=r(t) , 初期値z=K0, z'=K1とする。[b]L(y)=Y[/b], L(r)=Rとすると、[br] 1y''+ay'+by→1 L( y'')+aL(y')+bL(y)=1(s[sup]2[/sup][b]Y[/b]-sz-z' )+a(s[b]Y[/b]-z )+b [b]Y[/b]=(1 s[sup]2[/sup]+as+b) [b]Y - [/b](s+a)z- z' =R[br](ステップY)[br] Q= 1/ (s[sup]2[/sup]+as+b)とおくと、Y={(s+a)z+z'} Q+RQ。z=z'=0なら、Y=RQとなりQ=Y/R=L(出力)/L(入力)[br] Yを[b]部分分数展開[/b]などを使って、基本変換集のL(f)にあたる式の和に分解する。[br](ステップL-1)[br] 基本変換集のL(f)部分をf(t)に逆変換する。[br][color=#0000ff](例)[/color][br] y''-y=t , f(0)=zとするとき、z=1, z'=1。L(y)=Yとする。[br](Lする) Y(s[sup]2[/sup]-1)=sz+z'+1/s[sup]2[/sup]= s+1+1/s[sup]2[br][/sup](Yでとく)Q=1/(s[sup]2[/sup]-1)とおく、Y=(s+1+1/s[sup]2[/sup])Q=(s+1)/(s[sup]2[/sup]-1)+1/s[sup]2[/sup](s[sup]2[/sup]-1)=1/(s-1)+(1/(s[sup]2[/sup]-1)-1/s[sup]2[/sup])[br](逆Lする)f(t)=L[sup]-1[/sup](Y)=L[sup]-1[/sup][1/(s-1)]+L[sup]-1[/sup][1/(s[sup]2[/sup]-1)]-L[sup]-1[/sup][1/s[sup]2[/sup]]=e[sup]t[/sup]+sinh t -t =e[sup]t[/sup]+(e[sup]t[/sup]-e[sup]-t[/sup])/2-t=3/2e[sup]t[/sup]-1/2e[sup]-t [/sup]-t
[color=#0000ff](例)[br][/color]yはt秒後の自由落下の高さ[br]加速度はy'’=-g, 最初の高さf(0)=y0, 最初の速度f'(0)=0のとき[br]y''=-g , f(0)=zとするとき、z=y0, z'=0。L(y)=Yとする。[br](Lする) Ys[sup]2[/sup]-(sz+z')=-g/s Ys[sup]2[/sup]-sz=-g/s Ys[sup]2[/sup]=sz-g/s[sup][br][/sup](Yでとく)Q=1/s[sup]2[/sup]とおく、Y=(sz-g/s[sup][/sup])Q=z/s-g/s[sup]3[/sup][br](逆Lする)f(t)=L[sup]-1[/sup](Y)=zL[sup]-1[/sup][1/s]-gL[sup]-1[/sup][1/s[sup]3[/sup]]=z - 1/2g t[sup]2[/sup]=y0- 1/2g t[sup]2[/sup] [br][color=#0000ff](例)[br][/color]y'+2y=e[sup]-t[/sup], z=z'=0[br](Lする) Y(s+2)=z+1/(s+1)=1/(s+1)[sup][br][/sup](Yでとく)Q=1/(s+2)おく、Y=1/(s+1)(s+2)=(1/(s+1)- 1/(s+2))[br](逆Lする)f(t)=L[sup]-1[/sup](Y)=zL[sup]-1[/sup][1/(s+1)]-L[sup]-1[/sup][1/(s+2)[sup][/sup]]t=e [sup]-t[/sup]- e[sup]-2t[/sup] [br][color=#0000ff](例)[/color][br]y''+3y'+2y=e[sup]it[/sup], z=z'=0[br](Lする) Y(s[sup]2[/sup]+3s+2)=1/(s-i)[br] 特性方程式t[sup]2[/sup]+ 3s+2=0の解はp,q=-1,-2だから、Y(s+1)(s+2)=1/(s-i)[br](Yでとく) Y=1/(s-i)(s+1)(s+2)=[math]\frac{1}{s-i}\left\{\frac{1}{s+1}-\frac{1}{s+2}\right\}=\frac{1}{\left(s-i\right)\left(s+1\right)}-\frac{1}{\left(s-i\right)\left(s+2\right)}=\frac{1}{i+1}\left\{\frac{1}{s-i}-\frac{1}{s+1}\right\}-\frac{1}{i+2}\left\{\frac{1}{s-i}-\frac{1}{s+2}\right\}[/math] [br]=[math]\frac{1-i}{2}\left\{\frac{1}{s-i}-\frac{1}{s+1}\right\}-\frac{2-i}{5}\left\{\frac{1}{s-i}-\frac{1}{s+2}\right\}=\left(\frac{1-i}{2}-\frac{2-i}{5}\right)\cdot\frac{1}{s-i}-\frac{1-i}{2}\cdot\frac{1}{s+1}-\frac{2-i}{5}\cdot\frac{1}{s+2}[/math][br]=[math]\frac{1}{10}\left\{\frac{1-3i}{s-i}-\frac{5-5i}{s+1}+\frac{4-2i}{s+2}\right\}[/math][br](逆Lする)y(t)=[math]\frac{1}{10}\left(\left(1-3i\right)e^{it}-\left(5-5i\right)e^{-t}+\left(4-2i\right)e^{-2t}\right)[/math][br]=[math]\frac{1}{10}\left(\left(1-3i\right)\left(cost+isint\right)-\left(5-5i\right)e^{-t}+\left(4-2i\right)e^{-2t}\right)[/math][br] [math]=\frac{1}{10}\left\{\left(cost+3sint-\frac{5}{e^t}+\frac{4}{e^{2t}}\right)+\left(sint-3cost+\frac{5}{e^t}-\frac{2}{e^{2t}}\right)i\right\}[/math][br][color=#0000ff](例)[br][/color]y''+3y'+2y=sin t , z=z'=0 (A)[br]y''+3y'+2y=cos t , z=z'=0(B)[br]Aは上の解の虚部Im(y(t))、Bは上の解の実部Re(y(t))[br] [math]yA\left(t\right)=\frac{1}{10}\left(sint-3cost+\frac{5}{e^t}-\frac{2}{e^{2t}}\right)[/math][br] [math]yB\left(t\right)=\frac{1}{10}\left(cost+3sint-\frac{5}{e^t}+\frac{4}{e^{2t}}\right)[/math]
[color=#9900ff][u][b][size=150]質問:ラプラス変換をコードやるにはどうしたらよいでしょうか。[br][/size][/b][/u][/color][br][b][size=150]<Pytho with Sympy>[br][/size][/b]Sympyをつかいます。[br]残念ですがワンストップでは、計算できません。[br]1回目の入力で、ラプラス変換にはtの関数式f(t)を入力して渡します。[br]ラプラス変換の結果はsympy流の記述で等式が出力されます。[br]2回目の入力で、結果の等式を目で見ながら、通常に近い形の等式を入力します。[br]Yで解き、部分分数計算し、逆変換するを一気にできます。[br]青い部分が毎回変わる入力と出力です。他は、入力内容に応じて変わることがあるかもしれません。[br]Eq(左辺,右辺)で等式を作ります。[b]laplace_correspondenceで対応する関数{f: F}を等式に伝えましょう。[br]F = F.subs({z:1, dz:1})は等式Fに初期値を代入します。[/b][br]solve(F,Y)は等式FをY=にした右辺を答えます。[br]apart(s)はsの有理関数を部分分数に分解します。[br][br]手動のプロセスと比べてみよう。[br]y''-y=t , f(0)=zとするとき、z=1, z'=1。L(y)=Yとする。[br](Lする) Y(s[sup]2[/sup]-1)=sz+z'+1/s[sup]2[/sup]= s+1+1/s[sup]2[br][/sup](Yでとく)Q=1/(s[sup]2[/sup]-1)とおく、Y=(s+1+1/s[sup]2[/sup])Q=(s+1)/(s[sup]2[/sup]-1)+1/s[sup]2[/sup](s[sup]2[/sup]-1)=1/(s-1)+(1/(s[sup]2[/sup]-1)-1/s[sup]2[/sup])[br](逆Lする)f(t)=L[sup]-1[/sup](Y)=L[sup]-1[/sup][1/(s-1)]+L[sup]-1[/sup][1/(s[sup]2[/sup]-1)]-L[sup]-1[/sup][1/s[sup]2[/sup]]=e[sup]t[/sup]+sinh t -t =e[sup]t[/sup]+(e[sup]t[/sup]-e[sup]-t[/sup])/2-t=3/2e[sup]t[/sup]-1/2e[sup]-t [/sup]-t [br][br][b]#{IN]part1[br][/b]from sympy import *[br]from sympy.abc import t, s[br]y = Function('y')(t) #もとの関数[br]Y = Function('Y')(s) #ラプラス変換後の関数[br][b][color=#0000ff]ddy=Derivative(y, t, 2)[br]dy=Derivative(y, t)[br]f_left = ddy - y[br]f_right = t[br][/color][/b]ode = [b]Eq[/b](f_left, f_right)[br]print(f"もと: {ode}")[br]F_left = [b]laplace_transform[/b](f_left, t, s, noconds=True) [br]F_right = [b]laplace_transform[/b](f_right, t, s, noconds=True) [br]lap = [b]Eq[/b](F_left, F_right)[br]F = [b]laplace_correspondence[/b](lap, {y : Y})[br]print(f"Lのあと:{F}")[br][b]#[OUT][/b][b]もと: Eq(-y(t) + Derivative(y(t), (t, 2)), t)[br]Lのあと:Eq(s**2*LaplaceTransform(y(t), t, s) - s*y(0) - LaplaceTransform(y(t), t, s) - Subs(Derivative(y(t), t), t, 0), s**(-2))[br][/b][b]#[IN]part2[br][/b]#手動で、LaplaceTransform(y(t), t, s)をY、y(0)をz、Subs(Derivative(y(t), t), t, 0)をdzとして[br]#ラプラス変換後の式をかく。初期値z=y(0)=1,dz=y'(0)=1を後でsubstituteする。 [br][b][color=#0000ff]z, dz = symbols('z, dz')[br]F = Eq(s**2*Y - s*z - Y - dz, s**(-2))[br]F = F.subs({z:1, dz:1})[br][/color][/b]print(f"変換後の等式:{F}")[br][b]#Yで解いて、部分分数分解して、逆L[br][/b]Y_sol = [b]solve[/b](F, Y)[0] [br]Y_sol_apart = Y_sol.[b]apart[/b](s)[br]print(f"Y = {Y_sol_apart}")[br]y_t = [b]inverse_laplace_transform[/b](Y_sol_apart, s, t)[br]y_t = y_t.[b]subs[/b]({Heaviside(t):1})[br]print(f"y ={y_t}")[br][b]#[OUT][br][color=#0000ff]変換後の等式:Eq(s**2*Y(s) - s - Y(s) - 1, s**(-2))[br]Y = -1/(2*(s + 1)) + 3/(2*(s - 1)) - 1/s**2[br]y =-t + 3*exp(t)/2 - exp(-t)/2[/color][/b]
[color=#9900ff][u][b][size=150]質問:複素関数のラプラス変換をコードやるにはどうしたらよいでしょうか。[br][/size][/b][/u][/color][br][b][size=150]<Pytho with Sympy>[br][/size][/b]pythonでは虚数単位がiではなくjになることに気をつければ、実数と同様に作動します。[br][br]手動のプロセスと比べてみよう。[br]#y''+3y'+2y=e^it, z=z'=0[br]from sympy import *[br]from sympy.abc import t, s[br]y = Function('y')(t) #もとの関数[br]Y = Function('Y')(s) #ラプラス変換後の関数[br]ddy=Derivative(y, t, 2)[br]dy=Derivative(y, t)[br][b][color=#0000ff]f_left = ddy + 3*dy + 2*y[br]f_right = exp(j*t)[br][/color][/b]ode = Eq(f_left, f_right)[br]print(f"もと: {ode}")[br]F_left = laplace_transform(f_left, t, s, noconds=True) [br]F_right = laplace_transform(f_right, t, s, noconds=True) [br]lap = Eq(F_left, F_right)[br]F = laplace_correspondence(lap, {y : Y})[br]print(f"Lのあと:{F}")[br][color=#0000ff][b]#[OUT][br][/b][/color][b][color=#0000ff]もと: Eq(2*y(t) + 3*Derivative(y(t), t) + Derivative(y(t), (t, 2)), exp(j*t))[br]Lのあと:Eq(s**2*LaplaceTransform(y(t), t, s) + 3*s*LaplaceTransform(y(t), t, s) - s*y(0) + 2*LaplaceTransform(y(t), t, s) - 3*y(0) - Subs(Derivative(y(t), t), t, 0), 1/(-j + s))[br][/color][br][/b]#(Lする) Y(s2+3s+2)=1/(s-i)#手動で、LaplaceTransform(y(t), t, s)をY、y(0)をz、Subs(Derivative(y(t), t), t, 0)をdzとして[br]#ラプラス変換後の式をかく。初期値z=y(0)=0,dz=y'(0)=0を後でsubstituteする。 [br]z, dz = symbols('z, dz')[br]F = Eq(s**2*Y + 3*s*Y -s*z + 2*Y -3*z -dz, 1/(-j +s))[br]F = F.subs({z:0, dz:0})[br]print(f"変換後の等式:{F}")[br]#Yで解いて分数分解して逆L[br]Y_sol = solve(F, Y)[0][br]Y_sol_apart = Y_sol.apart(s)[br]print(f"Y = {Y_sol_apart}")[br]y_t = inverse_laplace_transform(Y_sol_apart, s, t)[br]y_t = y_t.subs({Heaviside(t):1})[br]print(f"y ={y_t}")[br][b]#[OUT][br][color=#0000ff]変換後の等式:Eq(s**2*Y(s) + 3*s*Y(s) + 2*Y(s), 1/(-j + s))[br]Y = 1/((j + 2)*(s + 2)) - 1/((j + 1)*(s + 1)) + 1/((-j + s)*(j + 1)*(j + 2))[br]y =exp(-2*t)/(j + 2) - exp(-t)/(j + 1) + exp(j*t)/((j + 1)*(j + 2))[br][br][/color][/b]複素数の計算を最低限しかやってくれないようなので、違う答えのように見えてしまいますが同じです。[br][br]y(t)=[math]\frac{1}{10}\left(\left(1-3i\right)e^{it}-\left(5-5i\right)e^{-t}+\left(4-2i\right)e^{-2t}\right)[/math][br][br][b][size=150]<geogebra>[br][/size][/b]geogebraでもラプラス変換、逆変換の関数はあるようですが、[br]基本関数が計算できるだけで、微分の変換ができないようです。だから、微分方程式を解く用途には[br]使えないですね。
[size=150][size=100]コンデンサーは文字通り、電気エネルギー(電荷)を蓄える。[br]さらに、直流を通さず、交流は周波数が高いほど通しやすい。[br]抵抗Rだけでなく、コンデンサーCをつないだ回路で、電流変化を考えてみます。[br][br][/size][b]<RC回路>[br][/b][/size]RC回路はRとCに電圧V[sub]0[/sub]をかけた回路です。[br]キルヒホフ第2法則(電池で高まった電圧はCとRで下がり、電池のマイナス側で0に戻る)を使います。[br]初期条件Q(0)=0でC,R,V0を順につないでスイッチを入れたときの回路に流れる電荷Qを求めよう。[br]Qは時間tの関数とします。[b]電流Iは電荷Qの時間変化なので、I=dQ/dtです[/b]。[br]コンデンサーの両端の電圧V=Q/Cで、抵抗Rの両端の電圧V=IR=RQ'[br]電圧の和がV0に等しいから、R Q' + Q/C = V[sub]0[/sub] 式変形 Q'+Q/CR=V[sub]0[/sub]/R[br][b][color=#0000ff]見やすさのために、一時的にQをyとして表記し、最後にQに戻すことにします。[br][/color][/b]y'+y/CR=V[sub]0[/sub]/R , y(0)=0とするとき、z=0。L(y)=Yとする。[br](Lする) Y(s+1/CR)=z+V[sub]0[/sub]/R/s=V[sub]0[/sub]/R/s[sup][br][/sup](Yでとく)Y=(V[sub]0[/sub]/R){1/s*1/(s+1/CR)}=(V[sub]0[/sub]/R){1/s-1/(s+1/CR)}CR=CV[sub]0[/sub]{1/s-1/(s+1/CR)}[br](逆Lする)f(t)=L[sup]-1[/sup](Y)=CV[sub]0[/sub]{L[sup]-1[/sup][1/s]-L[sup]-1[/sup][1/(s+1/CR]}=CV[sub]0[/sub](1 - e[sup]-t/CR[/sup])[br][b][size=150][color=#0000ff]Q=CV[sub]0[/sub](1 - e[sup]-t/CR[/sup]) [br][/color][/size][/b]e[sup]-t/CR[/sup]はt=0あたりではほぼ1でtが∞で0に近づく。Qは0から増えてCV0で飽和する。[br]また、I=dQ/dt=(CV[sub]0[/sub](1 - e[sup]-t/CR[/sup]))'=CV[sub]0 [/sub](-e[sup]-t/CR[/sup])'=CV[sub]0[/sub]/CR e[sup]-t/CR[/sup]=(V0/R) e[sup]-t/CR[br][/sup]e[sup]-t/CR[/sup]はt=0あたりではほぼ1でtが∞で0に近づく。Iは(V0/R)から減って0になる。
[b][size=150]<直流RLC回路>[br][/size][/b]Lはコイルです。RLC回路は、RC回路にさらにコイルをつないだものです。[br]回路といっても、計算上のものなので、数式の上での前提での[b]架空の推論[/b]をします。[br][b]工作をして実験するわけではありません。[br][/b][br]それぞれの電圧は、V=IR,V=L dI/dt, V=Q/Cとします。[br]定数はR=2、L=1, C=0.5, V0=10とし、[br]コンデンサーの電荷Q, 回路の電流Iは時間tの関数で、初期値Q(0)=I(0)=0としましょう。[br][b]L(I)=Yとします。[/b][b]電流Iは電荷Qの時間変化なので、I=dQ/dt=Q'です。 L(Q)= 1/s L(I)=1/s Y。[br]z=z'=0[br][/b]この回路に[b]減衰しながら振動する電流[/b]が流れるのはどんなときでしょう。[br][br]R,L,Qの電圧の和はR[b]I[/b]+L d[b]I[/b]/dt + Q/C=V[sub]0[br][/sub](Lする)RY+sLY+1/CsY=V0/s sR/LY+s[sup]2[/sup]Y+1/CL Y=V0/L[br](Yでとく) Y=V0/L(s[sup]2[/sup]+R/L s +1/CL)[br] 特性方程式 t[sup]2[/sup]+R/L s + 1/CL=0の判別式D=(R/L)[sup]2[/sup]-4/CL=(1/L)[sup]2[/sup](R[sup]2[/sup]-4L/C)が負のとき、[br] (つまり、Lに対してC,Rが小さいとき、C容量・抵抗RがインダクタンスLに対して小さいとき)[br] (t- α)[sup]2[/sup]+β[sup]2[/sup]=0と平方完成されるとして , V0/L=pとおくと、[br] Y=p *1/{(s- α)[sup]2[/sup]+β[sup]2[/sup]}=p/β *β/{(s- α)[sup]2[/sup]+β[sup]2[/sup]}[br](逆Lする)I=p/β *[b]e[sup]αt[/sup] sinβt [/b]α=-R/2L , β=√-(D/4)[br] 定数を入れると、p=10/1=10, R/L=2/1=2, 1/CL=1/(1*0.5)=2, [br] α=-R/2L=-2/2*1=-1,D=2^2-4*2=-4, β=√-(-4/4)=1[br] I(t) =10/ 1 * e[sup]-t [/sup]sin t = 10 e[sup]-t[/sup] sin t
[b][size=150]<交流RLC回路>[br][/size][/b][br]回路の電圧が時間とともに変化できる、つまり[color=#0000ff][b]交流電圧Vsinωt[/b][/color]をかけた回路ではどうでしょうか。[br]R[b]I[/b]+L d[b]I[/b]/dt + Q/C=[b]Vsinωtです。[/b][br]複素数の範囲で計算することで、複素数の範囲で解きましょう。[br]R,L,Qの電圧の和はR[b]I[/b]+L d[b]I[/b]/dt + Q/C=[b]Ve[/b][sup]iωt[/sup][b]とします。(iω=c, V/L=Aとおく)c[sup]2[/sup]=-ω[sup]2[/sup][br]L(Q)=Yとします。I=dQ/dt=Q'です。z=z'=0[br][/b]RI+L dI/dt + Q/C=Ve[sup]ct[/sup][br]L Q'' +R Q'+ 1/C Q=Ve[sup]ct[/sup][br](Lする)Y(s[sup]2[/sup]+R/L + 1/CL)=A/(s-c)[br](Yでとく) Y=A/(s[sup]2[/sup]+R/L s +1/CL)(s-c)[br] 特性方程式 t[sup]2[/sup]+R/L s + 1/CL=0の2解をa,bとすると、(t-a)(t-b)=0 [br] a=-R/2L +i√D/2L , b = -R/2L -i√D/2L [b]a,bの実部は負の数[br][/b] a+b=-R/L, ab=1/CL[br] Y=A/(s-a)(s-b)(s-c)=A*/(a-b){1/(s-a)(s-c)-1/(s-b)(s-c)}[br] =A/(a-b)[1/(a-c) * {1/(s-a)-1/(s-c)}- 1/(b-c)*{1/(s-b)-1/(s-c)}][br] = A/(a-b)[1/(a-c)(s-a)-1/(a-c)(s-c)-1/(b-c)(s-b)+1/(b-c)(s-c)][br] =-A/(a-b)[1/(c-a)(s-a)+1/(b-c)(s-b)+1/(s-c){1/(a-c)-1/(b-c}][br] =-A/(a-b)[1/(c-a)(s-a)+1/(b-c)(s-b)+1/(s-c){(b-a)/(a-c)(b-c)][br] =-A/(a-b)[(b-c)/(b-c)(c-a)(s-a)+(c-a)/(c-a)(b-c)(s-b)+(a-b)/(s-c)(c-a)(b-c)][br] =-A/{(a-b)(b-c)(c-a)[(b-c)/(s-a)+(c-a)/(s-b)+(a-b)/(s-c)] [br] =[math]\text{\frac{-A}{(a-b)(b-c)(c-a)}\cdot(\frac{b-c}{s-a}+\frac{c-a}{s-b}+\frac{a-b}{s-c})[br]}[/math][br] (逆Lする)Q(t)= [math]\frac{-A}{\left(a-b\right)\left(b-c\right)\left(c-a\right)}\cdot\left(\left(b-c\right)e^{at}+\left(c-a\right)e^{bt}+\left(a-b\right)e^{ct}\right)[/math] [br] a,bの実部は負の数だからtを十分大きくすると、e[sup]at[/sup], e[sup]bt[/sup]の極限は0になる。[br] Q(t)の極限は [math]\frac{-A}{\left(a-b\right)\left(b-c\right)\left(c-a\right)}\cdot\left(\left(a-b\right)e^{ct}\right)=\frac{A}{\left(c-a\right)\left(c-b\right)}e^{ct}=\frac{A}{c^2-\left(a+b\right)c+ab}e^{i\omega t}[/math] [br][br] =(V/L)/[-ω[sup]2[/sup]+R/L iω+1/CL] e[sup]iωt[br][/sup] = V/[-Lω[sup]2[/sup]+R iω+1/C] e[sup]iωt[br][/sup]ラプラス変換で求めた解から減衰する項を消し、交流定常状態の解を出すことができました。[br]この定常状態の虚部を求めてみよう。[br]分母Z(ω)=−Lω[sup]2[/sup]+Riω+1/C1=(1/C−Lω[sup]2[/sup])+i(Rω)=A+iB (複素インピーダンスという)と置くと、[br]|Z(ω)|=A[sup]2[/sup]+B[sup]2[/sup]で、その位相ϕ=arg(Z(ω))=arctan(B/A)です。[br]だから、Q(t)=Ve[sup]iωt[/sup]/(∣Z(ω)∣e[sup]iϕ[/sup])=Ve[sup]i(ωt−ϕ)[/sup]/|Z(ω)|=V/|Z(ω)| * (cos(ωt−ϕ)+i sin(ωt−ϕ))[br]この Q(t) の[b]虚部[/b]は、Im[Q(t)]=V sin(ωt−ϕ)/|Z(ω)|[br]