Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 195

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

сильно отличается от нормального, однако различие в формулах

(5.68), (5.69) невелико.

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

5.4. Ламповый генератор

Примером автоколебательной системы является генератор гар­ монических колебаний, одна из простейших схем которого пред­ ставлена на рис. 5.7. Колебательный контур включен в цепи анода. Величины сопротивлений конденсаторов Си С2 для переменной составляющей тока малы, а сопротивление R 1 велико, поэтому схема генератора может быть упрощена и представлена в виде схемы, изображенной на рис. 5.8. Эта схема составлена с учетом только пере­ менной составляющей тока.

Уравнения токов и напряжений имеют вид

 

 

J (/, +

1

L ^ r + R

I - ±

/др - I ) d t = ETi (/);

 

 

О

(5.70)

 

 

 

и = М

+

ДТ2 (0;

/ а = Su — Spu3,

+ 0

 

 

 

155'

где приняты следующие обозначения: L, R, С — параметры коле­ бательного контура; М — взаимная индуктивность катушек; / а — постоянный анодный ток; I — ток через индуктивность контура; i — дробовой шум; ETl — э. д. с. теплового шума сопротивления потерь R контура; Ег„— э. д. с. теплового шума активного сопро­ тивления катушки обратной связи; 5 — крутизна характеристики лампы; [} — коэффициент нелинейности; и — напряжение на сетке

лампы.

Разрешая систему уравнений (5.70) относительно напряжения, получаем нелинейное дифференциальное уравнение второго порядка

- g ------к - g - + + i - V2 4 г Ф («) = “ о т (5-71)

где

 

 

 

 

 

 

к =

= со5 (MS-Я С );

« 5 = -^ -;

у2 =

3co^MSp; (5.72)

 

м_

d2ET

9 di

d-ET 2

+

£ (0 =

- ^ L

dt-

 

d?

 

 

+

dET.

0)qEt

 

 

(5.73)

 

dt

 

 

Нелинейность имеет вид кубической параболы ср (и) = и3. Вследствие малости параметра 0 можно считать, что MS ^ RC,

поэтому приближенно индуктивность М = RC/S. С учетом этого случайное возмущение (5.73) можно представить в следующем виде:

Е (/) = —4

PC

d-ET

р

dinrt

LS

dt2

' SL

dt “г

ш0

 

 

 

 

+

 

 

 

(5.74)

Для выполнения расчетов возмущение, описываемое формулой (5.74), можно рассматривать как белый шум с постоянной спектраль­ ной плотностью SE (со0). Для вычисления данной величины примем, что

(со) = 2е/0Г2;

STi (со) =

4kTR\

ST! =

4kTRc,

(5.75)

где Rc — сопротивление цепи обратной связи; е — заряд электрона;

Г — коэффициент депрессии

дробового

шума; k — постоянная

Больцмана; Т — температура

резистора,

одинаковая для R и Rc.

В соответствии с формулой (5.74) спектральная плотность суммарного возмущения

SE (со) = bQ+ &2со2 + Ь4со4,

(5.76)

156


где

 

 

 

b0 =

4kTRc\

b-2

1

R2

-2сой

4A77?C+ 2е/0Г2 ( - ^ - ) 21;

9

L2

 

“о

 

(5.77)

4£7V?c + 4jfe77?(-g-)2' .

Полагая в уравнении (5.76) со = со0, получаем следующее выра­ жение для уровня спектральной плотности:

$е К ) —

г 4kTRcR*

 

4kTR

jj2Q2 _j_

qT^R^C

(5.78)

9 г 9

I

S-L2

S*L

 

“oL~

 

 

Приближенно можно считать, что отношение сопротивлений

Rc___RC

 

(5.79)

R ~ SL

-

 

Подставляя значение R c из данного равенства в первый член фор­ мулы (5.78), получаем

 

SEК )

4kTR

Г,

 

Q2

e!aT-Q- I

(5.80)

 

SZQ

L

5Z

2kTS J ’

 

 

 

где Q =

Lu>0/R — добротность

контура; Z = L/7?C — резонансное

сопротивление контура. При Q =

100, / 0 =

10 мА, S = 3 мА-В-1,

Г2 = 0,2,

Т = 300° К соотношение

членов

в квадратной скобке

выражения (5.80) составляет 1 : 333 : 1,3-10s. Следовательно, основ­ ной вклад в спектральную плотность вносит дробовой шум. Если учитывать только дробовой шум, то уровень спектральной плот­ ности составляет SE (со0) = 2eI0r 2R/S2Z ^ 7- 10~21В2-с.

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

ср (и) — k xii,

 

(5.81)

где k\ = За2 — коэффициент статистической

линеаризации

(при­

ложение 3, нелинейность 6); а~и — дисперсия напряжения.

 

Подставляя выражение (5.81) в уравнение

(5.71), приведем его

к следующему виду:

 

 

[р~ + (Т2°н — V.) р + щ] и =

®1Е,

(5.82)

где р = d/dt.

 

 

157


Вычислим дисперсию напряжения по известной формуле

а~„ — щБе

</й>

(5.83)

(у"0н —* ) ' ш + “ о |2'

| (im)" +

 

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

а" = "2у^^ ^

^ ^

(5.84)

где введены параметры

 

 

(.1 = 2соо/х2;

v — 2я5яу2.

(5.85)

Спектральная плотность переменного напряжения на сетке лампы

определяется формулой

 

Su (co) = |Ф M | 2Se,

(5.86)

где частотная характеристика линеаризованной системы

 

Ф (tco) = ___________ “ о___________

(5.87)

(ш)'2 + (у2а2 —к) /со + С05

 

Вычисляя квадрат модуля частотной характеристики и подстав­ ляя значение дисперсии напряжения в соответствии с формулой (5.86), получаем следующее выражение для спектральной плотности сеточ­ ного напряжения генератора:

 

2 я у 3 / с о \ 4 2/ ш

1

 

(5.88)

5 Ы(<*>)

1

(1 -1 Л + vp)2

+ 1

 

 

 

 

График этой зависимо­ сти при j . i = l и v = var представлен на рис. 5.9. Спектральная плотность имеет максимумы на ча­ стотах ± ыт. Дифферен­ цируя спектральную плот­ ность по частоте, прирав­ нивая производную к ну­ лю и решая уравнение от­ носительно частоты экст-

Рис. 5.9. График спектральной плот­ ности сеточного напряжения гене­ ратора

158


ремума, получаем

 

 

 

± Ш0 ] / 1—

(1 — / 1 + Vfj,)2

V < 4 (1 — У |р ) *

со,,

 

 

(5.89)

 

0

^

4 ( ' - у г ) '

Максимальная величина спектральной

плотности

Sa {®т)

V

 

(5.90)

2яу3

 

 

V 1+ vp)2 2

Естественную ширину спектральной линии генератора, обу­ словленную собственными шумами, можно определить, приравни­ вая спектральную плотность к определенной величине. Удобно в ка­ честве такой величины рассматривать половину максимального зна­ чения спектральной плотности. Приравнивая выражение (5.88) к половине максимального значения Su (со,,,), получаем биквадратное уравнение, определяющее четыре значения частоты, соответствую­ щие точкам пересечения кривой спектральной плотности с прямой, параллельной оси абсцисс:

г4 — 2z2

- 1

+

+ 2

= 0,

(5.91)

где z = co/co0.

Поскольку спектральная плотность симметрична относительно оси ординат, то для определения уширения спектральной линии

генератора достаточно

рассмотреть

только

разность

 

двух корней

 

 

AQ = z2=

 

 

 

 

= l / " 1 - ' i (l

+

v^ 2+ V

1

И 1 — Z

1 +

vn)2]'

— If

/ l

+ v|x)2— ] / ! — [■! —

(1— Z l

+

vji)2] 2.

 

 

 

 

 

 

 

(5.92)

При уменьшении интенсивности помехи резонансная частота спектральной плотности стремится к частоте со0, а высота пика воз­ растает. Спектральная плотность имеет вид узкого пика с пьедеста­ лом небольшого уровня. При увеличении интенсивности шума ре­ зонансная частота колебаний генератора отклоняется в сторону меньших частот от частоты <в0. Одновременно при этом увеличи­ вается ширина спектральной линии. При отсутствии шума спектраль­ ная плотность вырождается в 8-функцию на частотах ±со0.

159


5 .5 . С истем а ст а б и л и за ц и и угл а крена

Линейная модель системы стабилизации летательного аппарата по крену была рассмотрена в п. 3.5. В действительности углы откло­ нения элеронов всегда ограничены. Кроме того, вследствие конечной мощности привода ограничена также максимальная скорость вра­ щения элеронов. Наличие ограничений приводит к необходимости применения нелинейной модели, что особенно существенно при боль­ ших возмущающих моментах, действующих на летательный аппарат. Расчет системы стабилизации при использовании линейной модели в этом случае значительно отличается от результатов анализа нели­ нейной системы.

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

У + ахху

= — ахэ ср (б) + X (0,

(5.93)

где ф (б) — характеристика

нелинейного элемента,

учитывающая

ограничение по углу отклонения элеронов.

 

 

5 < ~ d ,

 

Ф(6)

| 6 | < d ,

(5.94)

/8 > d .

Вуравнении (5.93) ахх, ахъ — коэффициенты, характеризующие динамические свойства объекта; X (t) — внешнее возмущение в еди­ ницах углового ускорения, имеющее нулевое математическое ожи­ дание и постоянную спектральную плотность Sx.

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

Гб + б = /е6 (кгу + k 2y),

. (5.95)

где Г, /еб — параметры привода; k x, k 2— параметры

измерителей

и усилителей. Структурная схема системы стабилизации приведена на рис. 5.10.

Рис. 5.10. ьСхсма системы стабилизации крена

160