版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

tanzouqing

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 333.4
散金: 1
帖子: 70
在線: 14.4小時
蟲號: 747434
注冊: 2009-04-13
專業(yè): 固體力學

[求助] Matlab的遺傳算法求解多目標優(yōu)化問題

下面是經(jīng)典的基于Matlab的遺傳算法的多目標優(yōu)化程序,
我的疑問:
為何沒有為何沒有”重插入子代的新種群”這一步驟?
如reins重新插入函數(shù)

M文件func1.m
function f1=func1(x)
f1=x(:,1).*x(:,1)/4+x(:,2).*x(:,2)/4;    %第一目標函數(shù)

M文件func2.m
function f2=func2(x)
f2=x(:,1).*(1-x(:,2))+10;          %第二目標函數(shù)

M文件GA-func.m
function  GA()
clear;clc;close all
NIND=100;                         %個體數(shù)目
MAXGEN=50;                      %最大遺傳代數(shù)
NVAR=2;                         %變量個數(shù)
PRECI=20;                         %變量的二進制位數(shù)
GGAP=0.9;                         %代溝
trace1=[];trace2=[];trace3=[];          %性能跟蹤
%建立區(qū)域描述器
% rep([PRECI],[1,NVAR])
FieldD=[rep([PRECI],[1,NVAR]);rep([1;2],[1,NVAR]);rep([1;0;1;1],[1,NVAR])];
Chrom=crtbp(NIND,NVAR*PRECI);          %初始種群
v=bs2rv(Chrom,FieldD) ;                %初始種群十進制轉(zhuǎn)換
gen=1;
while gen [NIND,N]=size(Chrom);
M=fix(NIND/2);

ObjV1=func1(v(1:M,

); %分組后第一目標函數(shù)值
FitnV1=ranking(ObjV1); %分配適應度值
SelCh1=select('sus',Chrom(1:M,

,FitnV1,GGAP); %選擇

ObjV2=func2(v(M+1:NIND,

); %分組后第二目標函數(shù)值
FitnV2=ranking(ObjV2); %分配適應度值
SelCh2=select('sus',Chrom(M+1:NIND,

,FitnV2,GGAP);       %選擇

SelCh=[SelCh1;SelCh2];                %合并

SelCh=recombin('xovsp',SelCh,0.7);       %重組
Chrom=mut(SelCh);                      %變異
v=bs2rv(Chrom,FieldD);                %二進制轉(zhuǎn)化為十進制

trace1(gen,1)=min(func1(v));
trace1(gen,2)=sum(func1(v))/length(func1(v));
trace2(gen,1)=min(func2(v));
trace2(gen,2)=sum(func2(v))/length(func2(v));
trace3(gen,1)=min(func1(v)+func2(v));
trace3(gen,2)=sum(func1(v))/length(func1(v))+sum(func2(v))/length(func2(v));
gen=gen+1;
end
figure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid on;
legend('解的變化','種群均值的變化')
xlabel('迭代次數(shù)');ylabel('第一目標函數(shù)值');
figure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
plot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的變化','種群均值的變化');
xlabel('迭代次數(shù)');ylabel('第二目標函數(shù)值');
figure(3);clf;
plot(trace3(:,1));hold on;
plot(trace3(:,2),'-.');
plot(trace3(:,1),'.');
plot(trace3(:,2),'.');grid;
legend('解的變化','種群均值的變化');
xlabel('迭代次數(shù)');ylabel('目標函數(shù)值之和');
figure(4);clf;plot(func1(v));hold on;
plot(func2(v),'r-.');grid;

[v,func1(v),func2(v)] %查看數(shù)值結果

回復此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復
一志愿山東大學藥學學碩求調(diào)劑已經(jīng)有4人回復
07化學280分求調(diào)劑已經(jīng)有4人回復
298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑已經(jīng)有12人回復
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復
335求調(diào)劑已經(jīng)有5人回復
求調(diào)劑已經(jīng)有7人回復
一志愿吉大化學322求調(diào)劑已經(jīng)有4人回復
環(huán)境學碩288求調(diào)劑已經(jīng)有8人回復
341求調(diào)劑(一志愿湖南大學070300) 已經(jīng)有6人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

幫忙給看看這個matlab優(yōu)化函數(shù) 問題已經(jīng)有8人回復
由入門到精通--matlab求解高等數(shù)學問題課件已經(jīng)有100人回復
求粒子群優(yōu)化算法進行變量優(yōu)選的matlab程序已經(jīng)有20人回復
薛定宇《高等應用數(shù)學問題的MATLAB求解》已經(jīng)有52人回復
用matlab求解方程出問題，請幫忙看看已經(jīng)有3人回復
matlab線性矩陣不等式求解做優(yōu)化問題已經(jīng)有10人回復
求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼已經(jīng)有9人回復
求助有關Matlab有約束非線性最優(yōu)化問題已經(jīng)有9人回復
求matlab求解一個方程組的問題已經(jīng)有6人回復
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復
【求助】請問如何求解這個優(yōu)化問題已經(jīng)有7人回復
關于算法比較的問題，請大家?guī)兔?/a> 已經(jīng)有7人回復

【求助】matlab 二次規(guī)劃的優(yōu)化的問題已經(jīng)有4人回復

1樓 2012-10-15 17:07:41

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

djzh_7

2樓2012-10-31 19:20:13

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉(zhuǎn) 我要訂閱樓主 tanzouqing 的主題更新