24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 286  |  回复: 1

blesswj

金虫 (小有名气)

[求助] 非线性拟合问题求助!!!

不知道哪里出现问题了,总是运行不成功哦!有高人指点下吗?谢谢
代码如下:前两个函数都可以,后面两个总是不行!

clear all;
close all;
clc;
datas=[5.82        0.637
8.82        0.556
11.82        0.46
14.82        0.355
17.82        0.419
20.82        0.34
23.82        0.354
26.82        0.332
29.82        0.308
32.82        0.272
35.82        0.36
38.82        0.303
41.82        0.302
44.82        0.317
47.82        0.385
50.82        0.357
53.82        0.299
56.82        0.324
59.82        0.33];
x=datas(:,1);y=datas(:,2);
%4中函数形式分别为:

%y1=a*x^(-b)---Kostiakow公式
fun_1=inline('parm_1(1)*x.^(-parm_1(2))','parm_1','x');%inline用来定义函数
parm_1=nlinfit(x,y,fun_1,[0.01 0.01]);

%y2=fc+(f0-fc)*exp(-kx)---Horton公式
fun_2=@(parm_2,x) parm_2(1)+(parm_2(2)-parm_2(1))*exp(-parm_2(3)*x);%@用来定义函数
parm_2=lsqcurvefit(fun_2,[0.01 0.01 0.01],x,y);

%y3=k+k1/x^alpha---fang公式
fun_3=inline('parm_3(1)+parm_3(2)/x.^parm_3(3)','parm_3','x');
parm_3=nlinfit(x,y,fun_3,[0.01 0.01 0.01]);

%y4=fc+(f1-fc)/x^alpha---jiang公式
fun_4=@(parm_4,x) parm_4(1)+(parm_4(2)-parm_4(1))/x.^parm_4(3);
parm_4=lsqcurvefit(fun_4,[0.01 0.01 0.01],x,y);
回复此楼

» 猜你喜欢

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

多累的时间,想想未来,忍一忍,走下去......
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

blesswj

金虫 (小有名气)

问题已经解决,请斑竹关贴吧,谢谢!!!
多累的时间,想想未来,忍一忍,走下去......
2楼2013-01-06 09:11:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 blesswj 的主题更新
信息提示
请填处理意见