matlab如何求一个N*2的矩阵的任意两点间的距离?
假设有个矩阵C=[1304,2312;3639,1315;4177,2244;3712,1399;3488,1535;3326,1556;3238,1229;4196,1...
假设有个矩阵C=[1304,2312; 3639,1315; 4177,2244; 3712,1399; 3488,1535;3326,1556; 3238,1229; 4196,1004; 4312,790; 4386,570;3007,1970; 2562,1756; 2788,1491; 2381,1676; 1332, 695;3715,1678; 3918,2179; 4061,2370; 3780,2212; 3676,2578;4029,2838; 4263,2931; 3429,1908; 3507,2367; 3394,2643;3439,3201; 2935,3240; 3140,3550; 2545,2357; 2778,2826;2370,2975; ];
该矩阵N*2 每行是一个点的坐标 求matlab程序 计算出任意两点间的距离保存在矩阵C1(N*N)中,并求出最短的两点间的距离L。
求matlab程序 本人matlab菜鸟 只要结果 展开
该矩阵N*2 每行是一个点的坐标 求matlab程序 计算出任意两点间的距离保存在矩阵C1(N*N)中,并求出最短的两点间的距离L。
求matlab程序 本人matlab菜鸟 只要结果 展开
3个回答
展开全部
你这是求任意点的最短距,直接用欧拉公式,或叫欧氏距离求法,
问题矩阵C中有31个点,这样最短距矩阵是d(31*31),求解如下:
%1、求任意点距离
for i=1:31
for j=1:31
C1(i,j)=sqrt((C(i,1)-C(j,1))^2+(C(i,2)-C(j,2))^2);
end
end
%2、因为每个点到自身的距离为0,故求L之前要做处理,将每个点到自身的距离设为无限大:
for i=1:31
C1(i,i)=inf;
end
%3、求最短距L
[L,mini]=min(C1(:))
解得;
任意两点距离
C1(太多了,粘贴后不能发表,请复制上述问题自己求解。)
最短距L=
111.2879
mini = 35(表明在C1矩阵中标号为第2点和第4点之间)
另外,如果求带路径问题,如真实的道路路口距离可以用floyd算法,
function [D,path,min1,path1]=floyd(a,start,terminal)
D=a;n=size(D,1);path=zeros(n,n);
for i=1:n
for j=1:n
if D(i,j)~=inf
path(i,j)=j;
end, end, end
for k=1:n
for i=1:n
for j=1:n
if D(i,k)+D(k,j)<D(i,j)
D(i,j)=D(i,k)+D(k,j);
path(i,j)=path(i,k);
end, end, end,end
if nargin==3
min1=D(start,terminal);
m(1)=start;
i=1;
path1=[ ];
while path(m(i),terminal)~=terminal
k=i+1;
m(k)=path(m(i),terminal);
i=i+1;
end
m(i+1)=terminal;
path1=m;
end
问题矩阵C中有31个点,这样最短距矩阵是d(31*31),求解如下:
%1、求任意点距离
for i=1:31
for j=1:31
C1(i,j)=sqrt((C(i,1)-C(j,1))^2+(C(i,2)-C(j,2))^2);
end
end
%2、因为每个点到自身的距离为0,故求L之前要做处理,将每个点到自身的距离设为无限大:
for i=1:31
C1(i,i)=inf;
end
%3、求最短距L
[L,mini]=min(C1(:))
解得;
任意两点距离
C1(太多了,粘贴后不能发表,请复制上述问题自己求解。)
最短距L=
111.2879
mini = 35(表明在C1矩阵中标号为第2点和第4点之间)
另外,如果求带路径问题,如真实的道路路口距离可以用floyd算法,
function [D,path,min1,path1]=floyd(a,start,terminal)
D=a;n=size(D,1);path=zeros(n,n);
for i=1:n
for j=1:n
if D(i,j)~=inf
path(i,j)=j;
end, end, end
for k=1:n
for i=1:n
for j=1:n
if D(i,k)+D(k,j)<D(i,j)
D(i,j)=D(i,k)+D(k,j);
path(i,j)=path(i,k);
end, end, end,end
if nargin==3
min1=D(start,terminal);
m(1)=start;
i=1;
path1=[ ];
while path(m(i),terminal)~=terminal
k=i+1;
m(k)=path(m(i),terminal);
i=i+1;
end
m(i+1)=terminal;
path1=m;
end
展开全部
计算点之间的距离,距离矩阵distance(N*N),输出到表格中distance.xls。
%% Test_repmat
C=[1304,2312; 3639,1315; 4177,2244; 3712,1399; 3488,1535;3326,1556; 3238,1229;
4196,1004; 4312,790; 4386,570;3007,1970; 2562,1756; 2788,1491; 2381,1676;
1332, 695;3715,1678; 3918,2179; 4061,2370; 3780,2212; 3676,2578;4029,2838;
4263,2931; 3429,1908; 3507,2367; 3394,2643;3439,3201; 2935,3240; 3140,3550;
2545,2357; 2778,2826;2370,2975; ];
p = C;
k = size(p,1);
P_x = repmat(p(:,1),1,k);
P_y = repmat(p(:,2),1,k);
Pt_x = repmat(p(:,1)',k,1);
Pt_y = repmat(p(:,2)',k,1);
distance = sqrt((P_x - Pt_x).^2 + (P_y - Pt_y).^2);
xlswrite('distance.xls',distance)
%by huguohua_sysu@foxmail.com
%% Test_repmat
C=[1304,2312; 3639,1315; 4177,2244; 3712,1399; 3488,1535;3326,1556; 3238,1229;
4196,1004; 4312,790; 4386,570;3007,1970; 2562,1756; 2788,1491; 2381,1676;
1332, 695;3715,1678; 3918,2179; 4061,2370; 3780,2212; 3676,2578;4029,2838;
4263,2931; 3429,1908; 3507,2367; 3394,2643;3439,3201; 2935,3240; 3140,3550;
2545,2357; 2778,2826;2370,2975; ];
p = C;
k = size(p,1);
P_x = repmat(p(:,1),1,k);
P_y = repmat(p(:,2),1,k);
Pt_x = repmat(p(:,1)',k,1);
Pt_y = repmat(p(:,2)',k,1);
distance = sqrt((P_x - Pt_x).^2 + (P_y - Pt_y).^2);
xlswrite('distance.xls',distance)
%by huguohua_sysu@foxmail.com
本回答被网友采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
展开全部
建议你查一下Dijkstra算法和Floyd算法或者Kruskal或Prim算法,要想学好matlab还得多自学点东西~~~
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询