Главная страница
Навигация по странице:

  • Методы экспериментальных исследований АКС

  • Оценка эффективности испытаний на основе марковских моделей. ЛР. Отчет по лабораторной работе по курсу " Методы экспериментальных исследований акс " на тему Оценка эффективности испытаний на основе марковских моделей


    Скачать 169.58 Kb.
    НазваниеОтчет по лабораторной работе по курсу " Методы экспериментальных исследований акс " на тему Оценка эффективности испытаний на основе марковских моделей
    АнкорОценка эффективности испытаний на основе марковских моделей
    Дата22.02.2022
    Размер169.58 Kb.
    Формат файлаdocx
    Имя файлаЛР.docx
    ТипОтчет
    #370495

    МИНИСТЕРСТВО НАУКИ И ВЫСШЕГО ОБРАЗОВАНИЯ РОССИЙСКОЙ

    ФЕДЕРАЦИИ

    ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ

    УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ

    «МОСКОВСКИЙ АВИАЦИОННЫЙ ИНСТИТУТ

    (национальный исследовательский университет)»
    Институт №6 «Аэрокосмический»

    Кафедра 604 Системный анализ и управление

    ОТЧЕТ

    ПО ЛАБОРАТОРНОЙ РАБОТЕ

    по курсу "Методы экспериментальных исследований АКС"

    на тему: «Оценка эффективности испытаний на основе марковских моделей»

    Вариант 1
    Работу выполнили

    студенты гр. М60-401Б:

    Антонова А. А.

    Прокопьев М.Н.

    Работу проверил:

    Пельтихин А.В.

    Содержание


    1Постановка задачи 3

    2Результаты испытаний 5

    3Вывод 8


    1. Постановка задачи


    Целью данной лабораторной работы является получение оценок эффективности проведения испытаний сложных технических систем на основе ЛА с использованием марковской модели.

    Преимущество марковской модели заключается в том, что при выводе математической модели испытаний вводятся допущения:

    1. Предполагается, что все дефекты испытываемой системы одинаковы с точки зрения возможности их обнаружения;

    2. Предполагается, что дефекты являются независимыми.

    В лабораторной работе рассматривается 3 варианта соединения элементов системы: последовательное (рис. 1), параллельное (рис. 2) и смешанное (рис 3.) – для которых необходимо определить такие характеристики как вероятность не обнаружения дефекта и интенсивность обнаружения дефектов в системе



    Рис. 1.1 Последовательное соединение



    Рис. 1.2 Параллельное соединение



    Рис. 1.3 Смешанное соединение

    Для последовательного и параллельного соединений вероятности рассчитываются по соотношениям (1.1), (1.2), соответственно:



    (1.1)



    (1.2)

    Интенсивности рассчитываются по формулам



    (1.3)



    (1.4)

    где определяется из соотношения (1.2),

    L – число соединенных компонент.

    Задача заключается в построении графических зависимостей за время равное 500 с. По полученным результатам расчета вероятностей не обнаружения дефектов необходимо определить моменты времени , в которые функции попадают в интервал (0.05;0) для рассматриваемых вариантов построения систем.

    Исходные данные

    Числовые характеристики элементов при разных начальных вероятностях и интенсивностях компонент

    Таблица 2.1

    № элемента

    1

    2

    3

    4



    0.9

    0.8

    0.75

    0.7



    0.01

    0.015

    0.02

    0.025

    Числовые характеристики элементов при одинаковых начальных вероятностях и интенсивностях для всех компонент


    1. Результаты испытаний




    Рис. 3.4 Вероятность не обнаружения дефектов в системе с неоднородными элементами



    Рис. 3.5 Интенсивность обнаружения дефектов в системе с неоднородными элементами



    Рис. 3.6 Вероятность не обнаружения дефектов в системе с однородными элементами



    Рис. 3.7 Интенсивность обнаружения дефектов в системе с однородными элементами

    Программа аналитических расчетов на ПЭВМ представлена в приложении А.

    Результаты оценки продолжительности результатов представлены в таблице 3.1

    Таблица 3.2

    Оценка продолжительности результатов

    Однородные элементы

    Неоднородные элементы

    Последовательное

    соединение

    Параллельное

    соединение

    Смешанное соединение

    Последовательное

    соединение

    Параллельное

    соединение

    Смешанное соединение



    30

    236

    54

    29

    310

    46


    1. Вывод


    В результате проведенных испытаний было выявлено, что при любом виде системы (рис. 1.1-1.3) с течением времени вероятность обнаружения дефектов увеличивается, причиной чего может быть износ системы.

    Из полученных графиков (рис. 3.1, 3.3) выясняется, что вероятнее всего не обнаружить дефект в случае, если будем использовать систему со смешанным соединением как однородных элементов, так и неоднородных.

    С точки зрения трудности обнаружения дефектов в системе, то в системе с последовательным соединением элементов, в отличие от системы с параллельным соединением, обнаружить дефект легче. Это утверждение следует из графиков интенсивностей обнаружения дефекта (рис. 3.2, 3.4): интенсивность обнаружения дефекта при последовательном соединении выше, чем при параллельном.

    В таблице 3.1 приведены результаты расчёта момента времени, когда с вероятностью 0.95 обнаружены все дефекты. Из таблицы видно, что продолжительность испытаний выше при испытании системы с параллельным соединением в случае однородного и неоднородного соединения.

    Приложение А

    clc

    clear

    % Исходные данные

    P0=[0.9 0.8 0.75 0.7];% начальные вероятности необн. дефектов по элементам

    q=[0.01 0.015 0.02 0.025];% интенсивности обн. дефектов по элементам

    tk=500;% продолжительность моделирования

    Pm=mean(P0)% среднее значение P0

    qm=mean(q)% среднее значение q

    qsum=sum(q)% сумма интенсивностей

    Pmult=prod(P0)% произведение элементов вектора P0

    Tpoc_odn=0; Tpoc_raz=0; Tmix_odn=0; Tmix_raz=0; Tpar_odn=0; Tpar_raz=0;%начальное значение

    % оценки продолжительности испытаний

    for t=1:1:tk% цикл по времени моделирования

    Pcpocl(t)=Pmult*exp(-t*qsum);% Вер. необн. в системе с последоват. соед.

    %разнородных элементов

    if ((Pcpocl(t)<=0.05)&(Tpoc_raz==0))

    Tpoc_raz=t% оценка продолжительности испытаний

    end

    qpocl(t)=qsum;% интенсивность обн. в системе с последоват. соед.

    %разнородных элементов

    Pmpocl(t)=Pm^4*exp(-t*4*qm);% Вер. необн. в системе с последоват. соед.

    %однородных элементов

    if ((Pmpocl(t)<=0.05)&(Tpoc_odn==0))

    Tpoc_odn=t% оценка продолжительности испытаний

    end

    qmpocl(t)=4*qm;% интенсивность обн. в системе с последоват. соед.

    %однородных элементов

    for i=1:1:4

    Pobn(i)= 1-P0(i)*exp(-q(i)*t);

    end

    Pcpar(t)=1-prod(Pobn);% Вер. необн. в системе с параллельн. соед.

    %разнородных элементов

    if ((Pcpar(t)<=0.05)&(Tpar_raz==0))

    Tpar_raz=t% оценка продолжительности испытаний

    end

    Pmpar(t)=1-(1-Pm*exp(-qm*t))^4;% Вер. необн. в системе с параллельн. соед.

    %однородных элементов

    if ((Pmpar(t)<=0.05)&(Tpar_odn==0))

    Tpar_odn=t% оценка продолжительности испытаний

    end

    Pmmix(t)=(Pm^2*exp(-t*2*qm))*(1-(1-Pm*exp(-qm*t))^2);% Вер. необн.

    %в системе со смешанным соед.однородных элементов

    if ((Pmmix(t)<=0.05)&(Tmix_odn==0))

    Tmix_odn=t% оценка продолжительности испытаний

    end

    qmpar(t)=(Pm*4*qm*(1-Pm*exp(-qm*t))^3*exp(-qm*t))/(1-(1-Pm*exp(-qm*t))^4);% интенсивность обн. в системе с параллельн. соед.

    %однородных элементов

    qmmix(t)=2*qm+(Pm*2*qm*(1-Pm*exp(-qm*t))*exp(-qm*t))/(1-(1-Pm*exp(-qm*t))^2);% интенсивность обн. в системе со смешанным. соед.

    %однородных элементов

    Pcmix(t)=P0(3)*P0(4)*exp(-t*(q(3)+q(4)))*(1-Pobn(1)*Pobn(2));% Вер. необн.

    %в системе со смеш. соед. элементов (1 и 2 - паралл,3 и 4 - посл.)

    if ((Pcmix(t)<=0.05)&(Tmix_raz==0))

    Tmix_raz=t% оценка продолжительности испытаний

    end

    qpar(t)=-(P0(1)*q(1)*exp(-q(1)*t)*(P0(2)*exp(-q(2)*t) - 1)*(P0(3)*...

    exp(-q(3)*t) - 1)*(P0(4)*exp(-q(4)*t) - 1) + P0(2)*q(2)*exp(-q(2)*t)*...

    (P0(1)*exp(-q(1)*t) - 1)*(P0(3)*exp(-q(3)*t) - 1)*(P0(4)*exp(-q(4)*t)...

    - 1) + P0(3)*q(3)*exp(-q(3)*t)*(P0(1)*exp(-q(1)*t) - 1)*(P0(2)*exp(-q(2)*t)...

    - 1)*(P0(4)*exp(-q(4)*t) - 1) + P0(4)*q(4)*exp(-q(4)*t)*(P0(1)*exp(-q(1)*t)...

    - 1)*(P0(2)*exp(-q(2)*t) - 1)*(P0(3)*exp(-q(3)*t) - 1))/Pcpar(t);% интенсивность

    %обн. в системе с параллельн. соед. разнородных элементов

    qmix(t)=q(3)+q(4)-(P0(1)*q(1)*exp(-q(1)*t)*(P0(2)*exp(-q(2)*t) - 1) +...

    P0(2)*q(2)*exp(-q(2)*t)*(P0(1)*exp(-q(1)*t) - 1))/(1-Pobn(1)*Pobn(2));% интенсивность

    %обн. в системе со смешанным соед. элементов (1 и 2 - паралл,3 и 4 - посл.)

    end

    t=1:1:tk;

    figure(1)

    plot(t,Pcpocl, t,Pcpar, t,Pcmix),grid

    %title('Вероятность необнаружения дефектов в системе с неоднородными элементами');

    xlabel('время t'); ylabel('Вероятность P(t)');

    legend('Ppos(t)','Pmix(t)','Ppar(t)');

    figure(2)

    plot(t,qpocl, t,qpar, t,qmix),grid

    ylim([0 0.072]);

    %title('Интенсивность обнаружения дефектов в системе с неоднородными элементами');

    xlabel('время t'); ylabel('Интенсивность q(t)');

    legend('qpos(t)','qmix(t)','qpar(t)');

    figure(3)

    plot(t,Pmpocl, t,Pmpar, t,Pmmix),grid

    %title('Вероятность необнаружения дефектов в системе с однородными элементами');

    xlabel('время t'); ylabel('Вероятность Р(t)');

    legend('Ppos(t)','Pmix(t)','Ppar(t)');

    figure(4)

    plot(t,qmpocl, t,qmpar, t,qmmix),grid

    ylim([0 0.072]);

    %title('Интенсивность обнаружения дефектов в системе с однородными элементами');

    xlabel('время t'); ylabel('Интенсивность q(t)');

    legend('qpos(t)','qmix(t)','qpar(t)');

    x=['Tpoc_odn=',num2str(Tpoc_odn),' Tmix_odn=',num2str(Tmix_odn),' Tpar_odn=',num2str(Tpar_odn),...

    ' Tpoc_raz=',num2str(Tpoc_raz), ' Tmix_raz=',num2str(Tmix_raz), ' Tpar_raz=',num2str(Tpar_raz)];

    disp(x)

    Москва 2021 г.


    написать администратору сайта