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

he=[3 2 2 4 1 1 ]; // конечные узлы ветвей, по порядку, начиная с 1-й ветви
maxv=max(ta); // определение максимального номера узла
g = make_graph('cep1',1,maxv,ta,he); //создание направленного графа
g1= make_graph('cep2',0,maxv,ta,he); //создание ненаправленного графа
g('node_x') = [200 100 150 250 ]; //ввод х-координат узлов
g('node_y') = [200 200 100 100 ]; //ввод у-координат узлов
show_graph(g); //отображение графа
a=graph_2_mat(g,'node−arc'); //формирование ненулевых элементов узловой матрицы
A=full(a); //формирование полной формы узловой матрицы
node = node_number(g); // определение числа узлов графа
t=min_weight_tree(node,g1); //формирование дерева графа
arc = arc_number(g); //определение числа ветвей графа
at(1:node−1, 1:arc)=A(1:node−1, : ); // узловая часть расширенной матрицы Аt
hord=setdiff(1:arc, t); //формирование списка хорд
size_hord = size(hord); // определение числа хорд
for i=1:size_hord(2) at(node−1+i,hord(i))=1; end; //формирование полной матрицы Аt
at //вывод матрицы Аt в диалоговое окно Scilab
c=at^(-1) //вычисление и вывод матрицы С в диалоговое окно Scilab
Для введенного графа сформированные матрицы Аt и С получились следующими:
2
1
3
2
3
1
4
5
6
4
Аt=
! 0. 0. 1. 0. - 1. - 1. !
! 1. - 1. - 1. 0. 0. 0. !
! - 1. 1. 0. 1. 0. 0. !
! 0. 1. 0. 0. 0. 0. !
! 0. 0. 0. 1. 0. 0. !
! 0. 0. 0. 0. 0. 1. !
C=
! 0. 0. - 1. 1. 1. 0. !
! 0. 0. 0. 1. 0. 0. !
! 0. - 1. - 1. 0. 1. 0. !
! 0. 0. 0. 0. 1. 0 !
! - 1. - 1. - 1. 0. 1. - 1. !
! 0. 0. 0. 0. 0. 1. !
Как видно из полученных результатов, программа сформировала другой
состав матриц. Это связано с тем, что количество вариантов построения
дерева и соответственно контуров, может быть достаточно
большим. Выбор того или иного варианта определяется особенностями
алгоритма, в данном случае встроенной функцией построения дерева
пакета SciLab: min_weight_tree. В любом случае проверить правильность
сформированных матриц можно косвенно, результат умножения С на Аt
должен дать единичную матрицу: С·Аt=1.
Построить дерево можно отдельным алгоритмом, приведенным ниже.
Данный SciLab-алгоритм является модификацией фортрановской программы,
приведенной в [12]. Исходными данными являются два списка:
ta – список начальных узлов ветвей и he – список конечных узлов ветвей.
clear ;
ta= [2 3 1 3 4 4 ]; // начальные узлы ветвей, по порядку, начиная с 1-й ветви
he=[3 2 2 4 1 1 ]; // конечные узлы ветвей, по порядку, начиная с 1-й ветви
//определение размеров массивов
razmer = size( ta); JKV= razmer (2);
max1=max( ta); max2=max(he); KP=max(max1,max2);
IW=KP+2 *JKV; NEXT(IW)=0;
64