J2000toWGS84

function W = J2000_2_WGS84_matrix（JD_time）
% 计算J2000坐标系转换WGS84坐标系的旋转矩阵
% 目前暂未考虑地球极移矩阵
% 输入参数：
%       JD_time：    儒略日期
% 输出参数：
%       W：          WGS84坐标系与J2000.0坐标系的转换矩阵，
%% -------------------- 计算极移矩阵 A --------------------------
%-- load: 读取地球的地极坐标（xpyp）。
[nonRynon] = rotate_matrix（ -deg2radian（00-0.15） ）;
[Rxnonnon] = rotate_matrix（ -deg2radian（000.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（232621.448） - deg2radian（0046.8150）*T0 - deg2radian（000.00059）*T0^2 + deg2radian（000.001813）*T0^3;   %弧度角
DE_Blocks = DE405_load（JD_time）;           %读取DE405数据包
DE_angle =Planet_Position_J2000_Earth（JD_timeDE_Blocks.Blocks14）;   %获取对应时间的黄经章动和交角章动
delta_fai = DE_angle（1）;
delta_e   = DE_angle（2）;     
%--
% seta_G = （ 100.46061837 + 36000.770053608*Tu + 0.000387933*Tu^2 - Tu^3/38710000 ）/180*pi ...
%           + delta_fai *cos（e_M + delta_e） + 0.178401471734655;   %--  计算格林尼治恒星时[弧度];
seta_G =  （ deg2radia