24小时热门版块排行榜    

查看: 417  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 10/500 2026-09-06 20:08 by 阿美_Lml888
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 6/300 2026-09-06 18:41 by eObJPa4gWmp6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见