24小时热门版块排行榜    

查看: 469  |  回复: 2

坤少pop

新虫 (正式写手)

[求助] 1stopt运行不了,望大家赐教! 已有2人参与

Title "非线性拟合";Parameters b(1:10);SharedModel;Variable x,y1,y2,y3;Function y1 = b7-((1/b8+b5/b9)*((b1*b2*exp(b3*x)/(b2-b1+b1*exp(b3*x)))-b1))-(b2/b3)*(b6/b9+b10)*log((b2-b1+b1*exp(b3*x))/b2);y2 = b4+b5*((b1*b2*exp(b3*x)/(b2-b1+b1*exp(b3*x)))-b1)+(b6*b2/b3)*log((b2-b1+b1*exp(b3*x))/b2);y3 = b1*b2*exp(b3*x)/(b2-b1+b1*exp(b3*x));Data;//x,y1,y2,y30        28.21        0.39        1.156        27.65        0.92        2.212        26.96        1.97        4.118        22.88        2.41        7.2524        18.04        2.92        9.130        12.81        3.2        10.0536        8.71        3.42        10.6542        6.59        3.42        10.748        6.22        3.68        10.7554        5.86        3.43        1160        5.74        3.44        10.766        5.7        3.52        10.672        5.43        3.62        10.4578        5.26        3.52        10.184        4.14        3.41        10

@月只蓝 发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
参见以下链接中第5条:http://muchong.com/bbs/viewthread.php?tid=9176191
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2017-04-28 10:51:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by 月只蓝 at 2017-04-28 10:51:23
参见以下链接中第5条:http://muchong.com/bbs/viewthread.php?tid=9176191

仅供参考

迭代数: 21
计算用时(时:分:秒:微秒): 00:00:08:165
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判定标准
均方差(RMSE): 161.57197020968
残差平方和(SSE): 939798.056067761
相关系数(R): 0.616840194121187
相关系数之平方(R^2): 0.380491825083464
决定系数(DC): 0.380504389127482
F统计(F-Statistic): 6.141837032816

参数                  最佳估算
--------------------        -------------
b1                 2.33716162201893E-19
b2                 17.3958605321308
b3                 8.46776051292736
b4                 -2.0582990250142
b5                 21.1186527230704
b6                 -1.39858295194152
b7                 10.890468507053
b8                 -34615.4838100747
b9                 -279078.556769959
b10                0.0290034992138533

====== 输出结果 =====

No        实测y1        计算y1        实测y2        计算y2        实测y3        计算y3
1        28.21        10.890468507053        0.39        -2.0582990250142        1.156        2.33716162201893E-19
2        0.92        -0.333867754925933        2.212        -175.926090470188        26.96        17.3958605321308
3        4.12        10.890468507053        22.88        -2.05829902492817        2.41        4.10550916516905E-12
4        18.04        9.96054354576693        2.92        320.396750581853        9.13        17.3958577116967
5        3.2        7.15481184556644        10.0536        185.124490512736        8.71        17.3958605321308
6        10.65        10.8904685041219        6.59        -2.05828053229706        3.42        8.82560361355933E-7
7        6.22        8.19434553673991        3.68        235.243371700392        10.7554        17.3958605321308
8        3.43        10.6607023622293        1160        346.60961417123        5.74        17.0385627076581
9        10.77        10.890468503581        5.7        -2.05827711963598        3.52        1.04542876899124E-6
10        5.43        8.23471577717384        3.62        237.189736018359        10.4578        17.3958605321308
11        3.52        10.8754135740444        10.184        81.1938854301014        4.14        3.97743539181013
12        10        10.8904685043599        0        -2.05828203374925        0        8.10903926099291E-7
1stopt运行不了,望大家赐教!
y1.png


1stopt运行不了,望大家赐教!-1
y2.png


1stopt运行不了,望大家赐教!-2
y3.png

一团和气
3楼2017-05-02 10:56:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 坤少pop 的主题更新
信息提示
请填处理意见