
Литература / М.И.Финкельштейн Основы радиолокации
.pdfРассмотрим схему рис. 2.27. Для определения фазы в какой-либо момент времени t необходимо к начальному
t
значению фазы в момент / = 0 добавить интеграл 2л f b
В каком-либо периоде повторения, который назовем ус ловно первым, начинающемся при t = 0, фаза передатчика
• t |
фаза |
отраженного сигнала, за- |
<Рп = Фпо + 2л f /п |
||
о |
|
|
паздывающего на время /3, |
|
|
Фс = Фпо + 2л f |
fa(t)dt — срц, |
|
|
о |
|
t
фаза местного гетеродина (рмг = фмг0 + 2л J /мг (t)dt.
о
Фаза отраженного сигнала на ПЧ, действующего на фа зовый детектор,
|
* ^3 |
t |
Фспч = Фс Фмг — 2л |
fп (f) dt |
2л /мг (£) dt -f- |
|
о |
о |
Н"фпо |
ФмгО Фц« |
(5.6.23) |
В системе, показанной на рис. 2.27, |
на выходах первого |
и второго смесителей образуются разностные частоты, а следовательно, и разности фаз. Начальная фаза синхрони зируемого когерентного гетеродина, т. е. начальная фаза опорного сигнала, действующего на фазовый детектор, рав на разности фаз передатчика и местного гетеродина в мо мент окончания синхронизирующего импульса, когда про цесс синхронизации следует считать установившимся. Со ответственно начальная фаза когерентного гетеродина в
момент |
t = ти |
|
|
|
/• |
т |
\ / |
Т |
' |
Фыо=-(фпо + 2л \fn(t)dt ) — ( фмго + 2лрмг(0^ У (5.6.24)
\ о J \ о )
а фаза когерентного гетеродина в произвольный момент времени
t |
(5.6.25) |
фкг = Фкго-г 2л j* /кг (0 dt. |
320
Разность фаз когерентного гетеродина и сигнала на вхо
де фазового детектора |
в соответствии С (5.6.23)—(5.6.25) |
||
|
|
ти |
ти |
ср сркг Фспч ~ 2л j* |
fa (/) dt 2л J fмг (/) di |
||
|
|
b |
о |
/ |
ts |
|
/ |
= — 2л |
J |
fa (t) dt + 2л J fMr (/) dt + |
|
|
|
ти |
ти |
|
|
f |
|
Ч~ 2л fKV (/) dt ф- (рц.
Ь mvMCHr прихода отраженного сигнала t — t3 разность фаз на входе фазового детектора
ТИ |
*3 |
*3 |
Ф = 2л J fa (t) dt 4- 2л J fMP (t) dt + 2n J fKr (t) dt + q^,
0 |
ти |
ти |
(5.6.26)
Если частоты передатчика /п и гетеродинов fMr, fliT, ста бильны, то
ф = 2л/п ти + 2л/мг (t3 — т„) -t- 2л/кг (/3 — ти) 4- срц.
(5.6.27) Для одной и той же цели, расстояние до которой от пе риода к периоду повторения не меняется, т. е. t3 = const при фц ~ const, разность фаз в последующих периодах повторения остается неизменной (фх = ф2). Соответствен но амплитуда нескомпенсированного остатка на выходе
системы ЧПК
А(/ — UG cos ф! — Uc cos ф2 = О,
т.е. компенсация является полной.
Пусть частоты передатчика и когерентного гетеродина
стабильны, а частота местного гетеродина уходит по линей
ному закону: /мг = /мг0 + |
скорость ухода частоты |
fмг |
(5.6.28) |
kKr = |
где А/мг и &t — произвольные приращения частоты и вре мени.
321
Разность фаз в |
двух |
последующих |
периодах повторе |
ния О ...Тп и Тп |
...2ТП |
при /3»тп: |
|
|
|
Г*п+з |
|
ДфМг = Ф2—Ф1 = 2п J |
— |
||
|
^'г |
п |
2л/?мг Тп Тп |
ти |
~ 2л£мг t3 Тп. |
(5.6.29) |
|
|
Если разность фаз Дфмг достаточно мала, то разность амплитуд сигналов в двух последующих периодах повторе ния АС/ — Uс cos q?2 — Uс cos срх « (/сДфмг sin ср,
откуда коэффициент улучшения |
|
Л = 2/Д<р® г = 2/(2зт/гмг t3Tn)2. |
(5.6.30) |
Отсюда [см. (5.6.28)] допустимая скорость ухода часто ты местного гетеродина
^мгдоп ~ (^/мг/Д0доп ~ 1/4,441///1^Тп.
Если, например, принять допустимым 1Х ~ 33 дБ (ко эффициент компенсации /<к — —30 дБ), т. е. /1ДОП = 1995, то
(Д/мг/Д0доп ~ 1/200/зТп, |
(5.6.31) |
а если учесть, что /зтэх = 7П, то (Д/Мг/Д/) |
гаах |
«1/20075.
Точно так же для когерентного гетеродина
(Д/кг/ДОдоп ~ 1/200 4 7П; (А/кг/А0Доптах - 1/200 Г„,
Что же касается допустимой стабильности передатчика, то из (5.6.26) имеем
|
ти |
Лфп = ф2 — ф1 = 2л J |
(fno + ^n 0 dl—2л J (fno + ^nO#, |
Тп |
0 |
откуда Дфп = 2л/гптиТп и далее таким же способом, каким полу чено (5.6.31), найдем
(Д/п/Д0доп = 1 /200ти Та. |
(5.6.32 |
Так так в рассматриваемом случае фазирования по промежуточ ной частоте /кГ < /мг, то при равной абсолютной стабильности мест ного и когерентного гетеродина более жесткие требования предъ я в ляются к относительной стабильности местного гетеродина.
Пусть, например, при ти = 2 мкс, Тп = 2 мс, f3max = 1 мс <D = 150 км), /п — 3 ГГц, /кг = /пч = 30 МГц. Допустимые ско-
322
poctrt ухода частоты (А/кг/Д/)доп = (А/Мг/А/)дОП ~ 2,5 кГц/с; (Д/п'Д'ЛдОП ~ 1 ,25 МГц/с. Так как на работу системы СДЦ влияют уходы частоты, происходящие за период повторения Тп> то соот ветствующие допустимые нестабильности (А/мг)д0П (Д/кг)доп = = 5 Гц; (Д/П)доп = 2,5 кГц. Соответственно относительные неста бильности за период повторения:
(Д/мг//мг)дОП ~ I >7 • 10 |
9; |
(Д/кГ/7кг)дОП ~ 1,7-10 7 ; (А/д/^ц) ~ |
|
|
0,8-10-6. |
Требования к кратковременной относительной неста |
||
бильности 10-6 ...10-7 |
за |
период повторения для передат |
чика (например, магнетрона) и когерентного гетеродина сравнительно легко выполняются на практике без приме нения специальных мер. Относительная же стабильность частоты местного гетеродина должна быть примерно 10-9, поэтому требуются специальные меры стабилизации час тоты. к ним можно отнести использование местного гете родина, работающего на гармониках кварцевого генерато ра, применение отражательного клистрона, в котором ко лебательная система стабилизирована высокодобротным резонатором; применение отражательного клистрона с ав топодстройкой частоты и др.
При результирующей оценке следует иметь в виду, что изменение фаз колебаний передатчика и гетеродинов про исходят независимо. В наиболее неблагоприятном случае они суммируются арифметически, а .в благоприятном ком пенсируют друг друга. Среднее квадратическое значение результирующей фазы
*р< (ф) = К0йг (ф) + 0кг (Ф) + Оп (ф)- |
|
|
Полагая ап (ср) ж сгМР (ср) « аКР (ср) |
= а (ср), |
имеем |
ар (ф) ~ V3or (<р). Если подставить Дер = |
(ср) в |
(5.6.21), |
то при допустимом коэффициенте улучшения 33 дБ <тдоп (ср)=
- V2/1995 • |
3 = 0,018, откуда по |
аналогии с предыду |
щим сгМР (/) |
= огКР (/) = 0,018/2л4; |
оп (/) — 0,018/2лти. |
Если время задержки в системе ЧПК Т3 =£ Тп то, как нетрудно видеть, на выходе появляется остаток в виде двух< импульсов длительностью Дти = Та — Т3. Так как после подавителя имеется двухтактный детектор, изменяющий полярность отрицательных импульсов, то образуются два импульса общей длительностью 2Дти. Остальные цепи, яв ляющиеся для этих импульсов достаточно узкополосными, производят их усреднение, так что их амплитуда уменьша ется (при прямоугольной форме) до Д С7 = £/с2Дти/ти. От
323
сюда относительное допустимое изменение задержки [см. (5.6.20), (5.6.21)] (Дти/Тп)доп = 1/7/27^7, где q = = Тп1тп — скважность.
Например, при /1ДОП = 1995 и q — 103 получим, что (Ати/Гп)доп — 1,6 • 10~5, т. е. достаточно жесткое требо вание.
5.7. ЦИФРОВАЯ ФИЛЬТРАЦИЯ ПРИ ОБРАБОТКЕ РАДИОЛОКАЦИОННЫХ СИГНАЛОВ'
1. Общие сведения о цифровой фильтрации. Понятие цифрового фильтра (ЦФ) используется подобно тому, как это делается для аналоговых систем обработки, часто име нуемых фильтрами. Так же, как и аналоговый фильтр, ЦФ
Рис. 5.35. Структурная схема цифрового фильтра
преобразует сигнал по заданному алгоритму, например алгоритму СФ. Как известно из § 4.2, в аналоговой фильтра ции возможна обработка во временной области, описывае мая с помощью интеграла свертки и осуществляемая с по мощью элементов задержки сигналов, и обработка в час тотной области, описываемая интегралом Фурье и осущест вляемая с помощью фильтров частотной селекции. Такие же алгоритмы используются при цифровой фильтрации.
Структурная схема ЦФ изображена на рис. 5.35. Вход ной непрерывный сигнал х (t) с помощью импульсного эле мента ИЭ дискретизируется, т. е. превращается в последо вательность амплитудно-модулированных импульсов Хт (/)> следующих через интервал дискретизации Т. При этом длительность этих импульсов значительно меньше ин тервала Т. В аналого-цифровом преобразователе АЦП осу ществляется квантование (дискретизация по уровню) и переход от аналоговых величин к числовым хп, т. е. каждо му уровню присваивается кодовое слово, являющееся после довательностью бинарных единиц 0 и 1. Процессор (ариф метическое устройство—АУ) преобразует по определенному алгоритму последовательность хп в z/n; она с помощью циф ро-аналогового преобразователя ЦАП и устройства сгла живания УС превращается в выходной аналоговый сигнал
324
у (/); АУ является важнейшим элементом ЦФ, и поэтому
часто под цифровым фильтром понимают именно АУ, реа лизующее определенный алгоритм.
При переходе от аналогового прототипа системы обра ботки сигналов к его цифровому эквиваленту в простейшем случае исходят из того, что чем меньше интервал времен ной дискретизации и шаг квантования, тем ближе свойства цифрового эквивалента к свойствам аналогового прототи па. Однако такой подход вызывает усложнение аппарату ры, из-за роста разрядности чисел и емкости памяти запо минающего устройства (ЗУ) требуется высокое быстродей ствие узлов. Это имеет смысл, например, при построении подавителей систем СДЦ, но, как правило, не использует ся в устройствах обнаружения, так как последние облада ют высокой эффективностью при числе уровней квантования 2 ...4, когда потери в пороговом сигнале меньше 2 ...3 дБ.
Цифровые методы обработки сигналов отличаются ста бильностью и универсальностью, позволяют исключить дрейфы, обусловленные влиянием температуры и вибраций. Для получения требуемых Характеристик фильтров доста точно изменить весовые коэффициенты. Кроме того, появ ление больших интегральных схем (БИС) значительно об легчило практическую реализацию ЦФ.
2. Общие сведения о ^-преобразовании. Дискретизация сигнала, которая осуществляется на входе ЦФ, описывает ся с помощью решетчатых функций х (k) = х (kT), где k =
= 0, ±1, ±2, ... При анализе решетчатых функций вместо аппаратуры дифференциального и интегрального исчисле ния используется их дискретный аналог— метод конечных разностей и сумм. В общем случае линейного ЦФ с по стоянными параметрами входной и выходной сигналы свя заны разностным уравнением
У (■)* + М (k — 1) + b%y {k — 2) + ••• + bny (k — n) =
= а^х (k) + агх (k — 1) -J- ••• + отх (& —^)» (5.7.1
порядок которого равен наибольшему вычитаемому аргу менту у искомой функции у.
Реакция ЦФ у (k) в момент k Т определяется значе нием входного сигнала в этот момент, т. е. х (&), а также ли
нейной комбинацией предшествующих значений входного х (k — i) и выходного у (k — г) сигналов. Обратная связь
с задержанными входными сигналами характерна для так называемых рекурсивных фильтров. У нерекурсивных фильтров выходная величина зависит только от текущего
325
значения п конечного числа 'Предшествующих значений входного сигнала:
y(k) = 2 atx(k—i). |
(5.7.2) |
Л = 0 |
|
Для определения частотной характеристики и реакции |
|
системы целесообразно воспользоваться |
дискретным пре |
образованием Лапласа в виде ^-преобразования:
Рис. 5.36. Преобразование p-плоскости в г-плоскость
где z — ерг, р — о -{- jco, е_рг = а-1 соответствует за держке на интервал дискретизации в плоскости комплекс ного переменного р, a z-1 означает такую же задержку в
плоскости |
комплексного переменного z. Так как |
|
z = х + ]у = ерТ = еаТ ejb>T = eaT coscoT 4- je°TsincoT, |
||
то |
х = e°r coscoT; |
у — eaT sin о/Г, |
откуда |
|
|
г = | г | |
= ]/х2-р г/2 = еоГ; |
ср = argz = соГ + 2nk, (5.7.4) |
где k — целое число.
Отсюда следует, что левая часть p-плоскости в виде по лосы шириной 2л/7", параллельной вещественной оси о, превращается в круг единичного радиуса z-плоскости. Все параллельные полосы той же ширины соответствуют этому же кругу (рис. 5.36, а, б). Правая же часть p-плоскости со ответствует всей z-плоскости вне единичного круга (рис. 5.36, в, г).
326
Для z-преобразования сохраняются в силе основные тео ремы преобразования Лапласа:
теорема линейности
Z [AlX1 (k) |
+ Л2х2 (k)] = AXZ |
(6)1 + Л22 [x2 (6)1; (5.7.5) |
|
теорема |
смещения |
|
|
|
Z[x(k± n)] = z±«Z [x (6)1; |
(5.7.6) |
|
теорема |
свертки, согласно |
которой ‘ |
произведение |
Z [хх (k)\Z [х2 (6)1 представляется во временной области как
дискретная свертка:
п
f(n) = 2 М<)хг («-<). |
(5.7.7) |
1= о |
|
Связь между решетчатыми функциями и их г-преобразо- ванием дается с помощью специальных таблиц (см., напри
мер, табл. |
5.1). |
X (0 |
x (k) |
б (0 |
*«() |
1 (0 |
1 (k) |
t |
k |
V |
ft2 |
е“ае |
e—a.kT |
sin со t |
sinco k T |
cos co t |
cos co ft T |
Таблица 5.1
X (2)
1
г
г— 1
2
(г-1)2
z(z-M)
(z—I)3
2
2—e
______zsin co T______
z2—2 z cos co T-f-1
z2—z cos co *7
z2 — 2 z cos co T 4-1
327
Применяя (5.7.5) и (5.7.6) к (5.7.1), можно определить передаточную функцию системы:
~£ atz-i
K(Z) = X±L = -S-------------- |
|
|
X (2) |
1 + 2 М”* |
|
|
1=1 |
|
(2—Zbi) (z |
гоа)' * (2* —гот) |
(5.7>) |
(2—Р1)* (г— 2рз)- • • (2 — 2рп> ’
где zok — нули системы, являющиеся корнями полинома в
числителе, a |
zpk — полюсы |
системы, являющиеся корнями |
||||||||
|
|
|
|
полинома в знаменателе, при |
||||||
|
|
|
|
чем для выполнения усло |
||||||
|
|
|
|
вия |
физической |
реализуемо |
||||
|
|
|
|
сти |
требуется, |
чтобы |
степе |
|||
|
|
|
|
ни |
|
полинома |
|
числителя и |
||
|
|
|
|
знаменателя |
были |
связаны |
||||
|
|
|
|
соотношением |
т |
|
п |
(в про |
||
|
|
|
|
тивном случае в момент kT |
||||||
|
|
|
|
должно быть известно значе |
||||||
|
|
|
|
ние входного сигнала х (&+1)). |
||||||
|
|
|
|
|
Из (5.7.8) следует возмож |
|||||
|
|
|
|
ность наглядного |
представ |
|||||
Рис. 5.37. К построению АЧХ |
ления АЧХ как модуля К (г) |
|||||||||
по положению полюсов и нулей |
при |
z = е^г, |
когда |
изобра |
||||||
|
на 2-плоскости |
|
||||||||
|
|
жающая точка |
е* иГ |
движет |
||||||
ся по |
единичной |
|
||||||||
окружности |
z-плоскости [см. |
(5.7.4) |
||||||||
при о= 0] . |
При этом |
|
|
|
|
|
|
|
||
|
(со) — 1 е|<оГ~2о1 II е1оГ—гоа | |
• • -• |е1<оГ— zomI |
|
|||||||
|
|
|е^-2р1Це^-2р2[....|е^-2рП| |
|
|||||||
|
|
|
^01 ^02* |
' ’Г0ТП |
|
|
|
(5.7.9) |
||
|
|
|
ГРГГР2” ‘ГРП |
|
|
|
||||
|
|
|
|
|
|
|
||||
где г01, |
г02, |
•••» rQm, |
rpl, гРъ ..., |
грп~ расстояния от нулей |
и полюсов передаточной функции до произвольной точки на единичной окружности (рис. 5.37).
Так как для z — 0 расстояние rQp = 1, то АЧХ не зави
сит от нулей и полюсов, расположенных в точке z = 0. С помощью формулы (5.7.9) легко представить себе ход АЧХ
для достаточно сложных случаев. |
/ |
|
Если входное воздействие является единичной импульс |
||
ной решетчатой функцией, т. е. |
5 (k) = 1 при k |
0 и |
328
8 (k) — 0 при k < О, то согласно (5.7.3)
Z[6(fc)] = 2 6(Aj)z-ft = $(0)2°=l.
k—Q
Так как реакцией на единичную импульсную решетчатую функцию является решетчатая импульсная характеристика
4 г t/r„
а) |
S) |
|
Рис. 5.38. Система ЧПК.
а — структурная схема, б — положение нуля н полюса, в — импульсная ха рактеристика
g (k), то z-преобразованием импульсной характеристики
является передаточной функцией системы |
|
Z[g(k)] = 2 g(k)z~^X(z). |
(5.7.10) |
А = 0 |
|
Так как согласно (5.7.8) для нерекурсивного фильтра
т
K(z) = 2 fife2~ftfTO из (5.7.10) вытекает, что нерекурсивные
л= о фильтры обладают конечной импульсной характеристикой, для которой
g(k) = ah. |
(5.7.11) |
3. Подавители системы СДЦ как1 цифровые |
фильтры. |
Простейшим подавителем системы СДЦ, представляющим собой РГФ, является, система ЧПК. Ее можно представить, как показано на рис. 5.38, а. Разностное уравнение для
этой цепи |
(5.7.12) |
у (k) = х (k) —~x(k — 1). |
z-преобразование обеих частей дает Y (z) = X (z) — г~гХ (z),
откуда передаточная функция
К (z) = У (z)/X (z) « (z -- l)/z, |
(5.7.13) |
что соответствует нерекурсивному фильтру |
1-го порядка. |
Это иллюстрируется на z-плоскости (рис. 5.38, б) нулем
в точке z=l (полюс расположен в точке z — 0). В данном случае (в отличие от п. 2) оператор г-1 характеризует за
держку не на интервал дискретизации Т, а на период по
329