Определение характеристик оптимального обнаружения сигналов

МИНИСТЕРСТВО ОБЩЕГО И ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

ТАГАНРОГСКИЙ ГОСУДАРСТВЕННЫЙ РАДИОТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

КАФЕДРА ЭГА и МТ

КУРСОВАЯ РАБОТА

По курсу:

“Методы и системы обработки сигналов”

На тему:

«Определение характеристик оптимального обнаружения»

Выполнил: Озерин М.В. ­_______­______________

студент 4-го курса гр. Э-15 (дата) (подпись)

Руководитель: Черницер В.М. _____________________

(дата) (подпись)

Таганрог 1999

ВВЕДЕНИЕ

При проектировании гидроакустических систем (ГАС) различного функционального назначения на этапе проектирования решаются задачи оптимизации технических характеристик ГАС и выбора структуры приемного тракта, оптимизирующего отношение сигнал-помеха. Для таких систем основными параметрами является: дальность действия, пространственная разрешающая способность. В данном случае ищут компромиссное решение между этими параметрами, отдавая предпочтение одному или другому в зависимости от стоящей перед разработчиками задачи. Расчет ведется при фиксированной дальности и поэтому выбор оптимальной частоты определяется минимальной излучающей мощности. ЗАДАНИЕ

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

Таблица 1.1

Данные для расчеты

2. ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ

2.1. Характеристика обнаружения

Характеристиками обнаружения называются кривые, определяющие зависимость между вероятностью правильного обнаружения Р_о, вероятностью ложной тревоги P_лт и величиной сигнала выраженного в относительных еденицах:

(2.1)

Параметр q численно равен отношению сигнал-помеха (С/П) по напряжению на выходе согласованного фильтра (СФ). Если задано допустимое значение P_лт, то расписывается соответствующее значение P_o.

2.1.1 Случай полностью известного сигнала

Условные плотности вероятности корреляционного интеграла при отсутствии сигнала W(K/0) и при наличии сигнала W(K/1) определяются из выражения:

(2.2)

где величины K и E_s корреляционный интеграл и энергия сигнала. При сравнении значений К с порогом К_о вероятность ложной тревоги определяется отношением порогового уровня к среднеквадратичному значению s_к.

(2.3)

Вероятность правильного обнаружения зависит не только от отношения порога К_о среднеквадратичному значению s_к, но и от отношения

(2.4)

где q – параметр обнаружения.

(2.5)

В выражениях 2.3 и 2.5 в F(U) – интеграл вероятности:

(2.6)

Выражение 2.5 преобразуется к виду

(2.7)

где q_o=K_o/s_k.

Если интеграл вероятности определяется в виде

(2.6)

то выражение (2.3), (2.5), (2.6) приобретают вид

P_лт=1-Ф(q_o), (2.7)

P_лт=1-Ф(q-q_o), (2.8)

где q – параметр обнаружения.

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

Рис. 2.1 Кривые обнаружения

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

2.2. Случай сигнала со случайной начальной фазой

Условные плотности вероятности для корреляционного интеграла при наличии сигнала:

(2.9)

при отсутствии сигнала:

(2.10)

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

Максимально допустимая вероятность ложной тревоги

(2.11)

а пороговое значение отношение сигнал-помеха

(2.12)

Вероятность правильного обнаружения определяется, как

(2.13)

где S – переменная интегрирования.

Когда отношение сигнал-шум равен

формулы (2.9) и (2.13) упрощается, и расчет вероятности P_o можно вести по формуле

(2.14)

где Ф(U) – интеграл вероятности.

2.3. Случай со случайной амплитудой и начальной фазой

(2.15)

(2.16)

Вероятность ложной тревоги

(2.17)

Вероятность правильного обнаружения

(2.18)

Исключая q_o из (2.18), получим

(2.19)

В случае приема последовательности из n одинаковых когерентных импульсов энергетическое отношение сигнал/шум

(2.20)

где E_u/N_o – энергетическое отношение сигнал/шум, соответствующее одному импульсу последовательности.

По характеристикам обнаружения определяются значения q_n и пороговый сигнал, соответствующий полной энергии сигнала в пачке (E_S). Поэтому в случае когерентного обнаружения, энергия минимального порогового сигнала одного импульса должна быть – E_S/n. А в случае некогерентного обнаружения E_S/Ön. Выигрыш при когерентном приеме составляет Ön раз. Параметр обнаружения q может быть представлен как отношение максимального напряжения сигналаA_s к среднеквадратичного значения шума

(2.21)

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

для случая когерентного обнаружения

для случая некогерентного обнаружения

где W_и=A_s²/2 – импульсная мощность.

При n=1 различие между когерентным и некогерентным приемами отсутствует.

3. РАСЧЕТ ТЕХНИЧЕСКИХ ХАРАКТЕРИСТИК ГАС

Оптимальная частота работы ГАС

Оптимальную частоту выбираем из расчета, что сигнал будет иметь приемлемый шум и малое поглощение.

где r_max – дистанция до цели обнаружения (км).

Но так как можно перебирать частоту в некотором диапазоне, то выбираем частоту f_опт=39000, при этом получаем выигрыш в минимальном шуме, но имеем более сильное поглощение сигнала.

Полоса пропускания приемного тракта

Она складывается из доплеровского смещения частот и ширины спектра эхо-сигнала

Df=Df_д+Df_сп.

Найдем Df_д – доплеровское смещение частоты

где V_н – скорость носителя,

V_ц – скорость цели обнаружения,

с – скорость звука в среде.

Найдем Df_с – ширина спектра эхо-сигнала

Коэффициент 1,37 выбирается из того условия что отношение сигнал-шум является опртимальным для нашего случая.

где t_и=2×Dr/c=2×0,3/1483=0,67 (мс), где ×Dr – разрешающая способность по дальности. Тогда Df_сп=2032 (Гц).

Df=2032+2104=4136 (Гц).

Уровень шума, воздействующий на вход приемного тракта

Для расчета шума воспользуемся спектрально-энергетическими характеристиками шумов, в данном случай характеристикой для моря. Частота излученного сигнала равна 39000 Гц, тогда P_пр=2×10^-5 Па/Гц².

Уровень шумового давления на входе приемной антенны

P’_ш=U_ш/g,

где U_ш – уровень шумов на входе в приемный тракт и шум приемного тракта;

g - чувствительность антенны в режиме приема (мкВ/Па),

U_ш.эл – уровень шумов электронного тракта (мкВ).

тогда P’_ш=0,017 (Па).

Площадь антенны

S=a×b.

a=(50,5×с)/f_опт×Q_a=(50,5×1483)/39000×10=0,192 (м),

b=(50,5×с)/f_опт×Q_b=(50,5×1483)/39000×10=0,192 (м).

S=0,192 ×0,192 =0,037 (м²).

Где Q_a,Q_a - разрешающая способность по угловым координатам.

Интенсивность

I=| P’_ш /r×c |=0,017/10³×1483=1,127×10^-8,

где r - плотность среды распространения звука (вода),

с – скорость звука в среде.

Среднеквадратичное напряжение шума

W_ш=I×S=1,127×10^-8×0,037=4,157×10^-10.

Спектральная плотность мощности шумовой помехи

N_o= W_ш/Df=4,157×10^-10/4136=1,005×10^-13(Вт/Гц).

4. РАСЧЕТ ХАРАКТЕРИСТИК ОБНАРУЖЕНИЯ

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

Определим порог при двух заданных значениях вероятности ложной тревоги P_лт1=10^-3, P_лт2=10^-5 для трех случаев:

а) сигнал известен точно

Распределение помехи нормальное. При определении порога пользуемся таблицей интеграла вероятности

P_лт=1-Ф(q_o), тогда q_o=arg[Ф(1- P_лт)].

Из таблицы интеграла вероятности для:

P_лт1=10^-3, q_o=3,1;

P_лт2=10^-5, q_o=4,27.

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

P_лт=1-Ф(q_o-q).

Таблица 4.1

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

б) Сигнал со случайной начальной фазой

Распределение помехи релеевское, но при больших отношениях «сигнал-шум» распределение сводится к нормальному

q_o=Ö-2×ln(P_лт).

Тогда для P_лт1=10^-3, q_o=3,72;

P_лт2=10^-5, q_o=4,8.

Таблица 4.2

Точки построения кривой обнаружения для сигнала с неизвестной начальной фазой*

в) сигнал со случайной фазой и амплитудой

q_o=Ö-2×ln(P_лт).

Тогда для P_лт1=10^-3, q_o=3,72;

P_лт2=10^-5, q_o=4,8.

Расчет точек для кривой обнаружения.

Таблица 4.3

Точки построения кривой обнаружения для сигнала с неизвестной начальной фазой и амплитудой*

4.2. Расчет характеристик обнаружения

а) Находим энергию сигнала при P_omin=0,92

тогда

Данные наших расчетов приведены в приложении (рис.1) и (рис.2).

Таблица 3.4

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

б) энергия минимального сигнала при когерентном и некогерентном приеме.

Е_и=E_s/n –для когерентного приема.

Е_и=E_s/Ön – для некогерентного приема.

n=1 и n=20 – число сигналов принимаемой последовательности .

Для n=1 различие между когерентным и некогерентным приемами отсутствуют.

Таблица 4.5

Энергия минимального порогового сигнала

в) коэффициент распознавания

d=q_оп/Ön – для когерентного приема.

d=q_оп/⁴Ön – для когерентного приема.

Таблица 4.6

Коэффициент распознавания, d

г) импульсная мощность

W_и=E_s/t_и, для n=1;

W_и=E_и/t_и, для n=20.

Таблица 4.7

Импульсная мощность W_и, Вт

ВЫВОД

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

ПРИЛОЖЕНИЕ

* см. приложение (рис.1 и рис.2)

* см. приложение (рис.1 и рис.2)

* см. приложение (рис.1 и рис.2)