 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 795 | 回復(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))))); 兩個函數(shù)差不多 0.004《x1《0.006 4《x2《5 22《x3《25 x3表示角度 我想用遺傳求兩個函數(shù)同時取最小值時的x1，x2，x3的值并給出程序 求大神們指教，謝謝！ 根據(jù)一個題目編的的簡單程序： 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; 這個程序運行提示：Error using <= Not enough input arguments. 另外一個程序： 將下屬兩個目標函數(shù)分別保存在兩個m文件中 function f1=func1(x) %第一目標函數(shù) f1=x(:,1).*x(:,1)./4+x(:,2).*x(:,2)./4; function f2=func2(x) %第二目標函數(shù) f2=x(:,1).*(1-x(:,2))+10; 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<MAXGEN, [NIND,N]=size(Chrom) ; M=fix(NIND/2) ; ObjV1=func1(v(1:M,: )) ; %分組后第一目標函數(shù)值 FitnV1=ranking(ObjV1) ; %分配適應(yīng)度值 SelCh1=select('sus',Chrom(1:M,: ),FitnV1,GGAP) ; %選擇 ObjV2=func2(v(M+1:NIND,: )) ; %分組后第二目標函數(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('目標函數(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; 這個程序中我不明白：rep([1;0;1;1],[1,NVAR])]; 括號中的數(shù)字怎么定義的？ [ Last edited by jjdg on 2014-8-5 at 07:20 ] |
回復(fù)此樓» 猜你喜歡
290求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑
已經(jīng)有4人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有6人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有8人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有8人回復(fù)
-
317求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南昌大學(xué)，327分，材料與化工085600
已經(jīng)有5人回復(fù)
-
274求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有9人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有11人回復(fù)
elastic
金蟲 (小有名氣)
- 應(yīng)助: 43 (小學(xué)生)
- 金幣: 1426.4
- 紅花: 4
- 帖子: 207
- 在線: 54.1小時
- 蟲號: 3006535
- 注冊: 2014-03-01
- 性別: GG
- 專業(yè): 計算機硬件技術(shù)
2樓2014-08-03 17:54:10
3樓2014-08-04 22:36:17
4樓2014-08-11 02:49:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 木凡_123 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +5 | adaie 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +6 | 梁梁校校 2026-03-17 | 6/300 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +8 | 步川酷紫123 2026-03-13 | 8/400 |
|
|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿福大288有機化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研]
 0854，計算機類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|