disp_y_mid_graf(i,j)=disp_mid_y;
disp_eps_max_graf(i,j)=disp_eps_max;
disp_eps_mid_graf(i,j)=disp_mid_eps;
mat_y_max_graf(i,j)=mat_y_max;
mat_y_mid_graf(i,j)=mat_mid_y;
mat_eps_max_graf(i,j)=mat_eps_max;
mat_eps_mid_graf(i,j)=mat_mid_eps;
end
end
x =[100 200 400]
x=(meshgrid(x,1:4))'
subplot(4,1,1)
plot(x,disp_y_mid_graf((1:3),:)),YLABEL('Сред.дисп.ошибки по H'), grid on
subplot(4,1,2)
plot(x,disp_y_max_graf((1:3),:)),YLABEL('Мах.дисп.ошибки по H'), grid on
subplot(4,1,3)
plot(x,disp_eps_mid_graf((1:3),:)),YLABEL('Сред.дисп.бок.ошибки'), grid on
subplot(4,1,4)
plot(x,disp_eps_max_graf((1:3),:)),YLABEL('Мах.дисп.бок.ошибки'), grid on
legend('окно 20 м','окно 10 м','окно 7,5 м','окно 5 м')
figure(2)
subplot(4,1,1)
plot(x,mat_y_mid_graf((1:3),:)),YLABEL('Сред.ошибка по H'), grid on
subplot(4,1,2)
plot(x,mat_y_max_graf((1:3),:)),YLABEL('Мах.ошибка по H'), grid on
subplot(4,1,3)
plot(x,mat_eps_mid_graf((1:3),:)),YLABEL('Сред.бок.ошибка'), grid on
subplot(4,1,4)
plot(x,mat_eps_max_graf((1:3),:)),YLABEL('Мах.бок.ошибка'), grid on
legend('окно 20 м','окно 10 м','окно 7,5 м','окно 5 м')
figure(3)
subplot(4,1,1)
plot(x,disp_y_mid_graf((4:6),:)),YLABEL('Сред.дисп.ошибки по H'), grid on
title('прогноз 3 с')
subplot(4,1,2)
plot(x,disp_y_max_graf((4:6),:)),YLABEL('Мах.дисп.ошибки по H'), grid on
subplot(4,1,3)
plot(x,disp_eps_mid_graf((4:6),:)),YLABEL('Сред.дисп.бок.ошибки'), grid on
subplot(4,1,4)
plot(x,disp_eps_max_graf((4:6),:)),YLABEL('Мах.дисп.бок.ошибки'), grid on
legend('окно 20 м','окно 10 м','окно 7,5 м','окно 5 м')
figure(4)
subplot(4,1,1)
plot(x,mat_y_mid_graf((4:6),:)),YLABEL('Сред.ошибка по H'), grid on
title('прогноз 3 с')
subplot(4,1,2)
plot(x,mat_y_max_graf((4:6),:)),YLABEL('Мах.ошибка по H'), grid on
subplot(4,1,3)
plot(x,mat_eps_mid_graf((4:6),:)),YLABEL('Сред.бок.ошибка'), grid on
subplot(4,1,4)
plot(x,mat_eps_max_graf((4:6),:)),YLABEL('Мах.бок.ошибка'), grid on
legend('окно 20 м','окно 10 м','окно 7,5 м','окно 5 м')
33.1 Программа расчета дисперсии ошибки касания полосы
for i=1:10
link='C:\fly\crash\data_';
object_ip=['1234567890'];
load([link object_ip(i)]);
[n,m]=size(koord);
koord_matrix(:,i)=koord(:,m);
end
mat_crash_l_1=abs(sum(koord_matrix(4,:))/10+2012.5);
mat_crash_bok_1=abs(sum(koord_matrix(2,:))/10);
disp_crash_l_1=(1/9).*sum(((sum(koord_matrix(4,:))/10)-koord_matrix(4,:)).^2);
disp_crash_bok_1=(1/9).*sum((mat_crash_bok_1-koord_matrix(2,:)).^2);
for i=1:10
link='C:\fly\crash\data_';
type='_p3';
object_ip=['1234567890'];
load([link object_ip(i) type]);
[n,m]=size(koord);
koord_matrix(:,i)=koord(:,m);
end
mat_crash_l_2=abs(sum(koord_matrix(4,:))/10+2012.5);
mat_crash_bok_2=abs(sum(koord_matrix(2,:))/10);
disp_crash_l_2=(1/9).*sum(((sum(koord_matrix(4,:))/10)-koord_matrix(4,:)).^2);
disp_crash_bok_2=(1/9).*sum((mat_crash_bok_2-koord_matrix(2,:)).^2);
for i=1:10
link='C:\fly\crash\data_';
type='_p3';
type_1='_stop';
object_ip=['1234567890'];
load([link object_ip(i) type_1 type]);
[n,m]=size(koord);
koord_matrix(:,i)=koord(:,m);
end
mat_crash_l_3=abs(sum(koord_matrix(4,:))/10+2012.5);
mat_crash_bok_3=abs(sum(koord_matrix(2,:))/10);
disp_crash_l_3=(1/9).*sum(((sum(koord_matrix(4,:))/10)-koord_matrix(4,:)).^2);
disp_crash_bok_3=(1/9).*sum((mat_crash_bok_3-koord_matrix(2,:)).^2);
subplot(4,1,1),bar([mat_crash_l_1,mat_crash_l_2,mat_crash_l_3]),YLABEL('Мат.ожид.ошибки L')
subplot(4,1,2),bar([mat_crash_bok_1,mat_crash_bok_2,mat_crash_bok_3]),YLABEL('Мат.ожид.бок.ошибки')
subplot(4,1,3),bar([disp_crash_l_1,disp_crash_l_2,disp_crash_l_3]),YLABEL('Дисп.ошибки L')
subplot(4,1,4),bar([disp_crash_bok_1,disp_crash_bok_2,disp_crash_bok_3]),YLABEL('Дисп.бок.ошибки')
legend('1-без прогноза','2-с прогнозом','3-с указанием мнимой точки посадки'),
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.