应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】多元方程组求解及拟合
汕头大学海洋科学接受调剂
51/1返回列表
查看: 2117  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lilei119

新虫 (正式写手)

[交流] 【求助】多元方程组求解及拟合 已有2人参与

多元非线性方程组：
(x10-x1)/y=abx1/(1+bx1+cx2+dx3+ex4)
(x20-x2)y=acx2/(1+bx1+cx2+dx3+ex4)
(x30-x3)/y=adx3/(1+bx1+cx2+dx3+ex4)
(x40-x4)/y=aex4/(1+bx1+cx2+dx3+ex4)
其中x10,x20,x30,x40 都是已知的常数，abcde是待拟合的常数，现已知一系列的x1,x2,x3,x4和y值，要拟合出上述abcde5个参数，上面的方程组看样子是无法手工消元然后用一元非线性拟合了，类似的文献说用牛顿-拉弗逊迭代法，各位专家帮看看这个问题解决起来的难度？并告知具体或者大概步骤。不甚感激
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-03-12 19:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+2): 谢谢应助！ 2011-04-10 21:42:38
本帖内容被屏蔽
4楼2011-04-10 21:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

forcal

禁虫 (初入文坛)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+2): 鼓励讨论交流！ 2011-03-12 20:20:51
zzuwangshilei(金币+1): 鼓励交流 2011-03-12 20:26:58
本帖内容被屏蔽
2楼2011-03-12 20:19:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilei119

新虫 (正式写手)
★
xiegangmai(金币+1): 鼓励讨论交流！ 2011-04-10 21:42:25
引用回帖:
Originally posted by forcal at 2011-03-12 20:19:36:
如果楼主方便，可以给出数据，没有数据，就不知道难度有多大。
一般的拟合软件都可以给出解决此类问题的代码，这没有任何难度。
Forcal解法参考：http://www.forcal.net/sysm/forcal9/fchtm/fcopt.htm

谢谢回复，
问题是这样的，
有一个方程组：

(x10-x1)/y=a*b*x1/(1+b*x1+c*x2+d*x3+e*x4)
(x20-x2)y=a*c*x2/(1+b*x1+c*x2+d*x3+e*x4)
(x30-x3)/y=a*d*x3/(1+b*x1+c*x2+d*x3+e*x4)
(x40-x4)/y=a*e*x4/(1+b*x1+c*x2+d*x3+e*x4)

x10,x20,x30,x40,y是自变量，x1,x2,x3,x4是因变量，a,b,c,d,e是常数，现已给出26组关于x1,x2,x3,x4x10,x20,x30,x40,y数据，要拟合a,b,c,d,e 这5个
参数。
拟合时的约束条件是让SSD这个值最小：

SSD = Σ[ (x1 – x1 (model))2 / (x10)2+ (x2- x2(model))2 / (x20)2+ (x3- x3(model))2 /(x30)2 + (x4- x4(model))2 / (x40) 2]
这里：x1(model)预测值，

下面用化学语言重新表达一遍上述问题。

y是吸附剂，x1,x2,x3,x4是4个被吸附物质，每个物质的吸附都符合Langmuir 模型，

即单独吸附x1的时候，满足公式：

(x10-x1)/y=a*b*x1/(1+b*x1)
，x10是x1的初始浓度，x1是x1吸附平衡时的浓度，y是吸附剂的投量。a，b是吸附常数。

现把x1,x2,x3,x4和y同时放入反应器，简单计算后各物质平衡浓度会同时满足以下4个方程：

(x10-x1)/y=a*b*x1/(1+b*x1+c*x2+d*x3+e*x4)
(x20-x2)y=a*c*x2/(1+b*x1+c*x2+d*x3+e*x4)
(x30-x3)/y=a*d*x3/(1+b*x1+c*x2+d*x3+e*x4)
(x40-x4)/y=a*e*x4/(1+b*x1+c*x2+d*x3+e*x4)
即：x10,x20,x30,x40是各物质的初始浓度，x1,x2,x3,x4吸附平衡后各物质浓度，a,b,c,d,e是吸附常数，现有26组实验，已测各物质的初始浓度和吸附平衡之后的浓度，现要拟合a,b,c,d,e 这5个
参数。

拟合时的约束条件是让SSD这个值最小：

SSD = Σ[ (x1 – x1 (model))2 / (x10)2+ (x2- x2(model))2 / (x20)2+ (x3- x3(model))2 /(x30)2 + (x4- x4(model))2 / (x40) 2]
这里：x1(model)预测值，

数据如下
Y        X10        X1        X20        X2        X30        X3        X40        X4
0.098709898        0        0        0        0        0.001071429        1.2854E-05        0        0
0.034516927        0        0        0        0        0.001071429        2.27537E-05        0        0
0.023548457        0        0        0        0        0.001071429        4.57405E-05        0        0
0.009339838        0        0        0        0        0.001071429        0.000158461        0        0
0.001373626        0        0        0        0        0.001071429        0.000714286        0        0
0.000446429        0        0        0        0        0.001071429        0.000928571        0        0
2.67857E-05        3.74667E-05        0.0000324        0.000103226        9.03226E-05        0        0        0        0
5.35714E-05        3.81333E-05        2.86667E-05        0.000103226        7.41935E-05        0        0        0        0
8.92857E-05        3.69333E-05        0.000022        0.0001        5.51613E-05        0        0        0        0
0.000142857        0.0000364        0.0000136        9.83871E-05        2.99355E-05        0        0        0        0
0.00018125        0.0000368        8.66667E-06        0.000102903        1.67742E-05        0        0        0        0
0.000196429        0.000036        5.73333E-06        9.6129E-05        9.6129E-06        0        0        0        0
0.000221429        3.54667E-05        3.06667E-06        9.45161E-05        4.80645E-06        0        0        0        0
0.000267857        0.000034        2.53333E-07        9.35484E-05        2.77419E-06        0        0        0        0
0.000175        6.66667E-06        3.08889E-07        9.03226E-05        2.49336E-06        0        0        0        0
0.0001125        6.66667E-06        2.21111E-06        9.03226E-05        2.32827E-05        0        0        0        0
3.39286E-05        6.66667E-06        4.85556E-06        9.03226E-05        6.57495E-05        0        0        0        0
7.14286E-05        6.66667E-06        3.67778E-06        9.03226E-05        4.78178E-05        0        0        0        0
0.00011        6.66667E-06        2.28889E-06        9.03226E-05        2.57306E-05        0        0        0        0
0.000151786        6.66667E-06        0.0000012        9.03226E-05        1.05882E-05        0        0        0        0
0.000210714        6.66667E-06        0.00000011        9.03226E-05        1.0019E-06        0        0        0        0
8.03571E-05        2.26667E-06        1.5569E-06        9.67742E-05        4.96451E-05        0.001066071        0.001066071        0.0000044        4.30976E-06
0.000139286        0.0000034        1.13333E-06        9.67742E-05        1.22203E-05        0.00105        0.00105        3.26667E-06        3.06667E-06
0.000214286        4.06667E-06        4.46554E-07        9.67742E-05        1.06928E-06        0.001028571        0.001028571        0.0000026        2.22222E-06
0.000321429        4.53333E-06        1.51089E-07        9.67742E-05        2.6058E-07        0.000985714        0.000985714        2.13333E-06        1.3468E-06
0.000535714        4.99333E-06        4.66667E-08        9.67742E-05        9.67742E-08        0.000942857        0.000942857        1.67333E-06        6.06667E-07
赞一下(1人) 回复此楼
3楼2011-04-10 08:18:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilei119

新虫 (正式写手)
引用回帖:
Originally posted by forcal at 2011-04-10 21:30:21:
看不出楼主的目标函数有多少意义，以下目标函数不可以吗？
[(x10-x1)/y-a*b*x1/(1+b*x1+c*x2+d*x3+e*x4)]^2
+[(x20-x2)/y-a*c*x2/(1+b*x1+c*x2+d*x3+e*x4)]^2
+[(x30-x3)/y-a*d*x3/(1+b*x1+c*x2+d*x3+e*x4)]^2 ...

谢谢你，此外，关于目标函数，文献是用的这个，所以想对比一下我和他的结果如何。如果方便实现的话能否帮弄弄按我我的目标函数再走一下？此外，最后那个误差怎么理解？如何用它来说明拟合的精确度？非常感谢，您帮了我的大忙。
此外，我下载了forcal，按您的程序计算后结果不显示，说not registered。请教如何注册？

[ Last edited by lilei119 on 2011-4-11 at 07:13 ]
赞一下 回复此楼
5楼2011-04-11 06:27:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308求调剂 +6 相信必会光芒万� 2026-04-11 6/300 2026-04-14 08:57 by licg0208
[考研] 调剂求收留 +32 果然有我 2026-04-10 33/1650 2026-04-14 08:49 by 木木mumu～
[考研] 279求调剂 +11 张番茄不炒蛋 2026-04-11 11/550 2026-04-13 22:25 by zhq0425
[考研] 复试调剂 +17 积极向上； 2026-04-10 19/950 2026-04-13 19:53 by 闲人终南山
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 53/2650 2026-04-13 14:20 by 张zhihao
[考研] 291分调剂 +5 上岸小莹加油 2026-04-09 6/300 2026-04-11 21:06 by 逆水乘风
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 材料类284调剂 +40 想换手机不想解� 2026-04-08 48/2400 2026-04-10 23:28 by 314126402
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 计算机类求调剂，22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 一志愿矿大，材料工程专硕314分，0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 一志愿沪9，326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 277、学硕，求调剂 数一104， +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
信息提示
关闭
请填处理意见
关闭