program_5.m

用MATLAB编写关于飞机设计的matlab程序

%符号说明
%b1,b2,b3,b4为横航向小扰动运动方程的特征方程的系数
%m：飞机质量（kg）；V：飞行速度（m/s）；a：迎角(deg)；
%CY_b：侧力对侧滑角的导数
%Cl_b：滚转力矩系数对侧滑角的导数；Cl_R：滚转力矩系数对偏航角速度的导数
%Cl_P：滚转力矩系数对滚转角速度的导数；Cn_b：偏航力矩系数对侧滑角的导数
%Cn_P：偏航力矩系数对滚转角速度的导数；Cn_R：偏航力矩系数对偏航角速度的导数
%CY_b1、Cl_b1、Cl_P1、Cl_R1、Cn_b1、Cn_P1、Cn_R1只是中间符号，其含义如程序中表达式所示
%Ix：飞机相对于x轴的惯性矩；Iz：飞机相对于z轴的惯性矩；Ixz：飞机相对于xz轴的惯性积
%机体坐标系：x轴顺机身向后为正，y指向右翼为正，z轴与xy轴成右手系，原点为全机质心

%参照各参数的含义，输入对应值运行即可

m=2;V=15;Ix=;Iz=;Ixz=;g=9.81;a=3.39;
CY_b=-0.10948;Cl_b=0.034154;Cl_P=-0.54149;Cl_R=0.044063;
Cn_b=-0.037591;Cn_P=-0.050911;Cn_R=-0.034918;
CY_b1=CY_b/(m*V);
Cl_P=Cl_P/Ix;
Cl_R=Cl_R/Ix;
Cl_b=Cl_b/Ix;
Cn_R=Cn_R/Iz;
Cn_P=Cn_P/Iz;
Cn_b=Cn_b/Iz;
I=1-Ixz^2/(Ix*Iz);
Cl_P1=(Cl_P+Cn_P*Ixz/Ix)/I;
Cl_R1=(Cl_R+Cn_R*Ixz/Ix)/I;
Cl_b1=(Cl_b+Cn_b*Ixz/Ix)/I;
Cn_P1=(Cn_P+Cl_P*Ixz/Iz)/I;
Cn_R1=(Cn_R+Cl_R*Ixz/Iz)/I;
Cn_b1=(Cn_b+Cl_b*Ixz/Iz)/I;
b1=-(CY_b1+Cl_P1+Cn_R1)
b2=CY_b1*(Cn_R1+Cl_P1)+(Cn_R1*Cl_P1-Cn_P1*Cl_R1)-Cn_b1-(a*pi/180)*Cl_b1
b3=Cn_b1*Cl_P1-Cl_b1*(g/V+Cn_P1)+CY_b1*(Cn_P1*Cl_R1-Cl_P1*Cn_R1)+(a*pi/180)*(Cl_b1*Cn_R1-Cn_b1*Cl_R1)
b4=(g/V)*(Cl_b1*Cn_R1-Cn_b1*Cl_R1)
R=b1*b2*b3-b1^2*b4-b3^2
if b4>0&&R>0&&Cn_b<0&&Cl_b>0
    disp('横航向静安定且横航向动安定');
else if b4<=0&&R>0&&Cn_b<0&&Cl_b>0
            disp('横航向静安定且振荡安定（航向动安定），但螺旋不安定（横向动不安定）');
else if b4>0&&R<=0&&Cn_b<0&&Cl_b>0
            disp('横航向静安定且螺旋安定（横向动安定），但振荡不安定（航向动不安定）');
else if b4>0&&R>0&&Cn_b>=0&&Cl_b>0
            disp('横航向动安定且横向静安定，但航向静不安定');
else if b4>0&&R<0&&Cn_b<0&&Cl_b<=0
            disp('横航向动安定且航向静安定，但横向静不安定');
else if b4<=0&&R<=0&&Cn_b<0&&Cl_b>0
            disp('横航向静安定，但横航向动不安定');
else if b4<=0&&R>0&&Cn_b>=0&&Cl_b>0
            disp('横向静安定，但航向静不安定；振荡安定（航向动安定），但螺旋不安定（横向动不安定）');
else if b4<=0&&R>0&&Cn_b<0&&Cl_b<=0
            disp('航向静安定，但横向静不安定；振荡安定（航向动安定），但螺旋不安定（横向动不安定）');
else if b4>0&&R<=0&&Cn_b>=0&&Cl_b>0
            disp('横向静安定，但航向静不安定；螺旋安定（横向动安定），但振荡不安定（航向动不安定）');
else if b4>0&&R<=0&&Cn_b<0&&Cl_b<=0
            disp('航向静安定，但横向静不安定；螺旋安定（横向动安定），但振荡不安定（航向动不安定）');
else if b4>0&&R>0&&Cn_b>=0&&Cl_b<=0
            disp('横航向动安定，但横航向静不安定');
else if b4<=0&&R<=0&&Cn_b>=0&&Cl_b>0
            disp('横向静安定，但航向静不安定，且横航向动不安定');
else if b4<=0&&R<=0&&Cn_b<0&&Cl_b<=0
            disp('航向静安定，但横向静不安定；且横航向动不安定');
else if b4<=0&&R>0&&Cn_b>=0&&Cl_b<=0
            disp('横航向静不安定；振荡安定（航向动安定），但螺旋不安定（横向动不安定）');
else if b4>0&&R<=0&&Cn_b>=0&&Cl_b<=0
            disp('横航向静不安定；螺旋安定（横向动安定），但振荡不安定（航向动不安定）');
else if b4<=0&&R<=0&&Cn_b>=0&&Cl_b<=0
            disp('横航向静不安定且横航向动不安定');
        end
    end
end
end
end
end
end
end
end
end
end
end
end
end
end
end