版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

monixia

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2535.6
红花: 3
帖子: 80
在线: 24.2小时
虫号: 1817862
注册: 2012-05-15
性别: GG
专业: 化学反应工程

[求助] 求解决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调到很小了，可是还是没有满意的结果。
求助各位大牛，这个该怎么解？还是方程需要修改？

如果圆满解决，可以加金币~

回复此楼

» 猜你喜欢

070200求调剂，一志愿某211，288分已经有19人回复
法国博士后职位已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有188人回复
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。已经有0人回复
** 已经有1人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有3人回复
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？已经有2人回复
B区学生调剂-兰州交通大学材料科学与工程学院已经有8人回复
山西大同大学物理学专业还有调剂名额，欢迎调剂！已经有9人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab解二元方程组，错误提示，求解决方法已经有11人回复
MATLAB解非线性方程组已经有9人回复
matlab 解非线性偏微分方程组已经有6人回复
MATLAB求解五元非线性方程组，在线等已经有7人回复
非线性二阶微分方程组求解，matlab 已经有12人回复
Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！已经有4人回复
Matlab用牛顿法求解非线性方程组问题已经有8人回复
用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
matlab如何求解一个非线性微分方程组已经有8人回复
sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！已经有12人回复
matlab中利用LM算法解非线性方程组的问题已经有5人回复
MATLAB求解非线性方程组已经有5人回复
关于matlab求解非线性指数方程组出现问题已经有3人回复
【求助】一些A、B类错误求解已经有15人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
求助matlab---fsolve解非线性方程组已经有6人回复
求助！matlab用fsolve函数求解非线性方程组的问题！已经有19人回复
matlab报错，求解已经有5人回复
matlab求解非线性方程组已经有16人回复
matlab问题求助已经有16人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
请教Matlab高手解多元非线性方程组中的参数已经有6人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】求教matlab解非线性方程组已经有9人回复

1楼 2014-11-05 11:36:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +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

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2535.6
红花: 3
帖子: 80
在线: 24.2小时
虫号: 1817862
注册: 2012-05-15
性别: GG
专业: 化学反应工程

引用回帖:

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的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

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的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

比如：
让所有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

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2535.6
红花: 3
帖子: 80
在线: 24.2小时
虫号: 1817862
注册: 2012-05-15
性别: GG
专业: 化学反应工程

引用回帖:

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的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

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 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 一志愿沪9，326生物学求相关专业调剂 +3	刘墨墨 2026-04-09	3/150	2026-04-10 08:37 by 1753564080
[考研] 求调剂 +7	翩翩一书生 2026-04-09	7/350	2026-04-10 07:27 by bljnqdcc
[考研] 一志愿211电子信息347求调剂 +3	554916 2026-04-03	4/200	2026-04-09 22:24 by 554916
[考研] 0703化学调剂325分 +13	15771691647 2026-04-04	15/750	2026-04-09 16:55 by 15771691647
[考研] 一志愿郑州大学 22408 305分求调剂 +3	安小满zzz 2026-04-08	3/150	2026-04-09 16:16 by wp06
[考研] 调剂 +3	青灯不负 2026-04-09	3/150	2026-04-09 15:28 by 释放天性
[考研] 262求调剂 +10	天下第一文 2026-04-04	13/650	2026-04-09 15:16 by 探123
[硕博家园] 新一代电子信息294求调剂不挑学校 +5	Ytyt11 2026-04-09	6/300	2026-04-09 14:40 by Ytyt11
[考研] 086000调剂 +4	十七sa 2026-04-07	4/200	2026-04-09 09:27 by 壹往無前
[考研] 本科211，293分请求调剂 +12	莲菜就是藕吧 2026-04-03	13/650	2026-04-08 20:30 by 背对大海出发
[考研] 263分B区求调剂 +6	李nihao 2026-04-08	6/300	2026-04-08 09:38 by 南开小綦
[考研] 259求调剂 +5	就爱吃土豆呀呀 2026-04-07	5/250	2026-04-07 22:40 by JourneyLucky
[考研] 材料调剂 +11	一样YWY 2026-04-07	11/550	2026-04-07 15:13 by shdgaomin
[考研] 319分085702安全工程求调剂 +6	rious 2026-04-05	6/300	2026-04-07 09:42 by jp9609
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7	晚风不晚 2026-04-04	7/350	2026-04-06 14:06 by houyaoxu
[考研] 一志愿武汉理工大学080200机械工程308分，求调剂 +4	终不似从前 2026-04-05	4/200	2026-04-06 11:46 by 考研学校招点人
[考研] 材料专硕322分 +10	哈哈哈吼吼吼哈 2026-04-04	10/500	2026-04-05 21:22 by 学员8dgXkO
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11	adaie 2026-04-04	12/600	2026-04-05 19:04 by 蓝云思雨
[考研] 358求调剂 +7	秋gk 2026-04-04	7/350	2026-04-05 13:29 by huangmoli
[考研] 0854求调剂 +4	assdll 2026-04-04	4/200	2026-04-05 09:44 by zhq0425

24小时热门版块排行榜

monixia

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

月只蓝

【答案】应助回帖

monixia

月只蓝

【答案】应助回帖

月只蓝

【答案】应助回帖

monixia

月只蓝

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