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