我是matlab初学者,想请教一个求解三角函数方程组的问题.已知ClatClon循环计算中可以计算出xxyy待求解的变量是afadeta如何在循环体内求解方程组：xx=cos(Clat)*sin(Clon-afa)/(sin(Clat)*sin(deta)+cos

Question

我是matlab初学者,想请教一个求解三角函数方程组的问题.
已知 C_lat C_lon
循环计算中可以计算出 xx yy
待求解的变量是afa deta
如何在循环体内求解方程组：xx=cos(C_lat)*sin(C_lon-afa)/(sin(C_lat)*sin(deta)+cos(C_lat)*cos(deta)*cos(C_lon-afa))
yy=(sin(C_lat)*cos(deta)-cos(C_lat)*sin(deta)*cos(C_lon-afa))/(sin(C_lat)*sin(deta)+cos(C_lat)*cos(deta)*cos(C_lon-afa))

Answer 1

这两个方程有四个未知数,不知道是怎么算XX和YY的,得先把XX和YY求出
for C_lat=0:pi/10:pi
for C_lon=0:pi/10:pi
[afa,deta]=solve('xx=cos(C_lat)*sin(C_lon-afa)/(sin(C_lat)*sin(deta)+cos(C_lat)*cos(deta)*cos(C_lon-afa))',...
'yy=(sin(C_lat)*cos(deta)-cos(C_lat)*sin(deta)*cos(C_lon-afa))/(sin(C_lat)*sin(deta)+cos(C_lat)*cos(deta)*cos(C_lon-afa)) ')
end
end