Изучение курса ТЭЦ с использованием систем автоматизации инженерных расчетов. Регеда В.В - 88 стр.

                                                ⎡1    1⎤
где Δ – определитель системы, который равен ⎢       ⎥ – для случая
                                            ⎣p1 p 2 ⎦
                                                              ⎡ 1 0⎤
разных корней характеристического уравнения или               ⎢p 1⎥ –
                                                              ⎣ 1 ⎦
для случая кратных корней характеристического уравнения;
        ⎡A ⎤
    A = ⎢ 1 ⎥ – матрица-столбец постоянных интегрирования;
        ⎣A 2 ⎦
        ⎡ i L (0 + ) ⎤ ⎡0,1⎤
    С = ⎢ di L (0 + ) ⎥ = ⎢ ⎥ − матрица-столбец свободных членов.
        ⎢             ⎥    0
        ⎣ dt ⎦ ⎣ ⎦
    Решение системы (2.26) можно представить в матричной форме:
                             A = Δ−1 ⋅ C.
    Полученная система может быть легко решена с помощью про-
граммы MathCAD.
    На рис. 2.37 приведена программа, написанная в среде MathCAD,
реализующая описанный выше алгоритм для случая, когда сопро-
тивление R2 меньше критического сопротивления цепи в пять раз.
    В программе приняты следующие обозначения:
    – критическое сопротивление цепи – RKR;
    – ток через индуктивность до коммутации – iL0;
    – напряжение на емкости до коммутации – uC0;
    – принужденное значение тока через индуктивность iLpr.
    При вычислении главного определителя Δ системы (2.26) и
при задании аналитического выражения для iL(t) используются
программные блоки, которые создаются с помощью команды
Add Line.
    В первом программном блоке при выполнении условия р1 ≠ р2
                           ⎡1    1⎤             ⎡1   0⎤
главный определитель Δ = ⎢      ⎥ , иначе Δ = ⎢p 1⎥ .
                           p
                         ⎣ 1 p 2⎦             ⎣ 1 ⎦
    Во втором программном блоке при выполнении условия р1 ≠ р2
значения функции iL(t) в А вычисляются в соответствии с выраже-
нием (2.21), иначе – в соответствии с выражением (2.22). Значения
аргумента функции iL(t) в секундах задаются с помощью
ранжированной переменной t, которая изменяется в диапазоне от
                                87