求！Matlab的二元二次方程组解法

wwwaiting

我遇到一个难题，如何用Matlab程序求解下面的方程组，谢谢各位同学的指教！
" v1 r7 q" r9 S方程组：（x-x1)^2+(y-y1)^2=a
             (x-x2)^2+(y-y2)^2=c
             y=kx+b
& k- y% y! N9 \3 j& m. P, \' C其中x1,x2,y1,y2,a,b,c,k都是可输入的常数，求解x,y的值，请问如何编写Matlab程序？

karl_wang

这个是非线性方程，可以采用X=fsolve('fun',X0)求解。
以上是建议，供参考！

baiyh2008

Matlab解非线性方程组有多种解法,常用的有符号解法和数值解法，符号解法可以得到解析解，而数值方法需要给出最初的猜测值，然后得到数值解。你的方程可以直接用符号解法求得：
$ ?3 t! ~2 H: T, e% T4 I$ K( U
syms x y x1 x2 y1 y2 a b c
# V+ P" C, b. a% t$ }1 `
: L; o7 d' @- @2 ff = [(x-x1)^2+(y-y1)^2-a;
3 }( C2 D. L- o2 }7 Z    (x-x2)^2+(y-y2)^2-c]
[X,Y]=solve(f(1),f(2))

mnpfc

进来学习下

wwwaiting

2# karl_wang
2 `" u; u3 u4 E) e. p4 q) @谢谢哦~~ 是不是用solve也可以解决呢？如果在二次方程的等式右边再加个随机数，那个随机数怎么实现呢?

wwwaiting

3# baiyh2008
( f2 i, e6 l& D5 B% Y谢谢~~好像可以这个样子，我还有个问题同上面的问题。。。

karl_wang

2# karl_wang  
谢谢哦~~ 是不是用solve也可以解决呢？如果在二次方程的等式右边再加个随机数，那个随机数怎么实现呢?
wwwaiting 发表于 2009-3-26 15:42

# D- |6 n0 E* a0 ^% H& [你可以给一个数赋一个随机值，然后加在等式的右边，这样应该可以的。

baiyh2008

sorry! there is an error in orignal code, you have to define independent variables (x and y) in 'solve' command,

[X,Y]=solve(f(1), f(2), 'x','y')
2 O/ j3 B6 X0 a( {5 E0 L* cXx = simplify(X)
& P3 |: ]8 Y. v# [% wYy = simplify(Y)
& \( V, ]7 T) W$ `! ~: s2 x" i
) G9 G" T1 R! B' K& w& Z! oIn addition, describe more details for your second question. It will be good for solving.

jingyiloveyou

为啥solve的结果中含有复数，并且结果明显不对？？