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

24小时热门版块排行榜    

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

夜雨tju

金虫 (小有名气)

[求助] matlab求解非线性方程组 已有2人参与

方程组有36个，未知数为x1-x34,xa,xb，除了未知数，其他都已知，其中两个方程组是定积分，两个上限为未知量，分别为,xa,xb，下限皆为0.初值为
x0=[0.036 0.055 0.04 0.032 0.028 0.031 0.06 0.05 0.05 0.03 360 0.04 0.06 0.05 0.034 0.028 0.03 0.062 0.05 0.055 0.04 365 0.91 0.87 0.89 0.86 0.9 0.88 0.85 0.89 0.88 0.9 5900 20]；
xa0=0.01;
xb0=0.01;
开始做了不合理的简化，只解了前f(35),f(36),f(1)-f(11)，解释合理的，现在不做简化，四个组分一起列式子，方成数变多，方程组总是不能得到合理的解，请高手指点，如何解这个方程组，小女子不胜感激，由于积分显示不出来，我上传了附件和图片。谢谢各位大神
fff(35)=( K/k1)*∫_0^xa▒〖CA0/(K*CA0*CB0-CC0*CD0-(K*(CA0^2+CA0*CB0)+(CA0*CC0+CA0*CD0)).*x+(K-1)*CA0^2.*x.^2)dx〗 -Lc1/U
fff(36)= =( K/k1)*∫_0^xb▒〖x(11)/(K*x(11)*x(22)-x(33)*x(34)-(K*(x(11)^2+x(11)*x(22))+(x(11)*x(33)+x(11)*x(34))).*x+(K-1)*x(11)^2.*x.^2)dx〗 –Lc2/U

fff(1)=xa-(xa1-x(8))/xa1;
fff(2)=h1/HETP1-(x(4)-x(6))/(K1*x(7)-x(6));
fff(3)=L*(S1+S3+S5)/V*S-(x(4)-x(6))/(xa1-x(7));
fff(4)=x(2)-(x(7)*(S1+S3+S5)+x(8)*(S2+S4))/S;
fff(5)=h2/HETP2-mi*(x(6)-x(5))/(K1*x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(5));
fff(6)=L/V-mi*(x(6)-x(5))/(mi*x(2)-x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(7)=h3/HETP3-(HETP4*x(5)-h4*(K1*x(1)-ya1)-ya1*HETP4)/(K1*HETP4*x(9)-h4*(K1*x(1)-ya1)-ya1*HETP4);
fff(8)=L*(S2+S4)/V*S-(HETP4*x(5)-h4*(K1*x(1)-ya1)-ya1*HETP4)*mi/(HETP4*(xa1*M1+xb1*M2+xc1*M3+xd1*M4-mi*x(9)));
fff(9)=xb-(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(10))/(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(10)=x(3)-(x(10)*(S1+S3+S5)+x(9)*(S2+S4))/S;
fff(11)=L/V-(h4/HETP4)*(K1*x(1)-ya1)/(x(3)-x(1));
fff(12)=xa-(xb1-x(19))/xa1;
fff(13)=h1/HETP1-(x(15)-x(17))/(K2*x(18)-x(17));
fff(14)=L*(S1+S3+S5)/V*S-(x(15)-x(17))/(xb1-x(18));
fff(15)=x(13)-(x(18)*(S1+S3+S5)+x(19)*(S2+S4))/S;
fff(16)=h2/HETP2-mi*(x(17)-x(16))/(K2*x(22)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(16));
fff(17)=L/V-mi*(x(17)-x(16))/(mi*x(13)-x(22)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(18)=h3/HETP3-(HETP4*x(16)-h4*(K2*x(12)-yb1)-yb1*HETP4)/(K2*HETP4*x(20)-h4*(K2*x(12)-yb1)-yb1*HETP4);
fff(19)=L*(S2+S4)/V*S-(HETP4*x(16)-h4*(K2*x(12)-yb1)-yb1*HETP4)*mi/(HETP4*(xa1*M1+xb1*M2+xc1*M3+xd1*M4-mi*x(20)));
fff(20)=xb-(x(22)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(21))/(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(21)=x(14)-(x(21)*(S1+S3+S5)+x(20)*(S2+S4))/S;
fff(22)=L/V-(h4/HETP4)*(K2*x(12)-yb1)/(x(14)-x(12));
fff(23)=xa-(x(30)-xc1)/xa1;
fff(24)=h1/HETP1-(x(26)-x(28))/(K3*x(29)-x(28));
fff(25)=L*(S1+S3+S5)/V*S-(x(26)-x(28))/(xc1-x(29));
fff(26)=x(24)-(x(29)*(S1+S3+S5)+x(30)*(S2+S4))/S;
fff(27)=h2/HETP2-mi*(x(28)-x(27))/(K3*x(33)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-mi*x(27));
fff(28)=L/V-mi*(x(28)-x(27))/(mi*x(24)-x(33)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(29)=h3/HETP3-(HETP4*x(27)-h4*(K3*x(23)-yc1)-yc1*HETP4)/(K3*HETP4*x(31)-h4*(K3*x(23)-yc1)-yc1*HETP4);
fff(30)=L*(S2+S4)/V*S-(HETP4*x(27)-h4*(K3*x(23)-yc1)-yc1*HETP4)*mi/(HETP4*(xa1*M1+xb1*M2+xc1*M3+xd1*M4-mi*x(31)));
fff(31)=xb-(mi*x(32)-x(33)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4))/(x(11)*(xa1*M1+xb1*M2+xc1*M3+xd1*M4));
fff(32)=x(25)-(x(32)*(S1+S3+S5)+x(31)*(S2+S4))/S;
fff(33)=L/V-(h4/HETP4)*(K3*x(23)-yc1)/(x(25)-x(23));
fff(34)=(mi-(x(11)-x(22)-x(33))*(xa1*M1+xb1*M2+xc1*M3+xd1*M4))/(xa1*M1+xb1*M2+xc1*M3+xd1*M4)-x(34);matlab求解非线性方程组
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : fff.docx
    • 2015-08-16 11:38:06, 14.64 K

» 收录本帖的淘帖专辑推荐

matlab编程绘图

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-08-16 11:47:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
其它一直常数都一并给出来啊，比如：s1, s2..., M1, M2,..., Lc1, Lc2, U等等。
赞一下 回复此楼
8楼2015-08-16 22:38:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

hsdwlxhx

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
你是这是我遇到的用matlab解的最复杂问题，还是想想怎么简化吧，或者找个新的算法。fsolve这个函数估计是不够用了。

[ 发自小木虫客户端 ]
赞一下 回复此楼
BringLightonChip
2楼2015-08-16 14:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
随机优化中选择合适的，
差分演化，粒子群
再结合fmincon或fsolve里的方法进一步迭代

可能要反复尝试才能找到最合适的参数。
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
3楼2015-08-16 16:10:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜雨tju

金虫 (小有名气)
引用回帖:
3楼: Originally posted by cooooldog at 2015-08-16 16:10:17
随机优化中选择合适的，
差分演化，粒子群
再结合fmincon或fsolve里的方法进一步迭代

可能要反复尝试才能找到最合适的参数。

大神，我是化工的，对于matlab学的不是很精，你能帮我大体写一下程序么，我自己再摸索一下。这个程序已经困扰我好几个月了，谢谢，谢谢
赞一下 回复此楼
4楼2015-08-16 16:45:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308分求调剂（一志愿华东师大）接受跨专业 +8 相信必会光芒万� 2026-04-13 8/400 2026-04-13 20:26 by biomen
[考研] 一志愿沪9，326求生物学调剂 +9 刘墨墨 2026-04-13 9/450 2026-04-13 18:59 by lbsjt
[考研] 本科郑州大学，一志愿华东师范大学282求调剂 +28 熊哥xtk 2026-04-07 33/1650 2026-04-13 17:48 by 邹gv
[考研] 材料复试求调剂 +24 xhhdjdjsjks 2026-04-09 24/1200 2026-04-13 15:49 by 幸免 ..
[考研] 322求调剂 +6 123安康 2026-04-12 13/650 2026-04-12 15:51 by 123安康
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[基金申请] 山东省基金2026 +5 jerry681 2026-04-08 6/300 2026-04-12 08:33 by kudofaye
[考研] 化学工程调剂289 +44 yang婷 2026-04-07 50/2500 2026-04-12 02:36 by 秋豆菜芽
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 求调剂 +6 小聂爱学习 2026-04-11 9/450 2026-04-11 21:20 by 蓝云思雨
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 293调剂 +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by 柴小白
[考研] 282，电气工程专业，求调剂，不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
信息提示
关闭
请填处理意见
关闭