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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » MATLAB解非线性方程组
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 2129  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

crazyiyy

木虫 (小有名气)

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

本人有个非线性方程组，已用1stOpt解出，但需要将结果导入到Matlab中，十分麻烦。因此想用Matlab将此方程组解出来，肯定路过大神帮忙！题目如下
1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*(x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*(12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z=0;
         (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)=x;
         (1-(1-z)^8)/(2-(1-z)^9)=y;
注：其中x,y,z都是[0,1]数，一定需要用matlab解出来，感激不尽！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
用心做好每一件事
1楼 2014-07-14 15:43:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
crazyiyy: 金币+100, ★★★★★最佳答案, 衷心地感谢您对我的帮助！ 2014-07-14 16:20:17
CODE:
function jie_fun
u0=[0.5 0.5 0.5];
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
f(1)=1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*...
    (x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*...
    (12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z;
f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;
f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

计算结果：
u =

    0.8138    0.3049    0.0733


目标函数值f_val =

  1.0e-008 *

    0.0008    0.1291   -0.0216
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-07-14 16:03:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

以下程序第对2楼中程序的改进，提高解的小数点后的位数。
CODE:
function jie_fun
clear a;clc
format long
u0=[0.5 0.5 0.5];
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
f(1)=1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*...
    (x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*...
    (12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z;
f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;
f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

计算结果：
u =

   0.813803728151198   0.304867313866725   0.073289552279229


f_val =

  1.0e-008 *

   0.000793234922192   0.129125732328816  -0.021619173118381

已验证，上述结果与1stopt结果一致。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

crazyiyy
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2014-07-14 16:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crazyiyy

木虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by 月只蓝 at 2014-07-14 16:14:10
以下程序第对2楼中程序的改进，提高解的小数点后的位数。


function jie_fun
clear a;clc
format long
u0=;
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
...

您好！ 您的答案的确跟1stOpt一致，还有个问题希望您能为我解惑。如果我需要对方程组中的某些系数进行变化，进而得到一组解，我应该如何操作。
例如 我需要对方程1中的48*2^3 变为48*2^x  其中x=1,2,3,4 这样我应该得到4组解，
赞一下 回复此楼
用心做好每一件事
4楼2014-07-14 16:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
4楼: Originally posted by crazyiyy at 2014-07-14 16:39:24
您好！ 您的答案的确跟1stOpt一致，还有个问题希望您能为我解惑。如果我需要对方程组中的某些系数进行变化，进而得到一组解，我应该如何操作。
例如 我需要对方程1中的48*2^3 变为48*2^x  其中x=1,2,3,4 这样我应 ...

我的处理办法是把x定义为全局变量。
以前已经有过相关代码了，参见：
http://muchong.com/bbs/viewthread.php?tid=6267100&authorid=1122189
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

crazyiyy
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2014-07-14 16:45:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
4楼: Originally posted by crazyiyy at 2014-07-14 16:39:24
您好！ 您的答案的确跟1stOpt一致，还有个问题希望您能为我解惑。如果我需要对方程组中的某些系数进行变化，进而得到一组解，我应该如何操作。
例如 我需要对方程1中的48*2^3 变为48*2^x  其中x=1,2,3,4 这样我应 ...

以及：
http://muchong.com/bbs/viewthread.php?tid=6087640
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2014-07-14 16:47:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crazyiyy

木虫 (小有名气)
送红花一朵
引用回帖:
5楼: Originally posted by 月只蓝 at 2014-07-14 16:45:00
我的处理办法是把x定义为全局变量。
以前已经有过相关代码了，参见：
http://muchong.com/bbs/viewthread.php?tid=6267100&authorid=1122189...

您好！ 我还有最后一个问题，就是matlab种一定要给出初值， 请问初值的调试有没有什么方法或者经验啊，或者说如何判断当前初值得到的解是正解。
赞一下 回复此楼
用心做好每一件事
7楼2014-07-14 17:06:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
7楼: Originally posted by crazyiyy at 2014-07-14 17:06:19
您好！ 我还有最后一个问题，就是matlab种一定要给出初值， 请问初值的调试有没有什么方法或者经验啊，或者说如何判断当前初值得到的解是正解。...

把算出来的解代回原方程，看能有多接近于0，3楼正是给出了这些数值。
此外，对于多解的情况，最终解的如果有物理意义要求，举例说，比如未知数是开尔文温度T，是不可能小于0的，在解里面就应当取满足实际物理意义的。
在直接调用fsolve的情况下，初值基本靠自己尝试，大概判断一下。
究其本质而言，解方程组问题是最优化问题，调用具有全局优化能力的算法，比如说遗传算法，可以降低对初值的要求。
1stopt软件解决最优化问题的强大之处其一在于在操作者层面上不需要赋予初值，即可获得相当优秀的结果。
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
8楼2014-07-14 18:37:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mjjackey

新虫 (初入文坛)
引用回帖:
8楼: Originally posted by 月只蓝 at 2014-07-14 18:37:07
把算出来的解代回原方程，看能有多接近于0，3楼正是给出了这些数值。
此外，对于多解的情况，最终解的如果有物理意义要求，举例说，比如未知数是开尔文温度T，是不可能小于0的，在解里面就应当取满足实际物理意义 ...

用fsolve，貌似也只能设置三种算法，没办法设置更多的算法“调用具有全局优化能力的算法，比如说遗传算法，”
赞一下 回复此楼
9楼2015-03-02 16:40:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
9楼: Originally posted by mjjackey at 2015-03-02 16:40:35
用fsolve，貌似也只能设置三种算法，没办法设置更多的算法“调用具有全局优化能力的算法，比如说遗传算法，”...

“究其本质而言，解方程组问题是最优化问题，调用具有全局优化能力的算法，比如说遗传算法，可以降低对初值的要求。
”
先将原问题转化为最优化问题，调用遗传算法算出初值，该初值代入fsolve函数。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
10楼2015-03-02 17:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 crazyiyy 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学0703-一志愿211-338分求调剂 +7 vants 2026-04-05 7/350 2026-04-05 18:17 by cql1109
[考研] 081700学硕，323分，一志愿中国海洋大学求调剂学校 +16 披星河 2026-04-04 16/800 2026-04-05 11:27 by 猪会飞
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 323求调剂（计算机视觉和大模型项目经历） +3 chaoxiicy 2026-03-31 3/150 2026-04-05 10:33 by zhq0425
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 303求调剂 +9 DLkz1314. 2026-03-30 9/450 2026-04-03 18:34 by ls刘帅
[考研] 一志愿北京交通大学材料工程总分358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by 百灵童888
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 285求调剂 +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 285求调剂 +14 AZMK 2026-04-02 14/700 2026-04-02 15:54 by 上九天揽月（好�
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 08生物与医药专硕初试346找调剂 +6 dianeeee 2026-04-01 7/350 2026-04-02 08:23 by guoweigw
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 349求调剂 +6 吃的不少 2026-04-01 6/300 2026-04-01 17:55 by JYD2011
[考研] 358求调剂 +3 王向阳花 2026-03-31 3/150 2026-04-01 09:56 by zzchen2000
[考研] 调剂申请 +8 张张张张zy 2026-03-31 9/450 2026-04-01 08:29 by zjbkx
[考研] 085601一志愿西北工业大学初试346 +4 085601初试346 2026-03-30 4/200 2026-03-31 07:47 by jp9609
信息提示
关闭
请填处理意见
关闭