一、程序
/* compute dops ----------------------------------------------------------------
* compute DOP (dilution of precision)
* args : int ns I number of satellites
* double *azel I satellite azimuth/elevation angle (rad) - 方位角
* double elmin I elevation cutoff angle (rad) - 高程截止角(rad)
* double *dop O DOPs {GDOP,PDOP,HDOP,VDOP}
* return : none
* notes : dop[0]-[3] return 0 in case of dop computation error
*-----------------------------------------------------------------------------*/
#define SQRT(x) ((x)<0.0?0.0:sqrt(x))
extern void dops(int ns, const double *azel, double elmin, double *dop)
{
double H[4*MAXSAT],Q[16],