Курсовая / МЕТОДА
.pdfМИНОБРНАУКИ РОССИИ
_________________________________________
Санкт-Петербургский государственный электротехнический университет «ЛЭТИ»
___________________________________________
СОГЛАСОВАННАЯ ФИЛЬТРАЦИЯ СИГНАЛОВ
Методические указания к курсовой работе по дисциплине
«РАДИОТЕХНИЧЕСКИЕ ЦЕПИ И СИГНАЛЫ»
Санкт-Петербург 2022
|
ОГЛАВЛЕНИЕ |
|
Введение................................................................................................................... |
3 |
|
1. |
Задание............................................................................................................... |
4 |
2. |
Основные определения и комментарии.......................................................... |
6 |
3. |
Порядок выполнения курсовой работы........................................................ |
10 |
|
3.1. Расчет спектра входного сигнала ............................................................. |
10 |
|
3.2. Разработка структурной схемы согласованного фильтра...................... |
11 |
|
3.3. Расчет преобразования сигнала и шума согласованным фильтром ..... |
13 |
|
3.4. Расчет преобразования сигнала и шума квазиоптимальным RС- |
|
|
фильтром нижних частот.................................................................................. |
17 |
4. |
Содержание пояснительной записки ............................................................ |
20 |
5. |
Требования к оформлению курсовых расчетов по курсу РТЦС................ |
22 |
2
Введение
Вприемниках, предназначенных для обнаружения, а также разрешения
иоценки параметров сигналов, форма которых заранее известна (например, в приемниках импульсных радиолокационных станций), одной из важнейших составных частей является устройство, максимизирующее отношение сигнал/шум (с/ш) на выходе приемника. Такие устройства реализуются либо в виде корреляторов, либо в виде линейных фильтров, называемых согласованными [l], [2], коэффициент передачи которых определяется спектральными характеристиками обнаруживаемого входного сигнала, действующего на фоне шума.
Широкое применение в современных радиотехнических системах сложных сигналов и их корреляционной обработки требует от радиоинженера хорошего владения теорией согласованной фильтрации и практических навыков анализа и синтеза согласованных фильтров.
Цель настоящей курсовой работы — углубление представления студентов о методах и возможностях согласованной фильтрации сигналов, а также приобретение ими практических навыков в процессе решения конкретной задачи, включающей элементы синтеза и анализа согласованного фильтра.
3
1.Задание
1.В качестве входного сигнала в курсовой работе рассматривается видеоимпульс u(t) , заданный на интервале [0,T2 ] двумя отрезками прямых, как
показано на Рис. 1.1.
В каждом индивидуальном задании заданы три значения параметра T1. После ана-
лиза амплитудных спектров сигнала с разными T1 по заданному критерию выбирается
одно из значений.
2. В качестве шумового сигнала в курсовой работе рассматривается стационарный случайный процесс в виде квазибелого шума, спектральная плотность мощности W0 которого ограничена частотой fш.гр,
ω ≤ 2πfш.гр.
u(t) |
|
U4 |
|
U3 |
|
U2 |
|
U1 |
|
0 |
Т1 Т2 t |
|
Рис. 1.1 |
3. Задание на курсовую работу заключается в следующем:
а) рассчитав амплитудные спектральные функции для трех вариантов заданного входного сигнала, выбрать по указанному далее критерию один из них и, рассчитав для выбранного сигнала комплексную спектральную функцию, ее вещественную и мнимую части, а также фазовый спектр, определить комплексный коэффициент передачи согласованного с ним фильтра (СФ);
б) синтезировать структурную схему СФ и рассчитать форму напряжений в различных ее точках при подаче на вход фильтра импульса напряжения, имеющего вид дельта-функции δ(t ) ;
в) рассчитать форму выходного напряжения при воздействии на вход СФ выбранного сигнала, а также прямоугольного видеоимпульса, удельная энергия которого равна удельной энергии выбранного входного сигнала, а
длительность Т = 1 f0,1 , где |
f0,1 — частота, на которой уровень амплитудного |
||||||||
спектра |
|
Sɺ(ω) |
|
выбранного |
входного видеосигнала снижается до уровня |
||||
|
|
||||||||
0,1 |
|
Sɺ(ω) |
|
|
|
= 0,1S (0); |
|
||
|
|
|
|
|
|||||
|
|
|
|
max |
|
||||
|
|
|
|
|
|||||
|
|
г) рассчитать спектр мощности шума на выходе СФ; |
|||||||
д) проанализировать эффективность действия СФ, рассчитав выигрыш в отношении с/ш, обеспечиваемый им для сигналов, рассмотренных в п. 3в; е) исследовать прохождение выбранного сигнала и шума через квази-
оптимальный фильтр в виде RC-фильтра нижних частот, а именно:
- найти значение постоянной времени RC-фильтра τopt , при котором отношение с/ш на его выходе максимально (оптимальное значение);
4
-определить выигрыш в отношении с/ш, который дает квазиоптимальный фильтр, и сравнить его со значением, полученным для согласованного фильтра;
-рассчитать амплитудно-частотную и фазочастотную характеристики (АЧХ и ФЧХ) квазиоптимального фильтра с τopt , а также спектр мощности
шума на его выходе.
Задание может быть представлено в виде показанной здесь таблицы.
U1, В |
U2, В |
U3, В |
U4, В |
T1, мкс |
n |
m |
T2, мкс |
|
|
|
|
|
|
|
|
Параметры шумового сигнала определяют по формулам:
W = |
U 2 |
T |
ω |
|
= |
10π m |
|||
max 2 |
; |
|
|
|
. |
||||
|
ш.гр |
|
|||||||
0 |
n |
|
|
|
|
|
T2 |
||
|
|
|
|
|
|
|
|||
5
2. Основные определения и комментарии
При описании детерминированного сигнала важную роль играет спектральная функция Sɺ(ω) = Sɺ(ω) exp{jϕ(ω)}= S (ω)exp{jϕ(ω)}, связанная с сигна-
лом u(t ) парой преобразований Фурье:
. |
|
∞ |
|
|
|
|
S (ω ) = |
u |
(t )e− jωt dt , |
(2.1) |
|||
|
|
−∞ |
|
|
|
|
|
1 |
∞ . |
|
|
||
u (t ) = |
|
S (ω ) e jωt dω . |
(2.2) |
|||
2π |
||||||
|
−∞ |
|
|
|
||
|
|
|
|
|
||
Как следует из (2.1), спектральная функция имеет размерность В/Гц или В×с, если u(t ) — напряжение.
Модуль спектральной функции Sɺ(w) = S (w) и ее аргумент j(w) соответ-
ственно называют амплитудным и фазовым спектрами сигнала. С помощью спектральной функции можно определять энергетические характеристики сигналов. В теории сигналов под энергией сигнала подразумевается энергия, выделяемая сигналом на единичном сопротивлении, т. е. удельная энергия:
∞ |
(t ) dt = |
1 |
∞ |
S 2 (ω ) dω = |
1 |
∞ |
(ω ) dω . |
|
|
E = u 2 |
|
S 2 |
(2.3) |
||||||
2π |
π |
||||||||
−∞ |
|
−∞ |
|
0 |
|
|
|||
|
|
|
|
|
|
Функция S 2 (w) в (2.3) характеризует распределение удельной энергии сигнала по частоте и называется энергетическим спектром или спектральной плотностью энергии детерминированного сигнала, размерность которой есть
В2 × с/Гц.
Важную роль в описании детерминированного сигнала играет также его корреляционная функция (КФ), определяемая следующим образом:
|
|
|
∞ |
|
|
|
B (τ ) = u(t)u(t −τ )dt . |
(2.4) |
|||
|
|
|
−∞ |
|
|
Корреляционная функция связана с энергетическим спектром детерми- |
|||||
нированного сигнала парой преобразований Фурье: |
|
||||
|
|
|
|
∞ |
|
|
Sɺ(ω ) |
2 = S 2 (ω ) = B (τ ) e− jωτ dτ , |
(2.5) |
||
|
|
|
|
−∞ |
|
|
|
|
1 |
∞ |
|
|
B (τ ) = |
S 2 (ω ) e jωτ dω . |
(2.6) |
||
|
2π |
||||
|
|
|
−∞ |
|
|
|
|
|
|
|
|
КФ является четной функцией, B(τ) = B(− τ), и при t = 0 имеет максимум, причем в соответствии с (2.3) В(0) = Е. Как следует из (2.4) — (2.6), раз-
мерность КФ — В2 × c.
Для характеристики случайных сигналов в частотной области вводят понятие спектральной плотности мощности (СПМ) или спектра мощности:
6
|
W (ω ) = lim |
M {ST2 (ω )} |
, |
(2.7) |
|
|
|
||||
|
T →∞ |
T |
(w) усеченной (ограниченной |
||
где S |
(ω) — модуль спектральной функции Sɺ |
||||
T |
|
T |
|
|
|
по длительности интервалом времени Т) реализации случайного сигнала, M{..} – операция усреднения по ансамблю реализаций.
Как видно из выражения (2.7), спектр мощности (нередко, в связи с тем,
что размерность W (ω) есть В2 × c или В2 /Гц, т. е. удельная энергия, и здесь применительно к W (ω) используют термин «энергетический спектр», что не вполне корректно) случайного сигнала характеризует распределение по частоте его средней удельной мощности в отличие от энергетического спектра детерминированного сигнала (2.5).
Для случайных сигналов можно записать соотношение, подобное выражению (2.3), определяющее среднюю удельную мощность их переменной составляющей (дисперсию), следующим образом:
|
1 |
∞ |
|
1 |
∞ |
|
|
σ 2 = |
|
W (ω ) dω = |
W (ω ) dω . |
(2.8) |
|||
2π |
π |
||||||
|
−∞ |
|
0 |
|
|||
|
|
|
|
|
Случайный сигнал x(t) также характеризуется корреляционной функцией, которая для эргодических случайных сигналов с нулевым средним значением определяется следующим образом:
|
1 |
T / 2 |
|
|
R (τ ) = lim |
x (t ) x (t −τ ) dt . |
(2.9) |
||
T |
||||
T →∞ |
−T / 2 |
|
КФ случайного сигнала имеет размерность В2 . Спектр мощности (2.7) и корреляционная функция (2.9) связаны друг с другом парой преобразований Фурье:
|
∞ |
|
|
||
W (ω ) |
= R (τ ) e− jωτ dτ , |
(2.10) |
|||
|
−∞ |
|
|||
|
1 |
|
∞ |
|
|
R (τ ) = |
|
W (ω ) jωτ dω . |
(2.11) |
||
2π |
|||||
|
−∞ |
|
|||
|
|
|
|
||
Формулы (2.10) и (2.11) составляют содержание теоремы Винера—Хин- чина. КФ случайного сигнала, как и КФ сигнала детерминированного, обладает свойством четности, R(τ) = R(− τ), и имеет при t = 0 максимум, причем,
согласно (2.8) и (2.11), R(0) = s2.
Форма сигнала uвых (t ) , прошедшего через линейную цепь, определяется или с помощью интеграла Дюамеля (свертки):
∞ |
|
uвых (t ) = u(λ)h(t − λ )d λ , |
(2.12) |
−∞
где u(t ) — сигнал на входе, a h(t) — импульсная реакция цепи, или спектральным методом с помощью преобразования Фурье (2.2):
7
|
|
1 |
∞ . |
. |
|
||
uвых |
(t ) = |
|
S (ω) K (ω)e jωt dω . |
(2.13) |
|||
2π |
|||||||
|
|
|
−∞ |
|
|
|
|
Здесь Sɺ(ω) — спектральная функция входного детерминированного сигнала, Kɺ(ω) — комплексный коэффициент передачи цепи, связанный с ее импульсной реакцией преобразованием Фурье.
При прохождении случайного сигнала через линейную цепь энергетический спектр и КФ выходного сигнала определяются следующим образом:
W (ω) = W (ω) |
|
Kɺ |
(ω) |
|
2 , |
(2.14) |
|
|
|||||
вых |
|
|
|
|
|
|
∞ |
|
|
|
|
||
Rвых (τ ) = R(λ)Bh (τ − λ )d λ , |
(2.15) |
|||||
−∞
где W (w) и R(t) — соответствующие характеристики входного случайного сигнала; Bh — корреляционная функция импульсной реакции цепи:
∞
Bh (τ) = h(t)h(t − τ)dt . −∞
Частотно-избирательную систему, выполняющую обработку аддитивной смеси (суммы) сигнала и шума наилучшим образом в соответствии с тем или иным критерием, называют оптимальным линейным фильтром.
В задачах обнаружения известного по форме сигнала в качестве такого критерия используют критерий максимума отношения с/ш на выходе. Если помехой является стационарный случайный сигнал в виде белого шума, то комплексный коэффициент передачи фильтра, максимизирующего отношение с/ш, определяется следующим образом:
. |
|
K s (ω) = ASɺ (ω) e− jω t0 , |
(2.16) |
где Sɺ(ω) — спектральная функция входного сигнала; А — константа (коэффициент), имеющая размерность, обратную размерности спектральной плотности; t0 — момент времени, в который отношение с/ш на выходе фильтра до-
стигает максимума, * - знак комплексного сопряжения.
Фильтр, обладающий частотной характеристикой вида (2.16), называют согласованным. Необходимым (но недостаточным) условием физической реализуемости согласованного фильтра является неравенство t0 ³ T , где Т —
длительность входного детерминированного (полезного) сигнала. Используя (2.16), можно показать, что импульсная реакция согласованного фильтра имеет следующий вид:
hs (t) = Au(t0 − t ) , |
(2.17) |
где u(t0 - t ) есть зеркально отображенная относительно оси ординат и сдвинутая по оси времени на интервал t0 копия входного сигнала u(t ).
8
При подаче сигнала u(t ) на вход согласованного с ним фильтра на выходе фильтра появится определяемый из (2.5), (2.13) и (2.16) или из (2.12) и
(2.17) сигнал |
|
uвых (t ) = AB(t − t0 ) , |
(2.18) |
где B(τ) — корреляционная функция входного сигнала. Таким образом, отклик фильтра на согласованный с ним сигнал совпадает по форме с КФ последнего. Как следует из (2.18), вне зависимости от формы входного сигнала максимум сигнала на выходе uвых(t0 ) = AE , где E — энергия входного сигнала u(t ).
По формулам (2.13) и (2.16) или (2.12) и (2.17) можно определить форму сигнала на выходе согласованного фильтра и для произвольного входного сигнала:
uвых (t ) = ABv (t − t0 ) , |
(2.19) |
где |
|
∞ |
|
Bv (τ ) = uv (t)u(t −τ ) dt |
(2.20) |
−∞
— взаимная корреляционная функция (ВКФ) произвольного входного сигнала uv (t ) и сигнала u(t ), с которым фильтр согласован.
Спектр мощности шума (помехи) на выходе согласованного фильтра при воздействии на его вход белого шума со спектром мощности W0 может
быть найден из (2.14) или из (2.10) и (2.15): |
|
|
|
|
|
|||||||
Wвых (ω) = W0 |
|
Kɺs (ω ) |
|
2 = W0 Ks2 (ω ) . |
(2.21) |
|||||||
|
|
|||||||||||
Отношение с/ш на выходе согласованного фильтра в момент t0 из (2.8), |
||||||||||||
(2.18) и (2.21) определяется следующим образом: |
|
|||||||||||
|
|
uвых (t0 ) |
|
|
|
|
|
|
||||
(с/ш)вых |
= |
= |
|
E |
, |
(2.22) |
||||||
|
|
|||||||||||
|
|
|
|
σ вых |
W0 |
|
||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
σ |
2 |
|
|
= A2W E |
|
|
|
|
(2.23) |
|||
|
вых |
0 |
|
|
|
|
|
|||||
— дисперсия шума на выходе фильтра.
Соотношение (2.22) подчеркивает тот замечательный факт, что максимальное значение отношения с/ш на выходе фильтра, согласованного с сигналом, определяется только энергией последнего и спектром мощности помехи и не зависит от формы этого сигнала.
9
3. Порядок выполнения курсовой работы
При выполнении курсовой работы предусмотрено использование ПК как для вычислений, так и для проведения математического моделирования и исследований в соответствии с индивидуальным заданием. Выбор программы для расчетов остается за студентом (см. п.9 раздела 5).
3.1.Расчет спектра входного сигнала
Спомощью выбранной программы для ПК рассчитайте амплитудные спектры заданных вариантов входного сигнала. Расчет произведите в диапазоне, который ориентировочно можно оценить, исходя из того, что ширина главного лепестка спектра прямоугольного видеоимпульса длительности T2
равна 1/T2 ; поэтому целесообразно задать частотный интервал ~ 3/T2 . Резуль-
таты расчета сведите в таблицу и постройте соответствующие графики.
По графику определите верхние граничные частоты f0,1 спектров задан-
ных входных сигналов по критерию уменьшения их амплитудных спектров до уровня 0,1 от максимального значения амплитудного спектра S (0). Правее этой точки все значения амплитудного спектра меньше указанной величины. Для этого ориентировочно оцените частотные интервалы, в которых находятся интересующие вас частоты. Затем определите точные значения соответствующих граничных частот с помощью программы математического моделирования.
Выберите вариант сигнала, которому отвечает наибольшее (max), наименьшее (min) или среднее (mid) значение граничной частоты (в соответствии с индивидуальным заданием). В дальнейшем рассматривается только этот вариант сигнала.
Выведите аналитическое выражение для спектральной функции выбранного входного сигнала. Возможно прямое вычисление интеграла Фурье, но в данном случае целесообразно использовать известный метод дифференцирования (здесь — двукратного) кусочно-линейной функции u(t ). Получите расчетные соотношения для амплитудного спектра, вещественной и мнимой частей комплексного спектра и фазового спектра выбранного входного сигнала. Проведите машинный расчет этих функций в пределах f0,1 и постройте соот-
ветствующие графики. Используя полученные аналитические выражения, «вручную» рассчитайте значения амплитудного и фазового спектров в трех контрольных точках (на нулевой и двух произвольно выбранных в пределах f0,1 частотах), задав их таким образом, чтобы обеспечить возможность сопо-
ставления результатов с данными, полученными при машинном расчете. При
10
