3.1. Метод последовательного функционального анализа

3.2. Метод половинного разбиения

3.3. Метод «время-вероятность»

3.4. Метод на основе информационного критерия

3.5. Инженерный метод

3.6. Метод на основе иерархического принципа

3.7. Метод ветвей и границ

3.8. Метод диагностики на основе «белого шума»

3.9. Метод диагностики на основе ортогонального анализа отклика системы по базису гармонических функций

Выбор того или иного метода использования информации о техническом состоянии диагностируемой аппаратуры обусловлен структурой ОД и требуемой глубиной поиска неисправностей. Он накладывает определенные требования на принципы построения и структуру системы контроля и диагностики.

Вид алгоритма (программы) поиска неисправностей существенно влияет на эффективность процесса контроля и диагностики. При разработке алгоритма поиска обычно решают две задачи:

а) определяют наилучший набор контролируемых параметров;

б) получают наилучшую последовательность измерения контролируемых параметров;

Рассмотрим наиболее распространенные методы построения алгоритмов поиска неисправностей в РЭА.

3.1. Метод последовательного функционального анализа

Построение алгоритма поиска неисправностей методом последовательного функционального анализа, заключается в последовательном измерении параметров схемы, начиная «с конца», то есть с выхода устройства, где либо вовсе отсутствует сигнал, либо он находится вне допуска, до тех пор, пока не появится на каком – либо блоке (узле, элементе) сигнал, находящийся в поле допуска. Такой поиск можно проводить и «с начала», с входа устройства. Основные функции, характеризующие исправность ОД, например:

а) электрического питания;

б) приема и преобразования сигналов заданной частоты;

в) управления;

г) генерирования колебаний.

Рассмотрим составление алгоритма поиска неисправностей на примере канала звукового сопровождения телевизора типа УПИМЦТ (Рис. 4). Основной функцией этого канала является усиление и преобразование сигналов звукового сопровождения. Эта функция выполняется, если при входном сигнале с фиксированными параметрами на выходе будет наблюдаться вполне определенный сигнал. В процессе контроля принимается решение об исправности или неисправности канала. Последовательно контролируя сигналы на выходе каждого каскада, можно определить неисправный каскад. Получающуюся при этом схему поиска называют деревом функций (Рис. 5), а решения представляют обычно в виде матрицы:

Табл. 3

Обозначаются по схеме

Решение

Р0

Канал звука исправен

Р1

Неисправен УЧЗ

Р2

Неисправен РГ

Р3

Неисправен фильтр

Р4

Неисправна ИМС

Р5

Неисправна схема РГ

Р6

Неисправна ИСМ

Р7

Неисправен контур 6,5 МГц

Р8

Неисправен контур частного дискриминатора

Рис. 5 Схема поиска неисправностей

Данный метод прост, нагляден, требует минимум информации от ОД, однако он не оптимален ни по времени, ни по средним затратам.

3.2. Метод половинного разбиения

Этот метод часто используется в разработке алгоритмов поиска неисправностей в РЭА с последовательно соединенными элементами. Первым для контроля выбирается параметр, делящий всю схему пополам. При положительном результате контроля (сигнал на выходе элемента находится в допуске) следующим для контроля выбирается элемент, делящий неисправную часть схемы пополам и так далее до определения неисправного элемента (блока, узла). Такой алгоритм возможен в том случае, когда вероятности состояний P(Si) одинаковы для всех элементов, стоимости контроля выходных параметров Zi также одинаковы.

В том случае, когда вероятности состояний P(Si) для функциональных элементов неодинаковы, тогда вероятности необходимо первым контролировать такой параметр Zk , который делит ОД на части, вероятности состояния которых близки к 0,5.

Неопределенность состояния ОД до контроля оценивается величиной энтропии

N

H0 = - ∑ P(Si)*log2 P(Si) = log2 N (13)

i=1

Неопределенность состояния ОД при контроле параметра Zk будет:

H(Zk) = - (Pk *log2 Pk + (1- Pk) *log2 (1- Pk)) , (14)

N

где Pk = ∑ P(Si) , i = 1, 2, 3, …

i=1

Величина H(Zk) будет максимальна, если разность (Pk – 0,5) минимальна.

После контроля параметра (Zk) ОД будет разделен на две части: первая содержит K, а вторая (N - K) элементов. При выборе очередного параметра для контроля необходимо вероятности состояний в каждой из этих частей пронормировать, пересчитать по формулам:

k

P'(Si) = P(Si)/ ∑ P(Si) , i = 1, 2, 3, … k (15)

i = 1

N

P" (Si) = P(Si)/ ∑ P(Si) , i = k + 1, k + 2, … , N (16)

i=k+1

При этом

k N

∑ P' (Si) = 1 и ∑ P" (Si) = 1. (17)

i = 1 i=k+1

Тогда вторым параметром выбирается Zll , который делит одну из частей на две, вероятности которых

l

P" (Si) = 0,5 (18)

i=1

Такое деление продолжается до тех пор, пока состояние ОД не будет определено с заданной глубиной.

Метод половинного разбиения применим и для случая, когда в ОД неисправно несколько элементов.

3.3. Метод «время-вероятность»

Этот способ находит применение для РЭА, в которой функциональные элементы соединены произвольно и имеют разные вероятности P(Si) состояний и различные стоимости проведения контроля параметров С(Zi). Эффективность метода оценивается средним временем поиска неисправного элемента или средним временем контроля одного параметра. Последовательность контроля параметров устанавливается в порядке уменьшения величины:

P(S1)/ t1> P(S2)/ t2> …> P(SN)/ tN (19)

Располагая в порядке уменьшения величины P(Si)/ti , получим следующую последовательность для контроля параметров:

ZN → … → ZK → … → ZM → … → ZC

3.4. Метод на основе информационного критерия

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

Исходными данными являются функциональная модель и таблица неисправностей.

Предварительно ОД разделяются на N функциональных элементов, вероятности состояний, которых принимаем одинаковыми

P(Si) = P(S1) = P(S2) = …= P(SN) =1/N (20)

Неопределенность состояний ОД до контроля определяется оценивается величиной энтропии

H0 = log2 N (21)

Результат контроля к - го параметра ОД дает некоторое количество информации о его контроле:

IK = H0 - HK (22)

где HK - средняя условная энтропия ОД при условии контроля к - го параметра.

HK = P(Z'K) НZ'K + P(Z0K) НZ0K (23)

P(Z'K) = 1/N ; P(Z0K) = (N – m)/N, (24)

где m – количество единиц в к - ой строке.

HK = log2 m + log2 (N – m) (25)

Контроль к – го параметра дает следующее количество информации:

(26)

Последовательно вычисляем значения IK (где к = 1, N) и по убыванию IK определяем значимость параметра ZK . Первым контролируется параметр ZK , дающий максимальное количество информации.

После контроля первого параметра определяют количество информации, получаемое при контроле каждого n оставшегося параметра относительно состояния, характеризующегося энтропией НZK . Условная энтропия

H(zn/zK) = P(z'n/z'K)*Hz'n/z'K + P(z0n/z0K)*Hz0n/z0K +

+P(z0n/z'K)*Hz0n/z'K + P(z'n/z0K)*Hz'n/z0K (27)

где P(z'n/z'K) = m1/N – вероятность положительного решения при контроле параметра Zn в случае положительного решения при контроле параметра ZK; m1 – количество единиц в n-ой строке таблицы состояний относительно m единиц в к-ой строке; m2 – количество единиц в n-ой строке относительно (N – m) нулей к-ой строки.

(28)

(29)

(30)

Hz'n/z'K = log2 m1 ; (31)

Hz0n/z'K = log2 (m - m1 ) ; (32)

Hz'n/z0K = log2 m2 ; (33)

Hz0n/z0K = log2 (N - m - m2 ) ; (34)

I (zn / zK) = HK – H(zn / zK) . (35)

Выражение для вычисления количества условной информации имеет вид:

(36)

По максимуму условной информации выбирается второй контролируемый параметр. По такой же схеме выбираются все остальные параметры.

После всех расчетов строим схему поиска неисправностей.

Пример 3.4.1

Рис. 6 Функциональная модель ОД

Табл. 4

Zi

Si

S1

S2

S3

S4

S5

Z1

0

1

1

1

1

Z2

0

0

0

1

1

Z3

1

1

0

1

1

Z4

1

1

0

0

1

Z5

0

0

0

0

0

Из анализа табл. 4 находим, что контроль параметра Z5 для поиска неисправностей не дает никакой информации, поэтому его можно из дальнейшего рассмотрения исключить. Тогда энтропия (21) до контроля будет

H0 = log2 5 = 2,32

Количество информации (22) при контроле каждого параметра следующее:

I1 =

I2 =

I3 =

I4 =

Для контроля берем Z2 . После его контроля могут быть приняты два решения: значение параметра Z2 в допуске – функциональные элементы 1,2,3 исправны, а неисправность в элементе 4 или 5; (см. Рис. 6); значение параметра Z2 не в допуске - функциональные элементы 4 и 5 исправны, а не исправность в элементах 1, 2, 3 .

В соответствии с этим решением перестраиваем матрицу состояний (табл. 5)

Табл. 5

Zi

Si

S4

S5

S1

S2

S3

Z2

1

1

0

0

0

Z1

1

1

0

1

1

Z3

1

1

1

1

0

Z4

0

1

1

1

0

Теперь вычислим количество информации (36), которое дает контроль параметров Z1 , Z3 , Z4 при условии, что Z2 проконтролирован:

Следовательно, вторым для контроля выбираем Z4 . Исключим из табл. 5 строку Z2 (табл. 6)

Табл. 6

Zi

Si

S4

S5

S1

S2

S3

Z4

0

1

1

1

0

Z1

1

1

0

1

1

Z3

1

1

1

1

0

В результате построения алгоритма поиска неисправностей в заданном ОД получаем, что для поиска неисправностей достаточно контролировать последовательность из трех параметров (Z2 , Z4 , Z1) по определенной схеме (Рис. 7). Контроль параметра Z4 при условии, что Z2 = 1 дает два решения: если Z4=0, то неисправен элемент 4, если Z4=1, то неисправен элемент 5.

Если Z2 = 0 и Z4 = 0, то неисправен элемент 3 . Если Z2 = 0 и Z4 = 1, то надо контролировать параметр Z1 . Если Z1 = 0, то неисправен элемент 1, если Z1 = 1, то неисправен элемент 2 (Табл. 7).

Табл. 7

Zi

Si

S4

S5

S1

Z4

1

1

0

Z1

0

1

1

Z3

1

1

0

Рис.7 Схема поиска неисправностей

3.5. Инженерный метод

Этот метод построения алгоритмов диагностирования основан на вычислении функции предпочтения. При этом исходными данными являются функциональная модель ОД и таблица неисправностей.

Функция предпочтения выбирается в соответствии с решаемой задачей диагностики и исходными данными. При этом рассматриваются три случая определения перечня параметров:

а) для оценки работоспособности;

б) для поиска неисправностей;

в) для оценки работоспособности и поиска неисправностей.

Функция предпочтения при оценке работоспособности

N

W1= max Wi (Zi) = S0 (i j) , (37)

i N j=1

где S0 (i j) = 1, если состояние i j–го матричного элемента описывается нулем, и S0 (i j) = 0, если состояние i j–го матричного элемента описывается единицей.

Первым для контроля берут параметр Zi, у которого функция W1 (Zi) имеет максимальное значение. В результате контроля данного параметра матрица состояний делится на две части. В одну часть входят состояния, для которых результаты контроля выбранного параметра имеют положительные, а в другую отрицательные решения. В дальнейшем надо контролировать только первую часть, так как нужно фиксировать только факт исправности или отказа. Для нее аналогичным образом нужно вычислить значения Wi (Zi) и выбрать для контроля параметр по максимуму Wi (Zi).

Пример 3.5.2

Рис. 8 Функциональная модель ОД

Табл. 8

Zi

Si

Wi

S1

S2

S3

S4

S5

S6

Z1

0

1

1

1

1

1

1

Z2

0

0

1

1

1

1

2

Z3

0

0

0

1

1

1

3

Z4

0

0

1

0

1

0

4

Z5

0

0

0

0

0

0

6

Z6

0

1

1

1

1

0

2

Первым для контроля выбираем Z5 по максимуму W1 (Z5). В результате контроля в строке Z5 положительное решение имеет только Z3 , следовательно следующим для контроля следует выбрать Z3 .

Если для ОД известны вероятности состояний Р(Si), то функция предпочтения

W2 = max W2 (Zi) , (38)

i N

N

W2 (Zi) = P (Si)S0 (i j) , (39)

i =1

Если для ОД известны и стоимости контроля параметров С(Zi), то функция предпочтения

(40)

Рассмотрим составление алгоритма поиска неисправностей. Максимальное количество информации дает контроль параметра, который делит все возможные состояния на две равные части. Тогда функция предпочтения

N N

W4 (Zi) = min |S0 (i j) -S1 (i j)| (41)

j =1 j =1

где S1(i j) = 1, если состояние матричного элемента описывается единицей, и S1 (i j) = 0, если состояние матричного элемента описывается нулем, S0 (i j) = 1, если состояние матричного элемента описывается нулем, S0 (i j) = 0, если состояние матричного элемента описывается единицей.

Первым для контроля выбирается параметр, для контроля которого функция предпочтения W4 минимальна. Результаты контроля Zi разделяют матрицу состояний на две части. В первую часть входят состояния, для которых результат контроля отрицателен (на пересечении Zi – SJ стоит нуль), во вторую часть входят состояния, для которых результат контроля положителен. Для полученных частей вновь вычисляют функции предпочтения. По их минимальным значениям выбирают второй контролируемый параметр (при положительном исходе контроля первого выбранного параметра) или третий контролируемый параметр (при отрицательном результате контроля первого выбранного параметра). Совокупность параметров и последовательность их контроля выбирается до тех пор, пока множество N возможных состояний ОД не будет разделено на отдельные различимые состояния с требуемой глубиной поиска.

В случае, когда заданы вероятности состояний P(Si) диагностируемой аппаратуры и стоимости контроля параметров С(Zi) функциональных элементов, функции предпочтения примут вид

N N

W5 (Zi) = min |P(Si)S0 (i j) -P(Si)S1 (i j)| (42)

j =1 j =1

N N

W6 (Zi) = min С(Zi) |P(Si)S0 (i j) -P(Si)S1 (i j)| (43)

j =1 j =1

Пример 3.5.3

Рис. 9 Функциональная модель ОД

Табл. 9

Zi

Si

S1

S2

S3

S4

S5

S6

S7

S8

Z1

0

1

1

1

1

1

1

1

Z2

0

0

1

1

1

1

1

0

Z3

0

0

0

1

1

1

1

0

Z4

0

0

0

0

1

1

1

0

Z5

0

0

0

0

0

1

1

0

Z6

0

0

0

0

0

0

1

0

Z7

0

0

0

0

0

0

0

0

Z8

1

1

1

1

1

1

1

0

Для каждого параметра Zi табл. 9 по формуле (41) определяется значение W4 . Полученные значения записывают в табл. 10 и по минимальному значению выбирают для контроля параметр Z3 .

Табл. 10

В результате контроля параметра Z3 табл. 10 делится на две части: в первую часть входят функциональные элементы 4, 5, 6 и 7, для которых результат контроля Z3 положителен, в другую часть входят функциональные элементы 1, 2, 3 и 8, для которых результат контроля отрицателен, что означает неисправность одного из этих элементов. При поиске неисправностей контролю подлежат обе части матрицы, поэтому каждую часть из них представляем в виде табл. 11, а другую табл. 12. Для каждой из полученных таблиц определяем значения функций предпочтения W4. По минимуму W4 в табл. 11 выбираем для контроля параметр Z1 , в табл. 12 выбираем для контроля параметр Z5 . Аналогичные деления показаны в табл. 13 – 15. Схема поиска неисправностей показана на рис. 10.

Табл. 11

Табл. 12

Табл. 13

Табл. 14

Табл. 15

Рис. 10 Схема поиска неисправностей в диагностируемом устройстве

3.6. Метод на основе иерархического принципа

Построение алгоритмов диагностирования по иерархическому принципу целесообразно использовать для РЭА со встроенными устройствами контроля. При данном методе N первичных функциональных элементов ОД разбиваются к групп по N1 элементов в каждой группе. Выходные параметры первичных функциональных элементов объединяются в одной точке с измерительным устройством и индикатором неисправности. Таких индикаторов будет к штук. Последние еще разбиваются на r групп по N2 штук в каждой группе и т. д. В результате придем к одному индикатору неисправности.

В такой системе при выходе из строя функционального элемента последний индикатор покажет неисправность ОД. Для обнаружения неисправного функционального элемента просматриваются показания индикаторов первой ступени и при обнаружении индикатора, указывающего на неисправность, просматриваются индикаторы следующей ступени, соединенные только с этим индикатором. Проверки проверяются в указанной последовательности до тех пор, пока не будет обнаружен неисправный первичный элемент (рис. 11).

Поиск неисправного первичного функционального элемента по приведенной схеме позволяет значительно сократить время поиска по сравнению с поиском среди N элементов. Минимальное среднее время поиска неисправностей:

(44)

где

P(N) – вероятность пропуска неисправного элемента.

При и достигает максимума при z=4. Тогда количество ступеней иерархии

e = lnN /lnz (45)

Рис. 11 Схема поиска неисправности по иерархическому принципу

3.7. Метод ветвей и границ

Способ построения алгоритма поиска неисправностей методом ветвей и границ используется для синтеза алгоритмов поиска неисправностей в РЭА, функциональная модель, которой представляет собой произвольную структуру. Он позволяет определить наилучшую последовательность поиска среди возможных. Исходными данными являются функциональная модель ОД, таблица неисправностей с вероятностями различных состояний и стоимости контроля выходных параметров функциональных элементов. Целью данного метода является определить такую последовательность контроля параметров, которая будет обладать минимальной средней стоимостью при поиске любого неисправного элемента.

Средняя стоимость произвольной программы поиска неисправностей

или , (46)

где Ci – стоимость контроля i – параметра;

- сумма вероятностей состояний, которые рассматриваются при контроле i – параметра.

Расчет поиска начинается с любого i – параметра, который разбивает всё множество возможных состояний S на два подмножества: S0 (zi) и S1 (zi), соответствующие отрицательному и положительному результатам контроля параметра zi соответственно.

Последовательность контроля остальных параметров из приведенных подмножеств неизвестна, и определить значение средней стоимости алгоритма поиска невозможно. Поэтому значения средней стоимости заменяются их нижними границами Сн(S0) и Сн(S1) при контроле соответствующих параметров в подмножествах S0(zi) и S1(zi). Тогда нижняя граница средней стоимости всей программы поиска, которая начинается с контроля первого параметра, определяется как

(47)

Вычислив нижние границы стоимостей Сн(S0) и Сн(S1) для всех возможных алгоритмов поиска, выбирают первым такой параметр, контроль которого дает минимальную среднюю стоимость нижней границы алгоритма поиска.

Затем для подмножеств S0(zi) и S1(zi) вычисляют нижние границы стоимостей для всех возможных пар контролируемых параметров между первым и оставшимися для каждого подмножества.

Пусть для подмножества S0(zi) контролируется параметр zк , а для подмножества S1(zi) – параметр zj . Тогда нижняя граница средней стоимости программы поиска, начинающейся с контроля i – параметра, будет:

(48)

В общем виде записывается:

(49)

где (50)

(51)

где и - суммы вероятностей в подмножествах и соответственно.

Значения рассчитываются по формулам (46-48) для подмножеств

- подмножество состояний ОД, образующееся при отрицательном результате контроля параметров Ziи Zк;

- подмножество состояний ОД, образующееся при отрицательном результате контроля параметра Ziи положительном результате контроля параметра Zк;

- подмножество состояний ОД, образующееся при положительном результате контроля параметра Zк и отрицательном результате контроля параметра Zi;

- подмножество состояний ОД, образующееся при положительном результате контроля параметром Zк и Zi;

Вторым выбирается такой параметр, при контроле которого обеспечивается минимальная средняя стоимость нижней границы из всех возможностей алгоритмов поиска.

Аналогичным образом выбирают третий параметр и последующие параметры, пока получаемые пи контроле подмножества будут содержать более двух состояний ОД.

Рис. 12 Функциональная модель ОД

Матрица неисправностей и вероятности состояний

Таб.16

Si

Zi

Pi

Z1

Z2

Z3

Z4

Z5

Z6

S1

0

0

1

0

0

0

0,2

S2

1

0

1

0

0

0

0,05

S3

1

1

0

0

1

0

0,24

S4

1

1

1

0

1

0

0,15

S5

1

1

1

1

0

0

0,06

S6

1

1

1

1

1

0

0,3

Пусть первым контролируется параметр Z1. Тогда множество возможных состояний разбивается на два подмножества S0(Z1) и S1(Z1).Для подмножества S0(Z1) последовательность состояний вероятностей будет состоять из одного состояния Р(S1) =0,2 , а подмножества S1(Z1) составим первую последовательность: Р(S2) = 0,05; Р(S5) =0,06; Р(S4) =0,15; Р(S3) =0,24; Р(S6) =0,3. Затем определим сумму двух наименьших вероятностей: Р*1 = 0,05 + 0,06 = 0,1.

Составим вторую последовательность: Р*1 = 0,1; Р(S4) = 0,15; Р(S3) =0,24; Р(S6) =0,03. Затем определим Р*2 = 0,11 + 0,15 = 0,26.

Составим третью последовательность: Р(S3) = 0,24; Р*2 = 0,26; Р(S6) =0,3. Затем определим Р*3 = 0,24 + 0,26 = 0,5.

И наконец, составим четвертую последовательность: Р(S6) =0,3; Р*3 = 0,5 и определим Р*4 = 0,3 + 0,5 = 0,8.

При этом нижняя граница средней стоимости

N 0 6-1-1

Сн(Z1) = С [ ∑ Р(Si) + ∑ Рλ + ∑ Р*ν] = C (1 +0 +1,67) = 2,67 C

i=1 λ=1 ν=1

Следует указать, что если подмножество S0(Zi) содержит одно состояние, то нижняя граница стоимости Сн(S1) = 0.

Разделение множества возможных состояний на два подмножества S0(Z1) и S1(Z1) при контроле параметра Z1 как показано в табл. 16.

При контроле параметра Z2 множество возможных состояний разбивается на два подмножества: S0(Z2) и S1(Z2) (табл. 17).

Таб.17

Si

Zi

Pi

Z1

Z2

Z3

Z4

Z5

Z6

S1

0

0

1

0

0

0

S0 (Z2)

S2

1

0

1

0

0

0

S3

1

1

0

0

1

0

S0 (Z2)

S4

1

1

1

0

1

0

S5

1

1

1

1

0

0

S6

1

1

1

1

1

0

Для подмножества S0(Z2) последовательность двух значений вероятностей будет из двух значений: Р(S2) = 0,05 и Р(S1) = 0,2 , и тогда Р*1 = 0,05 + 0,2 = 0,25.

Для подмножества S1(Z2) первая последовательность будет иметь вид Р(S5) =0,06; Р(S4) =0,15; Р(S3) =0,24; Р(S6) = 0,3 , и тогда Р*1 = 0,21; Р*2 = 0,45; Р*3 = 0,75.

Следовательно, нижняя граница средней стоимости при контроле параметра

N 0 6-1-1

Сн(Z2) = С [ ∑ Р(Si) + ∑ Рλ + ∑ Р*ν] = C (1 +0,25 +1,41) = 2,66 C.

i=1 λ=1 ν=1

Поступая аналогичным образом при контроле параметров z3 , z4 и z5 , получим нижние границы стоимости при контроле соответствующих параметров: СН ( z3 ) = 2,59C; СН (z4) = 2,60С; СН (z5) = 2,50С.

Из результатов вычислений нижних границ средней стоимости следует, что минимальной стоимостью будет обладать алгоритм поиска, начинающийся с контроля параметра z5.

Посмотрим матрицу (табл. 18) состояний подмножества:

S0 (z5) и S1 (z5).

Анализируя табл. 18 , видим, что в случае отрицательного исхода при контроле параметра z5 следующими можно контролировать параметры z1 , z2 , z5. Параметр z3 контролировать не имеет смысла, так как подмножество S0(z5) не разбивает при его контроле на более мелкие подмножества. При положительном исходе можно контролировать параметры z3 , z4.

Таким образом, необходимо вычислить значения стоимостей нижних границ алгоритма поиска при контроле таких наборов параметров z5 , z1 , z3 , z5 , z1 , z5 , z2 , z3 , z5 , z2 , z4. Для этого воспользуемся формулой:

,

где

;

Таблица 18

S1

ZI

P(S1)

Z5

Z1

Z2

Z3

Z4

S1 0 0 0 1 0 0,2 S2 0 1 0 1 0 0,05 S3 0 1 1 1 1 0,06 S4 1 1 1 0 0 0,24 S5 1 1 1 1 0 0,15 S6 1 1 1 1 1 0,3

,

где ∑ Р(SK) и ∑ Р(SJ) суммы вероятностей Р(S) диагностируемого объекта, которые различаются контролем параметров zK и zJ в подмножествах и соответственно.

Вычислим значения стоимостей нижних границ при контроле параметров z5 , z1 , z3 . Контроль параметров z1 позволяет различать состояния S1 , S2 , S5 , следовательно,

∑ P(SK) = Р(S1) + Р(S2) + Р(S2) = 0,2 + 0,05 + 0,06 = 0,31

Также при контроле параметра zI выделяется из подмножества S0(zК) только одно состояние S001 , тогда CН (S001) = 0 , а значения CН (S101) вычисляется аналогичным способом , как и CН (S01) :

К

CН (S101) = С ∑Р*V = С [Р (S2) + Р (S5)] = C (0,05 + 0,06) = 0,11C

При контроле параметра z3 различаются состояния S3 , S4 , S6 , следовательно:

∑ P(SJ) = Р(S3) + Р(S4) + Р(S6) = 0,24 + 0,15 + 0,3 = 0,69

Кроме того, при отрицательном исходе из множеств S1(z5) выделяется только одно состояние, а при положительном – два , следовательно:

CН (S10J) = CН (S103) = 0

К

CН (S11J) = С ∑ Р*V = С [Р (S4) + Р (S6)] = C (0,15 + 0,3) = 0,45C

V

Тогда значение стоимости нижней границы при контроле параметров z5 , z1 , z3 будет равно

CН (z5 , z1 , z3,) = С(1 + 0,31 + 0,11 + 0,69 + 0,45) = 2,56С

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

CН (z5 , z1 , z4,) = 2,5С

CН (z5 , z2 , z4,) = 2,64С

CН (z5 , z2 , z3,) = 2,7С

Анализ результатов вычислений стоимостей нижних границ показывает, что минимальной стоимостью будет обладать последовательность параметров z5 , z1 , z4.

Таблица19

Преобразуем таблицу неисправностей в соответствии с результатами контроля последовательности параметров z5 , z1 , z4 (табл. 19) .

Из таблицы видно, что при отрицательном исходе контроля параметров z5 и z1 выделяется состояние S1 , а при отрицательном исходе контроля z5 и положительном z1 необходимо следующим контролировать параметр z2 (контроль параметра z3 не позволяет разделить подмножество S01 на более мелкие).

При положительном исходе контроля параметров z5 и z4 выделяется состояние S6 , а при отрицательном исходе контроля параметров z5 и отрицательном z4 необходимо контролировать параметр z3. Поскольку других вариантов нет, на этом заканчивается построение алгоритма поиска неисправностей, обладающего минимальной средней стоимостью (Рис. 13).

Таким образом, последовательность контролируемых параметров z5 , z1, z4, z2 , z3 дает среднюю стоимость ССР = 2,5С.

Рис. 13 Схема алгоритма поиска неисправностей, построенного методом ветвей и границ

3.8. Метод диагностики на основе «белого шума»

Интегральный метод диагностики на основе «белого шума» позволяет определить отклик диагностируемой системы во время ее нормальной работы, то есть без отключения. Известно, что для «белого шума» автокорреляционная функция равна нулю при всех значениях τ, кроме точки τ = 0, где случайная функция умножается сама на себя:

(51)

где Uш1 (t) – напряжение белого шума на входе ОД.

Связь между напряжениями входа и выхода диагностируемой системы через отклик h(t) определяется уравнением

(52)

Известно, что взаимная корреляционная функция входа и выхода системы выражается через автокорреляционную функцию

(53)

После перемены аргументов τ и t получим

(54)

Если на вход системы подать стимулирующее напряжение

δ(t) = U1(t) (55)

то уравнение (52) упрощается и напряжение на выходе системы становится равным отклику диагностируемой аппаратуры:

U2 (t) = h (t) (56)

Из сравнения выражений (52), (53) и (56) следует, что при подаче на вход ОД сигнала в виде «белого шума», корреляционная функция которого является дельта-функцией R11(t) = δ(t), взаимно корреляционная функция входа и выхода будет равна отклику системы:

R12(t) = h (t) (57)

Рис. 14 Схема для определения отклика диагностируемой системы.

Равенство (7) позволяет синтезировать схему (рис.1), которая с помощью коррелятора может определить отклик системы по напряжению «белого шума» на ее входе.

Если на вход подать напряжение Uс1 (t), то на ее выходе получим суммарный сигнал

U2 (t) = Uш2 (t) + Uс2 (t) (58)

где Uс2 (t) – напряжение, которое появляется на выходе контролируемой аппаратуры после прохождения через нее сигнала Uс1 (t), Uш2 (t) – напряжение на выходе аппаратуры.

3.9. Метод диагностики на основе ортогонального анализа отклика системы по базису гармонических функций

Известно, что динамические свойства любой радиотехнической системы можно описать её откликом h(t). Если отклик системы разложить в ряд Фурье и установить аналитическую связь между коэффициентами ряда Фурье для отклика и параметрами диагностируемой системы, то на этой основе можно проводить диагностирование. Отклик любой системы

< (59)

где М и С0 – постоянные положительные действительные числа:

h(t) = 0 при t < 0 и h(t) 0 при t > 0.

Для отклика системы, удовлетворяющего условию (59) можно применить преобразование Лапласа:

(60)

где функция h(t) является оригиналом, а функция Н(р) является изображением функции по Лапласу.

Учитывая условие (59) в формулах (60) можно положить С = 0 и р = .

Тогда получим преобразование Фурье:

(61)

Разложение функции h(t) в ряде Фурье в комплексной форме имеет вид:

(62)

где t0 – для отклика h(t);

Сn – комплексные коэффициенты.

(63)

Из сравнения (60) и (63) получаем формулы для определения коэффициентов Фурье:

, при (64)

, при (65)

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

S2 (P) = S1 (P)*K(p) (66)

Зная параметры элементов системы, можно вычислить её передаточную функцию К(р). Если параметры входного сигнала S1(P) будут неизменными, то любые отклонения параметров системы за допустимые пределы будут отражаться в спектре выходного сигнала.

Связь между параметрами системы и коэффициентами ряда Фурье разложения передаточной функции можно продемонстрировать на примере.

Пример 3.9.5 Передаточная функция видеосигнала.

(67)

где {х} = К, Т1, Т2 – диагностируемые параметры (вторичные).

Решение. Определить вещественную и мнимую часть передаточной функции:

;

;

;

ω = n*ω0 , где n – номер гармоники;

ω0 – частота основной гармоники.

Для определения текущих значений диагностируемых параметров, составим систему уравнений, используя действительную часть передаточной функции:

;

(68)

.

Используя мнимую часть передаточной функции, получим другую систему уравнений:

;

(69)

.

В дальнейшем для упрощения вычислений используем относительные коэффициенты разложения:

В1 = b1/a1 ; B2 = b2/a2 ; …; Bn = bn/a2 . (70)

Используя текущие значения в относительных коэффициентов, составим систему уравнений:

(71)

Систему уравнений приведем к виду:

(72)

Решая систему уравнений (71) и (72) относительно Т1 и Т2 , получим выражения для параметров диагностируемой системы:

(73)

Для определения коэффициента усиления К видеоусилителя воспользуемся уравнением:

. (74)

Отсюда:

. (75)

Полученные значения параметров сравниваются с их номинальными значениями. Параметр, вышедший из допуска, и определяет место неисправности.