电工杯B题小代码分享

第一问的遍历所有地点的最短路径

代码一：Lingo

model:
n=14; %总共个点14 A=zeros(n,n);
for i=1:nfor j=1:nif(i==j) A(i,j)=0;else A(i,j)=100000;endend
end
A(1,5)=54;A(1,7)=55;A(1,11)=26;A(2,3)=56;A(2,5)=18;A(3,5)=44;A(4,6)=28;
A(5,6)=51;A(5,7)=34;A(5,8)=56;A(5,9)=48;A(6,9)=27;A(6,10)=42;A(7,8)=36;
A(7,12)=38;A(8,9)=29;A(8,12)=33;A(9,10)=61;A(9,12)=29;A(9,13)=42;A(9,14)=36;
A(10,14)=25;A(11,12)=24;A(13,14)=47;
for j=1:n
for i=1:j-1 A(j,i)=A(i,j); %使对称 end
end
[m,n]=size(A);
B=zeros(m,n);
B=A;
%利用 Floyd 算法计算最短距离矩阵 for k=1:nfor i=1 :nfor j=1:nt=B(i,k)+B(k,j);if t<B(i,j) B(i,j)=t; endendend
end
%输出距离矩阵 fid=fopen('distance.txt','w');for i=1:nfor j=1:n
fprintf(fid,'%4d ',B(i,j));endfprintf(fid,'\n');endfclose(fid);
%初始数据%
D=[0 72 98 133 54 54 55 83 79 140 26 50 121 11572 0 56 97 18 69 52 74 66 111 98 90 108 10298 56 0 123 44 95 78 100 92 137 124 116 134 128133 97 123 0 79 28 113 84 55 70 108 84 97 9154 18 44 79 0 51 34 56 48 93 80 72 90 84105 69 95 28 51 0 85 56 27 42 80 56 69 6355 52 78 113 34 85 0 36 65 126 62 38 107 10183 74 100 84 56 56 36 0 29 90 84 33 71 6579 66 92 55 48 27 65 29 0 61 53 29 42 36140 111 137 70 93 42 126 90 61 0 114 90 72 2526 98 124 108 80 80 62 57 53 114 0 24 95 8950 90 116 84 72 56 38 33 29 90 24 0 71 65121 108 134 97 90 69 107 71 42 72 95 71 0 47115 102 128 91 84 63 101 65 36 25 89 65 47 0];
L=size(D,1); %求出有几行即有几个城市D c=[9 1 2 3 4 5:6 7 8 10 11 12 13 14 9]; %选择初始圈 %改良圈%
for k=1:L
flag=0; %退出标志 for i=1:L-2for j=i+2:Lif D(c(i),c(j))+D(c(i+1),c(j+1))<D(c(i),c(i+1))+D(c(j),c(j+1))
c(i+1:j)=c(j:-1:i+1); %替换圈原本的路径 flag=flag+1;
end %即i+1-之间的访问顺序没变j endendif flag==0
long=0;for n=1:Llong =long+D(c(n),c(n+1));
end
circle=c;
break; end
end end

代码二：MATLAB

function [xl,l]=ljpj(d,q,zx) %dqzx %q
if length(q(:,1))>1
q=q';
end
k=length(q); %qzx
for i=1:k
if q(i)==zx
q(i:k-1)=q(i+1:k);
k=k-1;
end
end
q=q(1:k);%
y=zx;
for j=1:k
sum0=inf; %zu for i=1:k-j+1
if d(y,q(i))<sum0
sum0=d(y,q(i));
z=i;
end
end
u(j)=q(z);
y=q(z); %qu
for i=1:k-j+1
if i==z
if i==1
q=q(2:end);
else
q=[q(1:i-1) q(i+1:end)];
end
end
end
end % xl=[zx u zx]';
i=1;
while xl(i+1)==zx xl=xl(2:end);
end
l=0;
for i=1:length(xl)-1 l=l+d(xl(i),xl(i+1));
end

上述两种求解电工杯B题第一问的代码，建议保存，以后求类似遍历所有地点的问题可以直接替换距离矩阵使用。

关注公众号，学习数学建模。

电工杯B题小代码分享相关推荐

2021电工杯B题建模思路代码完整版
2021电工杯B题思路和代码完整版 2021"中国电机工程学会杯"全国大学生电工数学建模竞赛B题思路代码完整版+参考文献,内容如下,可分享-
2023年电工杯B题详解
2023电工杯B题思路讲解,以及模型讲解,都会用视频he文件的方式发给你们! 资料主要包括模型部分及论文框架使用方法:模型由我们建立,大家根据视频讲解进行理解论文框架是论文的主体,文字的描述千变 ...
2021年电工杯B题光伏建筑一体化板块指数发展趋势分析及预测
电工杯B题光伏建筑一体化板块指数发展趋势分析及预测目录赛题简述具体分析赛后感想赛题简述电工杯这次有两题,a题是高铁牵引供电系统运行数据分析及等值建模,这题一看就是比较技术型的:然后这次比赛 ...
第一届华数杯A题完整思路分享
A题完整思路分享第一问: 要想问题做的好,美图肯定少不了.这道题,建模时不但要用数学公式把模型说清楚,还要一些图来帮助别人理解,让别人一眼就看懂你写的东西,话不多说,我们先放一张类似下面的图,记得说 ...
2019年‘泰迪杯’数据分析职业技能大赛A题——个人代码分享
目录题目任务 1 数据预处理与统计任务 2 数据分析与可视化代码展示任务一任务二题目任务 1 数据预处理与统计任务 1.1 对数据作必要的预处理,在报告中列出处理步骤,将处理后的结 ...
2022电工杯A题B题高质量保奖代码思路合集
目录 A题 B题 A题待研究系统包含火电.风电.储能和负荷,火电机组 3 台.装机容量 1050MW: 某日风电.负荷归一化功率(1.0p.u.风电对应其装机容量,1.0p.u.负荷对应最大负荷功率 ...
2022亚太杯C题思路代码分析
C题就是数学比较开放的题目了,属于一个数据分析类题目,跟前两年的华为杯差不多,考察的也是全球变暖问题.更多内容文末名片查看问题1.你同意有关全球气温的说法吗?使用2022_APMCM_C_Data. ...
吉首大学新星杯L题小李子的老年生活【容斥好题(只是问题比较裸)】
传送门 // 题意就是给定一个n, 求1 - n 中与n互质的数的平方和. // 思路: 问题我们可以进行转换, 就是我们得到1 - n的所有数的平方和, 然后我们要做的就是减去n的素数因子(不是素数 ...
【Python | Networks库详解】最佳实践（22年电工杯B题之路径可视化）
【2023年电工杯数学建模竞赛】选题分析+A题B题完整思路+代码分享
2023年电工杯B题(附带ChatGpt思路)思路已更新,请点击一下链接 [2023年电工杯数学建模竞赛B题人工智能对大学生学习影响的评价]完整思路分析+完整代码+(附带ChatGpt思路) 1.竞赛 ...

电工杯B题小代码分享

电工杯B题小代码分享相关推荐

最新文章

热门文章