книги / Радиолокационные измерители дальности и скорости. [Т.1]
.pdfк / -\ _ K0i(z - !) |
+K02TH(z -1) +K03TH |
(3.6.34) |
/ “ |
/ .\3 |
|
|
(z-i) |
|
и включает три интегратора и демпфирующие звенья, постоянные времени которых определяются соотношениями: Тф1=к02/коз;
Тф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