다음의 RLC회로를 미분방정식으로 풀기
위의 회로에서 입출력 및 초기조건은 다음과 같음.
input : $x(t) = 10e^{-3t}u(t)$ (voltage)
output : $y(t)$ (current)
initial condition : $y(0)=0$, $V_c(0)=5$
- $u(t)$ : unit step function.
1. Differential Equation
KVL에 의하여 다음이 성립
$$\begin{aligned}{ V }_{ L }+{ V }_{ R }+{ V }_{ C }=x\left( t \right) \\ 1\frac { dy\left( t \right) }{ dt } +3y\left( t \right) +{ \frac { 1 }{ \frac { 1 }{ 2 } } }\int { y\left( t \right) dt } =x\left( t \right) \\ \frac { { d }^{ 2 }y\left( t \right) }{ d{ t }^{ 2 } } +3\frac { dy\left( t \right) }{ dt } +2y\left( t \right) =\frac { dx\left( t \right) }{ dt } \\ { D }^{ 2 }y\left( t \right) +3Dy\left( t \right) +2y\left( t \right) =Dx\left( t \right)\end{aligned}$$
- $D$ : 미분 연산자.
2. Homogeneous solution
homogeneous solution을 구하기위해서 input signal을 0으로 하면 다음과 같음.
$$\begin{aligned} D^2 y \left( t \right) + 3D y \left( t \right) +2y \left( t \right) &=0 \\ \left\{ { D }^{ 2 }+3D+2 \right\} y\left( t \right) &=0 \end{aligned}$$
이를 characteristic equation으로 풀면 다음과 같음.
$${ D }^{ 2 }+3D+2=0 \\ \left( D+2 \right) \left( D+1 \right) =0 \\ \therefore { \lambda }_{ 1 }=-2,{ \lambda }_{ 2 }=-1$$
위 eigen value들로부터 구한 system mode로 이루어진 homogeneous sol.은 다음과 같음.
$$\begin{aligned}{ y }_{ h }\left( t \right) &= {C }_{ 1 }{ e }^{ { \lambda }_{ 1 }t }+{ C }_{ 2 }{ e }^{ { \lambda }_{ 2 }t } \\ &={C }_{ 1 }{ e }^{ -2t }+{ C }_{ 2 }{ e }^{ -t } \end{aligned}$$
3. Particular solution
input signal이 $10e^{−3t}u(t)$이므로 particular sol.은 다음과 같은 형태임.
$$
y_p\left(t\right)=Ce^{-3t}
$$
- $C$ : unknown constant.
이 particular sol.의 형태를 system의 미분방정식에 대입하여 풀면 다음과 같음.
$$\begin{aligned} \left\{ { D }^{ 2 }+3D+2 \right\} { y }_{ p }\left( t \right) &= \left\{ D \right\} 10{ e }^{ -3t } \\ \left[ 9C{ e }^{ -3t } \right] +3\left[ -3C{ e }^{ -3t } \right] +2\left[ C{ e }^{ -3t } \right] &= -3\cdot 10{ e }^{ -3t } \\ 2C{ e }^{ -3t } &= -30{ e }^{ -3t }\end{aligned} $$
$$ \therefore C=-15,\quad { y }_{ p }=-15{ e }^{ -3t } $$
4. Complete solution
complete soluton $i(t)$는 다음과 같이 homogeneous solution과 particular solution의 합임.
$$\begin{aligned} y\left( t \right) &={ y }_{ h }+{ y }_{ p } \\&={ C }_{ 1 }{ e }^{ -2t }+{ C }_{ 2 }{ e }^{ -t }-15{ e }^{ -3t }\end{aligned}$$
이제, $C_1$과 $C_2$ 와 같은 모르는 상수의 값을 초기조건을 통해 구해야 함.
일단, 초기값 조건과 complete solution의 형태를 미분방정식에 대입.
$$\begin{aligned}{ V }_{ L }\left( 0 \right) +{ V }_{ R }\left( 0 \right) +{ V }_{ C }\left( 0 \right) &= x\left( 0 \right) \\ 1\frac { dy\left( 0 \right) }{ dt } +3y\left( 0 \right) +{ V }_{ C }\left( 0 \right)&= 10e^{-3\cdot0}\\ D\left( { C }_{ 1 }{ e }^{ -2t }+{ C }_{ 2 }{ e }^{ -t }-15{ e }^{ -3t } \right) |_{t=0}+3\cdot 0+5&=10\\ \left( -2{ C }_{ 1 }{ e }^{ -2t }-{ C }_{ 2 }{ e }^{ -t }+45{ e }^{ -2t } \right)| _{ t=0 }+5&=10\\ -2{ C }_{ 1 }-{ C }_{ 2 }+45+5&=10\\-2{ C }_{ 1 }-{ C }_{ 2 }&=-40 \end{aligned}$$
또한 $y(0)=0$을 통해 다음이 성립.
$$\begin{aligned}y\left( 0 \right) &= \left( C_1e^{-2t}+C_2e^{-t}-15e^{-3t} \right ) \rvert _{t=0} \\0&=C_1+C_2-15\end{aligned}$$
위 두 식으로부터 $C_1$과 $C_2$ 를 구할 수 있음.
$$\begin{aligned}-2C_1-C_2&=-40 \\ C_1 +C_2&=15 \end{aligned}\\ \quad \\ \therefore C_1=25, C_2=-10$$
즉, complete solution은 다음과 같음.
$$y\left( t \right) =25{ e }^{ -2t }-10{ e }^{ -t }-15{ e }^{ -3t }$$
2023.10.03 - [.../Signals and Systems] - [SS] System and Differential Equation
[SS] System and Differential Equation
Differential Equation으로 System을 기술할 경우, 초기조건을 포함한 differential eq.은 system의 완전한 동작특성을 기술할 수 있음. Impulse response로 표기할 경우, system의 zero-state response만을 구할 수 있는 제
dsaint31.tistory.com
5. 여러 Response로 나누기
Complete solution은 다음과 같이 나눌 수 있음.
- Natural Response / Forced Response
- Zero-input Response / Zero-state Response
- Transient Response / Steady state Response
2023.10.02 - [.../Signals and Systems] - [SS] System Response
[SS] System Response
System의 Response는 다음과 같이 3가지 기준으로 분류할 수 있음. Zero-input response vs. Zero-state response "누가 response를 만드는가? (초기조건 vs. 외부입력)" 를 기준으로 분류하는 방식이며 다음과 같은 2
dsaint31.tistory.com
5-1. Natural / Forced Response
여기서, natural response는 다음과 같음.
$$y_{\text{natural}}\left( t \right) =25{ e }^{ -2t }-10{ e }^{ -t }$$
forced response는 다음임.
$$y_{\text{forced}}\left( t \right) =-15{ e }^{ -3t }$$
5-2. Zero-input / Zero-state Response
zero-input response/ zero-state response의 경우로 푼다면 다음과 같음.
zero-input response의 경우 input signal이 0이므로,
$$\begin{aligned}{ V }_{ L }\left( 0 \right) +{ V }_{ R }\left( 0 \right) +{ V }_{ C }\left( 0 \right) &= x\left( 0 \right)=0 \\ 1\frac { dy\left( 0 \right) }{ dt } +3y\left( 0 \right) +{ V }_{ C }\left( 0 \right) &= 0 \\ D\left( { C }_{ 1 }{ e }^{ -2t }+{ C }_{ 2 }{ e }^{ -t } \right) \rvert_{ t=0 }+3\cdot 0+5 &= 0 \\ \left( -2{ C }_{ 1 }{ e }^{ -2t }-{ C }_{ 2 }{ e }^{ -t }\right)\rvert_{ t=0 }+5 &= 0 \\ -2{ C }_{ 1 }-{ C }_{ 2 }+5 &= 0 \\ -2{ C }_{ 1 }-{ C }_{ 2 } &= -5\end{aligned}$$
이고,
$$\begin{aligned}y\left( 0 \right) &= \left(C_1e^{-2t}+C_2e^{-t} \right)\rvert_{t=0} \\ 0 &= C_1+C_2\end{aligned}$$
가 되어서
$C_1$과 $C_2$ 가 아래처럼 구해짐.
$$\begin{aligned}-2C_1-C_2&=-5\\ C_1 +C_2&=0\end{aligned} \\ \quad \\ \therefore C_1=5, C_2=-5$$
즉, zero-input response는 다음과 같음.
$$y_{\text{zero-input}}\left( t \right) =5{ e }^{ -2t }-5{ e }^{ -t }$$
zero-state response는 다음과 같음.
$$y_{\text{zero-state}}\left( t \right) =20{ e }^{ -2t }-5{ e }^{ -t }-15{ e }^{ -3t }$$
5-3. Transient / Steady-state Response
$$y\left( t \right) =25{ e }^{ -2t }-10{ e }^{ -t }-15{ e }^{ -3t }$$
위의 complete solution에서 모든 term들이 decay가 되고 있기 때문에 다음으로 나뉘어짐.
Transient Response
$$y_{\text{transient}} = 25r^{-2t}-10e^{-t}-15e^{-3t}$$
Steady-state Response
$$y_{\text{steady}}=0$$
더 읽어보면 좋은 자료
보통 미분방정식은 Laplace transform을 이용하여 구한다.
다음 URL은 input signal과 initial condition이 조금 다르긴 하지만,
미분방정식을 고전적 방법으로 푸는 방법과 Laplace transform으로 구하는 방법을 보여주고 있으니 참고할 것.
[SS] Differential Equation and Responses (zero-input, zero-state, natural, forced) w/o Laplace Transform
문제 다음 미분방정식 의 시스템이 있다고 하자 $$ \frac{d^2 y(t)}{dt^2} + 3\frac{dy(t)}{dt}+2y(t)=\frac{dx(t)}{dt} $$ 아래와 같은 입력과 초기조건에서 zero-state response (초기조건이 0.) zero-input response (input이 0.)
dsaint31.tistory.com
[SS] Differential Equation and Response (zero-input, zero-state, natural, forced) w/ Laplace Transform
Laplace Transform을 이용한 Differential Equation을 풀기. 문제 다음 미분방정식 의 시스템이 있다고 하자. $$ \frac{d^2 y(t)}{dt^2} + 3\frac{dy(t)}{dt}+2y(t)=\frac{dx(t)}{dt} $$ 아래와 같은 입력과 초기조건에서 zero-state
dsaint31.tistory.com
'... > Signals and Systems' 카테고리의 다른 글
[SS] Differential Equation and Response (zero-input, zero-state, natural, forced) w/ Laplace Transform (0) | 2023.08.22 |
---|---|
[SS] Differential Equation : 1st Canonical Form (0) | 2023.08.22 |
[SS] Differential Equation : 2nd Canonical Form (0) | 2023.08.22 |
[SS] Causal LTI System and Convolution (0) | 2023.08.22 |
[SS] Convolution Example : Pulse Function (0) | 2023.08.22 |