Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

МАОС без 1 главы

.pdf
Скачиваний:
47
Добавлен:
22.03.2016
Размер:
1.43 Mб
Скачать

3.3. Обнаружение сигналов с двоичным накоплением

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

x(t)

1

z(t)

2

Z(t)

y(t)

A*1

3

4

 

 

 

 

 

 

 

 

 

A*0

 

 

 

 

 

Z0

k0

Рис. 3.2

Однако при обнаружении сигналов в ряде случаев хорошие результаты дает и простейшее квантование на два уровня: бинарное или двоичное. Структурная схема такого обнаружителя приведена на рис. 3.2. Блок 1 – согласованный фильтр, блок 2 – амплитудный детектор. При превышении принимаемым сигналом Z(t) некоторого порогового уровня Z0 в пороговом устройстве 3 ему приписывается уровень 1 и вырабатывается стандартный импульс. В противном случае устанавливается 0. Далее производится счет стандартных импульсов. Если за заданное время (заданное число импульсов опроса или заданное число посылок l) число единиц достигает некоторого порогового значения k0 , то решающее устройство 4 принимает решение о наличии дефекта. На основании всего вышеизложенного правило принятия решения о наличии сигнала может быть записано в виде

l

 

1

k yi k0 ,

где

yi

i 1

 

0

Важная особенность обнаружителя пачки двоично-квантованных импульсов состоит в том, что у него кроме порога обнаружения пачки k0 имеется еще первый порог Z0 обнаружения отдельных импульсов. В соответствии с выражениями (2.20) и (2.22) вероятность превышения порога Z0 при отсутствии сигнала

 

 

 

2

 

 

 

 

 

Z0

 

,

2

PF exp

 

 

 

 

z

 

 

а вероятность превышения порога при наличии сигнала

63

 

 

 

Z

2

 

 

P

exp

 

2

 

0

 

 

D

 

E

2

 

 

z

s

 

 

 

 

 

.

Обозначим вероятность принятия правильного решения о наличии де-

фекта по результатам испытаний в l циклах P (l) , а вероятность ошибочного

 

 

D

решения (ложной тревоги) P

(l)

. Считая испытания в отдельных циклах неза-

 

F

 

висимыми, по формуле Бернулли можем записать, что вероятность превыше-

ния порога при наличии сигнала ровно в n циклах

n

C

n

D

)

n

(1

D

)

l n

,

 

 

 

 

 

 

 

 

P

l

(P

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где C

n

 

l!

– число сочетаний из l по n. Поскольку обнаружение про-

 

l

n!(l n)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

изводится при любых значениях n, не меньших k0 : k0 n l , то условная вероятность правильного обнаружения

P

(l)

 

 

D

 

l

n

 

C

(P

l

n k

D

 

 

0

 

 

) n (1 P D

)

l

 

n

.

(3.5)

Аналогично, условная вероятность ложной тревоги

P

(l)

 

 

F

 

l

n

 

C

(P

l

n k

F

 

 

0

 

 

) n (1 P F

)

l n

 

.

(3.6)

Можно показать [8], что при некоррелированных флуктуациях импульсов в пачке, когда справедливо приведенное рассмотрение, существует опти-

мальное значение второго порога

k0

opt

1,5

l , соответствующее минималь-

 

 

 

 

но необходимому отношению сигнал/помеха на входе накопителя. Проигрыш при этом по сравнению с идеальным некогерентным накопителем составляет не более 2 дБ.

При непрерывном сканировании обозреваемого пространства хранимая информация должна все время обновляться. Проще всего это выполнить, используя регистры сдвига. Схема решающего устройства с регистром сдвига изображена на рис. 3.3.

Устройство работает следующим образом. При сканировании на вход регистра сдвига 1 поступает с выхода порогового устройства 3 (рис. 3.2) сигнал y(t), состоящий из последовательности 0 и 1. Сигналы "0" и "1", поступающие на вход регистра 1, запоминаются в его первой ячейке 1.1 и затем продвигаются по регистру тактовыми импульсами, поступающими на регистр от специального генератора, срабатывающего с частотой посылок.

64

y(t)

1

 

1.11.2 . . . 1.(l-1)

A*1

2

A*0

Рис. 4.3

Сигналы с выходов ячеек регистра поступают на мажоритарную схему совпадений 2, имеющую l входов. Мажоритарная схема совпадений представляет собой устройство, сигнал "1" на выходе которого будет в случае, если "1" присутствует не менее чем на k0 входах. Примеры построения такого рода схем имеются в литературе [8], [9] . Таким образом, мажоритарная схема совпадений вырабатывает сигнал о наличии объекта лишь в случае, когда в последовательности из l посылок (импульсов опроса) имеется не менее k0 единиц (превышений принятым сигналом порогового уровня). Вероятности

правильного решения

P

(l)

и ошибочного решения

P

(l)

по-прежнему опреде-

 

 

 

D

 

F

 

ляются выражениями (3.5) и (3.6).

3.4. Последовательный обнаружитель

Рассмотренные ранее обнаружители пачек сигналов предполагают использование фиксированного времени наблюдения (время интегрирования Tи для обнаружителя с поcледетекторным интегрированием или максимальное количество импульсов k0 для обнаружителя с двоичным накоплением). Однако в некоторых случаях соотношение амплитуд сигнала и помехи позволяет надежно обнаруживать сигнал быстрее, чем в других циклах наблюдения. Поэтому, если заранее не фиксировать общую длительность наблюдения, можно за много циклов наблюдения получить значительную экономию времени. Такой метод анализа, названный последовательным, был предложен Вальдом [5]. При последовательном обнаружении производятся непрерывный анализ отношения правдоподобия (или монотонно связанной с ним величины) и сравнение его с двумя порогами l1 и l2:

65

l

 

1 P

l

 

P

 

D ,

D

1

 

1

P

2

 

P

 

 

 

 

 

 

 

F

 

 

F

, причем l2 l1.

Если выполняется условие l l1, то принимается решение о наличии только помехи. Если выполняется условие l l2 , то принимается решение о наличии полезного сигнала. И, наконец, если значение отношения правдоподобия лежит между порогами: l1 l l2 , то делается вывод, что имеющихся данных недостаточно для обоснованного заключения о наличии или отсутствии сигнала, и принимается решение «не знаю». При этом испытание продолжается до тех пор, пока не прекратятся ответы «не знаю».

Преимущество последовательного обнаружителя состоит в том, что можно независимо задавать вероятности PF и PD. Кроме того, он дает определенную экономию во времени принятия решения или в энергии излучаемого сигнала. Выигрыш по чувствительности (по значению отношения сигнал/помеха, необходимого для получения той же помехоустойчивости при том же среднем времени обнаружения) составляет от 3 до 20 дБ. Однако применение последовательного анализа требует значительного усложнения аппаратуры.

В настоящее время разработаны некоторые варианты последовательного анализа, упрощающие его реализацию. Например, в двухэтапном обнаружителе с первым порогом сравнивается результат накопления M импульсов. Если порог не превышен ни в одном элементе, то излучается еще N импульсов, производится сравнение со вторым порогом и выдается более уверенное решение (т. е. с меньшей вероятностью ложной тревоги) о наличии или отсутствии сигнала. Выигрыш по сравнению с однопороговой системой зависит от отношения сигнал/помеха, значений M и N и составляет не менее 2...3 дБ.

3.5. Обнаружение коррелированных сигналов

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

66

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

Рассмотрим пересечение стационарного дифференцируемого процесса x(t) с постоянным уровнем C . Согласно [6], среднее число выбросов этого процесса за уровень C :

N (C)

 

 

2

 

 

 

(C, x )dx

L

 

x p

 

0

 

 

,

где

x

dx

;

dt

 

 

p

(C, x )

2

 

p

(x, x )

x C

2

 

;

p

(x, x )

2

 

– двухмерная плотность веро-

ятностей ординаты случайного процесса x(t) и его производной; L – длина реализации. Интенсивность выбросов, т. е. количество выбросов на единицу длительности,

c x p2( C ,x )dx . 0

Для дальнейшего решения задачи необходимо задаться видом двухмерной плотности, т. е. видом закона распределения сигнала x(t). Рассмотрим случай ультразвукового теневого контроля стальных горячекатаных изделий. Как показывают эксперименты [12], [13], в этом случае закон распределения огибающей амплитуд U(t) можно считать гауссовским. Учитывая, что при теневом контроле на бездефектном участке среднее значение амплитуд должно быть выше порога, можно получить, что интенсивность ложной тревоги (перебраковки), т. е. среднее количество выбросов – ложных регистраций дефектов на единицу длины пути сканирования

п 1 2

 

 

 

1

mU C

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

 

.

r0

 

U

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

67

Если в изделии имеется протяженный дефект с акустической прозрачностью T , то на дефектном участке интенсивность пропуска цели (недобраковки), т. е. среднее количество пропусков в регистрации дефектов,

 

 

 

 

 

 

 

 

 

1

 

 

1 C m T

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

r0

exp

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь mU

– среднее значение огибающей амплитуд на бездефектном

участке контролируемого изделия;

2

дисперсия огибающей амплитуд;

U

 

d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

r

( )

 

 

– значение второй производной нормированной корреляци-

 

 

2

0

0

d

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

онной функции огибающей амплитуд rU ( )

в точке = 0.

 

Вводя в двух последних выражениях нормировку по уровню среднего

значения, получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

1

 

 

 

 

r exp

 

п

 

 

 

 

 

 

 

2

 

0

 

 

2

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

u /T

 

 

 

r exp

н

 

 

 

 

 

 

 

 

2

0

 

 

2

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

2

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

;

  

.

(3.7)

(3.8)

Здесь v U / mU – коэффициент вариации огибающей амплитуд. Формулы (3.7) и (3.8) позволяют получить выражения для относительного порога

регистрации u:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u 1 v

 

 

 

 

 

 

 

 

(3.9)

 

 

 

 

 

 

 

 

ln( r0 ) 2(ln п ln 2π) ,

 

 

 

u T 1 v

 

 

 

 

 

 

 

(3.10)

 

ln( r0 ) 2(ln νн ln 2π) .

 

 

Приравнивая правые части выражений (3.9) и (3.10), можно записать

выражение для максимальной акустической прозрачности Tmax

протяженного

дефекта,

который выявляется

с заданной

надежностью

(заданными

νп и νн

):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 v

 

 

п ln 2π)

 

 

 

 

Tmax

 

ln( r0 ) 2(ln

 

.

 

 

 

 

 

 

 

 

 

 

 

 

1 v

 

 

 

 

 

 

 

 

ln( r0 ) 2(ln νн ln 2π)

 

 

 

Взаимное расположение огибающей амплитуд и порогового уровня на бездефектном и дефектном участках показаны на рис. 3.4.

68

U

Без дефекта

С дефектом

mU

C mUT

t

Рис. 3.4

Проведенное рассмотрение справедливо для случая, когда необходимо учитывать лишь одномерные корреляционные связи огибающих амплитуд сигналов. Однако при контроле изделия по всей его площади при бесконечно малом шаге сканирования совокупность результатов контроля можно рассматривать как реализацию случайного поля амплитуд, что требует учета двухмерных корреляционных связей. Соответствующее рассмотрение на основе теории выбросов случайных полей проведено в [14] без учета и с учетом конечности значения шага сканирования в продольном и поперечном направлениях. Приведенные там же результаты расчетов показали, что в большинстве случаев для реальных параметров сканирования можно пренебречь корреляционными связями сигналов в одном направлении и рассчитывать чувствительность и надежность контроля по приведенным ранее формулам.

3.6. Достоверность результатов обнаружения

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

цели. При выполнении условия

1 u

2

выбросы гауссовского стационар-

v

 

 

 

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

69

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

цесса. Например, для гауссовского процесса при rU ( ) exp( 2 2 ) выбросы

независимы уже при

1 u

0.

v

 

Но если появление отдельных выбросов мож-

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

 

 

n

 

i

 

n

0

n

 

P e

 

 

 

i!

 

 

i 0

 

 

 

 

,

(3.11)

где n – среднее число выбросов на длине реализации L . Таким образом, задавшись доверительной вероятностью P получения не более n0 ложных регистраций сигнала на длине реализации L, можно определить среднее число регистраций n , и средняя интенсивность ложных регистраций п n / L .

Аналогично можно записать выражение для доверительной вероятности Pд того, что в записи сигнала длиною Lд будет не более nд пропусков:

 

 

n

 

n

д

n

j

д

 

 

 

 

e

 

пр

 

пр

P

 

 

 

 

 

 

 

j !

 

 

 

 

j 0

 

 

 

 

 

 

,

(3.12)

где nпр – среднее значение количества пропусков в записи сигнала длиною Lд

. Отсюда определяется средняя интенсивность пропусков (недобраковки):н nпр / Lд . Величины P, n0 , Pд , nд должны выбираться, исходя из необ-

ходимости обеспечения заданной надежности обнаружения и возможного снижения уровня ложной тревоги.

Точно так же могут быть определены вероятности ложной тревоги и пропуска цели и в случае, когда не требуется учета корреляционных связей сигналов в соседних точках. В этом случае задается доверительная вероятность P появления не более n0 ложных регистраций сигнала в N0 точках анализа. Среднее число ложных регистраций n при этом определится из условия (3.11), а вероятность ложной тревоги в каждой точке контроля PF n / N0 . Если минимальный подлежащий учету объект должен занимать не менее Mд элементов разрешения (Mд 1) и задана доверительная вероят-

70

ность PD получения не более nд пропусков регистрации такого объекта, то среднее число пропусков nпр определится из формулы (3.12) и вероятность

пропуска цели в каждой из точек наблюдения

P

ˆ

n

/M

Д

 

пр

 

D

 

 

 

. Полученные

значения PF и

P ˆ D

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

гов регистрации в рассмотренных схемах обнаружителей.

71

Список литературы

1.Харкевич А. А. Борьба с помехами. М.: Физматгиз, 1963.

2.Иванов М. Т., Сергиенко А. Б., Ушаков В. Н. Теоретические основы радиотехники. М.: Высш. шк., 2002.

3.Ольшевский В. В. Статистические методы в гидролокации. Л.: Судостроение, 1983.

4.Добротин Д. Д., Пигулевский Е. Д. Случайные сигналы и помехи в системах интроскопии: Учеб. пособие / ЛЭТИ. Л., 1990.

5.Тихонов В. И. Оптимальный прием сигналов. М.: Радио и связь, 1983.

6.Тихонов В. И. Статистическая радиотехника. М.: Радио и связь, 1982.

7.Гоноровский И. С. Радиотехнические цепи и сигналы. М.: Сов.радио,

1977.

8.Кузьмин С. З. Цифровая обработка радиолокационной информации. М.: Сов.радио, 1967.

9.Дымова А. И., Альбац М. Е., Бонч-Бруевич А. М. Радиотехнические системы. М.: Сов.радио, 1975.

10.Ермолов И. Н. Теория и практика ультразвукового контроля. М.: Машиностроение, 1981.

11.Горбунов В. И., Епифанцев Б. И. Автоматические устройства в радиационной дефектоскопии. М.: Атомиздат, 1979.

12.Голубев А. С., Добротин Д. Д., Паврос С. К. О выборе порога срабатывания теневых иммерсионных ультразвуковых дефектоскопов при контроле изделий с шероховатой поверхностью // Дефектоскопия. 1975. № 3. С.

71–77.

13.Добротин Д. Д. К вопросу о надежности ультразвукового теневого контроля горячекатаного металла // Электроакустика и ультразвуковая техника. - Л., 1977. С.47–52 (Изв. ЛЭТИ. Вып. 221).

14.Добротин Д. Д. Надежность сплошного ультразвукового контроля листов с шероховатой поверхностью// Электроакустика и ультразвуковая техника. Л., 1979. С.17–22 (Изв. ЛЭТИ. Вып. 252).

72