matlab 地惯坐标系转地固坐标系

根据提供的引用内容，可以得知Matlab可以完成大地坐标系、地心地固坐标系和站心坐标系之间的坐标变换。其中，地心地固坐标系到大地坐标系的坐标变换和地心地固坐标系到站心坐标系的坐标变换都可以实现地惯坐标系到地固坐标系的转换。下面是一个Matlab实现地惯坐标系到地固坐标系转换的例子：

% 地球椭球体参数
a = 6378137.0; % 长半轴
f = 1/298.257223563; % 扁率
b = a*(1-f); % 短半轴
e = sqrt(a^2-b^2)/a; % 第一偏心率
e2 = sqrt(a^2-b^2)/b; % 第二偏心率

% 地惯坐标系下的位置和速度
r = [1000; 2000; 3000]; % 位置矢量，单位：m
v = [10; 20; 30]; % 速度矢量，单位：m/s

% 地球自转角速度
w = 7.2921151467e-5; % 单位：rad/s

% 计算地球自转角
theta = w*100; % 假设时间为100s

% 地惯坐标系到地固坐标系的转换矩阵
C = [cos(theta), sin(theta), 0; -sin(theta), cos(theta), 0; 0, 0, 1];

% 地球上一点的经纬度和高程
lat = 39.9; % 纬度，单位：度
lon = 116.3; % 经度，单位：度
h = 100; % 高程，单位：m

% 计算地心地固坐标系到大地坐标系的转换矩阵
N = a/sqrt(1-e^2*sin(lat)^2);
X = (N+h)*cos(lat)*cos(lon);
Y = (N+h)*cos(lat)*sin(lon);
Z = (N*(1-e^2)+h)*sin(lat);
P = [X; Y; Z];
K = P/sqrt(P'*P);
I = [1 0 0; 0 1 0; 0 0 1];
J = cross(K, C*K);
J = J/sqrt(J'*J);
I = cross(J, K);
I = I/sqrt(I'*I);
T = [I'; J'; K'];

% 地惯坐标系到地固坐标系的转换
r_ECEF = T*r;
v_ECEF = T*v;