Рис. 9.3. Структурная схема формирующего фильтра
Из данной формулы, видно, что при GZN = 0 плотность вероят ности является симметричной, поскольку кумулянт третьего по рядка характеризует симметрию распределения вероятности.
Для реализации физической системы, соответствующей уравне нию (9.1), необходимо иметь два генератора шума ГШ1, ГШ2, блок произведения, усилители и интегратор. Функциональная схема реализации формирующего фильтра показана на рис. 9.3. Для созда ния коррелированных белых шумов сигналы генераторов белых шу мов суммируются в некоторой пропорции, определяемой коэффи циентом корреляции k 1. Генератор ГШ1 формирует мультиплика тивное возмущение Z (?). Генератор ГШ2 формирует белый шум V{t). Аддитивное возмущение формируется как сумма:
N (t) = k xZ (t) + V (t). |
(9.28) |
Интенсивность аддитивного белого шума определяется соотноше
нием |
|
GN = klGz + Gv, |
(9.29) |
где Gv — интенсивность шума генератора ГШ2. |
|
Взаимная интенсивность шумов N (t) и Z (() |
|
GZN = kxGz. |
(9.30) |
Таким образом, различная плотность вероятностей достигается |
изменением интенсивностей шумов генератора ГШ1 |
Gz , ГШ2 Gv |
и коэффициентом передачи kx. Параметры а и k также можно исполь зовать для изменения характеристик закона распределения выход
|
ной переменной |
фильтра. |
|
Для визуального контроля гра |
|
фика плотности вероятности мож |
|
но использовать устройство инди |
|
кации, блок-схема которого изо |
|
бражена на рис. 9.4. В схеме |
|
имеется два осциллографа 1 и 2, |
|
электрический мост с включенным |
|
в одно из плеч фоторезистором 3 и |
|
сглаживающая цепочка RC. Один |
Рис. 9.4. Блок-схема индикации |
из осциллографов |
(на рисунке ос |
циллограф 2) работает в режиме вынужденной пилообразной разверт ки по горизонтали. Этот же сигнал развертки подается на вертикаль ные пластины другого осциллографа. Сигнал развертки суммируется с наблюдаемым случайным процессом Y (t). На экране осциллографа 1 закрепляется фоторезистор, включенный в плечо моста. В резуль тате периодической развертки случайного процесса фоторезистор изменяет свое сопротивление пропорционально ординатам одно мерной плотности вероятности исследуемого случайного процесса. После сглаживания на интегрирующей цепочке RC напряжение, пропорциональное плотности вероятности, подается на вертикальные пластины осциллографа 2. Вследствие периодической горизонтальной развертки процесса на экране этого осциллографа воспроизводится график плотности вероятности процесса Y (t).
9.2. Система второго порядка
Математическое описание работы измерительных устройств при учете жесткости конструкции, маятниковых устройств при ускорен ных движениях точки подвеса, электрических цепей при флуктуа циях величин параметров во многих случаях приводит к дифферен циальному уравнению с аддитивными и мультипликативными воз мущениями вида
|
Y + (2|со0 + 2г (0) У + |
(со? + |
Zi (0) Y = N (t), |
(9.31) |
где |
Z2 (t), |
Zx (t) — параметрические |
возмущения; |
N (t) — адди |
тивное возмущение; g — коэффициент |
затухания; |
со0 — собствен |
ная |
частота |
системы. |
Z2 (t) |
' - ' |
|
коррели |
Случайные возмущения Z 2 (t), |
и N (t) являются |
рованными белыми шумами с нулевыми средними, интенсивностями
Gi, G2, Gn и взаимными интенсивностями GfN, G%n, G(2.
Составим уравнения для первых и вторых начальных моментов и дисперсии выходных переменных. Для этого предварительно вы
|
числим коэффициенты сноса |
и |
диффузии. Представим уравнение |
|
(9.31) в форме Коши: |
|
|
|
|
Yi |
= |
П; |
(9.32) |
|
Г2 = —(2Ы0+ Z2) Y2 - |
(w20+ Zi) Yi + N (t). |
|
|
Используя формулы (8.25), (8.26) коэффициентов сноса и диффу зии для линейных систем, получаем
А\ = У2', 'А2 = — ^2gcoo~— j-G2y«/2 —
/2 |
1r z\ „ |
1 r Z N . |
— ( ® . - т 0 , ) й - . 7 й , |
2BU = 0; |
2512 = |
2fJ21 = 0; |
2B22= GfУ? -(- 2Gi2!J\y2-f- G2y\ —
— 2 0 ? ^ — 2GlNy2+ Gn .
На основании уравнений (8.48) и (8.49) получаем следующие уравнения для математических ожиданий:
m i = m2; m2 = — ^ 2 gcoo ■— y G f ) m 2 — ^ <*>6------G f ^ m i ----------- |
G $ N |
(9.34)
идля вторых начальных моментов
«и = 2<х12>
|
Х12 = |
— ^2|С00---- i1-G f ) Й12— («о---- Y Gf:2^ ац — |
|
|
|
~2 |
|
|
|
|
|
|
|
|
---- Gfwmi |
агг', |
|
|
(9.35) |
|
|
агг = |
—2 ^2|соо---- Gf ) «22 — |
|
|
|
|
|
|
— |
2 ^соо----2" 0?2 ^ O12 — |
G i‘vni2 -f- Gfau -f- |
|
|
-(- 2G\2<xi2 -|- Gfo.22 — 2G\Ntill — 2G2Nm2 -)- G . |
|
|
В установившемся режиме |
|
|
|
|
|
|
m2 = а 1а = 0; mx= — |
GfN |
|
|
|
|
• Z N |
|
|
|
|
|
|
|
2o)0 — Gj2 |
|
|
|
a ll — |
GN (2o>q- |
Ofa) + |
G jN \2GfN + |
G jN ( 2jco0 - |
Of)] |
_ |
|
( 4 - |
°f2) [(2<b0S - °2 ) (2o)20 - Of2) - Gf] |
’ |
|
|
|
GN (2tog - |
Of2) + |
GfN [2Gfv + |
Of* (2Eo0 - |
Gf)] |
(9.36) |
|
a 22 |
2 |
[(2«Do| - |
G|) (2og - Cf2) - Gf] |
|
|
|
|
|
|
G f"
2 (2a»| — Gfa)
При отсутствии взаимной корреляции между параметрическими и аддитивными возмущениями моменты соответственно равны:
т, |
= т 9 = |
a 1 2 |
0; |
|
a n = |
2c0q(2co0| |
Gf) |
Gf ’ |
(9.37) |
|
|
|
|
gV
OC90—
(2ш0| — Gf) 2coq— Gf
Наконец, при отсутствии параметрических возмущений получаем
т1 — т2 « 1 2 = 0 ;- |
(9.38) |
gn |
«п — 6ц — 2t > «22 — ®2 |
4to0£ |
4a»Si |
|
Поскольку математическое ожидание производной Y в устано вившемся режиме равно нулю, то величина 022 есть дисперсия про изводной в установившемся режиме. Дисперсия координаты опре деляется вычитанием из второго начального момента квадрата ма тематического ожидания. Для установившегося режима
GN (2(og - Gfa) + 2Gf% 2ЛГ (2со* - Gf2) - GfO?"’
(9.39)
(2“ o — G f 2) 2 [(2co0£ — Of) (2(0q— Gf2) — Gf]
При отсутствии взаимной корреляции возмущений имеем
--N |
(9.40) |
Д ,= |
у 2со*(2со01 - G i) - |
Gf ’ . |
Вычислим теперь вторые начальные моменты. Эти моменты опре деляют по формулам (8.70):
|
|
k=I |
при ? |
> t, |
|
г „ ((, 0 |
= |
|
|
(9.41) |
|
|
1 |
при t |
> |
|
|
|
где gtk (t, О — весовые функции системы, а функции вычисляют по формулам
2 /'
F, (*', 0 = |
S |
J 8ik (Д Л) К (А ) + |
Gfn (h)] dh |
t > (, |
|
ft. i=i t |
|
|
|
2 |
t |
|
|
Fi (t, 0 = |
S |
J ^ (*. Л) К (Л) + |
G|h (A)] dh |
t > |
ft. /=i <'
Для рассматриваемого измерителя при (' > t начальные моменты соответственно равны:
Гц ((', 0 = |
g.n (/', t) 0ц (/) + |
£i2 (Д |
0 |
0U + |
m, (О Д (Г, /); |
|
г « ( Л |
о = £ я ((', |
о£ 2 20 и( Д(/)00+12 + |
mL{t) F 2 (t\ 0; |
(9.43) |
Г22 (/'. о = |
g2i ((', t) 021(О + |
£22 (/', |
0 |
022 + |
/па (О F2 (Г, t). |
|
Весовые функции определяются системой уравнений
(r>0 + 6 ( f - 0 ;
(9.44)
^ '■ „ 0 ^ g 22 ( /',/) ;
|
— ^2cOo£----g-G2^ gil (t , t)', |
|
ф ± = _ |
(ш*о _ 4 - |
(9.44) |
^ |
cf2) g a (t, t) - |
— |
(2co o £----- y |
Gf ) g 22 |
» О 4 - 6 — t ) . |
Дифференцируя четвертое уравнение системы (9.44) по f и исклю чая с помощью второго уравнения производную весовой функции
g 12 (t t ) , |
|
получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
d^ |
d% ' t] |
= |
- |
(со0 - 4 |
§ 2 2 {t, t) - |
|
|
|
|
|
|
- (2co<£---- TX |
Gf) |
+ |
6 (t |
- i). |
(9.45) |
|
Интегрируя это уравнение при t' > t и |
> Г и |
учитывая, что |
go2 |
(*', |
/) |
= |
g22 (t, |
|
t'), получаем |
|
|
|
|
|
|
|
g 22 (Л 0 = |
[4 (соI - |
4 |
0?2) - |
(2со0Е ---- y |
G?) |
] |
X |
|
|
|
XГ ^ |
|
^ |
0|)1 1 {- (2Ш,5 —i. Gf) X |
|
|
X |
sin |
j^-^- "|/ Г 4 |
( coq----- Gf2) — (2(o0| ; --------- ^ - G2 ) |
I * — * |
+ |
|
|
|
|
+ |
} ^ 4 |
(c o g - 4 |
Gf2) - |
(2co0| - |
4 |
G2 )2 X |
|
X |
cos |
4 |
] |
Л |
(co^----- G?2) |
— |
(2co0| |
--------- i - G 2 ) 2 1 / — t |
(9.46) |
или, |
вводя |
обозначение |
|
|
|
|
|
|
|
|
|
|
|
|
,= |
'|/ Л^ ( С0°-----g- Gf2) — (2соо^---- g-G2) |
, |
(9-47) |
будем |
иметь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
£ я (* '.0 = |
1 |
-4 -(2 ш0t - ±-в%)1 t-t' I |
|
|
|
|
|
|
4 е |
|
|
|
X |
|
|
|
X |
|
(2сооЕ ----- y |
G2 ) |
sin v U — ^ I + |
V COS V 1 1 — t |
(9.48) |
