版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲(chóng)號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

[求助] 1stOpt 求助，求高手指導(dǎo)�。。。�！已有1人參與

SharedModel;
Variable x,y(2);//y1,y2;
Parameter a,b,c,d;
Constant epsilon=0.3921,h0=6.3,R=125,Sf=0.03143,VR=1308.33,h1=3.82,Re=133,z0=134.91;
ConstStr z=z0-Sqrt(Re^2-x^2),Vx=h1*VR*(1+Sf)/(2*z),delta_f=15+9.2*epsilon^0.11,
p1=Exp(((a*Vx-VR*(1+Sf))*2*b*x)/(z^2)),p2=-(2/z)*Sqrt(Re^2-x^2),
p3=(4*b/(z^3))*(VR*(1+Sf)-a*Vx),p4=(-x/2)*Sqrt(Re^2-x^2)+((Re^2)/2)*(Arcsin(x/Re)),
p5=(1+4*Re/h1);
Function y1=(delta_f*p1)*(p2+p3*p4+p5);
y2=(delta_f*x)/(Sqrt(Re^2-x^2))+((2*d*y1)*(c*Vx-VR*(1+Sf)))/z;
Data;
//x,y1,y2
1.12891986062718e+000 6.40830449826990e+001 3.19723183391004e+001
2.25783972125436e+000 1.00069204152249e+002 3.47404844290657e+001
3.76306620209059e+000 8.56747404844291e+001 4.13840830449827e+001
5.01742160278746e+000 8.06920415224914e+001 2.58823529411765e+001
6.39721254355401e+000 7.62629757785467e+001 2.25605536332180e+001
7.77700348432056e+000 6.68512110726644e+001 1.75778546712803e+001
8.96864111498258e+000 5.96539792387543e+001 9.68858131487885e-001
1.02229965156794e+001 5.57785467128028e+001 -1.61937716262976e+001
1.16655052264808e+001 5.35640138408304e+001 -2.83737024221453e+001
1.28571428571429e+001 5.35640138408304e+001 -3.05882352941176e+001
1.42369337979094e+001 5.46712802768166e+001 -3.44636678200692e+001
1.54285714285714e+001 4.85813148788927e+001 -1.95155709342561e+001
1.67456445993031e+001 4.30449826989619e+001 -1.01038062283737e+001
1.75609756097561e+001 4.35986159169550e+001 -4.01384083044982e+000

運(yùn)行后，只有y1擬合圖像和a,b的估值， y2為一條等于零的直線，沒(méi)有c和d的估值，請(qǐng)問(wèn)這是怎么回事？是上面的代碼哪里有問(wèn)題嗎？？？跪求指導(dǎo)�。。。。�

回復(fù)此樓

» 猜你喜歡

211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) 已經(jīng)有7人回復(fù)
26申博已經(jīng)有3人回復(fù)
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專(zhuān)業(yè)都可以已經(jīng)有3人回復(fù)
材料專(zhuān)碩306英一數(shù)二已經(jīng)有7人回復(fù)
材料工程專(zhuān)碩274一志愿211求調(diào)劑已經(jīng)有6人回復(fù)
一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑已經(jīng)有4人回復(fù)
332求調(diào)劑已經(jīng)有5人回復(fù)
271求調(diào)劑已經(jīng)有14人回復(fù)
275求調(diào)劑已經(jīng)有4人回復(fù)
302求調(diào)劑已經(jīng)有8人回復(fù)

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

1stOpt 求助�。�！已經(jīng)有4人回復(fù)
用高版本1stOpt求解方程組已經(jīng)有6人回復(fù)
1stopt 語(yǔ)句修改已經(jīng)有23人回復(fù)
matlab高手進(jìn)！發(fā)現(xiàn)一個(gè)奇怪的現(xiàn)象…… 已經(jīng)有4人回復(fù)
求matlab大神指點(diǎn)，關(guān)于logistic擬合已經(jīng)有13人回復(fù)
求Matlab高手指導(dǎo) 已經(jīng)有12人回復(fù)
看了下1stOpt的介紹，感覺(jué)略顯浮夸已經(jīng)有14人回復(fù)
1stopt運(yùn)行報(bào)錯(cuò)“參數(shù)雖被定義，但是沒(méi)有在函數(shù)表示中出現(xiàn)”，順便求高手幫忙跑下程序已經(jīng)有3人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有6人回復(fù)
求大神跑下1stOpt，一個(gè)sharemodel，謝謝！已經(jīng)有7人回復(fù)
1stopt擬合的動(dòng)力學(xué)參數(shù)相關(guān)系數(shù)不高，有哪些原因已經(jīng)有3人回復(fù)
跪求高手幫忙看看這組數(shù)據(jù)如何進(jìn)行擬合已經(jīng)有7人回復(fù)
ELK簡(jiǎn)單計(jì)算失敗，求高手解答！已經(jīng)有4人回復(fù)
求高手給做一道題，用MATLAB或1stOpt，請(qǐng)把結(jié)果直接給出來(lái)，我的軟件就不行，謝謝已經(jīng)有5人回復(fù)
用1stOpt進(jìn)行復(fù)數(shù)曲線擬合時(shí)，結(jié)果出錯(cuò)。急求教！��！已經(jīng)有5人回復(fù)
求用高版本的1stOpt擬合方程已經(jīng)有8人回復(fù)
大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒(méi)有遇到過(guò)這種情況？已經(jīng)有8人回復(fù)
1stopt4階非線性常微分方程，幫小弟看一下，跪求。已經(jīng)有7人回復(fù)
求SWMM軟件高手指導(dǎo) 已經(jīng)有15人回復(fù)

數(shù)值計(jì)算

1樓 2015-06-18 16:36:59

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

月只藍(lán)

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

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

★
xzhdty: 金幣+1, 謝謝參與 2015-06-19 17:12:52

代碼沒(méi)有問(wèn)題。
1stopt版本原因。低版本不支持SharedModel。
你說(shuō)自己用的是5.0版本，這是不可能的，網(wǎng)上能找到的最高到1.5。
你的問(wèn)題可以用試用版解決：

Iterations: 20
Elapsed Time (Hr:Min:Sec:Msec): 00:00:06:342
Algorithms: Universal Global Optimation(UGO1)
Stop Reason: Convergence tolerance reached
Root of Mean Square Error (RMSE): 16.8998665796708
Residual Sum of Squares (SSE): 7996.95373149884
Correlation Coef. (R): 0.941605582652889
R-Square: 0.886621073283087
Determination Coef. (DC): 0.883913926992205
F-Statistic: 1.89550094274856
Parameters Best Estimate
-------------------- -------------
a 0.889811601928762
b 0.00162580722813606
c 1.15072355472944
d 0.00336596535381616
====== Output Results =====
File: Data file - 1
No Obs. y1 Cal. y1 Obs. y2 Cal. y2
1 64.083044982699 27.1022659715152 31.9723183391004 19.2079044663463
2 100.069204152249 40.8639262492391 34.7404844290657 27.199047
9997093
3 85.6747404844291 63.5888227024794 41.3840830449827 35.825427
767062
4 80.6920415224914 78.1728458807032 25.8823529411765 35.021544
9364615
5 76.2629757785467 84.8190532255592 22.560553633218 25.3497389960544
6 66.8512110726644 81.9846606375205 17.5778546712803 11.195662
6447939
7 59.6539792387543 74.6549102965342 0.968858131487885 -0.482941
881141063
8 55.7785467128028 65.0986508440345 -16.1937716262976 -9.972740
52598435
9 53.5640138408304 54.5596714437348 -28.3737024221453 -16.92488
18005384
10 53.5640138408304 47.3587938880552 -30.5882352941176 -20.20371
15454128
11 54.6712802768166 41.1038892035056 -34.4636678200692 -22.36185
86086396
12 48.5813148788927 37.4257598673193 -19.5155709342561 -23.58479
43182903
13 43.0449826989619 35.0029684495197 -10.1038062283737 -24.84909
53565655
14 43.598615916955 34.2715932747684 -4.01384083044982 -25.7687269913665

試用版：http://www.gaoyang168.com/bbs/viewthread.php?tid=8050245

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

3樓2015-06-19 17:08:40

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

查看全部 4 個(gè)回答

dingd

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

應(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
獨(dú)孤神宇: 金幣+5, ★★★★★最佳答案, 請(qǐng)問(wèn)怎么才算是高版本的，能否給個(gè)高版本的下載鏈接？我用的是5.0的，網(wǎng)上下的，應(yīng)該是盜版的了。。。 2015-06-18 20:47:34

盜版+版本低的原因：

均方差(RMSE):16.8998665796746
殘差平方和(SSE):7996.95373150243
相關(guān)系數(shù)(R): 0.64889751043271
相關(guān)系數(shù)之平方(R^2): 0.421067979045769
確定系數(shù)(DC): 0.0369393179593926
F統(tǒng)計(jì)(F-Statistic): 1.8955009427535

參數(shù) 最佳估算
-------------------- -------------
a 0.889811601926988
b 0.00162580722812602
c 1.15072355472983
d 0.00336596535380272

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No 目標(biāo) y1 計(jì)算 y1 目標(biāo) y2 計(jì)算 y2
1 64.083044982699 27.1022659714532 31.9723183391004 19.2079044662786
2 100.069204152249 40.8639262491168 34.7404844290657 27.1990479995981
3 85.6747404844291 63.588822702236 41.3840830449827 35.8254277668997
4 80.6920415224914 78.1728458803841 25.8823529411765 35.0215449363185
5 76.2629757785467 84.8190532252528 22.560553633218 25.3497389960066
6 66.8512110726644 81.9846606373267 17.5778546712803 11.1956626448538
7 59.6539792387543 74.6549102964756 0.968858131487885 -0.482941881027175
8 55.7785467128028 65.098650844115 -16.1937716262976 -9.97274052586747
9 53.5640138408304 54.5596714439447 -28.3737024221453 -16.9248818004721
10 53.5640138408304 47.3587938883414 -30.5882352941176 -20.203711545408
11 54.6712802768166 41.1038892038563 -34.4636678200692 -22.3618586087137
12 48.5813148788927 37.4257598677124 -19.5155709342561 -23.5847943184297
13 43.0449826989619 35.0029684499535 -10.1038062283737 -24.849095356772
14 43.598615916955 34.2715932752279 -4.01384083044982 -25.7687269916136

贊一下

回復(fù)此樓

2樓2015-06-18 20:28:41

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲(chóng)號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2015-06-19 17:08:40
代碼沒(méi)有問(wèn)題。
1stopt版本原因。低版本不支持SharedModel。
你說(shuō)自己用的是5.0版本，這是不可能的，網(wǎng)上能找到的最高到1.5。
你的問(wèn)題可以用試用版解決：

Iterations: 20
Elapsed Time (Hr:Min:Sec:Msec): ...

謝謝，網(wǎng)上下的，說(shuō)是5.0的，以前沒(méi)用過(guò)，也不知道是不是真的。。。
非常感謝提供的試用版本鏈接。。。

贊一下

回復(fù)此樓

數(shù)值計(jì)算

4樓2015-06-22 18:55:06

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

查看全部 4 個(gè)回答

普通表情龍兔虎貓

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

[考研] 268求調(diào)劑 +6	一定有學(xué)上- 2026-03-14	7/350	2026-03-17 10:14 by Sammy2
[考研] 一志愿，福州大學(xué)材料專(zhuān)碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 機(jī)械專(zhuān)碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[教師之家] 焦慮 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考研] 326求調(diào)劑 +4	上岸的小葡 2026-03-15	5/250	2026-03-16 08:39 by Linda Hu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 327求調(diào)劑 +4	Ffff03 2026-03-10	4/200	2026-03-14 00:17 by JourneyLucky
[考研] 材料工程,326分，求調(diào)劑 +6	KRSLSR 2026-03-10	6/300	2026-03-13 23:47 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 四川大學(xué)085601材料工程專(zhuān)碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[論文投稿] 投稿問(wèn)題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長(zhǎng)
[考研] 求調(diào)劑資源與環(huán)境 285 +3	未名考生 2026-03-10	3/150	2026-03-13 10:31 by houyaoxu
[考研] 290求調(diào)劑 +3	柯淮然 2026-03-10	8/400	2026-03-11 13:48 by 柯淮然