用solve解含小量的非线性方程组出现不同结果
本帖最后由 guohf 于 2014-4-21 17:53 编辑要解方程组
(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 mupadmex
Error in MuPAD command: Not a square matrix. [(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: Explicit solution could not be found.
> In solve at 179
x =[ empty sym ]
y = []
是我写错了，把第二个方程中的5e-5写成5e-6了，真马虎。
其实我要做一个循环，把这次得到的x, y值作为新的ax, ay值，再进行求解，第一种方法较实用一点了
页:
[1]