latlon2utm.m

GPS坐标转换UTM坐标的Matlab程序，转换结果与商业软件一致。转换的基本原理可以参考博客https://blog.csdn.net/Action_MYX/article/details/83421594 同样的其他语言的转换也可以自己写。

%% 输入经纬度 纬度 unit:degree

lat= 40.0691643333333;    
lon= 116.242161333333;

%%  %%%%%unit: km
a=6378.137; %e=0.0818192;
e=sqrt（0.00669438）;
k0=0.9996; E0=500; N0=0;
Zonenum=fix（lon/6）+31;
lamda0=（Zonenum-1）*6-180+3;   %%%%degree 
lamda0=lamda0*pi/180;    %%%%%radian
phi=lat*pi/180;   lamda=lon*pi/180;    %%%%%radian

v=1./sqrt（1-e^2*（sin（phi））.^2）;
A=（lamda-lamda0）.*cos（phi