Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO
3.2. Увеличение временного разрешения сигналовМы рассмотрели фильтрацию узкополосной помехи. Теперь на основе работы [23] рассмотрим фильтрацию уже не помехи, а сигнала. Цель такой фильтрации, примененной к сигналу, состоит в изменении его комплексного спектра, которое придает ему свойства, полезные для использования этого сигнала при локации в качестве зондирующего. Прежде всего мы рассмотрим задачу, в которой все параметры будут определяться нами численно. На этом этапе можно наиболее эффективно показать задачи и возможности их решения, а затем рассмотреть и эксперимент, который относится к локации на просвет. И здесь дадим другой подход к ее успешному решению, основанный не на фильтрации помехи, а на фильтрации зондирующего сигнала. На заре развития радиолокации в качестве зондирующего сигнала использовался прямоугольный импульс с монохроматическим заполнением. Разрешающая способность локатора по времени тогда определялась длительностью импульса: чем короче импульс, тем выше разрешающая способность локатора. На следующем этапе развития локации разрешение отраженных сигналов во времени стало получаться существенно меньше длительности импульса. Чтобы увеличить разрешение, зондирующий сигнал модулировался внутри импульса по частоте, что расширяло его спектр, увеличивая временное разрешение. Величина интервала временного разрешения лоцируемых объектов равна обратной величине ширины спектра зондирующего импульса [24]. Это справедливо как для импульса с монохроматическим заполнением, который для расширения полосы должен быть достаточно коротким, так и для импульса с широкополосным заполнением. Универсальное соотношение между интервалом временного разрешения
Соотношение (3.2.1) непосредственно следует из математического определения спектра на основе преобразования Фурье и поэтому нерушимо. В рамках этого соотношения возможна еще одна модификация зондирующего сигнала. Речь идет о зондирующем сигнале, имеющем вид затухающего синусоидального сигнала (ЗСС). Типичным примером ЗСС является отклик колебательного контура на импульсное воздействие. ЗСС справедливо считается узкополосным, так как медленная по сравнению с периодом высокой частоты модуляция амплитуды слабо расширяет его спектр. Если не иметь в виду специальную математическую обработку ЗСС, а ограничиться только теми видами обработки, которые применяются для сигналов с монохроматическим заполнением и широкополосным заполнением, то под шириной спектра такого сигнала Асов формуле (3.2.1) следует понимать эффективную ширину его спектра по уровню 3 дБ. Эта ширина определяется скоростью спада уровня сигнала и при достаточно медленном спаде высокого временного разрешения не обещает. Полная ширина спектра ЗСС по его самому низкому уровню теоретически бесконечна, а практически (как это следует из теоремы отсчетов) определяется той частотой, с которой этот сигнал проквантован при вводе в вычислительное устройство. Среди методов выделения сигналов существуют такие, для которых в соотношение (3.2.1) входит не эффективная ширина спектра зондирующего сигнала, а его полная ширина. В результате применения этих методов может быть получено высокое разрешение при использовании зондирующего сигнала, считающегося узкополосным. Ниже рассматриваются и сравниваются между собой два способа обработки принятого локационного сигнала, для которых в соотношении (3.2.1) работает полная ширина спектра зондирующего сигнала. Чтобы привязать известные методы обработки сигналов к задаче обработки локационных сигналов, введем ряд определений. Условимся зондирующий сигнал локатора считать откликом на импульсное воздействие некоторого фильтра. Пусть вид зондирующего сигнала будет
В соотношении (3.2.2) известны
где На рис. 3.8 представлен результат математического моделирования, выполненного по формуле (3.2.3). Вдоль горизонтальной оси отложено время в единицах квантования. По вертикали отложены для (0 - его уровень в децибелах, а для
Рис. 3.8. Отклик локатора. Вверху в линейном масштабе показан отклик локатора на узкополосный зондирующий импульс, внизу - результат обработки этого отклика путем расширения его спектра (дБ). линейная шкала. В качестве зондирующего сигнала была взята экспоненциально затухающая косинусоида - отклик резонансного контура на импульсное воздействие. На рисунке вверху показан сигнал Покажем, что зондирующий сигнал вида ЗСС является действительно узкополосным сигналом при существующих традиционных способах обработки принятого сигнала [24]. Будем это доказывать от противного. Пусть сигнал ЗСС в силу того, что он модулирован по амплитуде, является широкополосным сигналом. Основанием для такого допущения является то, что спектр ЗСС не обращается в нуль ни в одной точке всего рассматриваемого спектрального интервала. Обработаем его так, как принято обрабатывать широкополосные сигналы [24], а именно - путем определения взаимной корреляции зондирующего сигнала и принятого сигнала по формуле
(кликните для просмотра скана) На рис. 3.10 изображен результат обработки сигнала по формуле (3.2.4) в том случае, когда в качестве зондирующего использовался шумоподобный сигнал, имеющий равномерный спектр в пределах всего используемого спектрального интервала. Вдоль горизонтальной оси отложено время Соотношение (3.2.2) можно использовать иначе, не прибегая к операции деления спектра на спектр. Как принято поступать при комплексном кепстральном анализе [25], прологарифмируем (3.2.3):
Используемый здесь логарифм от комплексной функции
Рис. 3.10. Отклик локатора, обладающего длительным, широкополосным шумовым зондирующим сигналом, получен путем свертки сигнала, принятого локатором, с зондирующим сигналом. Моделировались две цели.
Рис. 3.11. Отклик локатора, обладающего узкополосным зондирующим сигналом, полученный путем свертки сигнала, принятого локатором, с зондирующим сигналом. Моделировались две цели.
Рис. 3.12. Отклик локатора, обладающего узкополосным зондирующим сигналом, полученный путем расширения спектра сигнала. Моделировались две цели. Далее, при кепстральном анализе полагается использовать различия в спектрах (кепстрах) двух слагаемых (3.2.5) с целью их разделения путем фильтрации [25]. В нашем случае разные спектры (кепстры) должны иметь логарифмы спектров отклика и принятого сигнала локатора. Это можно обеспечить специальными приемами. Основным из них является выбор такого значения аргумента в функции (3.2.6), которое обеспечивает постоянное нахождение на одном листе этой многолистной многозначной периодической функции [25]. Для этого надо в соответствующие моменты достижения этой функцией границ добавлять к ней или убавлять от нее постоянное слагаемое величиной в Сдвиг кепстра по частоте объясняется свойством преобразования Фурье превращать сдвиг функции в смещение ее частоты. Кепстр приобретает вид сигнала В результате кепстр суммы сигналов имеет вид, показанный на рис. 3.13, где приведены один под другим два кепстра, отличающиеся скоростью убывания со временем амплитуды зондирующего сигнала. Вдоль горизонтали отложена частота (сач-тота) (с нулем в центре рисунка). Вдоль вертикали отложен уровень кепстра в децибелах. Кепстр, имеющий более быстрый спад зондирующего сигнала, увеличен на 80 дБ. Низкочастотная часть кепстра содержит кепстр зондирующего сигнала, который в силу того, что аргумент превращен в однолистную функцию, максимально узок. Выше по частоте (сачтоте), начиная с частоты (сачтоты), соответствующей задержке сигнала, виден кепстр входного сигнала, представляющего спектр от его спектра, равный
Рис. 3.13. Кепстр отклика локатора, обладающего узкополосным зондирующим сигналом Две кривые отличаются шириной спектра зондирующего сигнала, одна из кривых смещена вверх
Рис. 3.14. Отклик локатора, обладающего узкополосным зондирующим сигналом Вверху - отклик без расширения спектра Внизу - отклик, кепстр которого выделен путем фильтрации из кепстра, показанного на рис 3.13 Полученный таким образом результат не отличается от результата, полученного путем фильтрации спектров ни по виду сигнала, ни по отношению сигнала к шуму Отличие заключается только в самой процедуре получения исходного локационного отклика. Принципиальным отличием использованных процедур является то, что деление спектра требует для своего осуществления отдельного канала, содержащего функцию Существует весьма близкая аналогия между рассмотренным способом нахождения исходного локационного сигнала и теми преобразованиями полей, которые осуществляются при формировании изображений [3]. Для описания изображения в монохроматическом волновом поле воспользуемся так называемой комплексной амплитудой волнового поля задать на плоскости
где
Здесь X - длина волны. Формула (3.2.7) полностью эквивалентна соотношению (3.2.2). В силу этого она может быть использована для нахождения одного поля по известному другому. Это облегчается тем, что частотная характеристика свободного пространства, определяемая (3.2.8), известна. Примеры применения алгоритма (3.2.3) для нахождения волновых полей приведены в [3, 6]. Особенностью, отличающей задачу распространения волн от локационной, является не только то, что частотная характеристика свободного пространства является заданной функцией и всегда известна, но и вид этой функции (3.2.8). У нее модуль равен единице для всех распространяющихся волн, и спектр ее расширяется за счет фазы. В ЗСС дело обстоит как раз наоборот: фаза практически постоянна (изменение фазы во всем спектре укладывается в интервал Очень существенным для использованного преобразования является изменение соотношения сигнала к шуму. Рассмотрим это подробнее. Эффект сжатия достигается благодаря выравниванию спектра сигнала по всей его ширине вплоть до частоты временного квантования. При этом длительность сжатого сигнала составляет один интервал временного квантования, в то время как длительность первоначального сигнала может составлять сотни и тысячи таких точек. При "выпрямлении" спектра области малых спектральных амплитуд, которые присутствуют в основном в спектрах рассматриваемых узкополосных сигналов, приходится умножать на относительно большие множители. Эта операция происходит в присутствии шумов, которые тоже умножаются на те же множители. Возникает вопрос о степени помехоустойчивости такой процедуры. Начнем с уточнения понятия помехоустойчивости. От этого зависит знак полученного результата. Целесообразно дать определение помехоустойчивости, специфическое для целей локации. Определение отношения сигнала к шуму оперирует со статистическими характеристиками как сигнала, так и шума. Основной статистической характеристикой является так называемое стандартное отклонение, определяемое как
где
Выражение (3.2.10) определяет обычное отношение сигнала к шуму. В локации отношение (3.2.10) много меньше единицы, так как узкий импульс, применяемый в локации, содержит очень мало энергии. Однако в этом случае сигнал отчетливо виден на экране локатора, если в наблюдаемом нами импульсе его максимальное значение
Отношение (3.2.11) положим в основу наших расчетов. Результат расчета в первую очередь зависит от того, на каком этапе обработки сигнала поступает шум, и от того, какой из этих этапов берется для сравнения. Есть зависимость и от формы импульса. Для начала рассмотрим такой случай, когда шум подается туда же, где помещается объект локации. На практике это соответствует так называемой реверберационной помехе, которая образуется вследствие рассеяния зондирующего импульса на случайных неоднородностях среды (каплях дождя, кустах, облаке пассивных помех и т. Процедура обработки сигнала с целью его сжатия заключается в том, что каждая мгновенная реализация комплексного фурье-спектра принятого сигнала вместе со спектром шума делится на комплексный спектр зондирующего сигнала Пусть форма зондирующего импульса до сжатия имеет вид затухающей косинусоиды. Такую форму принимает отклик колебательного контура на короткий импульс. Для нас существенно, что максимум отклика контура не будет изменяться в зависимости от добротности контура или от скорости затухания в нем свободных колебаний. Под действием короткого возбуждения контур будет совершать колебания от той амплитуды, которой он успел достичь за время толчка. Шум при сделанном нами предположении представляет собою толчки столь же малой длительности, но следующие хаотично один за другим. Примем, что средняя амплитуда одного толчка шума равна величине толчка бтах, который на контур оказывает сигнал. Это равносильно предположению, что отношение В, определяемое (3.2.11) в результате сжатия сигнала становится равным единице. Толчки от шума суммируются контуром в течение длительности его отклика Т. Условимся измерять Т числом интервалов временного квантования до спада амплитуды зондирующего сигнала на 2 дБ. Величина интервала Т, измеренная в точках отсчета, будет означать и величину сжатия узкополосного сигнала, поскольку он сжимается до одной точки отсчета. Стандартное отклонение колебаний в контуре под действием случайных толчков запишется следующим образом:
На основании (3.2.12) локационное отношение сигнала к шуму, если оно после сжатия сигнала равно 1, на начальном этапе до сжатия сигнала будет равно На рис. 3.15 представлен результат теоретического расчета и численного моделирования. По осям отложено: вдоль горизонтали - сжатие сигнала Т (в логарифмическом масштабе), а вдоль вертикали - увеличение отношения сигнала к шуму В, определяемое (3.2.11), которое произошло в результате сжатия сигнала. Сплошными линиями показаны результаты аппроксимации, выполненной в соответствии с расчетами, а точками отмечены результаты численного моделирования задачи при различных условиях. На рис. 3.15, а показан выигрыш, получаемый при зондирующем импульсе, имеющем форму затухающего косинусоидального сигнала, а на рис. 3.15, б - то же для импульса прямоугольной формы. Рассмотренный выше случай показан на рис. 3.15, а точками в виде косых крестиков. Аппроксимирующая прямая (пунктир) проведена по закону: Точками в виде прямых крестиков на рисунке а показан выигрыш, полученный в результате сжатия импульса после двукратного прохождения через фильтр, имеющий частотную характеристику колебательного контура. В этом случае выигрыш возрастает на 6 дБ. Увеличение выигрыша произошло в результате деформации формы отклика на одиночный начальный импульс. Отклик приобретает форму, при которой отношение максимума отклика к его стандартному отклонению убывает на 6 дБ. Рассмотрим случай, когда шум, имеющий равномерный спектр, поступает в локатор вместе с отраженными сигналами от целей. Если приемный тракт локатора широкополосный, то в процессе обработки сигнала с целью его сжатия принятый сигнал и шум делятся на одну общую для них функцию Рис. 3.15. (см. скан) Выигрыш в отношении сигнала к шуму, получаемый в результате сжатия узкополосного импульса экспоненциальной формы (а) и в результате прямоугольного узкополосного импульса (б): +++ - реверберационный шум, дважды подвергшийся фильтрации; Однако приемный тракт локатора не следует делать широкополосным. Различие в форме спектров сигнала и шума можно использовать для увеличения отношения сигнала к шуму, пропустив принятый локатором сигнал через фильтр с частотной рактеристикой На рис. 3.15, б линия, помеченная прямыми крестиками, проведенная для случая двукратного прохождения сигнала через фильтр, лежит ниже линии, показанной прямыми крестиками. Это объясняется тем, что для прямоугольного импульса, превращающегося в треугольный, отношение максимума к стандартному отклонению не убывает, а возрастает примерно на 5 дБ. Аппроксимирующие прямые на рис. 3.15, б проведены следующим образом (сверху вниз): Особого пояснения требует график, помеченный косыми квадратиками. Он приобретает вид, показанный на рис. 3.15, б, только в том случае, если форма зондирующего сигнала не совсем прямоугольная, т. е. импульс сформирован так, что не обращается в нуль сразу по истечении его длительности, а имеет крутой экспоненциальный спад. Если же оставить форму импульса в точности прямоугольной, то этот график приобретет форму, показанную на рис. 3.16. На этом рисунке вдоль горизонтали отложена длительность импульса в точках отсчета в линейном масштабе. По вертикали отложен выигрыш, получаемый при сжатии прямоугольного сигнала до одной отсчетной точки. Получающаяся на графике зависимость выигрыша от длительности импульса обязана преобразованию величины только шума. Сигнал сжимается без каких-либо особенностей. Причудливая регулярная зависимость выигрыша сжатия от длительности импульса объясняется тем, что в спектре импульса появляются нули, которые попадают (или не попадают) на дискретные точки отсчета спектра. Отметим, что это обстоятельство сказывается только на одной зависимости, все остальные получены при идеально прямоугольной форме импульса.
Рис. 3.16. Выигрыш в отношении сигнала к шуму при сжатии идеально прямоугольного импульса и белого шума на входе локатора. Процедуру сжатия узкополосного сигнала целесообразно использовать не только в том случае, когда получается выигрыш в отношении сигнала к шуму. Для локации несравненно более существенным оказывается выигрыш в разрешающей способности. Иллюстрацией этого может служить рис. 3.17, на котором приведены примеры картин на экране локатора, наблюдавшихся при численном моделировании. Вдоль горизонтали отложено время в точках отсчета, а по вертикали отклик локатора на сигнал и шум относительно максимального уровня сигнала локатора. Показаны отклики сигнала после сжатия при большом (А) и малом (Б) отношениях сигнала к шуму; отклики локатора до сжатия (В, Г) для импульса в виде затухающей косинусоиды при большой длительности (несколько сотен точек) (пунктир) и малой длительности (15 точек) (сплошная линия). Рис. 3.17. (см. скан) Осциллограммы сигналов локатора: А, Б - осциллограммы сжатых импульсов для большого и малого отношений сигнала к шуму; В, Г - осциллограммы до сжатия для тех же отношений сигнала к шуму. Сплошная линия - короткие импульсы, пунктир - длинные. Д, Е - то же, что и В, Г, но для прямоугольных импульсов. (кликните для просмотра скана) (кликните для просмотра скана) На рисунке В (см. рис. 3.17) отношение сигнала к шуму такое, как и на рисунке Заметим, что принятый локационный сигнал может быть обработан так, что применение рассмотренной процедуры сжатия будет крайне неэффективным. Например, можно с целью выравнивания формы спектров сигнала и шума применить фильтр, полностью обрезающий (превращающий в нуль) большую часть частотного диапазона. Такая процедура возможна, она увеличит обычное наблюдаемое отношение сигнала к шуму примерно также, как и рассмотренное умножение спектра на Существенно, что величина сжатия сигнала, достигаемая при процедуре деления спектров, не зависит от шума. Сигнал при увеличении шума может "утонуть" в нем, но он при этом не расширяется, сохраняя достигнутую малую длительность при увеличении шума. Соответствующая программа приведена на рис. 3.18.
|
1 |
Оглавление
|