数模论坛

 找回密码
 注-册-帐-号
搜索
热搜: 活动 交友 discuz
查看: 1746|回复: 5

高手来帮我改一下Matlab程序

[复制链接]
发表于 2004-8-27 06:25:33 | 显示全部楼层 |阅读模式
<> 帮我看一下,谢谢!</P>

<>clear all
global distMatrix
t=[1304 2312;3639 1315;4177 2244;3712 1399;3488 1535;3326 1556;...
      3238 1229;4196 1044;4312 790;4386 570;3007 1970;2562 1756;...
      2788 1491;2381 1676; 1332 695;3715 1678;3918 2179;4061 2376;...
      3780 2212;3676 2578;4029 2838;4263 2931;3429 1908;3507 2376;...
      3394 2643;3439 3201;2935 3240;3140 3550;2545 2357;2778 2826;2370 2975];...
sz=size(t,1);
distMatrix=distsz(t,t);
xFns='cyclicXover unniformXover partmapXover orderbasedXover'
xFns=[xFns,'singleptXover linerorderXover'];Opts=[2;2;2;2;2;2];%2;2;2;2;2;2;2];
mFns='inversionMutation adjswapMutation'
swapMutation threewapMutation';
mOpts=[2;2;2;2;2];
termFns='maxGeomSelect';
termOps=[100];
selectFn='normGeomSelect';
selectOps=[0.08];
evalFn='tspEval';
evalOps=[];
bounds=[sz];
gaOps=[le-6 1 1];
startPop=initalixeoga(80,bounds,'tspEval',[le=6 1]);
[x endPop bestPop trace]=ga(bounds,evalFn,evalOps,...
    startPop,gaPots,termFns,termOps,selectFn,...
    selectOps,xFns,xOpts,mFns,mOpts];
bestPop
trace
plot(trace(;,1),trace(;,2));
hold on
plot(trace(;,1),trace(;,3));
figue(2)
clf
A=ones(sz,sz);
A=xor(triu(A),tril(A));
[xg yg]=gplot(A,t);
clf
h=gca;
hold on
ap=x;
plot(t(x(1:sz),1),t(x(1:sz),2,),'r-')
plot(t([x(1):x(sz)],1),t([x(1):x(sz)],2),'r-')
plot(xg,yg,'b.','MarkerSize',24);</P>
发表于 2004-8-27 06:32:33 | 显示全部楼层
你这个程序想干什么呀?说一下比较好改
发表于 2004-8-27 07:15:06 | 显示全部楼层
<>??</P><>what do you want to do</P>
 楼主| 发表于 2004-8-27 07:53:09 | 显示全部楼层
<>是城市距离问题呀!</P>
发表于 2004-8-30 19:26:10 | 显示全部楼层
你不如给出题目,这样应该会好改点啊!
发表于 2004-8-30 20:42:12 | 显示全部楼层
<>clear all
global distMatrix
t=[1304 2312;3639 1315;4177 2244;3712 1399;3488 1535;3326 1556;...
      3238 1229;4196 1044;4312 790;4386 570;3007 1970;2562 1756;...
      2788 1491;2381 1676; 1332 695;3715 1678;3918 2179;4061 2376;...
      3780 2212;3676 2578;4029 2838;4263 2931;3429 1908;3507 2376;...
      3394 2643;3439 3201;2935 3240;3140 3550;2545 2357;2778 2826;2370 2975];...
sz=size(t,1);
<FONT color=#f73809>distMatrix=distsz(t,t);   %%distsz这个函数错了。</FONT>
xFns='cyclicXover unniformXover partmapXover orderbasedXover'
xFns=[xFns,'singleptXover linerorderXover'];Opts=[2;2;2;2;2;2];%2;2;2;2;2;2;2];
mFns='inversionMutation adjswapMutation'
swapMutation threewapMutation';
mOpts=[2;2;2;2;2];
termFns='maxGeomSelect';
termOps=[100];
selectFn='normGeomSelect';
selectOps=[0.08];
evalFn='tspEval';
evalOps=[];
bounds=[sz];
gaOps=[le-6 1 1];
startPop=initalixeoga(80,bounds,'tspEval',[le=6 1]);
[x endPop bestPop trace]=ga(bounds,evalFn,evalOps,...
    startPop,gaPots,termFns,termOps,selectFn,...
    selectOps,xFns,xOpts,mFns,mOpts];
bestPop
trace
plot(trace(;,1),trace(;,2));
hold on
plot(trace(;,1),trace(;,3));
figue(2)
clf
A=ones(sz,sz);
A=xor(triu(A),tril(A));
[xg yg]=gplot(A,t);
clf
h=gca;
hold on
ap=x;
plot(t(x(1:sz),1),t(x(1:sz),2,),'r-')
plot(t([x(1):x(sz)],1),t([x(1):x(sz)],2),'r-')
plot(xg,yg,'b.','MarkerSize',24);</P>
您需要登录后才可以回帖 登录 | 注-册-帐-号

本版积分规则

小黑屋|手机版|Archiver|数学建模网 ( 湘ICP备11011602号 )

GMT+8, 2024-11-28 21:37 , Processed in 0.068126 second(s), 18 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表