Расчёт переходных процессов методом синтетических схем, страница 3

Таблица 8. Данные полученные при шаге интегрирования

Номер ветви	I      А	U      В
1	1,37518812027379	1,37518812027379
2	1,37482322286120	4,12446966858360
3	1,37484215812389	4,12452647437168
4	1,66668074853078	8,33334423827692
5	1,66668074853078	1,66668074853078
6	6,02629689453238e-05	0,000602629689453238
7	1,37518812027378	0,375188120273784
8	1,93538444275276e-05	8,24962375945243
9	-1,89352626946970e-05	4,12515409086884
10	-1,37512743872311	-0,000627616497150374
11	1,66662048556184	-2,49868076971360e-05
12	0,000345543568160522	8,24899614295528

Таблица 8. Данные полученные в результате расчётов.

Номер ветви	I      А	U      В
1	1,375	1,375
2	1,375	4,125
3	1,375	4,125
4	1,667	8,33
5	1,667	1,375
6	0	0
7	1,375	0,375
8	0	4,125
9	0	4,125
10	1,667	8,25
11	0	0
12	0	8,25

Приложение 1. Расчёт после коммутационных процессов в MathCad.

Приложение2. Код программы для расчёта переходного процесса в MatLab.

function [Ukon]= perpr(NET);

R=NET(:,4);                 %Выделение матрицы сопротивлений.

L=NET(:,5);                 %Выделение матрицы индуктивности.

C=NET(:,6);                 %Выделение матрицы ёмкости.

E=NET(:,7);                 %Выделение матрицы ЕДС.

J=NET(:,8);                 %Выделение матрицы источников тока.

Nach=NET(:,9);              %Выделение матрицы начальных данных.

Vetvi=NET(:,1);             %Выделение матрицы ветвей.

p=max(Vetvi);               %Расчёт количества ветвей.

Uzli=NET(:,2);              %Выделение матрицы узлов.

q=max(Uzli);                %Расчёт количества узлов.

n=p-(q-1);                  %Расчёт числа  связей.

D1=length(C);           

C1=zeros(D1,1);

for k = 1:D1                %Обращение всех inf из матрицы С в ноль.

    if C(k,1)~=inf          %Необходимо для поиска максимального

        C1(k,1)=C(k,1);     %значения ёмкости.

    else

        C1(k,1)=0 ;

    end

end

C1max=max(C1);              % Максимальное значение ёмкости.

Rmax=max(R);                % Максимальное значение сопротивления.

Lmax=max(L);                % Максимальное значение индуктивности.

Cmin=min(C);                % Минимальное значение ёмкости.

R1=zeros(D1,1);