24小时热门版块排行榜    

查看: 520  |  回复: 1

fwmx851

铜虫 (初入文坛)

[求助] 求助1stopt 拟合求参数 已有1人参与

Variable: a,g
Parameter b,c,d,e,f,h,n[0,1],i
Function a=(50/c)+((50/d)*Sum(i=1:9)(((-e/d)^i)*((g/1)^n)*(1+i))/(n*(1+i)*Gamma((1+i)*n)))+(25/2*Sqrt((1/c+(e+d*(2*pi*1)^n*Cos(n*pi/2))/(e^2+d^2*(2*pi*1)^(2*n)+2*e*d*(2*pi*1)^(n)*Cos(n*pi/2))^2+d)^2+((d*(2*pi*d)^(n)*Sin(n*pi/2))/(e^2+d^2*(2*pi*1)^(2*n)+2*e*d*(2*pi*1)^(n)*Cos(n*pi/2)))^2))+(Sin(2*pi*g+Atan((1/c+(e+d*(2*pi*1)^n*Cos(n*pi/2))/(e^2+d^2*(2*pi*1)^(2*n)+2*e*d*(2*pi*1)^(n)*Cos(n*pi/2)))/(d*(2*pi*d)^(n)*Sin(n*pi/2))/(e^2+d^2*(2*pi*1)^(2*n)+2*e*d*(2*pi*1)^(n)*Cos(n*pi/2)))))
Data;
10 0.0034
30 0.0047
40 0.0053
50 0.0059
70 0.0069
100 0.0082
200 0.0113
300 0.013
500 0.0143
800 0.0148
1000 0.0149

参数超过了4个,求大神们帮帮忙,谢谢谢谢!n限制在0,1,中间用到Gamma函数,sum取到无穷,自己不知道多少比较好就取了个9,如果有什么问题希望大神们和我讲,谢谢谢谢!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
独孤神宇: 金币+1 2018-05-01 21:13:30
fwmx851: 金币+5, ★★★★★最佳答案, 嗯嗯,我再修正一下,麻烦您啦 2018-05-02 17:13:25
效果不好吧:

均方差(RMSE): 161.35267143553
残差平方和(SSR): 286381.530373202
相关系数(R): 0.868600219805926
相关系数之平方(R^2): 0.754466341846903
修正R平方(Adj. R^2): 0.693082927308629
确定系数(DC): 0.754466341846826
卡方系数(Chi-Square): 18916.8081807556
F统计(F-Statistic): 1.31689689297483

参数        最佳估算
----------        -------------
c        -0.00551746295993722
d        0.0100251834622464
e        -0.00841885439178249
n        0.999999999999998
i        0.666313686294838
2楼2018-05-01 11:58:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fwmx851 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 06:39 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 06:13 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 05:54 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:37 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 X5bHxOULxFqb 2026-10-30 27/1350 2026-11-02 05:28 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见