標(biāo)題: 分享下數(shù)模2015年A題的Matlab程序 [打印本頁(yè)]
作者: 51黑專(zhuān)家    時(shí)間: 2016-5-8 19:58
標(biāo)題: 分享下數(shù)模2015年A題的Matlab程序


全部資料（壓縮包）下載：

程序.rar (24.68 KB, 下載次數(shù): 12)

cumcm2012B附件6_可參考的相關(guān)概念.doc (77.5 KB, 下載次數(shù): 8)


哈哈！三天終于熬過(guò)頭了，今天也休息過(guò)來(lái)了，來(lái)做個(gè)分享！
國(guó)賽本科兩個(gè)題，我們考慮了自己的劣勢(shì)，對(duì)于創(chuàng)新類(lèi)的題做起來(lái)比較惱火，最終選擇了A題。

A題入手需要有一些相關(guān)的地理知識(shí)，例如方位角，時(shí)角等等。其實(shí)這些內(nèi)容在2012年B題的附件里都有解釋?zhuān)罴倨陂g老師讓我們做了題，當(dāng)時(shí)比較認(rèn)真的做。第一問(wèn)上午就解出來(lái)了，后來(lái)一直卡在后面，個(gè)人matlab編程水平不是很高超，所以隊(duì)友列出的非線性方程解得不是很好。拖了隊(duì)友的后腿。最后一問(wèn)在學(xué)長(zhǎng)的點(diǎn)撥下做了出來(lái)。論文尊重隊(duì)友意見(jiàn)，就不分享了。程序還是可以共享的。

做完了就是做完了！不考慮結(jié)果了。
上圖便于大家理解吧，空間想象你能給力不好的可以靠這個(gè)圖來(lái)思考的







下面是部分程序預(yù)覽（完整版本請(qǐng)下載本帖附件）：
  1. function [date,WE,NS]=fmain()  

  3. %參數(shù)賦初值
  4. maxgen=90;            %進(jìn)化代數(shù)，即迭代次數(shù)
  5. sizepop=200;           %種群規(guī)模
  6. pcross=0.6;            %交叉概率選擇
  7. pmutation=0.2;         %變異概率選擇
  8. lenchrom=[1 1 1];        %每個(gè)變量的字符長(zhǎng)度，如果是浮點(diǎn)，則長(zhǎng)度都變?yōu)?
  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];

  17. %附件4的影子的測(cè)量長(zhǎng)度
  18. LL1=textread('D.txt');


  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

  29. [bestfitness,bestindex]=min(fitness);
  30. zbest=chrom(bestindex,:);   %群體極值位置
  31. gbest=chrom;                %個(gè)體位置
  32. fitnessgbest=fitness;       %個(gè)體適應(yīng)度值
  33. fitnesszbest=bestfitness;   %群體極值
  34. %開(kāi)始迭代
  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算法進(jìn)行更新
  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.         %新粒子適應(yīng)度值
  78.         fitness(j)=fun(chrom(j,:),LL1);
  79.     end
  80.    
  81.     %對(duì)更新后群體進(jìn)行選擇和交叉，以生成剩下一半的個(gè)體；
  82.     chrom(ceil(sizepop/2)+1:sizepop,:)=Cross(lenchrom,chrom,fitness,(sizepop-ceil(sizepop/2)),bound);
  83.     %對(duì)其進(jìn)行變異
  84.     chrom(ceil(sizepop/2)+1:sizepop,:)=Mutation(pmutation,lenchrom,chrom(ceil(sizepop/2)+1:sizepop,:),(sizepop-ceil(sizepop/2)),[i maxgen],bound);
  85.     %計(jì)算變異后的適應(yīng)度
  86.     for j=ceil(sizepop/2)+1:sizepop
  87.         fitness(j)=fun(chrom(j,:),LL1);
  88.     end
  89.     %對(duì)個(gè)群體級(jí)值進(jìn)行更新
  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

  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; %求經(jīng)度附件4的基準(zhǔn)時(shí)間8.9
  109. NS=zbest(3);
復(fù)制代碼




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




歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1