Смоделируйте ряд схема 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 позволяет вам делать более простую блок-схему и избегать алгебраических циклов.

Смотрите также

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

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

Для просмотра документации необходимо авторизоваться на сайте