J2000.0坐标系转换wgs84坐标系的Matlab代码

J2000.0坐标系转换为WGS84坐标系的Matlab代码如下：

function [lat, lon, alt] = j2000_to_wgs84(x, y, z)
% J2000.0坐标系转换为WGS84坐标系
% x、y、z为J2000.0坐标系下的三维坐标
% lat、lon、alt为WGS84坐标系下的纬度、经度、高度

% J2000.0坐标系到WGS84坐标系的转换矩阵
T = [-0.000000095269014, 0.000000046898438, -0.000000008311758;...
     -0.000000046898438, -0.000000095269014, 0.000000010064163;...
     0.000000008311758, -0.000000010064163, -0.000000055350735];

% J2000.0坐标系到WGS84坐标系的平移向量
delta = [6378137.0; 6378137.0; 6356752.3142];

% 转换为WGS84坐标系下的三维坐标
xyz = T * [x; y; z] + delta;

% 计算经度、纬度和高度
lat = atan2(xyz(3), sqrt(xyz(1)^2 + xyz(2)^2));
lon = atan2(xyz(2), xyz(1));
N = calc_n(xyz(1), xyz(2));
alt = sqrt(xyz(1)^2 + xyz(2)^2)/cos(lat) - N;

% 将弧度转换为角度
lat = lat * 180/pi;
lon = lon * 180/pi;
end

function N = calc_n(x, y)
% 计算椭球体法线线长
a = 6378137.0; % WGS84椭球体长半轴
b = 6356752.3142; % WGS84椭球体短半轴
e = sqrt(1 - (b/a)^2); % WGS84椭球体第一偏心率
N = a^2/sqrt(a^2*(cos(y))^2+b^2*(sin(y))^2);
end

其中，`j2000_to_wgs84`函数实现了J2000.0坐标系到WGS84坐标系的转换，`calc_n`函数用于计算WGS84椭球体法线线长。在使用代码时，需要将J2000.0坐标系下的三维坐标作为输入参数传入，函数会返回WGS84坐标系下的纬度、经度和高度。