回复: [求助]这个方程怎么解啊?
首先说 我也是新手 我的做法也很业余
在这个方程上,用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
中间过程没有复制粘贴
祝你早日成为高手
|