Моделирование системы диагностирования, построенной с использованием структурных инвариантов, страница 2

Рис. 1.1.7. Графики сигналов в каналах при большом сдвиге (M = 3, D = 1).

Рис. 1.1.8. Графики оценок инвариантов при большом сдвиге (M = 3, D = 1).

Рис. 1.1.9. Графики сигналов в каналах при малом изменении уровня шумов

(M = 0, D = 1.5).

Рис. 1.1.10. Графики оценок инвариантов при малом изменении уровня шумов

(M = 0, D = 1.5).

Рис. 1.1.11. Графики сигналов в каналах при среднем изменении уровня шумов

(M = 0, D = 2).

Рис. 1.1.12. Графики оценок инвариантов при среднем изменении уровня шумов

(M = 0, D = 2).

Рис. 1.1.13. Графики сигналов в каналах при большом изменении уровня шумов

(M = 0, D = 3).

Рис. 1.1.14. Графики оценок инвариантов при большом изменении уровня шумов

(M = 0, D = 3).

1.2.  Разность между параметрами ИК и средним арифметическим.

clear all

count = 5000;   % число измерений

M = 0;          % новое значение МО

D = 1;          % новое значение Д

Dfm = 2500;     % момент возникновения неисправности

T=1:count;     

for i=1:count

    if (i < Dfm)

        base1(i) = normrnd(0,1);

        base2(i) = normrnd(0,1);

        base3(i) = normrnd(0,1);

    else

        base1(i) = normrnd(M,D);

        base2(i) = normrnd(0,1);

        base3(i) = normrnd(0,1);

    end

    % Инвариант

    delta1(i) = base1(i) - 1/3*(base1(i)+base2(i)+base3(i));

    delta2(i) = base2(i) - 1/3*(base1(i)+base2(i)+base3(i));

    delta3(i) = base3(i) - 1/3*(base1(i)+base2(i)+base3(i));

    % Оценка инварианта

    mu1(i) = mean(delta1);

    mu2(i) = mean(delta2);

    mu3(i) = mean(delta3);

    d1(i) = std(delta1);

    d2(i) = std(delta2);

    d3(i) = std(delta3);

end

figure('color', 'white')

subplot(3,1,1);

plot(T,base1,'g');

ylabel('y1');

subplot(3,1,2);

plot(T,base2,'b');

ylabel('y2');

subplot(3,1,3);

plot(T,base3,'r');

xlabel('t');

ylabel('y3');

figure('color', 'white')

% Оценка МО

subplot(2,1,1);

plot(T,mu1,'g',T,mu2,'b',T,mu3,'r');

ylabel('M');

% Оценка Д

subplot(2,1,2);

plot(T,d1,'g',T,d2,'b',T,d3,'r');

xlabel('t');

ylabel('D');

Результаты:

Рис. 1.2.1. Графики сигналов в каналах при отсутствии дефектов (M = 0, D = 1).