24小时热门版块排行榜    

查看: 1546  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wodaifei

银虫 (小有名气)

[求助] 非线性回归参数..急求帮助!!谢谢

各位学哥学姐:小弟刚开始学matlab(菜鸟级别),遇到问题求帮助:
公式:lny=a+b/x+clnx                求解a,b,c
x        y       
283.15        5.37997
293.15        4.79126
303.15        4.23979        ,,,,,,,
313.15        3.80361       
323.15        3.42896       
333.15        3.10245       
343.15        2.81782       
353.15        0.03688       
按非线性方式回归最后运行出来出现:
Matrix dimensions must agree.
对于这个问题要想得到参数该怎么解决啊,
小弟很着急,各位学哥学姐帮帮忙啊
万分感谢!!!!
回复此楼
要想拥有一切,就要让自己变得足够优秀!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

wodaifei(金币+1): 2011-06-30 19:20:44
如果原问题是非线性的,不建议线性化后求解,因为误差有可能非常大,以此题为例,线性化后公式及数据:
y1=a+b*x1+c*x2;
数据:
x1        x2        y1
0.00353169698039908        5.64597679255981        1.68268279794997
0.00341122292341805        5.68068442341022        1.56679342484846
0.0032986970146792        5.71422773251774        1.44451373971541
0.00319335781574325        5.74668230897142        1.33595061576793
0.00309453814018258        5.7781166117089        1.23225700820423
0.0030016509079994        5.80859283900855        1.13219212193131
0.00291417747340813        5.83816766935464        1.03596353634767
0.00283165793572137        5.86689289585493        -3.30008588029374

拟合结果:
a        2392.53187520471
b        -107971.570616336
c        -356.014799219908

直接非线性拟合:
公式: y=exp(a+b/x+c*ln(x))
数据:
x        y
283.15        5.37997
293.15        4.79126
303.15        4.23979
313.15        3.80361
323.15        3.42896
333.15        3.10245
343.15        2.81782
353.15        0.03688

结果:
a        392.19036548745
b        -16708.6338358188
c        -58.7228208089647

非线性直接误差平方和:3.26128834414432
线性化后求得参数代入原方程误差平方和: 30.7310566203494

误差是惊人的。
8楼2011-06-30 15:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

xiuyouxu

铁杆木虫 (职业作家)

【答案】应助回帖


soliton923(金币+1): 谢谢参与~~ 2011-06-29 20:17:53
wodaifei(金币+1): 2011-06-30 10:14:53
Matrix dimensions must agree.

矩阵维数不一致, 检查矩阵运算是否有问题.
忘记自己,忘记一切烦恼(欢迎访问我的网站兆字节:http://www.mathbeta.com/)
2楼2011-06-29 16:44:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


soliton923(金币+1): 谢谢参与~~ 2011-06-29 20:18:05
wodaifei(金币+1): 2011-06-30 08:39:53
你把方程“lny=a+b/x+clnx”改为:

y=exp(a+b/x+c*ln(x))

试试!
3楼2011-06-29 20:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pengyehui

木虫 (正式写手)

【答案】应助回帖

wodaifei(金币+2): 2011-06-30 10:14:38
检查数据的输入吧!
三楼的建议值得商榷,因为原来问题已经可以线性化了,按三楼的方法,把问题搞成非线性的,复杂了!商榷而已!
4楼2011-06-29 22:07:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
请填处理意见