24小时热门版块排行榜    

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

zhaoqian59

捐助贵宾 (小有名气)

这个我会,Matlab

[ 发自手机版 http://muchong.com/3g ]
爆炸力学
4楼2015-03-27 17:22:52
已阅   回复此楼   关注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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 8/400 2026-08-10 22:07 by demonstreets
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +3 Tide man 2026-08-10 3/150 2026-08-10 21:48 by GZUZM
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见