高斯正反算及换带计算matlab源代码 附截图.pdf

高斯投影坐标正、反算及相邻带的坐标换算MATLAB源代码 主程序截图 正算截图 换带计算截图 反算截图 以上是课本算例截图 精度0.001m 和0.001” 主程序 程序引用代码 MAIN 程序名称 主程序 高斯投影坐标正、反算及相邻带的坐标换算MATLAB源代码 disp(欢迎使用高斯投影正反算及相邻带的坐标换算程序); disp(1 ：高斯正算 2 ：高斯反算 3 ：换带计算); K=0; while (K1||K3) K=input(请根据上列选择计算类型K=); switch K case 1 GSZS; case 2 GSFS; case 3 HDJS; otherwise disp(K 值无效 （1-3 ）); end disp(程序作者：桂林理工大学 刘敬涛); disp(指导老师：王新桥 刘斌 卢献健); disp( 引用请注明出处); end 子程序1 程序引用代码 GSZS 程序名称 高斯正算 function GSZS %GSZS 是将大地坐标换算为高斯坐标的子函数 %此函数要调用DHH 和HHD 两个子函数 %此函数包含子午线收敛角的计算 disp(你选择的是高斯正算); B=input(输入大地坐标B=); L=input(输入大地坐标L=); L0=input(输入所用中央子午线L0=); B=DHH(B); L=DHH(L); L0=DHH(L0); disp(1:克拉索夫斯基椭球 2:1975 年国际椭球 3:WGS-84 椭球); T=0; while (T1||T2) T=input(请根据上列选择椭球模型T=); switch T case 1 a=6378245.0000000000; b=6356863.0187730473; X=111134.861*(B*180/pi)-16036.480*sin(2*B)+16.828*sin(4*B)-0.022*sin(6*B); case 2 高斯投影坐标正、反算及相邻带的坐标换算MATLAB源代码 a=6378140.0000000000; b=6356755.2881575287; X=111133.005*(B*180/pi)-16038.528*sin(2*B)+16.833*sin(4*B)-0.022*sin(6*B); otherwise disp(T 值无效 （1-2 ）); end end e=(sqrt(a^2-b^2))/a; e1=(sqrt(a^2-b^2))/b; V=sqrt(1+(e1^2)*(cos(B))^2); c=(a^2)/b; M=c/(V^3); N=c/V; t=tan(B); n=sqrt((e1^2)*(cos(B))^2); l=L-L0; xp1=X; xp2=(N*sin(B)*cos(B)*l^2)/2; xp3=(N*sin(B)*((cos(B))^3)*(5-t^2+9*n^2+4*n^4)*l^4)/24; xp4=(N*sin(B)*((cos(B))^5)*(61-58*t^2+t^4)*l^6)/720; x=xp1+xp2+xp3+xp4; yp1=N*cos(B)*l; yp2=N*(cos(B))^3*(1-t^2+n^2)*l^3/6; yp3=N*(cos(B))^5*(5-18*t^2+t^4+14*n^2-58*(n^2)*(t^2))