24小时热门版块排行榜    

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

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的回帖
查看全部 4 个回答

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的回帖

zhaoqian59

捐助贵宾 (小有名气)

这个我会,Matlab

[ 发自手机版 http://muchong.com/3g ]
爆炸力学
4楼2015-03-27 17:22:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 面上意见出来了 +10 黄鸟于飞Chao 2026-08-29 19/950 2026-08-31 08:01 by blueearth171
[考博] 找导师 +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 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +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 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +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
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见