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

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

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

к / -\ _ K0i(z - !)

+K02TH(z -1) +K03TH

(3.6.34)

/ “

/ .\3

 

(z-i)

 

и включает три интегратора и демпфирующие звенья, постоянные времени которых определяются соотношениями: Тф102/коз;

Тф2 - > 0 1 / К 03 •

В фильтре (3.6.83) формируются оценки трех компонент век­ тора состояния: координаты х г, производной (скорости измене­ ния) координаты х2 и второй производной (ускорения) координа­ ты х3.

Фильтр (3.6.33) в литературе называют а-Р-у фильтром. При этом полагают, что а=Коь р=к02Тн, у=к08Тд являются безразмер­ ными коэффициентами. Для линейной модели дискриминатора (3.6.2), следящая система описывается уравнениями (3.6.33), в которых йд(к) получают в соответствии с (3.6.3). Операторный ко­

эффициент передачи полученной замкнутой следящей системы ра­ вен

™ / \ _______________ (z-l) KQI +(z- + « 0 3 * 1 ^ _______________

(z - 1)3 +(z - 1)2(KQ! + КщТ;) +(z - 4 *02^ +Коз^) + 1^

(3.6.35) Из (3.5.35) следует, что коэффициент усиления разомкнутой сле­ дящей системы равен Ко~Коз/Т^. Приравнивая нулю знаменатель (3.6.35), получаем характеристическое уравнение системы

z8+z2( - 3 + Koi + КогТ;)+ я(з - 2коХKoai; + Коз^)+ Ко2Тн -1 = О.

(3.6.36)

Условия устойчивости следящей системы получаются из (3.6.36) и имеют вид

О < ^08^н ^ ^01^02» ^02 ^ ^ОЗ^н/^J Kgi ^ ^ОЗ^н/^* (3.6.37)

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

_________________ м ! ________________ . (3.6.38)

(г - ! ) 3 + (г-1 )2(к01 + КоД ) +(z ~

+ “ о Л ) + "оо^

При воздействии на входе следящей системы процесса х(к) с изо­ бражением x(z) изображение ошибки слежения определяется вы­ ражением

e(z) = КЙ1(z)x(z),

(3.6.39)

из которого, с учетом (3.6.38), следует, что, при постоянном, ли­ нейном и квадратичном воздействиях, ошибка слежения в устано­ вившемся режиме равна нулю, а при кубическом воздействии с

-/ \ uT|z(z + l)(z + 2)

изображением

x(z) = ------^ у с т а н о в и в ш а я с я

ошибка

равна

 

 

 

еуст = s (k) = lim(z - l)s(z) = — ^— ,

(3.6.40)

к—>оо

z_>1

К03/^н

 

где и - третья производная (скачок) фазовой координаты х(к). Так же, как и в следящей системе второго порядка, в следя­

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

Критический режим определяется условиями равенства кор­ ней характеристического уравнения (3.5.36), которые имеют вид

К01 + К02Тн = 3 л/йозТ?;

«огТе +

= 3 а/козТн

(3.6.41)

Корни характеристического уравнения равны

 

Z 1 = Z2= Z3=1+ YO;

Уо = - л/к0зтн •

(3.6.42)

Если уравнение (3.6.36) имеет три действительных корня, как минимум два из которых различны, то в системе возможны апе­ риодический или колебательный переходные процессы. При этом, установление того или иного процесса зависит от выполнения или невыполнения неравенства

£ 0 ,

(3.6.43)

Х = 4

27

где: ly = к01 + к02Тв; г2 = к02Тв + к08Тв ; г8 = к03Тв . Если х^О, то

имеет место апериодический режим, при х> 0 - колебательный. Различаются две разновидности апериодинекого режима: при

наличии двух и трёх различных корней характеристического уравнения (3.6.36). При имеются только два различных кор­ ня, определяемые выражениями

z1 = l +y1;

Z2 = Z8=1 +Y2,

(3.6.44)

где

 

 

Ti = -2 У Ф - гх/3;

у2 = ^ / 2 - гх/3,

(3.6.45)

а С = ^ /З )3 - ггг2/3 + г8.

При х<0 - все корни характеристического уравнения различ­ ны и равны

ZI = 1 + Y 3;

z2=i+y4;

Z3

= I +Y6,

(з.б.4б)

где

 

 

 

 

у3 = 2 ц - г 1/3 ;

Y4 = -Ц + Р>/3-г1/3 ;

у5 = -\i-pylB-rjB;

 

 

 

 

(3.6.47)

Ц+ i р = $J~ С/ 2

+ л/х - комплексное число.

 

Когда выполняется условие

 

переходной

процесс имеет

колебательный характер. Корни характеристического уравнения в этом случае определяются выражениями

Z i = l + Ye;

Z2 = 1 + Y 7;

 

Z3 = 1 + Y 8,

(3.6.48)

в которых следует положить

 

 

 

Y6 = В + Г - г 1/3 ;

Yr = - ^ ^

- - r 1/3 + i ^ ^ V 3 ;

В + Г

/

. В —Г

гт

 

(3.6.49)

У8 = ---- 2-----ri/3_1—

v8;

 

B = f-C/2 +Vx;

Г =V-C/2-Vx; вг =г2/9-г2/з.

Формулы, определяющие переходные процессы фильтра во време­ ни в различных режимах, приведены в табл. 3.6.3 - 3.6.6, где: у0 - у8, В и Г - определяются из (3.6.42), (3.6.45), (3.6.47), (3.6.49), а Хо, V, а, и - определяют тип входного воздействия - постоянного, линейного, квадратичного и кубичного, соответственно.

 

 

 

 

 

/

\

X0Z

Режим

 

 

Постоянное воздействие. X^ZJ -

z l

 

 

 

 

 

 

 

Критичес­

 

/

ч

,

2ky0

 

. у

кий

 

 

 

 

 

 

 

 

 

( о д и н

 

 

 

 

 

 

 

кореньYn)

 

 

 

 

 

 

 

Апериоди­

 

 

 

 

 

 

 

ческий (два

 

 

 

 

 

 

 

различных Е (к)=

 

 

jyf(l+Yl)k + f e +2f2) Y2 _Yl)_Y2 (1+Y2)k

корня: Yi.

 

\Y1 ~Y2)

l

L

 

 

Y2)

 

 

 

 

 

 

 

 

 

 

Апериоди­

 

 

 

 

 

 

 

ческий (три

hS

[

Yl(l+Y8)k

^(1+Y4)k

+

Yi(l+Ys)k j

различных

e(k)

Хо|(уз-у4)(у8-г 6) (Уз-Г4)(Г4-Гв)

(Ya-Ysb-Yoj j

корня:

Уз» Уа* Уд)

Колебате­ льный (три различных корня:

Ye Тт Та)

1 Ye(l+Ye)k

Y7(l+Y7)k

+ Y|(l+Ye)k ]

8U ” *°{(Ye-Yr)(Ye -Ye)

(Ye -Y7KY7 "Ye)

(Ye "YT)(Y7 - Ye) j

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

ную частоту следящей системы ю0 = ^/к08/Тн =

, а после этого

рассчитываются коэффициенты Ког> коъ которые в

основном оп­

ределяют условия демпфирования в системе

ф1 = к02/коз;

Тф2 = V K0 l / K03 )•

Шумовая полоса пропускания следящей системы, рассчиты­ ваемая по (3.6.11) с учетом (3.6.35), определяется следующим вы­ ражением

 

+ ^(^02 “ КохКоз)]~

 

+ 2KogTHj +

+

к озг^ н (к овг^ н + ^ ш 1^ ) ~

j /

( 3 . 6 . 5 0 )

/

“ КозТдЦв-^Со! - 2 KQ2TH + К озТ ^ ||.

Режим

 

Линейноевоздействие: x(V) =

^

>hZ

 

 

 

 

(z -l)2

 

Критичес­

 

 

 

 

 

 

кий

 

 

 

 

 

 

Апериоди­

 

 

 

 

 

 

ческий

 

 

 

 

 

 

(два раз­

 

 

 

 

 

 

личных

 

 

 

 

 

 

корня)

 

 

 

 

 

 

Апериоди­

 

 

 

 

 

 

ческий (три

 

Y8(l+Ya)k

Y4(l+Y4)k

H

Y5(l+Y5)k

]

различных

8(h)= Vt.

корня)

 

(Y3-Y4)(Y3-Y5)

(Y3-Y4)(Y4-Y5)

(Y8-Y6XY4~Y5) J

Колеба­

 

Ye(1+Ye)k

Y7(l+Y7)k

(

Ys(l+Ye)k

1

тельный

e(k)=vi;<

1

 

(Ye -Y7)(Ye "Ye)

[ УбУт^ т- Ув )

(Ye-YeXYT-Ys) }

 

 

 

 

 

 

При Тн->0 (3.6.50) переходит в известное выражение для не­ прерывной следящей системы с тремя интеграторами и демпфи­ рующими звеньями.

Дисперсия Бф оценки информационного процесса, обусловлен­ ная воздействием аддитивного дискретного белого шума на входе следящей системы, определяется общим выражением (3.6.13) при подстановке в него AF3, вычисленным в соответствии с (3.6.50).

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

2

D =Бф+8уСТ= 2D„THAFa + VK03 /т„ »

для чего целесообразно использовать численные методы нахожде­ ния экстремума функции по двум переменным.

Нетрудно показать, что рассмотренная следящая система с тремя дискретными интеграторами оптимальна для наблюдаемого процесса, описываемого уравнениями

Режим

 

 

 

 

it f M z 2 + 4z + 1)

 

Кубичное воздействие: X^zj = -------- ^--------- -------- -

 

 

 

 

 

( z - l )

 

 

 

 

 

\

Крити­

, ,

V flff

[

^ o jr o - 6)

b(k-l)Yo(y§ +6y0 +б)

 

 

 

 

(i+Yo)

ческий

' l ) '

Г *Д

[

e ( l + y0) +

u ( l + r ,J»

/

Аперио­

дический

(два

различ­

ных

корня)

Аперио­

дический (три раз­ личных корня)

Колеба­

тельный

e ( k ) = ^ { -

 

 

 

k

e , + T? + e , i + e ( l . T l)

W

6

[

V iy2

Y l

*

l l >

I i [ ^ f 2 + h l 2 + 12l 2 - &h )

k(y| + ву2 + б)

 

 

 

 

 

 

( i - r 2)k

[

Y2(YI - Y 2)2

 

' Y2(1+Y2)(YI - Y 2)

E ( k ) ' «

{ '

Y , Y , 7 . + Y a( 7 » - Y 1X Y .'- Y s) ( l n *) *

% . a : - Y . ) ( i + ^

+ Yi a

;

e- Y , ) (i - - » 1

S ( k ) * e

f

Y „ Y ,Y .'+ Y . ( Y . - 7 , t Y r - Y . ) (l + T ,) *

хх(к) = хэ1(к) +кх(к^2(к) - хэ1(к)],

х1(0)=х1О;;

xai(k) = ^i(k - 1)+ Тах2(к -1);

 

i 2(k) = хэ2(к -

1) + K2(k)[z(k) - хэ1(к)],

x2(0)=x2o; (3.6.52)

хэг(к) = x2(k -

1) + THx8(k -1 );

 

x3(k) = x3(k -

1) + K8(k)[z(k) - хэ1(к)],

х8(0)=х30;

 

Ki(k)

"1

т

0'

К(к) =

К2(к) ; Ф =

А

Ан

т ; К(к) = D(k)HTD-1;

0

1

_Кз(к).

0

0

1

1>и=1>и;

Dx==D x;

 

 

 

D(k) = (I - K(k)H)D3(k);

D3(k) = <M>(k - 1)ФТ+Dx.

Коэффициенты усиления оптимального фильтра меняются во времени в соответствии с уравнениями (3.6.53). В установившемся режиме коэффициенты усиления постоянны и находятся из реше­ ния системы уравнений (3.6.53) при к->оо.

3.6.5. А нализ следящ их систем с накоплением в дискриминаторе

Рассмотренные в предыдущем параграфе характеристики дис­ кретных следящих систем с различными фильтрами в следящем контуре основаны на представлении процесса на выходе дискри­ минатора в виде (3.6.3), которое является хорошим приближением при формировании выходного отсчета дискриминатора по каждо­ му импульсу зондирующего сигнала. Однако, в реальных дискри­ минаторах, выходные отсчеты формируются с существенно мень­ шим темпом (Тн=0,005-0,1 с) и получаются в результате накопле­ ния соответствующих отсчетов на данном интервале. Процедура накопления в дискриминаторе изменяет динамические свойства следящей системы, что особенно проявляется при больших време­ нах накопления. Поэтому этот факт необходимо учитывать при более детальном анализе характеристик следящих систем.

Для учета промежуточных N отсчетов, формируемых в дис­ криминаторе, введем соответствующие моменты времени tkji=tk+iTd; i = 1>N, Td = TB/N (tk>N=tk+1). В существующих дискриминаторах используется равновесное суммирование промежу­ точных отсчетов

(3.6.54)

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

копления (см. §1.5). Для пояснения сути происходящих при нако­ плении процессов в дальнейшем в качестве входного сигнала дис­ криминатора будем рассматривать простейшую модель (3.6.54). Полагая, что работа происходит на линейном участке дискрими­ национной характеристики представим с учётом (3 .6.2 1 ),

в виде

®л(*к-1. i) = К *к-1, i) " Xal(tk-1, i);

(3 .6 .66)

5(tk-i, j) = xftk-i, i) + S „(tk_i, i).

где ^H(tk-i,i) - дискретный белый гауссовский шум единичного из­ мерения с нулевым математическим ожиданием и дисперсией DH.

В (3.6.21) входит xsl(tk. lf экстраполированная на момент времени tk_u = tk_x + iTd оценка координаты, которая определяет­

ся принципом построения дискриминатора. Поясним это на при­ мере радиодальномера, в котором оценка xal(tk -1 определяется

положением опорного сигнала

uroe\"k-l, i) ^ »

что в реальных схемах соответствует управлению положением расщепленных стробирующих импульсов. Если на всем интервале накопления Тв опорный сигнал неподвижен, то x8i(tk-i, i) = Xi(tk-i). Если опорный сигнал постоянно (или от им­

пульса к импульсу) перемещается по линейному закону, то

*.i(*k-i. i) = X i(tk -i) + iTdx 2 (tk -i)-

(3.6.56)

Такой режим работы достаточно сложно реализуется практически, однако более оптимален. В дальнейшем будем рассматривать ре­ жим экстраполяции (3.6.56).

Рассмотрим обобщенную следящую систему, которая в п- мерном пространстве состояний описывается векторным уравнени­

ем

 

х(к) = Фк(к -1) + Кйд(к),

(3.6.57)

где йд(к) = ufl(tk) определяется соотношениями (3.6.64Н 3.6.56). Экстраполированную оценку (3.6.56) запишем в виде

 

1

Td

О

о

 

о

1

Td

о

где H«|l О О ...|; Ф

0

0

1

о - матрицы размера

 

0

0

0

1

1хп и пхп соответственно. Подставляя (3.6.54), (3.6.55), (3.6.58) в (3.6.57), получаем

S(k)=<Цк- 1)+К± - 1 * ^ - 1))=

= (ф -КИi £ф '}х(к-1) + K i i 2(tk-j,i) (8.6.59)

В (3.6.59) второе слагаемое представляет собой внешнее воз­ мущение, действующее на систему фильтрации, а первое слагае­ мое определяет динамические свойства системы, такие как: устой­ чивость, переходные процессы, полоса пропускания и др. Вид пер­ вого слагаемого (выражение, стоящее в скобках) в (3.6.59) отлича­ ется от того, какое имеет место в системах фильтрации без накоп­ ления (Ф-КНФ), что и определяет изменение динамических свойств. Конкретизируем (3.6.59) для следящих систем второго и третьего порядков, рассмотренных в п.п. 3.6.2, 3.6.3.

Для системы второго порядка

ф1 = '1 iTd‘ 0 1 .

* •

rH1!

II е

NiTi

•н

 

'1 0,5(N + l)Td'

.0 1

Введем параметр b=(N +l)/2N . Тогда, расчет операторного коэф­ фициента передачи дает следующее выражение

К » (z) = ---------

( z - l K + к у Г .----------------

.

,3.6.60)

(z -1)

+ (z- 1)(KOI + K02THb)+ К02ТН

 

При отсутствии накопления Ь=1 и данное выражение переходит в (3.6.22). При N>>1, что характерно для реальных систем, Ы),б.

Устойчивость следящей системы с коэффициентом передачи (3.6.60) определяется соотношениями (3.6.24), в которых первое условие заменяется на

260