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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求解决matlab解非线性方程组出现的问题
71/1返回列表
查看: 1684  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

monixia

木虫 (小有名气)

[求助] 求解决matlab解非线性方程组出现的问题 已有1人参与

各位matlab大神，我计算的方程组如下，
>> m=0.00040;
f=@(x)([3.14159*m^3+1.9635e-7*x(1)-0.002*3.14159*0.0005^2;
6*(m+x(2))*1.9635e-7-3.14159*m^3+1.9635e-7*x(3)+1.9635e-7*x(6)*(m+x(2))/x(4)-0.002*3.14159*0.0005^2;
3.14159*m^3*x(4)/(m+x(2))+1.9635e-7*x(1)*x(5)/(x(1)+x(3))-1/6*1e-6;
(6*1.9635e-7*(m+x(2))-3.14159*m^3)*x(4)/(m+x(2))+1.9635e-7*x(3)*x(5)/(x(1)+x(3))+1.9635e-7*x(6)-1/6*1e-6;
fix(0.05/(x(1)+x(3)))*(32*0.0010087*x(5)/0.0005^2*x(3)+7.16*0.07275*(3*(1.0087*x(5)/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x(5)/72.75)^(2/3))))-6455.68*x(6)
0.05/(m+x(2))*(32*0.0010087*x(4)/0.0005^2*x(2)+7.16*0.07275*(3*(1.0087*x(4)/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x(4)/72.75)^(2/3))))-6455.68*x(6)]);
options = optimset('MaxIter',1e30,'MaxFunEvals',1e50,'TolFun',1e-50);
x=fsolve(f,[0.005 0.0005 0.001 0.15 0.9 0.3],options)

No solution found.

fsolve stopped because the relative size of the current step is less than the
default value of the step size tolerance squared, but the vector of function values
is not near zero as measured by the selected value of the function tolerance.

<stopping criteria details>


x =

    0.0070   -0.0001    0.0030    0.0344    1.0563    0.4352


fsolve stopped because the relative norm of the current step, 2.848645e-15, is less than
max(options.TolX^2,eps) = 1.000000e-12. However, the sum of squared function values,
r = 3.833839e-15, exceeds sqrt(options.TolFun) = 1.000000e-25.

Optimization Metric                                               Options
relative norm(step) =   2.85e-15                max(TolX^2,eps) =   1e-12 (default)
r =   3.83e-15                                    sqrt(TolFun) =  1.0e-25 (selected)

如上，初值已经取的很接近实验值了，还是没有满意的结果。将matlab算的带入六个方程，
ans =
   4.7168e-12
ans =
  -8.4160e-11
ans =
   1.5716e-09
ans =
  -1.5222e-09
ans =
   15.7955
ans =
  232.7567
后两项收敛不是很好。我已经尝试把optimset调到很小了，可是还是没有满意的结果。
求助各位大牛，这个该怎么解？还是方程需要修改？

如果圆满解决，可以加金币~
回复此楼

» 猜你喜欢

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

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
1stOpt代码（注意不是MATLAB软件代码）:
CODE:
Constant m=0.00040;
Parameters x(1:6);
Function 3.14159*m^3+1.9635e-7*x1-0.002*3.14159*0.0005^2;
6*(m+x2)*1.9635e-7-3.14159*m^3+1.9635e-7*x3+1.9635e-7*x6*(m+x2)/x4-0.002*3.14159*0.0005^2;
3.14159*m^3*x4/(m+x2)+1.9635e-7*x1*x5/(x1+x3)-1/6*1e-6;
(6*1.9635e-7*(m+x2)-3.14159*m^3)*x4/(m+x2)+1.9635e-7*x3*x5/(x1+x3)+1.9635e-7*x6-1/6*1e-6;
(0.05/(x1+x3))*(32*0.0010087*x5/0.0005^2*x3+7.16*0.07275*(3*(1.0087*x5/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x5/72.75)^(2/3))))-6455.68*x6;
0.05/(m+x2)*(32*0.0010087*x4/0.0005^2*x2+7.16*0.07275*(3*(1.0087*x4/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x4/72.75)^(2/3))))-6455.68*x6;

计算结果：
CODE:
函数表达式 1: 3.14159*0.0004^3+1.9635e-7*x1-0.002*3.14159*0.0005^2
         2: 6*(0.0004+x2)*1.9635e-7-3.14159*0.0004^3+1.9635e-7*x3+1.9635e-7*x6*(0.0004+x2)/x4-0.002*3.14159*0.0005^
            2
         3: 3.14159*0.0004^3*x4/(0.0004+x2)+1.9635e-7*x1*x5/(x1+x3)-1/6*1e-6
         4: (6*1.9635e-7*(0.0004+x2)-3.14159*0.0004^3)*x4/(0.0004+x2)+1.9635e-7*x3*x5/(x1+x3)+1.9635e-7*x6-1/6*1e
            -6
         5: (0.05/(x1+x3))*(32*0.0010087*x5/0.0005^2*x3+7.16*0.07275*(3*(1.0087*x5/72.75)^(2/3)/0.0005/(1+3.34*(1.0087
            *x5/72.75)^(2/3))))-6455.68*x6
         6: 0.05/(0.0004+x2)*(32*0.0010087*x4/0.0005^2*x2+7.16*0.07275*(3*(1.0087*x4/72.75)^(2/3)/0.0005/(1+3.34
            *(1.0087*x4/72.75)^(2/3))))-6455.68*x6
目标函数值: 8.44256812628831E-7
x1: 1.38116851215046
x2: 0.213584946434081
x3: 0.106474657467166
x4: 0.0176943028978422
x5: 0.2480594437683
x6: 0.0180988189842738
6个方程目标函数值：
ans =

  2.6982e-007


ans =

  3.1421e-007


ans =

-1.2143e-007


ans =

-1.3880e-007


ans =

  1.4211e-013


ans =

-8.5265e-014

1.5版本计算结果不唯一，以下再例举6组解：
CODE:
目标函数值: 8.81354127610376E-9
x1: 0.0270233957846693
x2: -4.08878886642551E-5
x3: 0.0126231922603729
x4: 0.0323838969507018
x5: 1.10773562695781
x6: 0.383966561657027


目标函数值: 5.88138553710961E-7
x1: 2.08840573452721
x2: -0.0811262956407749
x3: 0.416992754705403
x4: 0.19145905588163
x5: 1.1205806122223
x6: 0.187004870338678


目标函数值: 6.80817986394196E-7
x1: 2.10149925984767
x2: -0.218796314718131
x3: 0.718745532551135
x4: 0.109052962294469
x5: 0.42221573528721
x6: 0.10785410685195

目标函数值: 1.10547660307994E-6
x1: 2.53429273796369
x2: -0.218841716843157
x3: 1.13544861896696
x4: 0.431086688290742
x5: 1.38395310620881
x6: 0.42858810590596

目标函数值: 1.32464424554872E-6
x1: 4.35394653830248
x2: -0.286471171270745
x3: 2.04705166768062
x4: 0.401695771492555
x5: 1.24966197787108
x6: 0.39985126552927

目标函数值: 7.67093610988707E-7
x1: 1.6997802237498
x2: -0.00709436222448859
x3: -0.505554887761384
x4: 0.00623112985186669
x5: 0.00101373791179563
x6: -0.000417366283335469
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-11-05 13:40:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

monixia

木虫 (小有名气)
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-11-05 13:40:03
1stOpt代码（注意不是MATLAB软件代码）:
Constant m=0.00040;
Parameters x(1:6);
Function 3.14159*m^3+1.9635e-7*x1-0.002*3.14159*0.0005^2;
6*(m+x2)*1.9635e-7-3.14159*m^3+1.9635e-7*x3+1.9635e-7*x6*(m+ ...

非常感谢你的解答，1stOpt用起来真方便。不过，我有个问题就是，我希望找到初值附近的零点，但是几乎每次算出来的值都差很多（至少应该没有负数出现）。是不是不太容易实现在初值附近寻找呢？
赞一下 回复此楼
3楼2014-11-05 19:09:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by monixia at 2014-11-05 19:09:09
非常感谢你的解答，1stOpt用起来真方便。不过，我有个问题就是，我希望找到初值附近的零点，但是几乎每次算出来的值都差很多（至少应该没有负数出现）。是不是不太容易实现在初值附近寻找呢？...

可以设定x的取值范围
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2014-11-05 20:52:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by monixia at 2014-11-05 19:09:09
非常感谢你的解答，1stOpt用起来真方便。不过，我有个问题就是，我希望找到初值附近的零点，但是几乎每次算出来的值都差很多（至少应该没有负数出现）。是不是不太容易实现在初值附近寻找呢？...

比如：
让所有x都大于0：
CODE:
Constant m=0.00040;
Parameters x(1:6)[0,];
Function 3.14159*m^3+1.9635e-7*x1-0.002*3.14159*0.0005^2;
6*(m+x2)*1.9635e-7-3.14159*m^3+1.9635e-7*x3+1.9635e-7*x6*(m+x2)/x4-0.002*3.14159*0.0005^2;
3.14159*m^3*x4/(m+x2)+1.9635e-7*x1*x5/(x1+x3)-1/6*1e-6;
(6*1.9635e-7*(m+x2)-3.14159*m^3)*x4/(m+x2)+1.9635e-7*x3*x5/(x1+x3)+1.9635e-7*x6-1/6*1e-6;
(0.05/(x1+x3))*(32*0.0010087*x5/0.0005^2*x3+7.16*0.07275*(3*(1.0087*x5/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x5/72.75)^(2/3))))-6455.68*x6;
0.05/(m+x2)*(32*0.0010087*x4/0.0005^2*x2+7.16*0.07275*(3*(1.0087*x4/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x4/72.75)^(2/3))))-6455.68*x6;

计算结果：
目标函数值: 2.37290707692599E-8
x1: 0.0105848644277172
x2: 0.0161346934843435
x3: 0.000420183576622565
x4: 0.115861612534128
x5: 0.876247555934376
x6: 0.132259599167984
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2014-11-05 20:57:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

monixia

木虫 (小有名气)
引用回帖:
5楼: Originally posted by 月只蓝 at 2014-11-05 20:57:57
比如：
让所有x都大于0：

Constant m=0.00040;
Parameters x(1:6);
Function 3.14159*m^3+1.9635e-7*x1-0.002*3.14159*0.0005^2;
6*(m+x2)*1.9635e-7-3.14159*m^3+1.9635e-7*x3+1.9635e-7*x6*(m+x2)/x4-0.0 ...

十分感谢您的回答！本打算用matlab解，现在就用1stopt吧。我有个小问题问一下，像这种多解的方程，我可以取接近实验值的一组解作为方程的解来验证方程的合理性吗？1stopt得到的解是否可靠呢？
赞一下 回复此楼
6楼2014-11-08 10:42:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
6楼: Originally posted by monixia at 2014-11-08 10:42:38
十分感谢您的回答！本打算用matlab解，现在就用1stopt吧。我有个小问题问一下，像这种多解的方程，我可以取接近实验值的一组解作为方程的解来验证方程的合理性吗？1stopt得到的解是否可靠呢？...

1stopt得到的解是否可靠，你把求出来的解代回原方程就知道能达到多少精度了。
多解取接近实验值有物理意义。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
7楼2014-11-08 11:18:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 monixia 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕297已过四六级求调剂推荐 +5 adaie 2026-03-19 5/250 2026-03-20 10:44 by EBSD
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +4 @taotao 2026-03-20 4/200 2026-03-20 10:37 by EBSD
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考研] 329求调剂 +6 想上学吖吖 2026-03-19 6/300 2026-03-20 09:00 by 每天只摆一小会
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +3 葵梓卫队 2026-03-18 5/250 2026-03-19 19:35 by 给你你注意休息
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
信息提示
关闭
请填处理意见
关闭