temp.m

来自「对一个数组进行操作」· M 代码 · 共 37 行

clc;
clear;

% b = 106;
% c = 185.35;
% 
% beita = pi*( 29 + 39/60 )/180;
% garma = pi*( 120 + 7/60 )/180;
% 
% a1 = sin(garma)/sin(beita)
% a2 = b*cos(garma)/sin(beita)
% a3 = b*sin(garma)*cos(beita)/sin(beita)/sin(beita)
% 
% A = [a1 a2 a3];
% 
% dPar = [0.06*0.06               0                 0;
%             0       (120/206265)*(120/206265)     0;
%             0                   0      (60/206265)*(60/206265)];
% 
% Res = A*dPar*A'

a0 = 30461e-9;  %起始历元对应的a0，单位：s
a1 = 1.0e-10;   %a1，单位：s/s
a2 = 1.1e-18;   %a2，单位：s/s2

RefTime = datenum([2006,5,21,18,47,35]) ;

ObsTime(1,1) = datenum([2006,5,21,23,10,0]) ;
ObsTime(2,1) = datenum([2006,5,21,20,5,0]) ;
ObsTime(3,1) = datenum([2006,5,21,19,15,0]) ;

for iEph=1:3
    dt = (ObsTime(iEph,1) - RefTime) * 86400 ;
    Bias(iEph,1) = a0 + a1*dt + 0.5*a2*dt*dt;
end

Bias = Bias * 299792458 ;