Расчет переходного процеса осуществляется оператором

Z2 = rkfixed(x0,0,t_п, 150,D2)

где t_п – время переходного процесса.

Для построения графика необходимо из матрицы Z2 выделить столбцы

t = Z2,^<1>, i₂ = Z2^<2>, u_C = Z^<3>.

Расчет токов и напряжений

Пусть по заданию кроме тока в индуктивности и напряжения на емкости требуется найти токи i₈, i₅, i₃. Образуем вектор

Y = [i₈ i₅ i₃]^T.

Из уравнений Кирхгофа находим

- i₈ = - i₂ – u₈/R₂ + e/R₂ + i

- i₅ = - u₈/R₂ + e/R₂

- i₃ = - i₂ + i

или в матричной форме

i₈ - 1 1/R₂ i₂ 1/R₂ 1 e(t)

Y = i₅ = 0 -1/R₂ u₈ + 1/R₂ 0 i(t)

I₃ - 1 0 0 1

Отсюда следует, что

- 1 - 1/R₂ - 1/R₂ 1

C = 0 - 1/R₂ , D = 1/R₂ 0

- 1 0 0 1

Для расчета искомых токов составляем программу

k = 1..150

t_k = Z2_k,1

i_k = Z2_k,2

u_Ck = Z2_k,3

W_{k,1 =} C₁₁*i_k + C₁₂ * u_Ck + D₁₁*E_msin(2512t + y_e) + D₁₂*I_msin(2512t +y_i)

W_k,2 = C₂₁*i_k + C₂₂ * u_Ck + D₂₁*E_msin(2512t + y_e) + D₂₂*I_msin(2512t +y_i)

W_k,3 = C₃₁*i_k + C₃₁* u_Ck + D₃₁* E_msin(2512t + y_e) + D₃₂*I_msin(2512t +y_i)

W_k,4 = t_k

Для построения графиков необходимо образовать ранжированные вектора с помощью следующих операторов

i₈ = W^<1>, i₅ = W^<2>, i₃ = W^<3>, t = W^<4>.

Приложение 1