EDA365电子论坛网

标题: 求解复杂方程组 [打印本页]

作者: 大小的小    时间: 2020-12-7 15:18
标题: 求解复杂方程组
>> syms x y z a b c d h
% E0 ^! S  ]( {4 v* Z* x% N  x$ seq1= 1/3*1/[((1/6)[(x - y)^2 + (y - z)^2 + (x - z)^2])^1.5 + c (((2 x - y - z)*(2 y - x - z)*(2 z - x - y))/27)]^(2/3)*[1.5*Sqrt[(1/6)[(x - y)^2 + (y - z)^2 + (x - z)^2]]*(2 x - y - z)/3 +c*(4 y*z - 2 x*z - 2 x*y + 2*x^2 - y^2 - z^2)/9]-a/h;# T9 ~% e" z+ K% Y
eq2= 1/3*1/[((1/6)[(x - y)^2 + (y - z)^2 + (x - z)^2])^1.5 +c (((2 x - y - z)*(2 y - x - z)*(2 z - x - y))/27)]^(2/3)*[1.5*Sqrt[(1/6)[(x - y)^2 + (y - z)^2 + (x - z)^2]]*(2 y - x - z)/3 +c*(4 x*z - 2 y*z - 2 x*y + 2*y^2 - x^2 - z^2)/9]-b/h;
9 i- m2 a8 c* Q' Oeq3= 1/3*1/[((1/6)[(x - y)^2 + (y - z)^2 + (x - z)^2])^1.5 + c (((2 x - y - z)*(2 y - x - z)*(2 z - x - y))/27)]^(2/3)*[1.5*Sqrt[(1/6)[(x - y)^2 + (y - z)^2 + (x - z)^2]]*(2 z - y - x)/3 + c*(4 y*x - 2 x*z - 2 z*y + 2*z^2 - y^2 - x^2)/9]-d/h;7 J$ _# D4 |, B3 \4 r
[x,y,z]=solve(eq1,eq2,eq3,a,b,c,d,h)
/ y' {: L4 U1 F6 X各位大神 是因为我输入的方程括号太多吗,这个解不出来,用牛顿迭代能解吗?
, q9 i/ B" X" n5 n8 U
作者: Uifhjvv    时间: 2020-12-7 16:27
你这式子连着犯至少四个语法错误
作者: IBB-EUT    时间: 2020-12-7 18:52

作者: kekek    时间: 2020-12-8 10:24





欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2