查看单个帖子
旧 2008-12-20, 09:22   #4
snower1895
初级会员
 
注册日期: 2008-12-20
年龄: 40
帖子: 11
声望力: 17
snower1895 正向着好的方向发展
默认 回复: [求助]这个方程怎么解啊?

首先说 我也是新手 我的做法也很业余
在这个方程上,用solve是最不明智的选择,因为solve解决线性方程还可以,解决非线性方程就差远了,这也是新手容易犯的错误 碰到方程就solve,解非线性方程还是fsolve比较好。
下面是我的结果,同样期待高手批阅。

m文件
function y=fun(x)
y=[168/334-acos((x(1)^2+272.5674^2-(x(2)+63.56)^2)/(2*x(1)*272.5674))-acos((x(1)^2+272.5674^2-(x(2)-63.56)^2)/(2*x(1)*272.5674)),0.2461-acos((-x(1)^2+272.5674^2+(x(2)+63.56)^2)/(2*(x(2)+63.56)*272.5674))-acos((-x(1)^2+272.5674^2+(x(2)-63.56)^2)/(2*(x(2)-63.56)*272.5674))]


command 窗口
x0=[1.0 1.0]
fsolve(@fun,x0)

ans =

48.6881 + 8.1128i 0.9401 - 0.7512i





中间过程没有复制粘贴
祝你早日成为高手
snower1895 当前离线   回复时引用此帖