下载后可任意编辑测量平差程序设计1.角度（度分秒）到弧度 AngleToRadian#define PI 3.14159265double AngleToRadian（double angle）{int D，M；double S,radian,degree, angle,MS；D=int（angle+0.3）；MS=angle-D；M=int（（MS）*100+0.3）；S=（MS*100-M）*100;degree=D+M/60.0+S/3600.0;radian=degree*PI/180.0;return radian;}注意:防止数据溢出,要加个微小量,例如 0.3.2.弧度换角度(度分秒) RadianToAngle#define PI 3.14159265double RadianToAngle(double radian){int D,M;double S,radian,degree,MS,angle;degree=radian*180/PI;D=int(degree);MS=degree-D;M=int(MS*60);S=(MS*60-M)*60;angle=D+M/100.0+S/10000.0;return angle;}3.已知两点求坐标方位角 Azimuth#include double Azimuth(double xi,double yi,double xj,double yj){double Dx,Dy,S,T;Dx=xj-xi;Dy=yj-yi;S=sqrt(Dx*Dx+Dy*Dy);if(S<1e-10) return 0;T=asin(Dy/S);if(Dx<0) T=PI-T;if(Dx>0&&(Dy<0)||T<0) T=2*PI+T;1下载后可任意编辑return T;}4.开辟二维数组的动态空间的宏#include #define NewArray2D(type,A,i,n,m){A=(type**)malloc(n*sizeof(type*));\ for(i=0;ivoid countermatrix(double **T, double **s, double **r, double **Q,double **N, double **rt,int n){for(i=0;i=0;i++){r[i][i]=1/T[i][i];for(j=i...
