在simulink中使用函数编写一段由滑跑、加速拉起、爬高、改平、平直飞行、倾斜预转弯、转弯、改平、平直飞行、加速拉起、加速爬高、改平、平直飞行、减速飞行、倾斜转弯、转弯、改平、平直飞行、低头、俯冲、改平、倾斜预转弯、转弯、改平、平直飞行、加速飞行、低头、俯冲等各种飞行机动动作构成的飞行轨迹的代码示例

以下是一个简单的示例，展示如何在Simulink中使用函数编写一个由滑跑、加速拉起、爬高、改平、平直飞行、倾斜预转弯、转弯、改平、平直飞行、加速拉起、加速爬高、改平、平直飞行、减速飞行、倾斜转弯、转弯、改平、平直飞行、低头、俯冲、改平、倾斜预转弯、转弯、改平、平直飞行、加速飞行、低头、俯冲等各种飞行机动动作构成的飞行轨迹的代码：

function out = flight_trajectory(in)

% 该函数定义了一段由滑跑、加速拉起、爬高、改平、平直飞行、倾斜预转弯、转弯、改平、平直飞行、加速拉起、加速爬高、改平、平直飞行、减速飞行、倾斜转弯、转弯、改平、平直飞行、低头、俯冲、改平、倾斜预转弯、转弯、改平、平直飞行、加速飞行、低头、俯冲等各种飞行机动动作构成的飞行轨迹
% 输入参数in是一个n*1的矩阵，其中n是轨迹中的总时间步数，in(i)表示第i个时间步长
% 输出参数out是一个n*6的矩阵，其中n是轨迹中的总时间步数，out(i,:)表示第i个时间步长的机体姿态(quaternion)和位置(x,y,z)

% 初始化变量
t = linspace(0, sum(in), length(in));
q = [1; 0; 0; 0]; % 初始机体姿态为四元数[1 0 0 0]
pos = [0; 0; 0]; % 初始位置为[0 0 0]

% 定义转弯半径和速度等参数
R = 1000; % 转弯半径
V = 100; % 飞行速度

% 循环计算每个时间步长的机体姿态和位置
for i = 1:length(in)
    dt = t(i+1) - t(i); % 计算时间步长
    if i == 1 % 第一个时间步长为滑跑
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 2 % 第二个时间步长为加速拉起
        pitch = 10; % 俯仰角为10度
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 3 % 第三个时间步长为爬高
        pitch = 20; % 俯仰角为20度
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 4 % 第四个时间步长为改平
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 5 % 第五个时间步长为平直飞行
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 6 % 第六个时间步长为倾斜预转弯
        pitch = 0; % 俯仰角为0
        roll = 5; % 横滚角为5度
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 7 % 第七个时间步长为转弯
        pitch = 0; % 俯仰角为0
        roll = 5; % 横滚角为5度
        yaw = atan2d(in(i)*sind(roll), R+in(i)*cosd(roll)); % 计算偏航角
        q = quatmultiply(q, angle2quat(deg2rad(yaw), [0 0 1])); % 更新机体姿态
        pos = pos + [R*(sind(yaw)-sind(yaw-dt*V/R)); -R*(cosd(yaw)-cosd(yaw-dt*V/R)); 0]; % 更新位置
    elseif i == 8 % 第八个时间步长为改平
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 9 % 第九个时间步长为平直飞行
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 10 % 第十个时间步长为加速拉起
        pitch = 10; % 俯仰角为10度
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 11 % 第十一个时间步长为加速爬高
        pitch = 20; % 俯仰角为20度
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 12 % 第十二个时间步长为改平
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == 13 % 第十三个时间步长为平直飞行
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif i == length(in) % 最后一个时间步长为减速飞行
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif mod(i-14, 4) == 0 % 倾斜转弯时间步长
        pitch = 0; % 俯仰角为0
        roll = 60; % 横滚角为60度
        yaw = 0; % 偏航角为0
        q = quatmultiply(q, angle2quat(deg2rad(yaw), [0 0 1])); % 更新机体姿态
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif mod(i-14, 4) == 1 % 转弯时间步长
        pitch = 0; % 俯仰角为0
        roll = 60; % 横滚角为60度
        yaw = atan2d(in(i)*sind(roll), R+in(i)*cosd(roll)); % 计算偏航角
        q = quatmultiply(q, angle2quat(deg2rad(yaw), [0 0 1])); % 更新机体姿态
        pos = pos + [R*(sind(yaw)-sind(yaw-dt*V/R)); -R*(cosd(yaw)-cosd(yaw-dt*V/R)); 0]; % 更新位置
    elseif mod(i-14, 4) == 2 % 改平时间步长
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    elseif mod(i-14, 4) == 3 % 平直飞行时间步长
        pitch = 0; % 俯仰角为0
        roll = 0; % 横滚角为0
        yaw = 0; % 偏航角为0
        pos = pos + [in(i)*cosd(pitch)*cosd(yaw); in(i)*cosd(pitch)*sind(yaw); -in(i)*sind(pitch)]; % 更新位置
    end
    out(i,:) = quaternion2euler(q, 'XYZ')'; % 将四元数机体姿态转换为欧拉角机体姿态
    out(i,4:6) = pos'; % 更新位置
end

以上代码只是一个简单的示例，实际飞行轨迹的编写需要根据具体需求进行调整和修改。同时，还需要在Simulink中定义输入和输出端口，并将该函数作为MATLAB函数模块嵌入到Simulink模型中。