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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】多元方程组求解及拟合
51/1返回列表
查看: 2478  |  回复: 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的回帖

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 老师的回复，约束方程稍微改了该，能否再耽误您时间帮忙看看？

(x10-x1)/y=a*b*x1/(1+b*x1+c*x2+d*x3+e*x4)    (1)
(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)     (3)
(x40-x4)/y=a*e*x4/(1+b*x1+c*x2+d*x3+e*x4)     (4)

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 + (x2- x2(model))2 + (x3- x3(model))2  + (x4- x4(model))2 ]
SSD = Σ[ (x1 – x1 (model))2 + (x2- x2(model))2 + (x3- x3(model))2  + (x4- x4(model))2 ]
这里：x1,x2,x3,x4是实验测出来的值（数据见下），x1(model)，x2(model)，x3(model)，x4(model)是指上面方程组(方程（1）（2）（3）（4））根据给出的数据（数据见下）y, x10,x20,x30,x40和5个参数求出的解。此外，求解该方程组会有多组解，所以还需要限制条件（即 0 不知道这个是否能实现，非常感谢您的帮助！！




数据如下
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

[ Last edited by lilei119 on 2011-5-25 at 10:16 ]
赞一下 回复此楼
8楼2011-05-25 09:59:18
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭