Exponenta Banner

Моделируйте последовательную схему RLC

Откройте модель

Физические системы могут быть описаны как ряд дифференциальных уравнений в неявной форме, $F\left(t,x,\dot{\left\lbrace x\right\rbrace } \right)=0$или в форме неявное пространство состояний $E\dot{x} =A\;x+B\;u\;$

Если$E$ он несингулярен, то система может быть легко преобразована в систему обыкновенных дифференциальных уравнений (ОДУ) и решена как таковая:

$$\dot{x} =\left(E^{-1} A\right)x+\left(E^{-1} B\right)u$$

Много раз состояния системы появляются без прямого отношения к их производным, обычно представляющим физические законы сохранения. Для примера:

$$\begin{array}{l}\dot{x_1 } =x_2 \\0\;=x_1 +x_2 \end{array}$$

В этом случае$E$ сингулярна и не может быть инвертирована. Этот класс систем обычно называют дескрипторными системами, а уравнения - дифференциально-алгебраическими уравнениями (ДАУ).

Последовательная схема RLC

Рассмотрим схему RLC простой серии.

Из Закона о Напряжении Кирхгофа, падение напряжения через схему равен сумме падения напряжения через каждый из его элементов:

$$V_{AC} = V_R+V_L+V_C$$

Из закона тока Кирхгофа:

$$I_{AC}=I_R=I_L=I_C$$

где нижние индексы, $R$$L$и$C$ обозначают сопротивление, индуктивность и емкость соответственно.

$V_{R} =I\left(t\right)R$

$V_L =L\dot{I_L }$ или $\dot{I_L } = \frac{1}{L}V_L$

$V_C =V_{AC} \left(0\right)+\int_0^t I_C \left(\tau \right)d\tau$ или $\dot{V_c } =\frac{1}{C}I_c$

В неявной форме пространства состояний

Моделируйте систему в Simulink с,$R=10\;\Omega$, $L=1\times {10}^{-6} \;H$чтобы$C=1\times {10}^{-4} F$ найти напряжение на резисторе. Чтобы $V_R$использовать блок Descriptor State-Space, система может быть записана в неявную, или дескриптор, форму пространства состояний, как$E\dot{x}=Ax+Bu$ показано ниже.

$$\left\lbrack \begin{array}{ccccc}1 & 0 & 0 & 0 & 0\\0 & 0 & 0 & 0 & 0\\0
& 0 & 0 & 0 & 0\\0 & 0 & 0 & 1 & 0\\0 & 0 & 0 & 0 & 0\end{array}\right\rbrack
\left\lbrack \begin{array}{c}\dot{V_C } \\\dot{V_L } \\\dot{V_R } \\\dot{I_L
} \\\dot{I_{AC} } \end{array}\right\rbrack =\left\lbrack \begin{array}{ccccc}0
& 0 & 0 & \frac{1}{C} & 0\\1 & 1 & 1 & 0 & 0\\0 & 0 & -1 & R & 0\\0 & \frac{1}{L}
& 0 & 0 & 0\\0 & 0 & 0 & 1 & -1\end{array}\right\rbrack \left\lbrack \begin{array}{c}V_C
\\V_L \\V_R \\I_L \\I_{AC} \end{array}\right\rbrack +\left\lbrack \begin{array}{c}0\\-1\\0\\0\\0\end{array}\right\rbrack
V_{AC}$$

где$x = {\left\lbrack \begin{array}{ccccc}V_C &V_L& V_R& I_L& I_{AC}\end{array}\right\rbrack}^T$ является вектором состояния.

Установите$C=\left\lbrack \begin{array}{ccccc}0&0& 1& 0&0\end{array}\right\rbrack$, так как напряжение на резисторе измеряется.

Сравните это с моделированием системы алгебраическим циклом, чтобы найти.$V_R$

Симуляция обеих моделей дает одинаковые результаты. Однако блок Descriptor State-Space позволяет вам сделать более простую блок-схему и избежать алгебраических циклов.

См. также

Концепции Алгебраических циклов

Моделируйте дифференциальные алгебраические уравнения

Документация Simulink

Поддержка