24小时热门版块排行榜    

查看: 742  |  回复: 7

m452193160

新虫 (小有名气)

[求助] 数据拟合求助 已有2人参与

x=1:6;
y=[0.105935612        0.042331078        0.024863921        0.028627497        0.033000851        0.063137188];
f=inline('a(1)./(a(2)+x)-a(3)./(a(4)+x)+a(5)','a','x');
[a,res]=lsqcurvefit(f,[1,1,1,1,1],x,y);


为什么matlab一直错误。。。求助。。。

发自小木虫Android客户端
回复此楼

» 猜你喜欢

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

m452193160

新虫 (小有名气)

如果模型是  a(1)x/a(2)^x也是可以的,但是为什么都不能进行拟合。。。

发自小木虫Android客户端
2楼2016-12-28 22:08:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luoshengfeng

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
我这里是可以的 (matlab2013a) 以下是结果
>> x=1:6;
>> y=[0.105935612        0.042331078        0.024863921        0.028627497        0.033000851        0.063137188];
>> f=inline('a(1)./(a(2)+x)-a(3)./(a(4)+x)+a(5)','a','x');
>> [a,res]=lsqcurvefit(f,[1,1,1,1,1],x,y)

Local minimum possible.

lsqcurvefit stopped because the final change in the sum of squares relative to
its initial value is less than the default value of the function tolerance.

<stopping criteria details>


a =

    0.6771    0.8328    0.6623    1.1782    0.0258


res =

    0.0018

>>

» 本帖已获得的红花(最新10朵)

知行合一
3楼2016-12-29 08:33:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

m452193160

新虫 (小有名气)

送红花一朵
引用回帖:
3楼: Originally posted by luoshengfeng at 2016-12-29 08:33:42
我这里是可以的 (matlab2013a) 以下是结果
>> x=1:6;
>> y=;
>> f=inline('a(1)./(a(2)+x)-a(3)./(a(4)+x)+a(5)','a','x');
>> =lsqcurvefit(f,,x,y)

Local minimum possible.

...

是跟版本有问题吗我是2010r

发自小木虫Android客户端
4楼2016-12-29 08:35:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小圆盾

铜虫 (小有名气)

没有问题,是楼上的结果

发自小木虫IOS客户端
5楼2016-12-29 08:56:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
用1stOpt吧,Matlab得到的明显是不对,是局部最优:

均方差(RMSE): 0.0016219348788169
残差平方和(SSE): 1.57840365067368E-5
相关系数(R): 0.998332583373811
相关系数之平方(R^2): 0.996667947025828
修正R平方(Adj. R^2): 0.99444657837638
确定系数(DC): 0.996667947025828
卡方系数(Chi-Square): 0.000286448253875885
F统计(F-Statistic): 74.7788190961433

参数        最佳估算
----------        -------------
a1        0.0924403935311779
a2        -0.228191932872472
a3        0.11076326893216
a4        -7.41895318645034
a5        -0.0310503432998637
6楼2016-12-29 09:23:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

m452193160

新虫 (小有名气)

引用回帖:
6楼: Originally posted by dingd at 2016-12-29 09:23:03
用1stOpt吧,Matlab得到的明显是不对,是局部最优:

均方差(RMSE): 0.0016219348788169
残差平方和(SSE): 1.57840365067368E-5
相关系数(R): 0.998332583373811
相关系数之平方(R^2): 0.996667947025828
修正 ...

不知是否好用,从来没用过。。。

发自小木虫Android客户端
7楼2016-12-29 09:29:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

拟合问题估计目前还没有比得上1stOpt的,不论是效果还是易用性。
8楼2016-12-29 09:52:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 m452193160 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 05:13 by gy1nBQXYQJqL
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 04:18 by gy1nBQXYQJqL
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 04:06 by gy1nBQXYQJqL
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 3/150 2026-08-29 03:23 by gy1nBQXYQJqL
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见