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

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (4876)
>考研 (898)
>導師招生 (747)
>文獻求助 (278)
>蟲友互識 (178)
>考博 (170)
>博后之家 (127)
>碩博家園 (106)
>休閑灌水 (90)
>基金申請 (77)
>論文投稿 (75)
>公派出國 (43)
>招聘信息布告欄 (42)
>教師之家 (30)
>找工作 (30)
>SciFinder/Reaxys (23)

返回列表

M.bySong

鐵蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 419.1
散金: 26
帖子: 145
在線: 10.3小時
蟲號: 4403056
注冊: 2016-02-15
專業(yè): 零件成形制造

[求助] 1stopt高版本代跑，很急，先謝謝大神們了！已有1人參與

大神上午好，最近碰到一個多目標優(yōu)化的問題，想用理想點法與平方和加權法求多目標最優(yōu)。三個目標函數(shù)，分別跑出最優(yōu)值以后，加權求解。希望大神能夠幫忙代跑一下，不勝感激。就這點金幣了，完成后加v有大紅包，現(xiàn)在很急，先謝謝大神們了�。�！

求第一個目標目標函數(shù)最優(yōu)值op1：
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
minfunction 6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3));
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;
——————————————————————————————————————————————————————————————————————————————————第2個目標目標函數(shù)最優(yōu)值op2：
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
minfunction sqrt((x1-x2)^2+(-h3+h2)^2);
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;       //收攏長度
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;          //H
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;             //展開長度
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;    //導向比
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;       //導向行程
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;       //導向行程
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;                         //拉桿長度
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;

——————————————————————————————————————————————————————————————————————————————————第3個目標目標函數(shù)最優(yōu)值op3：
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
minfunction sqrt((h3-269)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6422)^2);
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;       //收攏長度
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;          //H
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;             //展開長度
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;    //導向比
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;       //導向行程
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;       //導向行程
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;                         //拉桿長度
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;

————————————————————————————————————————————————————————————————————————
將下面代碼中的op1，op2，op3替換為上面計算出的結果，求多目標加權結果
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
ConstStr f1=6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3));
conststr f2=minfunction sqrt((x1-x2)^2+(-h3+h2)^2);
conststr f3=sqrt((h3-269)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6422)^2);
minfunction ((f1-op1)/op1)^2+((f2-op2)/op2)^2+((f3-op3)/op3)^2;
//minfunction 6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3));
//minfunction sqrt((x1-x2)^2+(-h3+h2)^2);
//minfunction sqrt((h3-269)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6422)^2);
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;       //收攏長度
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;          //H
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;             //展開長度
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;    //導向比
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;       //導向行程
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;       //導向行程
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;                         //拉桿長度
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;

回復此樓

» 猜你喜歡

幫我的英語口語老師找學生已經(jīng)有4人回復
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有277人回復
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復
哈爾濱理工大學物理系招收物理學考研調劑已經(jīng)有2人回復
0702一志愿吉大B區(qū)求調劑已經(jīng)有2人回復
求調劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復

1樓 2019-07-11 11:09:33

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

獨孤神宇

版主 (知名作家)

應助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學
管轄: 計算模擬

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應助指數(shù) +1
小紅豆: 金幣+20, 應助指數(shù)+1, 回帖應助 2019-07-11 14:39:04
M.bySong: 金幣+38 2019-07-11 14:53:51

目標函數(shù)值(最小): 1.87232999181561
x1: 6422
x2: 2282.96873670372
x3: 4485.76040752828
h1: 2279.99154310903
h2: 1600.03227522495
h3: 1241.37911564834
h: -1006.24586353898

贊一下(2人)

回復此樓

數(shù)值計算

2樓2019-07-11 14:11:41

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

相關版塊跳轉我要訂閱樓主 M.bySong 的主題更新