用MATLAB编写ICRF 惯性系到地球固连系 ECF 的坐标转换矩阵，采用IAU2006 标准。

在MATLAB中编写ICRF（国际天文学联合会定义的地球指向的惯性参考框架）到ECF（地球固定坐标系）的坐标转换矩阵，可以遵循IAU（国际天文学联合会）2006年的定义。IAU2006定义的转换包括几个步骤，主要包含从ICRF到地心地固坐标系（ECEF，Earth-Centered Earth-Fixed）的转换，以及从ECEF到ECF（通常是根据UT1和极移来调整）的转换。

以下是一个简化的示例MATLAB代码，它展示了如何根据IAU2006定义进行坐标转换：

function R_eci_ecf = icrf_to_ecf(jdUt1)
    % 计算从ICRF到ECEF的转换矩阵，jdUt1为UT1儒略日时间
    % 此函数使用IAU-2006/2000预cession-nutation模型和地球自转模型
    
    % 计算儒略世纪数
    T = (jdUt1 - 2451545.0) / 36525.0;
    
    % 计算恒星时间
    GMST = 24110.54841 + (8640184.812866 + (0.093104 - 6.2E-6 * T) * T) * T / 3600;
    % 将恒星时间转换为恒星角度
    GMST = mod(GMST, 86400);
    % 为了与IAU定义一致，需要对恒星时间进行调整
    GAST = mod(GMST + (0.00264 * sin(2*pi*(jdUt1 - 2451545.0)/36525) + 0.000063 * sin(2*pi*jdUt1/36525)), 86400);
    
    % 极移参数
    x_pole = 0; % 以弧秒为单位的极点在X轴上的移动
    y_pole = 0; % 以弧秒为单位的极点在Y轴上的移动
    % 将极移转换为弧度
    x_pole = x_pole * pi / (180 * 3600);
    y_pole = y_pole * pi / (180 * 3600);
    
    % 转换矩阵构建
    % 首先计算从ICRF到ECEF的旋转矩阵
    R = [ cosd(GAST) sind(GAST) 0; -sind(GAST) cosd(GAST) 0; 0 0 1];
    
    % 然后考虑极移
    R = R * [1 0 0; 0 cos(y_pole) sin(x_pole); 0 -sin(y_pole) cos(x_pole)];
    
    % 最后转换到ECF
    R_eci_ecf = R;
end

需要注意的是，该代码只是一个示例，实际应用中可能需要更精确的极移参数和恒星时间计算，以及考虑岁差和章动的影响。因此，在实际使用中，你应该引入最新的天文常数和精确的岁差章动模型。通常，天文学软件包会提供这些功能。

请记住，使用这些算法时，需要确保输入的儒略时间是正确的，并且你有访问准确的天文参数和模型的途径。