用solve解含小量的非线性方程组出现不同结果
本帖最后由 guohf 于 2014-4-21 18:08 编辑要解方程组 (x-6e-5)^2+(y-5e-5)^2=(8e-6)^2
x^2+(y-5e-5)^2=3.6e-9
方程有实数解，可采用两种方法，得不到想要的结果，且会报错，请高手给看下，先谢了法一：>> ax=6e-5;ay=5e-5;>> syms x y;>>=solve((ax-x)^2+(ay-y)^2-(8e-6)^2==0,x^2+(y-5e-6)^2-3.6e-9==0) x, y得到很怪的复数解。 将ax, ay的具体值代入，法二：=solve((x-6e-5)^2+(y-5e-5)^2-(8e-6)^2==0,x^2+(y-5e-6)^2-3.6e-9==0) Error using mupadmexError in MuPAD command: Not a squarematrix. [(Dom::Matrix(Dom::ExpressionField()))::_power] Error in sym/privBinaryOp (line 1693) Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:}); Error in sym/mpower (line 189) B = privBinaryOp(A, p, 'symobj::mpower');出现上面错误，不太理解这些“错误”指的是什么？在网上查了一下，添加了点：=solve((x-6e-5).^2+(y-5e-5).^2-(8e-6)^2==0,x.^2+(y-5e-6).^2-3.6e-9==0) Warning: Explicitsolution could not be found. > In solve at 179 x = [ empty sym ]y = [] 是我写错了，把第二个方程中的5e-5写成5e-6了，真马虎。
其实我要做一个循环，把这次得到的x, y值作为新的ax, ay值，再进行求解，第一种方法较实用一点了
页:
[1]