cy_fct.m

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

function c_y_alpha = cy_fct(M)
% Interpolation of c_y_alpha 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_y_alpha
%           1    2     3     4     5     6     7     8     9     10
cy_1   = [0.35  0.35  0.35  0.35  0.35  0.35  0.35  0.35  0.35  0.35 ];
cy_2   = [0.35  0.35  0.35  0.35  0.35  0.35  0.35  0.35  0.35  0.35 ];
cy_3   = [0.35  0.35  0.35  0.35  0.355 0.36  0.375 0.41  0.43  0.45 ];
cy_4   = [0.46  0.48  0.49  0.50  0.51  0.52  0.53  0.54  0.545 0.55 ];
cy_5   = [0.555 0.56  0.57  0.568 0.565 0.561 0.555 0.55  0.54  0.525];
cy_6   = [0.51  0.50  0.49  0.48  0.47  0.46  0.455 0.45  0.45  0.45 ];
cy_7   = [0.45  0.45  0.45  0.45  0.45  0.45  0.45  0.45  0.45  0.45 ];
cy_8   = [0.45  0.45  0.45  0.45  0.45  0.45  0.45  0.45  0.45];
cy_val = [cy_1  cy_2  cy_3  cy_4  cy_5  cy_6  cy_7  cy_8];
%
if M <= 5
   c_y_alpha = interp1(M_val,cy_val,M);
else
   c_y_alpha = 0.45;
end