查看: 21331|回复: 8

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

字体大小: 正常放大

4 主题	4 听众	15 积分

升级 10.53%

该用户从未签到

电梯直达

1^#

发表于 2009-3-25 14:41 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

我遇到一个难题，如何用Matlab程序求解下面的方程组，谢谢各位同学的指教！
方程组：（x-x1)^2+(y-y1)^2=a
(x-x2)^2+(y-y2)^2=c
y=kx+b
其中x1,x2,y1,y2,a,b,c,k都是可输入的常数，求解x,y的值，请问如何编写Matlab程序？

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

karl_wang

10 主题	4 听众	422 积分

升级 40.67%

TA的每日心情

	开心 2014-11-27 19:08

签到天数: 21 天

[LV.4]偶尔看看III

群组: 第六届国赛赛前冲刺培

群组: 国赛讨论

群组: 学术交流A

群组: 学术交流B

群组: 2014数学建模国赛备战

2^#

发表于 2009-3-25 16:38 |只看该作者 |招呼Ta 关注Ta

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

我也说一句

发表

使用道具举报

baiyh2008

0 主题	0 听众	2 积分

升级 40%

该用户从未签到

3^#

发表于 2009-3-26 01:27 |只看该作者 |招呼Ta 关注Ta

Matlab解非线性方程组有多种解法,常用的有符号解法和数值解法，符号解法可以得到解析解，而数值方法需要给出最初的猜测值，然后得到数值解。你的方程可以直接用符号解法求得：

syms x y x1 x2 y1 y2 a b c

f = [(x-x1)^2+(y-y1)^2-a;
(x-x2)^2+(y-y2)^2-c]
[X,Y]=solve(f(1),f(2))

我也说一句

发表

使用道具举报

mnpfc

131 主题	38 听众	1万积分

升级 0%

TA的每日心情

	开心 2018-12-4 08:49

签到天数: 282 天

[LV.8]以坛为家I

群组: 2010MCM

群组: 数学建模

群组: 中国矿业大学数学建模协会

群组: 华中师大数模协会

群组: Mathematica研究小组

4^#

发表于 2009-3-26 08:00 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

进来学习下

我也说一句

发表

使用道具举报

wwwaiting

4 主题	4 听众	15 积分

升级 10.53%

该用户从未签到

5^#

发表于 2009-3-26 15:42 |只看该作者 |招呼Ta 关注Ta

2# karl_wang
谢谢哦~~ 是不是用solve也可以解决呢？如果在二次方程的等式右边再加个随机数，那个随机数怎么实现呢?

我也说一句

发表

使用道具举报

wwwaiting

4 主题	4 听众	15 积分

升级 10.53%

该用户从未签到

6^#

发表于 2009-3-26 15:46 |只看该作者 |招呼Ta 关注Ta

3# baiyh2008
谢谢~~好像可以这个样子，我还有个问题同上面的问题。。。

我也说一句

发表

使用道具举报

karl_wang

10 主题	4 听众	422 积分

升级 40.67%

TA的每日心情

	开心 2014-11-27 19:08

签到天数: 21 天

[LV.4]偶尔看看III

群组: 第六届国赛赛前冲刺培

群组: 国赛讨论

群组: 学术交流A

群组: 学术交流B

群组: 2014数学建模国赛备战

7^#

发表于 2009-3-26 16:54 |只看该作者 |招呼Ta 关注Ta

2# karl_wang ! n2 `8 p: j/ m: p E; h
谢谢哦~~ 是不是用solve也可以解决呢？如果在二次方程的等式右边再加个随机数，那个随机数怎么实现呢?' \; P6 W/ v! F/ d% s
wwwaiting 发表于 2009-3-26 15:42

你可以给一个数赋一个随机值，然后加在等式的右边，这样应该可以的。

我也说一句

发表

使用道具举报

baiyh2008

0 主题	0 听众	2 积分

升级 40%

该用户从未签到

8^#

发表于 2009-3-26 22:00 |只看该作者 |招呼Ta 关注Ta

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')
Xx = simplify(X)
Yy = simplify(Y)

In addition, describe more details for your second question. It will be good for solving.