9-当地重力加速度

1.已知某点位置经纬度(deg)
g=9.7803(1+0.0053024sin²ψ-0.000005sin²2ψ)   //ψ为物体所在处的地理纬度

2.已知某点位置的经纬度(deg)
% 当地重力加速度 matlab
latitude = latitude/180*pi;
g=9.7803*(1+0.0053024*((sin(latitude))^(2))-0.000005*((sin(2*latitude))^(2)))