应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 1stOpt 求助，求高手指导！！！！！
41/1返回列表
查看: 610  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

独孤神宇

版主 (知名作家)

[求助] 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的估值，请问这是怎么回事？ 是上面的代码哪里有问题吗？？？跪求指导！！！！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
数值计算
1楼2015-06-18 16:36:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
★
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
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2015-06-19 17:08:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)
引用回帖:
3楼: Originally posted by 月只蓝 at 2015-06-19 17:08:40
代码没有问题。
1stopt版本原因。低版本不支持SharedModel。
你说自己用的是5.0版本，这是不可能的，网上能找到的最高到1.5。
你的问题可以用试用版解决：

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

谢谢，网上下的，说是5.0的，以前没用过，也不知道是不是真的。。。
非常感谢提供的试用版本链接。。。
赞一下 回复此楼
数值计算
4楼2015-06-22 18:55:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 独孤神宇 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭