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

2.3.6. Алгоритм формирования и расчета уравнений резистивной
схемы.
Далее приводится алгоритм, позволяющий автоматически сформировать
необходимые векторы и матрицы по методу контурных токов и
рассчитать токи, напряжения и мощность ветвей резистивной схемы.
Исходными данными является список с информацией о каждой ветви.
Каждая ветвь должна быть представлена в виде обобщенной ветви с последовательно
включенным источником напряжения и параллельно
включенным источником тока. Список оформляется оператором list по
правилам пакета SciLab. В состав списка входит имя ветви, номер начального
и номер конечного узла ветви, номинал ветви, значение источника
напряжения ветви, значение источника тока ветви. Имя в данном
алгоритме не обрабатывается и вводится для удобства ввода. Фактически
ветви идентифицируются номерами, соответствующими порядку их
ввода. Узлы должны быть пронумерованы последовательно, начиная с
единицы. В алгоритме наибольший номер узла принимается за общий
узел.
Алгоритм переформировывает исходные данные, представленные
списком соединений в массивы векторов источников тока и напряжения
ветвей и матрицы сопротивлений ветвей. Далее формируются матрицы
соединений и затем матрицы, необходимые для расчета токов и напряжений
ветвей (в соответствии с формулами (2.3.5)-(2.3.8) ).
В алгоритме применены следующие встроенные функции SciLab:
size — определение размера списка или массива, zeros — обнуление
массива, max — определение максимального числа из списка чисел,
make_graph — создание графа, node_number — определение числа узлов
графа , show _graph — отображение графа на экране, graph_2_mat — формирование
упакованной матрицы инциденций графа, full — преобразование
упакованной матрицы в полную матрицу, min _weight _tree — формирование
дерева графа, setdiff — вычитание множеств.
clear ( ); //очистка данных.
// Ввод списка соединений:
// имя нач. кон. парам. e J
newshema= list ('r1', [2 3], 2, 0, 0,... //многоточие означает перенос строки
'r2', [3 2], 3, 0, 0.33333,... //в столбце парам. задаются номиналы
'r3', [1 2], 4, 0, 0,... // ветвей
'r4', [3 4], 5, 0, 0,...
'r5', [4 1], 6, 0, 0,...
'r6', [4 1], 7, 0, 0);
//определение числа ветвей - arc, объявление массива сопротивлений - Z
arc=size(newshema)/5; Z=zeros(arc,arc);
//Извлечение информации из списка соединений и формирование новой информации:
//ta-массив начальных узлов, he-массив конечных узлов, Z-матрица сопротивлений,
// e-вектор источников напряжений ветвей, J-вектор источников тока ветвей
k=2;
for i=1:arc
ta(i)=newshema(k)(1); he(i)=newshema(k)(2);
76