重力加速度计算公式:
g=9.7803(1+0.0053024sin²ψ-0.000005sin²2ψ)
ψ为物体所在处的地理纬度
本以为将上述公式转换为matlab程序语言即可,但是转换完算了几个值发现都不太对,问了同学才知道,matlab里面三角函数公式用的弧度! 直接代入纬度值是不对的,要先转换为弧度才可。
修改后代码如下:
% 当地重力加速度
latitude = latitude/180*pi;
g=9.7803*(1+0.0053024*((sin(latitude))^(2))-0.000005*