Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Радиолокационные измерители дальности и скорости. [Т.1]

.pdf
Скачиваний:
14
Добавлен:
12.11.2023
Размер:
21.62 Mб
Скачать

Д= V + КфпАй! + кф12Д22,

Д(0)=До;

 

А

V(0)=V0;

(6.3.12)

V = ан 4-ац +кф21Azx +ЦдаAz2,

 

£ц = -сйц + КфзхАг! + Кфз2А22,

ац(0)=апо;

 

Azi гд-кдД—идд-кд(Д-Ду),

 

 

AZ2==ZV“KvV“-иду"Ку(V-Vy))

 

 

где были использованы представления (6.8.6).

Следует отметить, что коэффициенты Кф^, Кфу1> КФУ2’ K4»U (i—1.2. j-1,2,3), используемые в алгоритмах (6.3.8и[б.3.13), рас-

считываются по формулам (1.4.4), (1.4.5).

Соотношения (6.1.39) и (6.3.8)-(6.3.13) определяют алгоритм функционирования ИДС с декомпозированным фильтром. Струк­ турная схема дальномерного канала БРЛС, соответствующая этому алгоритму, приведена на рис. 6.3.1.

Анализ этого алгоритма позволяет отметить следующие осо­ бенности.

Для его реализации необходимо решать

Nfl=Na4-Ny+NOT=<l-K),6-l*2)+{2+0>6-2*3>f(3-rt),6*34)=s16 (6.3.14)

уравнений, где Na=2 определяет количество уравнений фильтра ускорений, Ny=5 - уравнений фильтра управителя, а N^=9 - фильтра отслеживаемых координат. Сравнение (в.3.14) с (6.1.59) свидетельствует о существенном упрощении процедур формирова­ ния оценок в декомпозированном фильтре. Следует однако отме­ тить, что в результате декомпозиции исходного вектора состояния в полученных фильтрах были утрачены взаимные связи и коррек­ тирующие поправки, обусловленные невязками, вычисленными в других фильтрах. Несмотря на это, в полученном дальномере со­ хранились все типы контуров, имевших место в дальномере, схема которого приведена рис. 6.1.1. В связи с этим, полученный алго­ ритм сохраняет все достоинства, присущие многоконтурным сле­ дящим системам, которые были рассмотрены в п. 6.1.6. Как и ра­ нее полученный алгоритм фильтрации основан на использовании коэффициентов усиления невязки, вычисленных для вполне опре­ делённой статистики. При изменении условий функционирования, связанных с типом сопровождаемой цели, дальностями до неё и видами её манёвров, принятая статистика не будет соответствовать той, которая была заложена при синтезе фильтра.

Исследования полученного алгоритма проводились путём со­ вместного моделирования отслеживаемого процесса (6.1.4), задан­ ной части (6.1.5), результатов измерений (6.1.19) и алгоритмов функционирования фильтров (6.3.8)-(6.3.13) и регулятора (6.1.39), (6.1.52), (6.1.53).

В процессе этих исследований проводился анализ потенциаль­ ной точности, переходных процессов, реальной точности и време­ ни памяти при тех же условиях, что и для дальномера с нерасщеплённым дальномером (см. §6.2)

Потенциальные дисперсии ошибок фильтрации вычислялись в процессе решения нелинейных дифференциальных уравнений Риккати (1.4.5) при соблюдении допущений, изложенных в п.6.2.1. Анализ полученных результатов позволяет сделать сле­ дующие выводы.

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

Сравнение с результатами, полученными в п.6.2Л, позволяет придти к заключению, что потенциальная точность декомпозиро­ ванного фильтра превышает точность оценивания аналогичных координат в нерасщеплённом фильтре. Причиной этого является разрыв взаимных связей, по которым в каналы нерасщеплённого фильтра поступали шумы измерений от всех измерителей.

Ошибки захвата устраняются несколько быстрее и их значе­ ния меньше, чем в нерасщеплённом фильтре, во всём возможном диапазоне их значений при любом сочетании по знаку. Динамиче­ ские ошибки в установившемся режиме практически равны нулю.

Постоянные времени Тд* Tv, Тду, Туу, Тан отработки первона­ чальных ошибок оценивания составляют тысячные доли секунды и несколько меньше, чем в нерасщеплённом фильтре. В тоже вре­ мя постоянные времени отработки ошибок АДр=Д-Ду и AVp=V-Vy составляющие десятые доли секунды, такие же как и в дальноме­ ре, рассмотренном в §§ 6.1, 6.2. Улучшение быстродействия даль­ номера с декомпозированным фильтром по сравнению с фильт­ ром, исследованным в п. 6.2.1 можно объяснить снижением инер­ ционности оптимальных фильтров, которая обусловлена уменьше­ нием числа используемых невязок и, соответственно, уменьшени­ ем количества взаимных связей.

Фильтры устойчиво формируют оценки всех фазовых коорди­ нат во всех условиях, включая и не соответствующие тем, на ко­ торые они были оптимизированы. Изменения начальных дально­ стей и скоростей практически не влияют на точность оценивания.

Реальные точности несколько хуже потенциальных, однако в целом, точность оценивания всех фазовых координат несколько лучше, чем в дальномере с нерасщеплённым фильтром.

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

Результаты, полученные в процессе исследования функциони­ рования декомпозированного измерителя при наличии не только возмущений, но и ошибок захвата ДДрДО^Д^-ДуДО) и AVpCO^VCOJ-V^O), в основном совпадают с теми, которые анализи­ ровались в п.6.2.2, с учётом небольших количественных отклоне­ ний за счёт влияния конкретных реализаций возмущений. Полу­ ченные результаты свидетельствуют о высокой робастности и точ­

ности синтезированного дальномера в процессе его функциониро­ вания в условиях влияния ошибок захвата и возмущений.

Время памяти в ИДС с декомпозированным фильтром на 2-3 секунды больше, чем при использовании обобщённого фильтра. Это обусловлено меньшим уровнем реальных ошибок оценивания на момент пропадания входных сигналов.

Подсчёт количества операций, используемых в алгоритмах функционирования ИДС с декомпозированным фильтром показал, что для их реализации требуется примерно в 4-5 раз меньшее бы­ стродействие ЦВМ, чем при использовании нерасщеплённого фильтра.

6 .3 .3 . А в т о м а т и ч е с к и й и з м е р и т е л ь д а л ь н о с т и и с к о ро с ти со

СНИЖЕННОЙ ИЗБЫТОЧНОСТЬЮ ИНФОРМАЦИИ

Анализ алгоритмов функционирования ИДС с декомпозиро­ ванным фильтром и результатов его исследования, проведённые в п.6.3.2, позволяют придти к заключению, что в нём по-прежнему сохраняется высокая степень избыточном информации. В частно­ сти, информация о состоянии заданной части дублируется фильт­ ром управителя. Эта особенность позволяет сделать вывод о воз­ можности дальнейшего упрощения процедур функционирования ИДС за счёт исключения из его состава фильтра управляемого процесса, используя в сопрягаемых алгоритмах вместо оценок Ду

и Vy их естественные значения Ду и Vy, если это позволяет струк­

тура заданной части.

Кроме того, анализ точности функционирования фильтра соб­ ственного ускорения свидетельствует о возможности использова­ ния в алгоритмах работы ИДС непосредственно результатов изме­ рений, формируемых акселерометром, отмасштабировав их соот­ ветствующим образом.

С учётом этих замечаний алгоритм функционирования упро­ щенного дальномера будет определяться следующими уравнения­ ми:

заданной части

Ду =

Ду(0) = Дуо>

 

(6.3.15)

Vy =bvuv+^vy,

Vy(0) = vy0,

измерителей

 

384

идд"кд(Д-Ду)+4да;

u ^ = Kv( v - v y) + u ;

(6.3.16)

Z* ~

4аи»

 

регулятора

 

 

иу=кдДД+куДУ,

(6.3.17)

где

 

 

ДД = Д -Д у, ДV = V - Vy;

(6.3.18)

фильтра отслеживаемых координат

 

Д = V + кф11Дгг +K$12AZ2 ,

Д(0)=До;

А

 

 

V = za /к а +ац +K^ XAZX+кф22Дг2,

V(0)=V0; (6.3.19)

а„ = - о а ц + КфзгД21 + « « 32^ 2 »

£ц(0)=ацо

Дг1=идд-кд(Д-Ду);

(6.3.20)

Az2~Ujy“Ky(V -Vy)•

Структурная схема упрощенного дальномера, соответствующая этим уравнениям, приведена на рис. 6.3.2.

Исследования упрощенного ИДС проводилась путём совмест­ ного моделирования отслеживаемого процесса (6.1.4), заданной части (6.3.15), результатов измерений (6.3.16), алгоритма функ­ ционирования фильтра отслеживаемого процесса (6.3.19), (6.3.20) и регулятора (6.3.17), (6.3.18) в условиях идентичных тем, в ко­ торых исследовались рассмотренные ранее, более сложные даль­ номеры.

Анализ полученных результатов позволяет придти к следую­ щим заключениям.

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

Рис. 6.3.2.

Упрощённый дальномер обладает несколько большей скоро­ стью отработки ошибок захвата, однако худшей (на 7-40%) точно­ стью оценивания дальности, скорости сближения и ускорений по сравнению с дальномером с декомпозированным фильтром. Для реализации рассмотренного алгоритма требуется на 20-25% меньше быстродействия ЦВМ, чем при использовании декомпози­ рованного фильтра.

Проведённые исследования дают возможность утверждать, что

вупрощённом дальномере, как и в декомпозированном, можно использовать неадаптивный фильтр, оптимизированный на работу

вусловиях средних интенсивностей возмущений*

6.4. ИЗМЕРИТЕЛЬ СКОРОСТИ, ОПТИМАЛЬНЫЙ ПО ЛОКАЛЬНОМУ КРИТЕРИЮ

Рассматриваемый измеритель скорости (ИС) предназначен для формирования оценок скорости и её производных в процессе со­ провождения интенсивно маневрирующих целей в РЛС с непре­ рывным или квазинепрерывным зондирующим сигналом. Сопро­ вождение таких целей можно потребовать оценивания не только скорости и ускорения, как это было принято в рассмотренных ра­ нее измерителях (§6.1, п.п. 6.3.2, 6.3.3), но и производной уско­ рения. Знание этой производной необходимо для повышения точ­

ности и устойчивости сопровождения и для констатации факта начала интенсивного маневрирования цели.

Синтез ИС будет осуществляться по алгоритмам СТОУ опти­ мальным по локальному критерию (см. §1.11). Эти алгоритмы по­ зволяют для заданной части (управителя)

Vy=ay,

Vy(0) = Vy0;

 

ау = baua + £ау

ау(0) = ауо,

(6.4.1)

 

предназначенной для отслеживания процесса

 

V = а = ац +ан,

V(0) = V0 ;

 

ан = с н.

ан( ° ) = ан о;

 

а ц = ^ц»

ац(0) - ад о ;

(6.4.2)

Сн(0) = Си0 ;

 

С„ = ^н .

 

 

“CtдСц + ^Д >

Сц(0) = Сцо ,

 

при наличии наблюдений

 

 

Иду

Ky(V - Vy) +

 

 

Zyy

КууУу+^ууд;

 

(6.4.3)

^*ан — ^анан ^ани

 

 

сформировать сигнал управления иа, оптимальный по минимуму локального функционала качества

1—►

Я

II

 

>

а

-----1 >* 1

-ау

Qvii <lvi2

_4.v21 Чу22.

>

1

1

 

 

1

 

+

а-ау

t

Jufkadt[,

0

В соотношениях (6.4.1)-(6.4.4) ау “ управляемое (отслеженное) ус­ корение сближения; Ьа - коэффициент усиления сигнала управле­ ния и*; ац, ан и Сц, Сн - соответственно проекции ускорений цели и ОУ и их производных на линию визирования; ац - коэффициент учитывающий маневренные свойства цели; uflv, и zM - изме­ ренные сигналы на выходе частотного дискриминатора, датчика управляемой скорости (управляемого гетеродина) и акселерометра; ку, КууИКадкоэффициенты передачи датчиков; £ау, £н, £ц и ^ £ууи* ^аын ~ центрированные белые шумы состояния и измерений с

387

известными односторонними спектральными плотностями; qvn, qvi2=(lv2i и 4v22 ~ коэффициенты штрафов за точность слежения по скорости и ускорению; ка - коэффициент штрафа за величину управляющего сигнала.

Состав вектора управляемых координат (6.4.1) выбирают, ис­ ходя из необходимости формировать сигнал управления с учетом ошибок как по скорости, так и по ускорению, и с учетом обеспе­ чения памяти по ускорению при кратковременных пропаданиях радиосигналов. При выборе состава вектора отслеживаемых коор­ динат (6.4.2) принималось во внимание требование устойчивой се­ лекции по доплеровской частоте сигналов, отражённых от интен­ сивно маневрирующих целей, и необходимость формирования оценок V, а, ан, ац, Сн и Сц. Оценки V могут быть использова­

ны в современных и перспективных алгоритмах наведения и помехозащиты, информация о а, ан, ац, Сн и Сц необходима для

повышения точности и устойчивости оценивания скорости при ин­ тенсивном маневрировании ОУ и цели, а также при экстраполя­ ции её пространственного положения. Кроме того, эти оценки также используются в алгоритмах помехозащиты, особенно от уводящих по дальности и скорости помех [31], и для констатации факта начала маневра цели.

Поскольку все исходные модели линейные, шумы гауссовские, функционал качества квадратичный, то в соответствии с выводами теоремы разделения (п.1.9.3) синтез оптимального регулятора и фильтра будет выполняться раздельно. При синтезе регулятора для уравнивания размерностей векторов управляемых и отслежи­ ваемых координат будем полагать, что в (6.4.2) ан=0 и ап=0 при

анО= 0 , ацо=*0. Тогда, сопоставляя (6.4.1) с (1.9.1), (6.4.2) с (1.9.2) и (6.4.4) с (1.11.9), получим

*, = [V а г.

X

II

Ат = А

1

0

*

У “

1

 

 

0

 

F = F

0

1 '

у

 

 

хт ху

0

0

 

 

 

а у ]Т’

 

11==ua.

к = к а,

 

"

0

 

(6.4.5)

 

 

 

9

 

В У = А.

 

 

Q

Qvii

<lvl2

 

=

 

 

 

 

 

, Q

V

2 1

Q V 2 2 _

 

Используя (6.4.5) в (1.11.10) получим алгоритм функциониро­ вания оптимального регулятора

ua =ЬаЧуггк!1^ -Vy) +b^q^k;1^ -ay) = KvAV +каДа, (6.4.6)

где AV= V -Vy и Aa=a-ay - ошибки слежения по скорости и уско­

рению, a Kv = baqv21k“1, ка * baqv22ka1 - коэффициенты передачи регулятора по ошибкам слежения. Бели использовать методику оптимизации коэффициентов ку и ка, изложенную в п. 6.1.4, то можно найти их значения:

RV _

Ьа'ЦцЦудоп "

. ка _ ЬаТ^опЦудОД+АУ0

^

 

^а^доп(^0 + ДафТудод)

ЬаТуд0П(ДУэ+ДаоТудоп)

 

обеспечивающие минимальные ошибки слежения в установившем­ ся режиме при заданных ограничениях сигнала управления иа<иудоп и постоянную времени Т^Т^од отработки максимально возможных ошибок захвата по скорости ДУ0 и ускорению Да<).

Анализ (6.4.7) позволяет сделать следующие заключения. Значения коэффициентов KV и ка зависят не только от параметров заданной части (Ъа) и накладываемых на регулятор ограничений (ЦудоЕр ТУдОП), но и от точности устройства поиска и захвата сигна­ ла цели по доплеровской частоте, предопределяющего величину ошибок AVQ и Дао* Следует отметить, что в полученном регуляторе можно реализовать лишь постоянные времени, удовлетворяющие условию

ТУДоп-|Дао + +ЬаиудопДУ0)у/ьаиудоп.

Для упрощения процедуры формирования оценок всех требуе­ мых фазовых координат обобщенного вектора состояния (6.4.1), (6.4.2)

х[V ан ац Сн Сц Vy ayJ

используем принцип декомпозиции, рассмотренный в п. 6.3.2. В результате будут отдельно синтезированы:

фильтр собственного ускорения на основе уравнений

ан = Сн, ав(0) = ан0 ;

СН= 4Н, Сн(0) = Сн0; (6.4.8)

^ан ^анан ~^~^эани»

(6*4.9)

фильтр управителя, использующего модели

Vy = ау, Уу(0) - Ууо; ау -

baua + ^ау, ау(0) = а^; (6.4.10)

Z^K^Vy+SvyH;

(б -4 -11)

фильтр отслеживаемых координат, базирующийся на соотно­ шениях:

V = ац +ан

II О >

 

 

йц = Сц,

ац(0) = ац0;

(6.4.12)

Сц= ■ ац^ц + *

Сц(0) =

Сц0;

 

zv ЧдУ "I* KvVy

KvV 4*

.

(6.4.13)

В моделях (6.4.12) и (6.4.13): ан и Vy - оценки собственного

ускорения и отслеженной скорости, полученные в первых двух фильтрах; шум измерений отличается от шума (6.4.3) на величину добавок, вызванных ошибками оценивания Vy.

Использование для выбранных трех типов исходных моделей алгоритма оптимальной линейной фильтрации (1.4.3)-(1.4.б) дает возможность сформировать оценки:

в фильтре собственного ускорения

= СН+ K HlAzft,

а н( ° ) = г ан(0) / к ан;

=

*

 

Сн(0) = Сн0;

 

ZaH- Канан »

 

 

в фильтре управителя

 

 

Vy

&у 4 KyyiAZyy,

«!>*

> II о

ау —baua 4 Куу2 А2 уу,

ау(0) = ауо;

Azуу

Л

 

 

Zyy - КууVy,

 

 

ua вычисляются по закону (6.4.6);

в фильтре отслеживаемых координат

V =

a 4 + a H + K vl Azv,

V(0) = Vo;

Л

аР> О II О

Нц

Сц 4 Kv2 AZy ,

(6.4.14)

(6.4.15)

(6.4.16)

(6.4.17)

(6.4.18)