cx_fct.m

导弹控制系统的鲁棒控制设计实例

function c_x = cx_fct(M)
% Interpolation of c_x for a given value of Mach number
%
% Values of the Mach number
%          1      2      3      4      5      6      7      8      9      10       
M_1    = [0.00   0.10   0.20   0.30   0.40   0.50   0.60   0.62   0.64   0.66];
M_2    = [0.68   0.70   0.72   0.74   0.76   0.78   0.80   0.82   0.84   0.86];
M_3    = [0.88   0.90   0.92   0.94   0.96   0.98   1.00   1.02   1.04   1.06];
M_4    = [1.08   1.10   1.12   1.14   1.16   1.18   1.20   1.22   1.24   1.26];
M_5    = [1.28   1.30   1.40   1.50   1.60   1.70   1.80   1.90   2.00   2.10];
M_6    = [2.20   2.30   2.40   2.50   2.60   2.70   2.80   2.90   3.00   3.10];
M_7    = [3.20   3.30   3.40   3.50   3.60   3.70   3.80   3.90   4.00   4.10];
M_8    = [4.20   4.30   4.40   4.50   4.60   4.70   4.80   4.90   5.00];
M_val  = [M_1    M_2    M_3    M_4    M_5    M_6    M_7    M_8];
%
% Values of c_x
%          1      2      3      4      5      6      7      8      9      10
cx_1   = [0.1020 0.1010 0.1000 0.0990 0.0980 0.0970 0.0960 0.0958 0.0956 0.0954];
cx_2   = [0.0952 0.0950 0.0940 0.0946 0.0944 0.0948 0.0954 0.0960 0.0968 0.0976];
cx_3   = [0.0986 0.1000 0.1018 0.1038 0.1060 0.1130 0.1450 0.1780 0.2000 0.2320];
cx_4   = [0.2660 0.2730 0.2760 0.2770 0.2764 0.2760 0.2744 0.2720 0.2660 0.2620];
cx_5   = [0.2570 0.2530 0.2370 0.2250 0.2150 0.2070 0.2010 0.1950 0.1900 0.1845];
cx_6   = [0.1800 0.1750 0.1710 0.1665 0.1620 0.1585 0.1545 0.1510 0.1495 0.1465];
cx_7   = [0.1435 0.1420 0.1400 0.1385 0.1360 0.1340 0.1324 0.1310 0.1294 0.1278];
cx_8   = [0.1262 0.1250 0.1240 0.1230 0.1220 0.1210 0.1200 0.1190 0.1180];
cx_val = [cx_1 cx_2 cx_3 cx_4 cx_5 cx_6 cx_7 cx_8];
%
if M <= 5
   c_x = interp1(M_val,cx_val,M);
else
   c_x = 0.11;
end  
%
c_form = 2.0;           %        form coefficient
%
c_x = c_form*c_x;