24小时热门版块排行榜    

查看: 461  |  回复: 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文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:24 by jHcrHMhbKCFV
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:23 by jHcrHMhbKCFV
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:21 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[基金申请] 面上提前没消息,有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 8/400 2026-08-02 15:01 by jHcrHMhbKCFV
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 13:34 by jHcrHMhbKCFV
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
信息提示
请填处理意见