地心地固直角坐标系转换站心直角坐标系matlab代码

地心地固直角坐标系是以地球质心为原点的直角坐标系，而站心直角坐标系是以观测站点为原点的直角坐标系。要将地心地固直角坐标系转换为站心直角坐标系，可以使用matlab代码实现。下面是一个简单的示例matlab代码：

function [x_s, y_s, z_s] = ECEFtoENU(x_e, y_e, z_e, lat, lon)
    a = 6378137.0;  % 地球长半轴
    b = 6356752.3;  % 地球短半轴
    e = sqrt(1 - (b/a)^2);  % 地球偏心率
    N = a / sqrt(1 - e^2 * sind(lat)^2);  % 卯酉圈半径
    x_p = (N + z_e) * cosd(lat) * cosd(lon);  % 转换矩阵
    y_p = (N + z_e) * cosd(lat) * sind(lon);
    z_p = (N * (1 - e^2) + z_e) * sind(lat);
    x_s = -sind(lon) * x_p - sind(lat) * cosd(lon) * y_p;  % 转换为站心直角坐标系
    y_s = cosd(lon) * x_p - sind(lat) * sind(lon) * y_p;
    z_s = cosd(lat) * y_p - sind(lat) * z_p;
end

这段代码实现了地心地固直角坐标系到站心直角坐标系的转换。其中，x_e, y_e, z_e是地心地固直角坐标系下的坐标，lat和lon分别是站点的纬度和经度，x_s, y_s, z_s即为转换后的站心直角坐标系下的坐标。可以使用这个函数来实现坐标系的转换。