| 查看: 658 | 回复: 3 | ||||
独孤神宇版主 (知名作家)
|
[求助]
1stOpt 求助，求高手指导！！！！！ 已有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 运行后，只有y1拟合图像和a,b的估值， y2为一条等于零的直线，没有c和d的估值，请问这是怎么回事？ 是上面的代码哪里有问题吗？？？跪求指导！！！！！ |
回复此楼» 猜你喜欢
291分调剂
已经有9人回复
-
调剂求收留
已经有34人回复
-
291 求调剂
已经有38人回复
-
22408 312求调剂
已经有17人回复
-
一志愿华中农业071010，320求调剂
已经有6人回复
-
290调剂生物0860
已经有41人回复
-
291求调剂
已经有3人回复
-
211本科材料化工求调剂
已经有23人回复
-
山东省基金2026
已经有9人回复
-
药学求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stOpt 求助！！！
已经有4人回复
- 用高版本1stOpt求解方程组
已经有6人回复
- 1stopt 语句修改
已经有23人回复
- matlab高手进！发现一个奇怪的现象……
已经有4人回复
- 求matlab大神指点，关于logistic拟合
已经有13人回复
- 求Matlab高手指导
已经有12人回复
- 看了下1stOpt的介绍，感觉略显浮夸
已经有14人回复
- 1stopt运行报错“参数虽被定义，但是没有在函数表示中出现”，顺便求高手帮忙跑下程序
已经有3人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有6人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 1stopt拟合的动力学参数相关系数不高，有哪些原因
已经有3人回复
- 跪求高手帮忙看看这组数据如何进行拟合
已经有7人回复
- ELK简单计算失败，求高手解答！
已经有4人回复
- 求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢
已经有5人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 求用高版本的1stOpt拟合方程
已经有8人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- 求SWMM软件高手指导
已经有15人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
独孤神宇: 金币+5, ★★★★★最佳答案, 请问怎么才算是高版本的，能否给个高版本的下载链接？ 我用的是5.0的，网上下的，应该是盗版的了。。。 2015-06-18 20:47:34
感谢参与，应助指数 +1
独孤神宇: 金币+5, ★★★★★最佳答案, 请问怎么才算是高版本的，能否给个高版本的下载链接？ 我用的是5.0的，网上下的，应该是盗版的了。。。 2015-06-18 20:47:34
|
盗版+版本低的原因： 均方差(RMSE):16.8998665796746 残差平方和(SSE):7996.95373150243 相关系数(R): 0.64889751043271 相关系数之平方(R^2): 0.421067979045769 确定系数(DC): 0.0369393179593926 F统计(F-Statistic): 1.8955009427535 参数 最佳估算 -------------------- ------------- a 0.889811601926988 b 0.00162580722812602 c 1.15072355472983 d 0.00336596535380272 ====== 结果输出 ====== 文件: 数据文件-1 No 目标 y1 计算 y1 目标 y2 计算 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 |
2楼2015-06-18 20:28:41
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
★
xzhdty: 金币+1, 谢谢参与 2015-06-19 17:12:52
xzhdty: 金币+1, 谢谢参与 2015-06-19 17:12:52
|
代码没有问题。 1stopt版本原因。低版本不支持SharedModel。 你说自己用的是5.0版本，这是不可能的，网上能找到的最高到1.5。 你的问题可以用试用版解决： 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 试用版：https://muchong.com/bbs/viewthread.php?tid=8050245 |
3楼2015-06-19 17:08:40
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
4楼2015-06-22 18:55:06
5