Курсовая работа3

5

Федеральное агентство по образованию РФ

Санкт – Петербургский Государственный электротехнический

университет «ЛЭТИ» им. В. И. Ульянова (Ленина)

Кафедра ЭУТ

Пояснительная записка

к курсовой работе по дисциплине:

«Методы анализа и обработки сигналов»

Вариант№13

Выполнила студентка: Николенко С.

Группа № 5586

Преподаватель: Коновалов

Санкт – Петербург

2008

Содержание

3

4

4

4

7

8

10

10

13

18

20

21

22

Задание…………………………………………………………………………………….

1. Оценка статистических характеристик случайного процесса…………………………

1.1. Нахожднние математического ожидания и дисперсии …………………………

1.2. Нахождение корреляционной матрицы ………………………………………….

1.3. Проверка стационарности случайного процесса

в широком смысле…………………………………………………………….

1.4. Нахождение нормированной корреляционной функции………………….

2. Определение структуры согласованного и квазиоптимального

фильтра………………………………………………………………………………

2.1. Построение согласованного фильтра……………………………………..

2.2. Нахождение параметров квазиоптимального фильтра…………………………

3. Определение характеристик обнаружения………………………………………….

Выводы………………………………………………………………………………….

Литература……………………………………………………………………………..

Приложения…………………………………………………………………………….

ЗАДАНИЕ N 13

1.Случайная функция X(t) задана 12 реализациями x_i(t) в 21 сечении. Значения реализаций с шагом 0,1 сек заданы в файле в виде матрицы. Вычислить математическое ожидание, дисперсию и корреляционную матрицу случайной функции, проверить, является ли функция X(t) стационарной, и в последнем случае определить ее нормированную корреляционную функцию.

2.На вход приемного устройства поступает сигнал

x(t)=s(t)+n(t), где

s(t) = A cos²(t/2_) cos(_t+_), t  [-_; _]

A - случайная амплитуда, распределенная по закону Рэлея:

,

₀ = 2 мсек.; ₀- случайная начальная фаза, распределенная по закону:

n(t) - квазибелый гауссовский шум, имеющий спектральную плотность:

S()=N₀/2

в полосе || = ₂ - ₁, полностью перекрывающей спектр сигнала.

₀=2f₀; f₀ = 10⁵ Гц; || = 2*2*10⁵

Требуется определить:

А.Структуру согласованного фильтра и параметры (ширина полосы пропускания и изменение отношения сигнал/помеха на выходе по сравнению со входом) квазиоптимального фильтра, состоящего из 3 колебательных контуров.

В.Зависимость P_D, где^_s²/_n² на входе приемного тракта, если обнаружитель выполнен по схеме согласованный фильтр - линейный детектор - пороговое устройство. При этом с доверительной вероятностью P=0.9 должно быть не более n₀ = 0 ложного срабатывания регистратора за 8 часов непрерывной работы при частоте посылок 10 имп/с.

1. Оценка статистических характеристик случайного процесса.

1.1. Нахождение математического ожидания и дисперсии.

Математическое ожидание и дисперсию случайного процесса определим по следующим формулам:

, где n – количество реализаций, k- номер реализаций,

i, j –номера сечений, - оценка математического ожидания в i – ом сечений, - оценка дисперсии в i - ом сечений.

1.2. Нахождение корреляционной матрицы.

Корреляционная матрица случайной функции может быть вычислена по формуле:

, где - оценка корреляционной функций для сечений и . Результаты вычисления корреляционной функции представлены в таблице 1.

Таблица 1

1.3. Проверка стационарности случайного процесса в широком смысле.

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

1. Проверим по оценкам :

- условие стационарности выполняется.

2. Проверим по оценке :

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

1.4 Нахождение нормированной корреляционной функции.

Нормированная корреляционная функция вычисляется по формуле: . Матрица представлена в таблице 2.

Таблица2

2. Определение структуры согласованного и квазиоптимального фильтра.

2.1. Построение согласованного фильтра.

Найдём выражение для спектра заданного сигнала S(t).( =0)

В общем виде:

Где

Получаем:

Частотная характеристика согласованного фильтра определяется выражением:

,

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

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

Вычисление произведём в Mathcad:

Интеграл сходится, следовательно, фильтр физически реализуем.

Построим структурную схему фильтра:

2.2. Нахождение параметров квазиоптимального фильтра.

Зачастую построение оптимального (согласованного) фильтра на практике невозможно, в таком случае используют квазиоптимальные фильтры, т.е. практически осуществимые, на выходе которых отношение сигнал\помеха лишь немного меньше значения, рассчитываемого для согласованного фильтра. Тогда ухудшение отношения сигнал\помеха на выходе квазиоптимального фильтра по сравнению с этим отношением на выходе оптимального фильтра равно:

.

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

, где n=3 –число контуров. Путём замены , получаем

, где Q – добротность колебательного контура.

Преобразуя функцию r, получим числитель:

Верхний предел интегрирования y выберем из условия, что в полосе частот содержится 90% энергии сигнала:

, где Δω – полоса частот шума.

Расчет произведем в пакете MathCad 2001.

Значение полосы пропускания фильтра, состоящего из трёх контуров: ,

где .

Найдём (полуширину спектра сигнала) и (полуширину полосы пропускания фильтра), используя следующие выражения:

Отношение сигнал-помеха на выходе квазиоптимального фильтра выражается через отношение сигнал-помеха на его входе:

, где .

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

Обнаружитель выполнен по схеме согласованный фильтр - линейный детектор - пороговое устройство. Вероятность ложной тревоги определяется:, где в - среднее количество ложных регистраций, N₀ – число независимых точек контроля.

находится из уравнения . В данном случае n₀ = 0. Вероятность превышения порога для отдельного импульса при наличии сигнала:. Проведём расчёт в MathCAD и построим кривую обнаружения.

Вывод:

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

Литература:

1. Д.Д. Добротин, С.И. Коновалов Методические указания к выполнению курсовой работы по дисциплине “Методы анализа и обработки сигналов”. СПб.: Изд-во СПбГЭТУ “ЛЭТИ”, 2008.

2. Добротин Д.Д., Паврос С.К. Методические указания “Методы анализа случайных сигналов”. СПб.: Изд-во СПбГЭТУ “ЛЭТИ” 2006.

3. Добротин Д.Д., Пигулевский Е.Д. Случайные сигналы и помехи в системах интроскопии Мет. указания. ЛЭТИ, Ленинград, 1990. 80с.

4. Гоноровский И.С., Демин М.П. Радиотехнические цепи и сигналы. – Учеб. Пособие М.: Радио и связь, 1994.

Приложение