 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 800 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
求程序指導(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ù)此樓» 猜你喜歡
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有13人回復(fù)
-
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有6人回復(fù)
-
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷
已經(jīng)有10人回復(fù)
-
298求調(diào)劑
已經(jīng)有6人回復(fù)
-
310求調(diào)劑
已經(jīng)有4人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
307求調(diào)劑
已經(jīng)有11人回復(fù)
-
材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑
已經(jīng)有5人回復(fù)
-
315分，誠(chéng)求調(diào)劑，材料與化工085600
已經(jīng)有3人回復(fù)
-
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò)
已經(jīng)有10人回復(fù)
xajdwangpeng
新蟲(chóng) (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 3.5
- 帖子: 4
- 在線: 2.4小時(shí)
- 蟲(chóng)號(hào): 3350539
- 注冊(cè): 2014-08-03
3樓2014-08-04 22:36:17
elastic
金蟲(chóng) (小有名氣)
- 應(yīng)助: 43 (小學(xué)生)
- 金幣: 1426.4
- 紅花: 4
- 帖子: 207
- 在線: 54.1小時(shí)
- 蟲(chóng)號(hào): 3006535
- 注冊(cè): 2014-03-01
- 性別: GG
- 專業(yè): 計(jì)算機(jī)硬件技術(shù)
2樓2014-08-03 17:54:10
4樓2014-08-11 02:49:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|