 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 797 | 回復(fù): 3 | ||
[求助]
求程序指導(dǎo) 已有2人參與
|
|
題目如下： minf(1) = ((7240.6*x(1)+11.3)*sqrt(9.81*x(2)*cos(x(3))))/(300.91*pi*sqrt(x(1))*x(3)*(1+22*tan(x(3))*sqrt(tan(x(3))))); minf(2) = ((12908.88*x(1)+11.1877)*sqrt(10.791*x(2)*cos(0.87*x(3))))/(290.89*pi*sqrt(1.8*x(1))*x(3)*(1+22*tan(0.87*x(3))*sqrt(tan(0.87*x(3))))); 兩個(gè)函數(shù)差不多 0.004《x1《0.006 4《x2《5 22《x3《25 x3表示角度 我想用遺傳求兩個(gè)函數(shù)同時(shí)取最小值時(shí)的x1，x2，x3的值并給出程序 求大神們指教，謝謝！ 根據(jù)一個(gè)題目編的的簡(jiǎn)單程序： f1=inline('((7240.6*x(1)+11.3)*sqrt(9.81*x(2)*cos(x(3))))/(300.91*pi*sqrt(x(1))*x(3)*(1+22*tan(x(3))*sqrt(tan(x(3)))))'); f2=@(x) ((12908.88*x(1)+11.1877)*sqrt(10.791*x(2)*cos(0.87*x(3))))/(290.89*pi*sqrt(1.8*x(1))*x(3)*(1+22*tan(0.87*x(3))*sqrt(tan(0.87*x(3))))); fun=@(x) [f1(x),f2(x)]; Lb=[0.004,4,22]; ub=[0.006,5,25]; options=gaoptimset('TolFun',le-4); [x,fval]=gamultiobj(fun,3,[],[],[],Lb,ub) x,fval; 這個(gè)程序運(yùn)行提示：Error using <= Not enough input arguments. 另外一個(gè)程序： 將下屬兩個(gè)目標(biāo)函數(shù)分別保存在兩個(gè)m文件中 function f1=func1(x) %第一目標(biāo)函數(shù) f1=x(:,1).*x(:,1)./4+x(:,2).*x(:,2)./4; function f2=func2(x) %第二目標(biāo)函數(shù) f2=x(:,1).*(1-x(:,2))+10; function GA() clear;clc;close all NIND=100; %個(gè)體數(shù)目 MAXGEN=50; %最大遺傳代數(shù) NVAR=2; %變量個(gè)數(shù) PRECI=20; %變量的二進(jìn)制位數(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) ; %初始種群十進(jìn)制轉(zhuǎn)換 gen=1; while gen<MAXGEN, [NIND,N]=size(Chrom) ; M=fix(NIND/2) ; ObjV1=func1(v(1:M,: )) ; %分組后第一目標(biāo)函數(shù)值 FitnV1=ranking(ObjV1) ; %分配適應(yīng)度值 SelCh1=select('sus',Chrom(1:M,: ),FitnV1,GGAP) ; %選擇 ObjV2=func2(v(M+1:NIND,: )) ; %分組后第二目標(biāo)函數(shù)值 FitnV2=ranking(ObjV2) ; %分配適應(yīng)度值 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) ; 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;l egend('解的變化','種群均值的變化') xlabel('迭代次數(shù)') ;ylabel('目標(biāo)函數(shù)值') ; figure(2);clf;plot(trace2(:,1)) ;hold on; plot(trace2(:,2),'-.') ; plot(trace2(:,1),'.') ; plot(trace2(:,2),'.') ;grid; legend('解的變化','種群均值的變化') ; xlabel('迭代次數(shù)') ;ylabel('目標(biāo)函數(shù)值') ; figure(3) ;clf; plot(trace3(:,1)) ;hold on; plot(trace3(:,2),'-.') ; plot(trace3(:,1),'.') ; plot(trace3(:,2),'.') ;grid; legend('解的變化','種群均值的變化') ; xlabel('迭代次數(shù)');ylabel('目標(biāo)函數(shù)值') ; figure(4);clf; plot(func1(v));hold on; plot(func2(v),'r-.');grid; 這個(gè)程序中我不明白：rep([1;0;1;1],[1,NVAR])]; 括號(hào)中的數(shù)字怎么定義的？ [ Last edited by jjdg on 2014-8-5 at 07:20 ] |
回復(fù)此樓» 猜你喜歡
287求調(diào)劑
已經(jīng)有5人回復(fù)
-
22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù)
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有10人回復(fù)
-
廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑）
已經(jīng)有3人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有8人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有25人回復(fù)
-
085600材料與化工
已經(jīng)有6人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
elastic
金蟲 (小有名氣)
- 應(yīng)助: 43 (小學(xué)生)
- 金幣: 1426.4
- 紅花: 4
- 帖子: 207
- 在線: 54.1小時(shí)
- 蟲號(hào): 3006535
- 注冊(cè): 2014-03-01
- 性別: GG
- 專業(yè): 計(jì)算機(jī)硬件技術(shù)
2樓2014-08-03 17:54:10
xajdwangpeng
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3.5
- 帖子: 4
- 在線: 2.4小時(shí)
- 蟲號(hào): 3350539
- 注冊(cè): 2014-08-03
3樓2014-08-04 22:36:17
4樓2014-08-11 02:49:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 木凡_123 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +4 | Mqqqqqq 2026-03-19 | 4/200 |
|
|
|
[考研] 材料學(xué)碩297已過四六級(jí)求調(diào)劑推薦 +6 | adaie 2026-03-19 | 6/300 |
|
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|