24小时热门版块排行榜    

查看: 1580  |  回复: 3

不曾走远

新虫 (初入文坛)

[求助] 跪求: matlab三元非线性回归,请指教。。。

跪求: matlab三元非线性回归,请指教。。。
有一组数据,进行三元非线性回归分析,Matlab怎么编程,有识之虫请指教,多谢!
y=a(1)*x1.^(a(2)).*(1+a(3)*ln(x2)).*exp(-a(4)*x3))
x1、x2、x3,为三个自变量,a(1)、a(2)、a(3)、a(4)为待定系数。
数据取了一部分,如下:
  y       x1    x2     x3
644.8 0.78 23795 143
633.1 0.80 16673 140
678.0 0.80 47284 149
677.5 0.80 79812 152
621.7 0.81 10752 137
658.5 0.82 24781 150
629.0 0.83 5126  137
641.9 0.83 14673 143
683.1 0.83 32167 152
632.6 0.86 9025  141

跪求用matlab将其4个待定系数求出来!多谢大侠相助!
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
不曾走远: 金币+5, ★★★很有帮助 2012-10-15 22:26:58
不曾走远: 金币+3, ★★★很有帮助 2012-10-15 22:27:51
拟合问题用1stOpt更简单些,不用猜初值:
CODE:
Variable y,x1,x2,x3;
Function y=a1*x1^(a2)*(1+a3*ln(x2))*exp(-a4*x3);
Data;
//y       x1    x2     x3
644.8 0.78 23795 143
633.1 0.80 16673 140
678.0 0.80 47284 149
677.5 0.80 79812 152
621.7 0.81 10752 137
658.5 0.82 24781 150
629.0 0.83 5126  137
641.9 0.83 14673 143
683.1 0.83 32167 152
632.6 0.86 9025  141

均方差(RMSE): 5.82535725444792
残差平方和(SSE): 339.34787141949
相关系数(R): 0.96266480547321
相关系数之平方(R^2): 0.926723527696774
决定系数(DC): 0.926722803069674
卡方系数(Chi-Square): 0.257128283384663
F统计(F-Statistic): 27.2936204410447

参数        最佳估算
----------        -------------
a1        293.462825876497
a2        -0.0349876564552783
a3        0.00423477464786676
a4        -0.00517038892366596

jv1.jpg

2楼2012-10-15 16:33:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不曾走远

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2012-10-15 16:33:36
拟合问题用1stOpt更简单些,不用猜初值:

Variable y,x1,x2,x3;
Function y=a1*x1^(a2)*(1+a3*ln(x2))*exp(-a4*x3);
Data;
//y       x1    x2     x3
644.8 0.78 23795 143
633.1 0.80 16673 140
678.0 0 ...

谢谢啊!  这种方法确实很方便。(matlab也应该能做出来,可能要麻烦些)。多谢了!
3楼2012-10-15 18:23:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不曾走远

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2012-10-15 16:33:36
拟合问题用1stOpt更简单些,不用猜初值:

Variable y,x1,x2,x3;
Function y=a1*x1^(a2)*(1+a3*ln(x2))*exp(-a4*x3);
Data;
//y       x1    x2     x3
644.8 0.78 23795 143
633.1 0.80 16673 140
678.0 0 ...

谢谢啊!  这种方法确实很方便。(matlab也应该能做出来,可能要麻烦些)。多谢了!
4楼2012-10-15 22:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 不曾走远 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见