houchuyue

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

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

1.我的目標(biāo)是求X1,X2,X3,且是正數(shù)的解，如何修改，應(yīng)該如何添加約束條件。2。為何每次運(yùn)行，結(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
計(jì)算用時(shí)(時(shí):分:秒:毫秒): 00:00:31:104
計(jì)算中止原因: 達(dá)到收斂判定標(biāo)準(zhǔn)
優(yōu)化算法: 準(zhǔn)牛頓法(BFGS) + 通用全局優(yōu)化法
函數(shù)表達(dá)式 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))
目標(biāo)函數(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

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

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

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有101人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(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ù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有54人回復(fù)

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

求助1stOpt作圖已經(jīng)有6人回復(fù)
求助1stopt含約束非線性擬合已經(jīng)有5人回復(fù)
1stopt求最小值已經(jīng)有37人回復(fù)
誰(shuí)用1stopt這個(gè)軟件啊，離散數(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ù)
哪位朋友有注冊(cè)版1stOpt軟件？急需已經(jīng)有5人回復(fù)
1stOPT1.5版每次仿真結(jié)果不一樣，求助！急！已經(jīng)有3人回復(fù)
運(yùn)算過(guò)程有復(fù)數(shù)出現(xiàn)，求高版本1stOpt 已經(jīng)有6人回復(fù)
1stopt運(yùn)行程序，總是出現(xiàn)字符函數(shù)錯(cuò)誤的提示，請(qǐng)高手幫助已經(jīng)有5人回復(fù)
請(qǐng)問(wèn)matlab擬合二元方程的問(wèn)題已經(jīng)有3人回復(fù)
麻煩高人指點(diǎn)一下1stopt和spss的非線性回歸中weibull的兩參數(shù)的公式應(yīng)該怎么表達(dá)�。� 已經(jīng)有9人回復(fù)
MATLAB擬合工具如何引入誤差棒？已經(jīng)有11人回復(fù)
學(xué)1stopt要多久？已經(jīng)有12人回復(fù)
用1stopt進(jìn)行數(shù)據(jù)擬合時(shí)，無(wú)法運(yùn)行sharedmodel ,求助高手已經(jīng)有7人回復(fù)

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

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

月只藍(lán)

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

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(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小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

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

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

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

目標(biāo)函數(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

銀蟲(chóng) (初入文壇)

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

引用回帖:

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

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

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

贊一下

回復(fù)此樓

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

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

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

1.5還是破解版，原因就說(shuō)不清了。

贊一下

回復(fù)此樓

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

houchuyue

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

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

houchuyue

銀蟲(chóng) (初入文壇)

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

送紅花一朵

內(nèi)容已刪除

回復(fù)此樓

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

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

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

[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4	Ymlll 2026-03-18	6/300	2026-03-18 22:15 by li123456789.
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 266求調(diào)劑 +5	陽(yáng)陽(yáng)哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 085601材料工程專碩求調(diào)劑 +6	慕寒mio 2026-03-16	6/300	2026-03-18 14:26 by 007_lilei
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛(ài)生物生物愛(ài)?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂(lè) 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu