4楼: Originally posted by sinyfiny at 2015-09-12 08:10:54
能不能帮我也看看，我没得权限，折腾一晚上了...

function as
format long
clear;clc
u0=[0 0 0];
[u fval]=fsolve(@eqs,u0);
[u ;fval]'


function f=eqs(u)
x1=u(1);
x2=u(2);
x3=u(3);

f(1)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((15*12.7-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;
f(2)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((15*12.75-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;
f(3)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((15*12.8-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;

3楼: Originally posted by 月只蓝 at 2015-09-12 07:52:01
代码如下：
function as
format long
clear;clc
u0=;
=fsolve(@eqs,u0);
'


function f=eqs(u)
x1=u(1);
x2=u(2);
x3=u(3);

f(1)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((1 ...