24小时热门版块排行榜    

查看: 437  |  回复: 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的回帖

zhi123chao

铜虫 (小有名气)

WWJ1991

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2013-09-30 20:56:15

dingd

铁杆木虫 (职业作家)

相关版块跳转 我要订阅楼主 zhi123chao 的主题更新
信息提示
请填处理意见