查看: 1616  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 不曾走远 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 7/350 2026-10-12 14:47 by chuans0923
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见