?????????????? ?????? ? ????????? ?????????? ??????? ?????

формируется массив номеров ветвей, соответствующих емкостным и
индуктивным ветвям.
В третьем блоке выполняется формирование топологических матриц.
Для этой цели используется функция formc, написанная на SciLab
и приведенная далее за основным алгоритмом. Особенность функции
состоит в том, что она может сформировать топологические матрицы
для несвязных схем, то есть схем, в состав которых входят топологически
изолированные подсхемы.
В четвертом блоке выполняется вычисление матриц сопротивлений в
координатах схемы, то есть в контурном и узловом базисах. Поскольку
расчет схемы соответствует методу контурных сопротивлений, то в этом
блоке вычисляется матрица контурных проводимостей.
В пятом блоке формируется цикл расчета переходного процесса. Так
как схема линейная, то этот цикл организуется после вычисления матриц
сопротивлений, проводимостей и топологических матриц. При формировании
цикла производится проверка на условия окончания цикла.
Если параметры цикла заданы неверно, например, число расчетных точек
оказалось меньше или равно начальному значению переменной цикла,
то цикл выполнятся не будет.
В шестом, седьмом и восьмом блоках, выполняется расчет схемы по
соотношениям, приведенным в разделе 2.3.3. В восьмом блоке для инерционных
ветвей выполняется перерасчет источников тока и напряжения
в соответствии с методом Эйлера 1-го порядка.
clear ();
//Ввод расчетного шага dt и времени расчета T
dt=0.000001;T=0.0002;
//Вычиcление числа расчетных точек NN
NN=round (T/dt); //вычисляет наименьшее целое от деления
//задание массива времени
t=[0:dt:T]'; //одинарная кавычка означает транспонирование
/ /Ввод списка соединений
// имя нач. кон. парам. e J
newshema=list ( ' r1', [2 3], 2, 1, 0,...
'c1', [3 2], 10.0e-6, -0.5, 0.,...
'r3', [1 2], 4, 0, 0,...
'L1', [3 4], 0.10e-3, -0.00, -1.0,...
'r5', [4 1], 6, 0, 0,...
'r6', [4 2], 7, 0, 0);
//определение числа ветвей-arc, объявл массива сопротивлений Z
arc=size(newshema)/5; Z=zeros(arc,arc);
/ / Извлечение информации из списка соединений
k=2; kc=0; kL=0;
for i=1:arc
ta(i)=newshema(k)(1); he(i)=newshema(k)(2);
sim(i)= convstr (newshema(k-1),'u');
a=stripblanks (sim(i)); //удаление пробелов
a=part ([a],[1]); //определение 1-го символа в имени ветви
select a
80