分享下数模2015年A题的Matlab程序

鑫玺

全部资料（压缩包）下载：

程序.rar (24.68 KB, 下载次数: 0)

2017-4-21 16:09 上传

点击文件名下载附件

cumcm2012B附件6_可参考的相关概念.doc (77.5 KB, 下载次数: 0)

2017-4-21 16:09 上传

点击文件名下载附件

哈哈！三天终于熬过头了，今天也休息过来了，来做个分享！
国赛本科两个题，我们考虑了自己的劣势，对于创新类的题做起来比较恼火，最终选择了A题。

A题入手需要有一些相关的地理知识，例如方位角，时角等等。其实这些内容在2012年B题的附件里都有解释，暑假期间老师让我们做了题，当时比较认真的做。第一问上午就解出来了，后来一直卡在后面，个人matlab编程水平不是很高超，所以队友列出的非线性方程解得不是很好。拖了队友的后腿。最后一问在学长的点拨下做了出来。论文尊重队友意见，就不分享了。程序还是可以共享的。

做完了就是做完了！不考虑结果了。
上图便于大家理解吧，空间想象你能给力不好的可以靠这个图来思考的







下面是部分程序预览（完整版本请下载本帖附件）：

1. function [date,WE,NS]=fmain()  
   
2. 
   
3. %参数赋初值
   
4. maxgen=90;            %进化代数，即迭代次数
   
5. sizepop=200;           %种群规模
   
6. pcross=0.6;            %交叉概率选择
   
7. pmutation=0.2;         %变异概率选择
   
8. lenchrom=[1 1 1];        %每个变量的字符长度，如果是浮点，则长度都变为1
   
9. bound=[1 365;2 10;0 90];
   
10. c1=1.5;
    
11. c2=1.5;
    
12. popmax=[365;10;90];
    
13. popmin=[1;2;0];
    
14. Vmax=[60,2,15];
    
15. Vmin=[-60,-2,-15];
    
16. 
    
17. %附件4的影子的测量长度
    
18. LL1=textread('D.txt');
    
19. 
    
20. 
    
21. %种群初始化
    
22. for i=1:sizepop
    
23.     chrom(i,:)=Code(lenchrom,bound);
    
24.     V(i,:)=5*rand(1,length(lenchrom));
    
25.     x=chrom(i,:);
    
26.     fitness(i,1)=fun(x,LL1);
    
27. end
    
28. 
    
29. [bestfitness,bestindex]=min(fitness);
    
30. zbest=chrom(bestindex,:);   %群体极值位置
    
31. gbest=chrom;                %个体位置
    
32. fitnessgbest=fitness;       %个体适应度值
    
33. fitnesszbest=bestfitness;   %群体极值
    
34. %开始迭代
    
35. for i=1:maxgen
    
36.     A=[fitness,chrom];
    
37.     A=sortrows(A,1);         %按fitness小到大排序
    
38.     A=A(1:ceil(sizepop/2),:);%保留前面一半
    
39.     fitness=A(:,1);
    
40.     chrom=A(:,2:end);
    
41.     gbest=chrom;
    
42.     %采用pso算法进行更新
    
43.     for j=1:ceil(sizepop/2)
    
44.         %速度更新
    
45.         V(j,:)=V(j,:)+c1*rand*(gbest(j,:)-chrom(j,:))+c2*rand*(zbest-chrom(j,:));
    
46.         if  V(j,1)>Vmax(1)
    
47.             V(j,1)=Vmax(1);
    
48.         end
    
49.         if  V(j,2)>Vmax(2)
    
50.             V(j,2)=Vmax(2);
    
51.         end
    
52.         if  V(j,3)>Vmax(3)
    
53.             V(j,3)=Vmax(3);
    
54.         end
    
55.         if  V(j,1)<vmin(1)
    
56.             V(j,1)=Vmin(1);
    
57.         end
    
58.          if  V(j,2)<vmin(2)
    
59.              V(j,2)=Vmin(2);
    
60.          end
    
61.          if  V(j,3)<vmin(3)
    
62.              V(j,3)=Vmin(3);
    
63.         end
    
64.      
    
65.         %粒子更新
    
66.         chrom(j,:)=chrom(j,:)+1*V(j,:);
    
67.         if chrom(j,1)>popmax(1)
    
68.            chrom(j,1)=popmax(1);
    
69.         end
    
70.         if chrom(j,2)>popmax(2)
    
71.            chrom(j,2)=popmax(2);
    
72.         end
    
73.         if chrom(j,3)>popmax(3)
    
74.            chrom(j,3)=popmax(3);
    
75.         end
    
76.         chrom(j,find(chrom(j,:)<popmin(3)))=popmin(3);
    
77.         %新粒子适应度值
    
78.         fitness(j)=fun(chrom(j,:),LL1);
    
79.     end
    
80.    
    
81.     %对更新后群体进行选择和交叉，以生成剩下一半的个体；
    
82.     chrom(ceil(sizepop/2)+1:sizepop,:)=Cross(lenchrom,chrom,fitness,(sizepop-ceil(sizepop/2)),bound);
    
83.     %对其进行变异
    
84.     chrom(ceil(sizepop/2)+1:sizepop,:)=Mutation(pmutation,lenchrom,chrom(ceil(sizepop/2)+1:sizepop,:),(sizepop-ceil(sizepop/2)),[i maxgen],bound);
    
85.     %计算变异后的适应度
    
86.     for j=ceil(sizepop/2)+1:sizepop
    
87.         fitness(j)=fun(chrom(j,:),LL1);
    
88.     end
    
89.     %对个群体级值进行更新
    
90.     [newbestfitness,newbestindex]=min(fitness);
    
91.     if fitnesszbest>newbestfitness
    
92.        fitnesszbest=newbestfitness;
    
93.        zbest=chrom(newbestindex,:);
    
94.     end
    
95.     trace(i)=fitnesszbest;
    
96. end
    
97. fitnesszbest
    
98. zbest
    
99. 
    
100. time=[31 28 31 30 31 30 31 31 30 31 30 31];
     
101. ssum=0;
     
102. i=0;
     
103. while ssum<round(zbest(1))
     
104.       i=i+1;
     
105.       ssum=ssum+time(i);
     
106. end
     
107. date=[i,round(zbest(1))-(ssum-time(i))];
     
108. WE=120+(zbest(2)-8.9)*15; %求经度附件4的基准时间8.9
     
109. NS=zbest(3);
     

复制代码

</round(zbest(1))
</popmin(3)))=popmin(3);
</vmin(3)
</vmin(2)
</vmin(1)