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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】多元方程组求解及拟合
51/1返回列表
查看: 2110  |  回复: 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):给个红包，谢谢回帖
臭水沟(金币+1): 谢谢参与！ 2011-05-26 22:59:42
本帖内容被屏蔽
9楼2011-05-26 18:21:58
已阅   回复此楼   关注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的回帖

forcal

禁虫 (初入文坛)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+2): 谢谢应助！ 2011-04-10 21:42:38
本帖内容被屏蔽
4楼2011-04-10 21:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 344 材料专业 求调剂211 无地域要求 +5 hualkop 2026-04-11 5/250 2026-04-11 23:13 by 852137818
[考研] 一志愿厦大0856，306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 266求调剂，一志愿哈工程电子信息，本科获多项国奖和省奖 +8 lumine1 2026-04-06 8/400 2026-04-11 18:35 by 逆水乘风
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 一志愿北理工298英一数二已上岸，感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 08600生物与医药-327 +10 18755400796 2026-04-05 10/500 2026-04-10 08:14 by kangsm
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 生物学308分求调剂（一志愿华东师大） +13 相信必会光芒万� 2026-04-06 16/800 2026-04-09 13:54 by 徐良白眉大侠
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 计算机408｜在校多次国家级竞赛获奖｜申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 求助 +3 卡卡东88 2026-04-06 4/200 2026-04-06 15:28 by going home
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
信息提示
关闭
请填处理意见
关闭