MATLAB计算地磁代码

以下是一个简单的MATLAB代码，用于计算地球磁场在某一点的磁场强度和方向：

% 地球磁场计算代码
% 作者: XXX
% 日期: XXX

% 输入该点的地理坐标（经度、纬度、高度）
longitude = 120; % 经度，单位：度
latitude = 30; % 纬度，单位：度
altitude = 0; % 高度，单位：米

% 计算该点的地磁场（单位：uT）
[Bx, By, Bz, H, F, D, I] = wrldmagm(altitude, latitude, longitude, decyear(now));

% 显示计算结果
fprintf('在经度 %f°，纬度 %f°，高度 %f 米处的地磁场为：\n',longitude,latitude,altitude);
fprintf('水平分量 H = %f uT\n', H);
fprintf('总强度 F = %f uT\n', F);
fprintf('倾斜角 D = %f°\n', D);
fprintf('磁偏角 I = %f°\n', I);
fprintf('三个分量分别为 Bx = %f uT，By = %f uT，Bz = %f uT\n', Bx, By, Bz);

该代码利用了MATLAB的`wrldmagm`函数来获取地球磁场数据，并使用输入的地理坐标计算该点的磁场强度和方向。其中，`altitude`表示海拔高度，`latitude`和`longitude`表示该点的纬度和经度，`decyear(now)`表示当前的日期。`Bx`、`By`和`Bz`表示地磁场在该点的三个方向上的分量，`H`表示水平分量，`F`表示总强度，`D`和`I`分别表示倾斜角和磁偏角。

需要注意的是，该代码仅适用于低纬度地区（赤道附近），在高纬度地区可能存在误差。此外，`wrldmagm`函数需要下载额外的数据文件，首次运行可能需要较长时间。