j2000 经度_j2000坐标系转换WGS84坐标，看不太懂

functionW=J2000_2_WGS84_matrix(JD_time)%计算J2000坐标系转换WGS84坐标系的旋转矩阵%目前暂未考虑地球极移矩阵%输入参数：%JD_time：儒略日期%输出参数：%W：WGS84坐标系与J2000.0坐...

function W = J2000_2_WGS84_matrix(JD_time)

% 计算J2000坐标系转换WGS84坐标系的旋转矩阵

% 目前暂未考虑地球极移矩阵

% 输入参数：

% JD_time： 儒略日期

% 输出参数：

% W： WGS84坐标系与J2000.0坐标系的转换矩阵，

%% -------------------- 计算极移矩阵 A --------------------------

%-- load: 读取地球的地极坐标(xp,yp)。

[non,Ry,non] = rotate_matrix( -deg2radian(0,0,-0.15) );

[Rx,non,non] = rotate_matrix( -deg2radian(0,0,0.3) );

%-- A = Ry(-xp)*Rx(-yp);

A = Ry*Rx;

%% ------------------- 计算自转矩阵 B----------------------------

T0 = ((JD_time - 2400000.5) - 51544.5)/36525.0; %-- 从J2000起到t时刻的儒略世纪数

Tu = (JD_time - 2451545.0)/36525.0;

%--

e_M = deg2radian(23,26,2