%%---------極移矩陣計算---%
function?polarmotion=EP（MJD）
%這里的MJD是UTC的簡約儒略日
%t為自J2000.0算起至t的儒略世紀數
%?206264.806247096?=?3600*DEG???用于秒化弧度

%MJD是約簡的儒略世紀數
%MJD=JD-2400000.5
%t=（MJD-51544.0333）/36525.0;
%eop---地球定向文件?在LERS可以查到??注：盡量將文件補充完整
eop=importdata（‘地球定向參數_EOP.txt‘）;
%%
if?MJD>eop（length（eop）1）||MJD%?時間超出表中時間范圍
????xp?=?0.0;
????yp?=?0.0;
else
%?地極坐標?插值得到UTC時的XpYp
????xp?=?interp1（eop（:1）eop（:2）MJD）/3600*pi/180;
????yp?=?interp1（eop（:1）eop（:3）MJD）/3600*pi/180;
end
%%?極移矩陣
polarmotion=rotationy（-xp）*rotationx（-yp）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????650??2015-03-18?20:17??EP.m

?????文件???????1407??2015-03-18?15:46??ER.m

?????文件???????1274??2015-03-18?15:40??NR.m

?????文件????????526??2015-03-17?09:56??PR.m

?????文件??????93131??2015-03-18?09:04??地球定向參數_EOP.txt

?????文件???????2317??2015-03-17?14:00??黃經和傾角章動序列表.txt

-----------?---------??----------?-----??----

????????????????99305????????????????????6