求欧拉距离 matlab,欧拉回路的求解matlab.doc

欧拉回路的求解matlab

欧拉回路的求解

左图是一个井田图，由于2、3、5、8、9、12、14、15几个点都是奇数连线，故不存在欧拉回路。而右图增加几条连线后，该图就存在欧拉回路。

假设点1和点2 之间的连线消失，

建立数学模型把右图的拓扑关系(并考虑a中连线消失的因素)表达出来

理解fleury算法，并计算一条欧拉迹，使得该欧拉迹从点1出发，经过b中的每一条边，最终达到点2

使用plot命令把该欧拉迹显示出来。这个动画过程可以用一个for循环语句实现，如下。其中 pos是个2x16的矩阵，2行分别代表x/y轴坐标，每一列表示每个点的坐标，共16个点；另外，T是个2xN的矩阵，每一列表示一条边从 T(1，i) 点 到 T(2，i)点。fleury算法的目的就是要产生这样一个 T 矩阵。

for i …

draw_arrow(pos(:,T(1,i))',pos(:,T(2,i))',0.5)

pause;

end

以下是这段动画的其中几个截图

clear all

hold off

A=zeros(16);

for i=1:16

%A(i,i)=1/2;

if i+1<=16 && mod(i,4)~=0

A(i,i+1)=1;

end

if i+4<=16

A(i,i+4)=1;

end

end

A(2,5)=1;

A(3,8)=1;

A(9,14)=1;

A(12,15)=1;

A=A+A';

[T3 c3]=fleury3(A);

pos3(1:2,1)=0;

for i=1:16

if i==1, pos3(1:2,i)=0;

else

if mod(i-1,4)~=0

pos3(1,i)=pos3(1,i-1)+1;

pos3(2,i)=pos3(2,i-1);

else

pos3(1,i)=pos3(1,i-4);

pos3(2,i)=pos3(2,i-4)-1;

end

end

end

figure (1), title('3rd Question')

for i=1:16

for j=i:16

if A(i,j)==1,

plot([pos3(1,i),pos3(1,j)],[pos3(2,i),pos3(2,j)]); hold on,

end

end

end

for i=2:size(T3,2)

draw_arrow(pos3(:,T3(1,i))',pos3(:,T3(2,i))',0.5)

x = mean(pos3(1,T3(:,i)));

y = mean(pos3(2,T3(:,i)));

text(x,y,num2str(i),'FontSize',18);

pause;

end

pause off;

hold off

function [T,sleds]=fleury3(A)

[m,n]=size(A);

if m~=n

fprintf('there is sth wrong.\n');

return;

end

temp=sum(A,1);

tteds=sum(temp); % total nos edges

sleds=0; %selected edges

mtr = A;

startp = 1;

eulerPath = startp;

while sleds~=tteds % 当sleds = tteds 时，表示所有的边已经进入到 eulerPath里面了

listNp = find(mtr(startp,:)==1); %列出与 startp 相通的点

nosNgbr = length(listNp); %计算 上述 点的个数

if nosNgbr ==1, %没有其他选择情况下

nextp = listNp(1);

else

for i=1:length(listNp) %依次判断是否是割边

flag = isGeBian(mtr,startp,listNp(i));