24小时热门版块排行榜    

查看: 1558  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 不曾走远 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见