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

Z(i,i)=newshema(k+1); e(i)=newshema(k+2);
J(i)=newshema(k+3); k=k+5;
end;
//определение числа узлов. Число узлов должно быть равно наибольшему номеру узла
//наибольший номер узла принимается за общий узел
yzlov1=max (ta);yzlov2=max (he);yzlov=max (yzlov1,yzlov2);
//создание направленного графа g по введенной информации
g = make_graph ('cep1', 1, yzlov, ta', he');
//создание ненаправленного графа g0 по введенной информации
g0 = make_graph ('cep1', 0, yzlov, ta', he');
g('node_x')=[200 100 150 250 ]; //ввод х-координат узлов
g('node_y')=[200 200 100 100 ]; //ввод у-координат узлов
// show _ graph (g);
yzlov = node_number (g); //определение числа узлов введенного графа
//формирование матрицы инциденций "узел-ветвь" для введенного графа g
A=full (graph_2_mat (g,'node-arc'));
//формирование дерева графа
derevo= min _ weight _ tree(yzlov,g0);
//формирование расширенной узловой матрицы
at=zeros(arc,arc); //объявление матрицы At;
at(1:yzlov-1, 1:arc)=A(1:yzlov-1, 1:arc); //заполнение узловой части At;
hord=setdiff (1:arc,derevo); //формирование массива хорд;
size_hord=size(hord); //определения числа хорд;
for ii=1:size_hord(2) //заполнение контурной
at(yzlov-1+ii, hord(ii))=1; //части расширенной узловой
end; //матрицы;
c=at^(- 1); //вычисление контурной матрицы;
er=c'*e; Jr=at*J; //вычисление векторов источников схемы;
Zr=c'*Z*c; //вычисление матрицы сопротивлений схемы;
Yk=Zr(yzlov:arc, yzlov:arc)^(- 1); //вычисление матрицы контурных проводимостей
//вычисление матриц взаимных сопротивлений узловых и контурных координат
zok=Zr(1:yzlov-1,yzlov:arc); zko=Zr(yzlov:arc,1:yzlov-1);
//вычисление узловых матриц сопротивлений
zo =Zr(1:yzlov-1,1:yzlov-1); zor=zo-zok*Yk*zko;
//вычисление контурных и узловых источников напряжений
ek=er(yzlov:arc); eo=er(1:yzlov-1);
//вычисление контурных и узловых источников тока
Jk=Jr(yzlov:arc); Jo=Jr(1:yzlov-1);
//вычисление контурных токов и узловых потенциалов
ik=Yk*ek−Yk*zko*Jo-Jk; Uo=zor*Jo+zok*Yk*ek−eo;
//вычисление токов и напряжений ветвей
iv= c(1:arc,yzlov:arc)* ik; Uv= at(1:yzlov−1,1:arc)'* Uo;
//вычисление суммарной мощности источников схемы
P= J'* Uv // P= e'* iv
Для того, что бы вывести в диалоговое окно нужный вектор или матрицу,
достаточно заменить в конце соответствующего выражения точку
с запятой на запятую.
В приведенном алгоритме можно выделить следующие части:
(1) Ввод данных;
(2) Перевод введенных данных во внутреннюю форму;
(3) Вычисления;
77