houchuyue

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 248
帖子: 20
在線: 12.6小時
蟲號: 2690455
注冊: 2013-09-29
專業(yè): 冶金物理化學與冶金原理

[求助] 求指點，1stopt如何是所求值為正值已有2人參與

1.我的目標是求X1,X2,X3,且是正數(shù)的解，如何修改，應(yīng)該如何添加約束條件。2。為何每次運行，結(jié)果都不同，如何處理？
Constant a=0.01195,b=0.0256,k=0.0105,B3=0.0027;
Parameters x1, x2, x3;
Function
a+b=x1+x2;
x3+2*b=x1+2*x2;
x2*x3*y/x1=k;
I=0.5*(x1+4*x2+x3+4*b);
f=-2*(1.303*ln(1+1.2*I^(0.5))+1.303*I*(1/(1+1.2*I^(0.5))*0.5*1.2*I^(-0.5)));
B1=0.2103+(2*0.4711/(1.4^(2)*I))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))));
B2=-2*0.4711/(1.4^(2)*I^(2))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))))-2*0.4711/(1.4^(2)*I)*(0.5*1.4*I^(-0.5)*exp(-1.4*I^(0.5))+(1+1.4*I^(0.5))*exp(-1.4*I^(0.5))*(-0.5*1.4*I^(-0.5)));
B4=0.2+(2*2.7/(1.4^(2)*I))*(1-(1+1.4*I^(0.5))*exp(-1.4*I^(0.5)))+(2*(-30.7)/(12^(2)*I))*(1-(1+12*I^(0.5))*exp(-12*I^(0.5)));
B5=-2*2.7/(1.4^(2)*I^(2))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))))-2*2.7/(1.4^(2)*I)*((0.5*1.4*I^(-0.5))*exp(-1.4*I^(0.5))+(1+1.4*I^(0.5))*exp(-1.4*I^(0.5))*(-0.5*1.4*I^(-0.5)))-2*(-30.7)/(12^(2)*I^(2))*(1-(1+12*I^(0.5)*exp(-12*I^(0.5))))-2*(-30.7)/(12^(2)*I)*((0.5*12*I^(-0.5))*exp(-12*I^(0.5))+(1+12*I^(0.5))*exp(-12*I^(0.5))*(-0.5*12*I^(-0.5)));
B6=0.375+(2*1.5/(1.4^(2)*I))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))));
B7=-2*1.5/(1.4^(2)*I^(2))*(1-(1+1.4*I^(0.5)*exp(-1.4*I^(0.5))))-2*1.5/(1.4^(2)*I)*(0.5*1.4*I^(-0.5)*exp(-1.4*I^(0.5))+(1+1.4*I^(0.5))*exp(-1.4*I^(0.5))*(-0.5*1.4*I^(-0.5)));
ln(Y)=f+2*(x1-x3)*B1+4*x3*x1*B2+2*(x2+x3)*B3+2*((x2+x3)*(x3+b)+x3*x2)*0.0416+2*b*(2*x2*B4+B5+2*x1*B6-B7);

迭代數(shù): 72
計算用時(時:分:秒:毫秒): 00:00:31:104
計算中止原因: 達到收斂判定標準
優(yōu)化算法: 準牛頓法(BFGS) + 通用全局優(yōu)化法
函數(shù)表達式 1: 0.01195+0.0256-(x1+x2)
      2: x3+2*0.0256-(x1+2*x2)
      3: x2*x3*y/x1-(0.0105)
      4: i-(0.5*(x1+4*x2+x3+4*0.0256))
      5: f-(-2*(1.303*ln(1+1.2*i^(0.5))+1.303*i*(1/(1+1.2*i^(0.5))*0.5*1.2*i^(-0.5))))
      6: b1-(0.2103+(2*0.4711/(1.4^(2)*i))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5)))))
      7: b2-(-2*0.4711/(1.4^(2)*i^(2))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5))))-2*0.4711/(1.4^(2)*i)*(0.5*1.4*i^(
         -0.5)*exp(-1.4*i^(0.5))+(1+1.4*i^(0.5))*exp(-1.4*i^(0.5))*(-0.5*1.4*i^(-0.5))))
      8: b4-(0.2+(2*2.7/(1.4^(2)*i))*(1-(1+1.4*i^(0.5))*exp(-1.4*i^(0.5)))+(2*(-30.7)/(12^(2)*i))*(1-(1+12*i^(0.5))
         *exp(-12*i^(0.5))))
      9: b5-(-2*2.7/(1.4^(2)*i^(2))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5))))-2*2.7/(1.4^(2)*i)*((0.5*1.4*i^(-0.5))
         *exp(-1.4*i^(0.5))+(1+1.4*i^(0.5))*exp(-1.4*i^(0.5))*(-0.5*1.4*i^(-0.5)))-2*(-30.7)/(12^(2)*i^(2))*(1
         -(1+12*i^(0.5)*exp(-12*i^(0.5))))-2*(-30.7)/(12^(2)*i)*((0.5*12*i^(-0.5))*exp(-12*i^(0.5))+(1+12*i^(0.5))
         *exp(-12*i^(0.5))*(-0.5*12*i^(-0.5))))
      10: b6-(0.375+(2*1.5/(1.4^(2)*i))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5)))))
      11: b7-(-2*1.5/(1.4^(2)*i^(2))*(1-(1+1.4*i^(0.5)*exp(-1.4*i^(0.5))))-2*1.5/(1.4^(2)*i)*(0.5*1.4*i^(-0.5)
         *exp(-1.4*i^(0.5))+(1+1.4*i^(0.5))*exp(-1.4*i^(0.5))*(-0.5*1.4*i^(-0.5))))
      12: ln(y)-(f+2*(x1-x3)*b1+4*x3*x1*b2+2*(x2+x3)*0.0027+2*((x2+x3)*(x3+0.0256)+x3*x2)*0.0416+2*0.0256*(2*x2
         *b4+b5+2*x1*b6-b7))
目標函數(shù)值: 0.0346057342149256
x1: -0.613628828746281
x2: 0.651181234844569
x3: 0.637530987966978
y: 0.0356211237987828
i: 1.36550919200013
f: -3.04451434010513
b1: 0.0981303646848507
b2: 0.149337201828133
b4: 0.869550118362684
b5: 0.855878503182631
b6: 0.0178474428783668
b7: 0.475498544863951

====== 計算結(jié)束 ======

[ Last edited by houchuyue on 2014-3-12 at 16:18 ]

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有217人回復(fù)
哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學工程學碩調(diào)劑信息已經(jīng)有27人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有45人回復(fù)

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

求助1stOpt作圖已經(jīng)有6人回復(fù)
求助1stopt含約束非線性擬合已經(jīng)有5人回復(fù)
1stopt求最小值已經(jīng)有37人回復(fù)
誰用1stopt這個軟件啊，離散數(shù)據(jù)擬合方程，方程未知已經(jīng)有5人回復(fù)
求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有13人回復(fù)
求助，用1stopt解方程組，兩次得到的解不唯一？什么原因，大神幫忙已經(jīng)有6人回復(fù)
求高版本1stopt擬合，已經(jīng)有11人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有9人回復(fù)
求大俠幫助跑一下1stopt 已經(jīng)有12人回復(fù)
哪位朋友有注冊版1stOpt軟件？急需已經(jīng)有5人回復(fù)
1stOPT1.5版每次仿真結(jié)果不一樣，求助！急！已經(jīng)有3人回復(fù)
運算過程有復(fù)數(shù)出現(xiàn)，求高版本1stOpt 已經(jīng)有6人回復(fù)
1stopt運行程序，總是出現(xiàn)字符函數(shù)錯誤的提示，請高手幫助已經(jīng)有5人回復(fù)
請問matlab擬合二元方程的問題已經(jīng)有3人回復(fù)
麻煩高人指點一下1stopt和spss的非線性回歸中weibull的兩參數(shù)的公式應(yīng)該怎么表達啊？已經(jīng)有9人回復(fù)
MATLAB擬合工具如何引入誤差棒？已經(jīng)有11人回復(fù)
學1stopt要多久？已經(jīng)有12人回復(fù)
用1stopt進行數(shù)據(jù)擬合時，無法運行sharedmodel ,求助高手已經(jīng)有7人回復(fù)

1樓 2014-03-12 16:11:48

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
houchuyue: 金幣+5, ★★★很有幫助 2014-03-12 16:33:27

Parameters x1[0,], x2[0,], x3[0,];

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2014-03-12 16:25:04

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

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
houchuyue: 金幣+2, ★★★很有幫助 2014-03-12 17:08:20

不加大于0的范圍約束求得的x1、x2、x3也都是大于0的啊，且穩(wěn)定唯一：

目標函數(shù)值(最小): 2.44738091577421E-22
x1: 0.0174024208040297
x2: 0.0201476131995002
x3: 0.00649756064023006
y: 1.39581797264244
i: 0.103445344352546
f: -1.21342501037417
b1: -1.12354414258258
b2: 15.7971962194402
b4: -1.4883351519721
b5: 81.0412352437473
b6: -3.87200958150146
b7: 50.2988629360697

贊一下

回復(fù)此樓

3樓2014-03-12 16:46:46

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

houchuyue

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 248
帖子: 20
在線: 12.6小時
蟲號: 2690455
注冊: 2013-09-29
專業(yè): 冶金物理化學與冶金原理

引用回帖:

3樓: Originally posted by dingd at 2014-03-12 16:46:46
不加大于0的范圍約束求得的x1、x2、x3也都是大于0的啊，且穩(wěn)定唯一：

目標函數(shù)值(最小): 2.44738091577421E-22
x1: 0.0174024208040297
x2: 0.0201476131995002
x3: 0.00649756064023006
y: 1.3958179726424 ...

我用的是1.5破解版。算不出滿意的結(jié)果，一直變動。請問是版本原因？還是我需要在算法設(shè)置上要進行調(diào)整？

贊一下

回復(fù)此樓

4樓2014-03-12 17:07:29

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

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

1.5還是破解版，原因就說不清了。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

houchuyue

5樓2014-03-12 17:11:57

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

houchuyue

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 248
帖子: 20
在線: 12.6小時
蟲號: 2690455
注冊: 2013-09-29
專業(yè): 冶金物理化學與冶金原理

送紅花一朵

內(nèi)容已刪除

回復(fù)此樓

6樓2014-03-12 17:13:15

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 材料與化工求調(diào)劑 +6	為學666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以 +4	小丁想進步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 0703化學調(diào)劑 +6	妮妮ninicgb 2026-03-15	9/450	2026-03-16 16:40 by houyaoxu
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 266求調(diào)劑 +4	學員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 求調(diào)劑 +5	一定有學上- 2026-03-12	5/250	2026-03-13 18:31 by ms629