24小时热门版块排行榜    

查看: 440  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhi123chao

铜虫 (小有名气)

[求助] 非线性拟合结果不好,求大神给修正意见

需通过非线性拟合出a和k的值。
其中,给出一组t和x相对应的值。
x=1-(1+(a-1)*k*14000^(a-1)*t)^(1/(1-a))
上式是已知的方程式,t和x值如下:
t=[3,6,9,15];
y=[0.838,0.8627,0.8688,0.8699];

我写的编程如下:
t=[3,6,9,15];
y=[0.838,0.8627,0.8688,0.8699];
myfunc=inline('1-(1+(beta(1)-1).*beta(2).*14000.^(beta(1)-1).*t).^(1/(1-beta(1)))','beta','t');
beta=nlinfit(t,y,myfunc,[0,0]);
a=beta(1),k=beta(2)
%test the model
tt=min(t):max(t);
yy=1-(1+(a-1)*k*14000^(a-1)*t)^(1/(1-a));
plot(t,y,'o',t,yy,'r')

但是拟合结果是非常不好,拟合的a和k值都出现虚数,弄好久了,不知怎么解决,求大神了。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

1stOpt求解:
CODE:
Variable t,x;
Function x=1-(1+(a-1)*k*14000^(a-1)*t)^(1/(1-a));
Data;
t=[3,6,9,15];
y=[0.838,0.8627,0.8688,0.8699];

均方差(RMSE):0.00459568570854431
残差平方和(SSE):8.44813085268739E-5
相关系数(R): 0.934583731677735
相关系数之平方(R^2): 0.873446751516681
确定系数(DC): 0.873274869081417
卡方系数(Chi-Square): 4.90332655183441E-5
F统计(F-Statistic): 15.782189258774

参数        最佳估算
----------        -------------
a        7.79553088593713
k        9.05778386695621E-25

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

No        实测值x        计算值x
1        0.838        0.841783499215192
2        0.8627        0.857125849902977
3        0.8688        0.865401290924269
4        0.8699        0.875148244423607
4楼2014-04-03 23:48:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

WWJ1991

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2013-09-30 20:56:15
CFTOOL工具箱不可以做吗?
3楼2013-09-30 16:39:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见