| 查看: 857 | 回复: 4 | ||||
[求助]
复杂方程组求解 已有1人参与
|
|
syms x1 x2 x3 x4 x5 x6 x7 x8 eqn1 = (x1-x2)*(1-x1-x3+x2+x4+x6+x7)-1.63734*10^(-13)*(1-2*x1-x3)^2; eqn2 = (x2-x3)*(x2+x4+x7)-5752818*(x1-x2)*(1-x1-x3+x2+x4+x6+x7); eqn3 = (x3-x4)*(1-x1-x3+x2+x4+x6+x7)-2.17528*10^(-12)*(x2-x3)*(1-2*x1-x3); eqn4 = (x2+x4+x7)*(x4-x5)-285869.4*(x3-x4)*(1-x1-x3+x2+x4+x6+x7); eqn5 = x5-x6-x7-8.26*10^(-5)*(x4-x5); eqn6 = (x6-2*x8)*x6-1.10948*10^15*(x5-x6-x7)*(1-x1-x3+x2+x4+x6+x7); eqn7 = x7*(x2+x4+x7)-4254857432*(x5-x6-x7)*(1-x1-x3+x2+x4+x6+x7); eqn8 = x8^2-0.000750773*(x6-2*x8)^2; eqn1=matlabFunction(eqn1) eqn2=matlabFunction(eqn2) eqn3=matlabFunction(eqn3) eqn4=matlabFunction(eqn4) eqn5=matlabFunction(eqn5) eqn6=matlabFunction(eqn6) eqn7=matlabFunction(eqn7) eqn8=matlabFunction(eqn8) f=@(x)[eqn1(x(1),x(2),x(3),x(4),x(6),x(7)); eqn2(x(1),x(2),x(3),x(4),x(6),x(7)); eqn3(x(1),x(2),x(3),x(4),x(6),x(7)); eqn4(x(1),x(2),x(3),x(4),x(5),x(6),x(7)); eqn5(x(4),x(5),x(6),x(7)); eqn6(x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8)); eqn7(x(1),x(2),x(3),x(4),x(5),x(6),x(7)); eqn8(x(6),x(8))] options = optimoptions('fsolve','algorithm','levenberg-marquardt') [x,exitflag]=fsolve(f,[1 1 1 1 1 1 1 1],options) 这是程序，但是猜测的初始值不同解出来的结果就不同。我并不知道初始值的大概是多少。 有没有不用猜测初始值的方法。或者说怎么加一些不等式限制结果？ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有92人回复
-
时间戳变了，能看出什么问题？
已经有5人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助，给予经济报酬。关于微分方程组求解。
已经有3人回复
- 史上最复杂的多元微分方程组的差分求解
已经有9人回复
- 求助Matlab大神帮忙求解复杂的参数方程组
已经有3人回复
- 线性方程组matlab求解，遇到难题。。。。
已经有3人回复
- MATLAB中采用同fsolve求解复杂非线性方程组的问题
已经有1人回复
- 较复杂非线性方程组的求解
已经有13人回复
- 牛顿迭代法接非线性方程组求解
已经有1人回复
- 求助解复杂非线性方程组的好的方法
已经有24人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 求助matlab解一个复杂非线性方程组
已经有2人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 拟牛顿法求解复杂（带根号的）非线性方程组
已经有17人回复
- 请大家帮我看看这样复杂的方程组用什么方法求解比较好？
已经有3人回复
- 求解待定系数线性方程组
已经有6人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- fslove命令如何求解这个复杂的四元一次方程组
已经有0人回复
- fslove命令如何求解这个复杂的四元一次方程组
已经有0人回复
- 【求助】非线性方程组求解
已经有3人回复
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
2楼2017-07-23 19:30:15
3楼2017-07-23 19:36:21
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
【答案】应助回帖
感谢参与，应助指数 +1
|
转化为极值问题之后，这类问题常有多个局部极值点，所以，迭代法收敛到局部极值点或者鞍点很正常 {x1=0.6355106711326112, x2=-0.1030334431563012, x3=0.35402353260730507, x4=0.3456964184083561, x5=-0.2531286421850123, x6=-0.46081051722303384, x7=0.2076817467914361, x8=-0.42714029617040405} {x1=0.9028736908415601, x2=0.48571680999216704, x3=-0.43455809512663685, x4=-0.9741912923385155, x5=-0.043209954815936744, x6=-0.03908048783795981, x7=-0.004129453990300291, x8=0.43950983982809383} |
4楼2017-07-23 19:38:26
5楼2017-07-24 08:34:28
