24小时热门版块排行榜    

查看: 830  |  回复: 3

zhaoqian59

捐助贵宾 (小有名气)

[求助] matlab非线性拟合模型参数 已有2人参与

matlab拟合,求参数x1~x8
y=den*(1-x1*ω)(x2*den+x3* a+x4* b+x5* c+x6* d+x7 *ω+x8* hf)
den        a        b        c        d        e        hf        w        y
1.673        3        6        6        6        0        70.71        -48        29.39
1.763        1.42        3.04        2.16        2.16        0.56        26.34        -50.8        23.91
1.865        1.22        2.23        1.77        1.77        1.11        37.41        -55.29        22.21
1.816        3        6        6        6        0        70.71        -48        32.6
1.667        3.28        2.81        1.26        2.51        0        55.08        -87.28        23.31
1.72        2.97        2.59        2.25        2.25        0.37        48.01        -88.64        22.58
1.788        1.37        2.7        1.98        1.98        1.74        40.94        -51.52        22.04
1.853        1.23        2.44        1.72        1.72        1.11        33.87        -52.24        20.79
1.891        4        8        8        8        0        75.02        -64        41.91
1.93        1.22        2.43        2.43        2.43        0.37        67.52        -28.48        38.29
1.971        1.08        2.16        2.16        2.16        0.72        60.02        -34.56        34.52
2.009        0.94        1.89        1.89        1.89        1.11        52.51        -41.6        30.9
2.047        0.81        1.62        1.62        1.62        1.48        45.01        -48.48        26.97
1.787        3        6        6        6        0        70.71        -48        34.35
1.83        1.22        2.43        2.43        2.43        0.37        63.64        -28.48        31.42
1.872        1.08        2.16        2.16        2.16        0.72        56.57        -34.56        28.53
1.921        0.95        1.89        1.89        1.89        1.11        49.5        -41.92        25.49
1.975        0.81        1.62        1.62        1.62        1.48        42.43        -48.48        22.42
1.128        1        3        1        2        0        -112.6        -24        12.83
1.177        1.47        4.42        1.47        2.95        0.37        -101.34        -44.08        12.08
1.226        1.31        3.93        1.31        2.62        0.74        -90.08        -49.2        11.36
1.286        1.15        3.44        1.15        2.29        1.11        -78.82        -54.32        10.49
1.344        0.98        2.95        0.98        1.97        1.48        -67.56        -58.96        9.52
1.63        7        5        3        6        0        -65.5        -168        20.56
1.652        2.77        1.98        1.19        2.38        0.37        -58.95        -75.28        19.37
1.675        2.46        1.76        1.06        2.11        0.74        -52.4        -76.8        18.08
1.699        2.16        1.54        0.92        1.84        1.11        -45.85        -78.64        16.67
1.726        1.84        1.23        0.79        1.58        1.48        -39.3        -78.96        15.1
回复此楼

» 猜你喜欢

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

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

zhangzhiguang

至尊木虫 (文坛精英)

【答案】应助回帖

感谢参与,应助指数 +1
楼主给的数据有的不是数值格式。我把程序给你,你运算一下。或者你弄个xls格式的数据也可以。另外,好像所给数据中没有用到e。参考这个帖子http://muchong.com/bbs/viewthread.php?tid=8239371&authorid=492270
先编写函数
function y=myfunc(beta,x)
x1=beta(1);
x2=beta(2);
x3=beta(3);
x4=beta(4);
x5=beta(5);
x6=beta(6);
x7=beta(7);
x8=beta(8);
den=x(:,1);
a=x(:,2);
b=x(:,3);
c=x(:,4);
d=x(:,5);
hf=x(:,6);
w=x(:,7);
y=den.*(1-x1.*w).*(x2.*den+x3.*a+x4.*b+x5.*c+x6.*d+x7.*w+x8.*hf);
再在命令窗口中输入
>> x=xlsread('你的数据,每行为一个自变量.xls');
>> y=数据;
>> beta0=[1,1,1,1,1,1,1]';
>> [beta,r,j]=nlinfit(x,y,'myfunc',beta0)
你自己先试试吧,不行再说。不过matlab的非线性拟合依赖于初值的选择。
好好学习,天天向上!
2楼2015-03-27 14:55:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1stOpt计算:

Root of Mean Square Error (RMSE):1.53722146173097
Sum of Square Error:66.1653950273761
Correlation Coef. (R): 0.98344679661229
R-Square: 0.967167601766974
Determination Coef. (DC): 0.967165369235768
Chi-Square: 1.2611968147024
F-Statistic: 87.0161409049629

Parameter        Best Estimate
----------        -------------
x1        -0.00146362899420343
x2        9.85830995434919
x3        5.05749953752554
x4        1.04663297028699
x5        -1.528254895988
x6        0.0460108206492992
x7        0.209167180906584
x8        0.011540122741908
matlab非线性拟合模型参数
c1.jpg

3楼2015-03-27 16:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoqian59

捐助贵宾 (小有名气)

这个我会,Matlab

[ 发自手机版 http://muchong.com/3g ]
爆炸力学
4楼2015-03-27 17:22:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhaoqian59 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 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 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (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
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见