24小时热门版块排行榜    

查看: 391  |  回复: 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 的主题更新
信息提示
请填处理意见