飞机系统matlab建模

%飞机系统测试模块可设置参数如下： %*************************Aircraft model*********************************** M = 12000; % 飞机质量 V0 = 72; % 飞机着陆初速度 Dis_b = 6.727; % 前轮与重心之间的水平距离 Dis_a = 1.076; % 主轮与重心之间的水平距离 H_wei = 2.178; % 飞机重心到地面的高度 To = 4260; % 剩余推力 Cx = 0.1027; % 飞机滑跑阻力系数 Cy = 0.6 ; % 飞机滑跑的升力系数 Ss = 20 ; % 阻力伞的迎风面积 Hs = 0.67; % 阻力伞作用点离飞机重心的距离 Ht = 0.1; % 飞机发动机离重心的距离 J_air = 4000; % 机体的转动惯量 Kv = 1; % 发动机推力系数 Air_den = 0.125; % 空气密度 S_wing = 50.88; % 机翼面积 Cs = 0.75; % 阻力伞的阻力系数 %*************************landing gear parameter*************************** K1=42529; % 主起落架刚度系数 K2=2500; % 前起落架刚度系数 （参数值不可超过5000） C2=800; % 前起落架阻尼系数 C1=800; % 主起落架阻尼系数 K0=536000; % 动态刚度系数 Wn=60; % 固有频率 Ksn=0.2; % 阻尼比 %内部固定参数参数（包括测试产生信号模块参数） %***************************Servo_value_m********************************** kx=0.36; % 伺服阀的相对阻尼系数 Wgn=17.74; % 伺服阀的固有频率 C0 = 4; % 电流门限 C_fixv=0; % 电流索定值 K_c=0.25; % 电流比例系数 C_cons=10; % 电流常数 P_rat=10^6; % 压力比例系数 P_max = 9*10^6; % 最大刹车压力 P_up=10^7; % 防滑压力最大极限 P_low=0; % 防滑压力最小极限 P_ini = 0; % 初始刹车压力 P_loss = 10^6; % 损失刹车压力 %***************************Servo_value_m********************************** %********************************Break_System***************************** Msm=5000; % 最大刹车力矩 Pm=10*10^6; % 最大刹车压力 P0=1*10^6; % 刹车压力损失 Px=8.3*10^6; % 大迟滞刹车压力 u_mc=0.23; % 摩擦系数 N_mc=4; % 摩擦面面数 R_effe=0.142; % 有效摩擦半径 %*******************************Break_System******************************* %**********************************Aircraft_wheel************************** V0 = 72; % 飞机着陆时初速度 V_aup = 100; % 飞机速度上限 V_alow = 0; % 飞机速度下限 V_wcon = 0; % 机轮速度最小值 J_turn = 5*0.377; % 机轮转动惯量 %**********************************Aircraft_wheel************************** %*************************landing gear parameter*************************** Fro_lim_up = 0.5; % 前起落架抖动饱和输出上界 Fro_lim_low = -0.5; % 前起落架抖动饱和输出下界 Ang_rate_lim_up = 1; % 俯仰角变化率饱和输出上界 Ang_rate_lim_low = -1; % 俯仰角变化率饱和输出下界 Main_lim_up = 0.5; % 前起落架抖动饱和输出上界 Main_lim_low = -0.5; % 前起落架抖动饱和输出下界 %*************************landing gear parameter*************************** %***************************Combined with torque*************************** Num_tire=4; % 主起落架轮胎个数 %****************************Combined with torque************************** %*************************Aircraft model*********************************** g = 9.8; % 重力加速度 %um = 0.3 ; % 主轮与地面的结合系数 %uf = 0.33; % 前轮与地面的结合系数 %*************************Aircraft model************************************ %***************************Tire parameter******************************** Ro_ini = 0.4; % 轮胎初始半径 Pres_fac = 1.07*10^-5; % 轮胎压缩系数 Num_tire = 2; % 轮胎个数 Min_radius = 0.345; % 机轮最小压缩半径 %***************************Tire parameter******************************** %**************************Runway parameter****************************** %Select_runway = [1,2,3]; % 跑道选择 Select_runway = 1; % 跑道选择 D_dry = 0.65; % 干跑道倔强因子 C_dry = 1.5344; % 干跑道波形因子 B_dry = 14.5; % 干跑道峰值因子 D_wet = 0.35; % 湿跑道倔强因子 C_wet = 2.0; % 湿跑道波形因子 B_wet = 8.2; % 湿跑道峰值因子 D_snow = 0.15; % 积雪跑道倔强因子） C_snow = 2.0875; % 积雪道波形因子 B_snow = 10; % 积雪跑道峰值因子 %**************************Runway parameter****************************** %***************************Controller parameter*************************** Scale_fac = 1; % 速度频率和电压比例系数 Dec_rate = 2.6; % 减速率 Delta_v = 0; % 速度差门限 Amp_fac = 5; % 速度差放大系数 Kp = 1; % 瞬时级放大系数 Delta_vp = 1; % 瞬时级门限 Ins_const = 0; % 瞬时级常数 KI1 = 3.5; % 升压速度 KI2 = 1.8; % 放电加速度 Delta_vt1 = 0.4; % PBM级门限（下） Delta_vt2 = 1.7; % PBM级门限（上） Vpbm_up = 50; % PBM饱和输出值(上) Vpbm_low = 0; % PBM饱和输出值(下) Delta_vd = 1; % 微分级门限 Diff_const = 0; % 微分级常数 Kps = 0.3; % 瞬时级平衡系数 Kis = 0.4; % PBM级平衡系数 Kds = 0.3; % 微分级平衡系数 Kc = 1; % 电流放大倍数 Vsm_up = inf; % 综合级饱和输出值(上) Vsm_low = -inf; % 综合级饱和输出值（下） %***************************Controller parameter***************************