24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1557  |  回复: 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的回帖

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的回帖

wodaifei

银虫 (小有名气)

引用回帖:
Originally posted by pengyehui at 2011-06-29 22:07:59:
检查数据的输入吧!
三楼的建议值得商榷,因为原来问题已经可以线性化了,按三楼的方法,把问题搞成非线性的,复杂了!商榷而已!

这应该不是线性吧,否则应该能算出来,我就按照正确步骤一步一步来的,最后出现以上问题,我也知道是矩阵的问题,可小弟就是不知道该怎么改。。。
要想拥有一切,就要让自己变得足够优秀!!!
5楼2011-06-30 08:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pengyehui

木虫 (正式写手)

【答案】应助回帖


lovibond(金币+1): 鼓励交流 2011-07-01 08:32:14
令Y=lny; X1=1/x; X2=lnx
不就线性了吗
6楼2011-06-30 09:38:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodaifei

银虫 (小有名气)

引用回帖:
Originally posted by pengyehui at 2011-06-30 09:38:33:
令Y=lny; X1=1/x; X2=lnx
不就线性了吗

非常感谢啊,呵呵,小弟愚昧。。。
要想拥有一切,就要让自己变得足够优秀!!!
7楼2011-06-30 10:14:27
已阅   回复此楼   关注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的回帖

wodaifei

银虫 (小有名气)

引用回帖:
Originally posted by dingd at 2011-06-30 15:48:50:
如果原问题是非线性的,不建议线性化后求解,因为误差有可能非常大,以此题为例,线性化后公式及数据:
y1=a+b*x1+c*x2;
数据:
x1        x2        y1
0.00353169698039908        5.64597679255981        1.68268279794997
0.0034112 ...

哥,真的很感谢啊。。。你说的对,我做的时候也感觉有问题。。同学给的数据确实存在问题,不过方法是对的。。。再次感谢!!!
要想拥有一切,就要让自己变得足够优秀!!!
9楼2011-06-30 19:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wodaifei 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856,材料与化工321分求调剂 +7 大馋小子 2026-03-27 8/400 2026-03-27 22:31 by 无际的草原
[考研] 322求调剂 +4 宋明欣 2026-03-27 4/200 2026-03-27 22:02 by 催化大白
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 求调剂推荐 材料 304 +14 荷包蛋hyj 2026-03-26 14/700 2026-03-27 17:49 by kiokin
[考研] 085600,材料与化工321分调剂 +4 大馋小子 2026-03-27 6/300 2026-03-27 14:11 by 松花缸1201
[考研] 334求调剂 +3 雨清天晴 2026-03-21 3/150 2026-03-27 11:30 by 不吃魚的貓
[考研] 求调剂323材料与化工 +7 1124361 2026-03-24 7/350 2026-03-27 10:22 by wangjy2002
[考研] 材料学硕,求调剂 6+5 糖葫芦888ll 2026-03-22 10/500 2026-03-27 08:18 by hypershenger
[考研] 0703化学求调剂 +3 丹青奶盖 2026-03-26 5/250 2026-03-26 20:11 by macy2011
[考研] 化学工程085602 305分求调剂 +17 RichLi_ 2026-03-25 17/850 2026-03-26 19:44 by plmuchong
[考研] 085600 材料与化工 329分求调剂 +9 Mr. Z 2026-03-25 9/450 2026-03-26 10:36 by baoball
[考研] 0854AI CV方向招收调剂 +4 章小鱼567 2026-03-23 4/200 2026-03-25 17:04 by CoderLoser
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 化学调剂 +6 yzysaa 2026-03-21 6/300 2026-03-25 09:27 by aa331100
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +4 我爱学电池 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
信息提示
请填处理意见