matlab方程组求解

Question

想求出x,y的值，可是程序不能运行，大神给瞅瞅
a=(10+43/60+13/3600)*pi/180;
syms x y
[x,y]=solve('sin(x)*sin(a)+cos(x)*cos(a)*cos(pi*6)=sin(atan(y/sqrt(1.1383^2+0.5142^2))','sin(x)*sin(a)+cos(x)*cos(a)*cos(pi*5/24)=sin(atan(y/sqrt(1.516^2+0.5715^2))');

Answer 1

有时候，用solve（）函数不一定能得到解析值或数值解。特别是含有较复杂的三角函数方程往往是得不到你想要的值。所以求解这类方程可以考虑数值方法。一般可以用fsolve（）函数。

求解步骤：

1、用ezplot(eq1,eq2)绘出两函数方程的图形，来判断两函数曲线的交点的位置。

从图形来看，很容易判断x、y初值为[-12 0]或[-8 0]或[4 0]。。。。。。时，有交点。

2、将两函数方程自定义为一个变量{x}的函数，即f（x1,x2）

fun=@(x)[sin(x(1)).*sin(a)+cos(x(1)).*cos(a)*cos(pi*6)-sin(atan(x(2)/sqrt(1.1383^2+0.5142^2))),sin(x(1))*sin(a)+cos(x(1))*cos(a)*cos(pi*5/24)-sin(atan(x(2)/sqrt(1.516^2+0.5715^2)))]

3、用fsolve（）函数求解。

[x,f,h]=fsolve(fun,[-12 0])

[x,f,h]=fsolve(fun,[-8 0])

[x,f,h]=fsolve(fun,[4 0])

。。。。。。

x=x(1)，y=x(2)

f1≈0，f2≈0说明所求的x、y是方程的解。