Sg (ω)= Wопт(jω)Sg+n(ω),
откуда
Wопт (jω)= |
Sg (ω) |
|
|
|||||
Sg +n (ω) |
|
|||||||
|
|
|
|
|||||
или в случае независимости сигнала и шума |
|
|
|
|
||||
Wопт (jω)= |
|
|
Sg (ω) |
|
|
|
. |
|
|
|
|
|
|
|
|||
Sg |
(ω)+ Sn |
(ω) |
||||||
|
|
|||||||
(3.2.18)
(3.2.19)
Применив к обеим частям выражения (3.2.19) преобразование Фурье, получим
w |
(t)= |
1 |
∞W |
(jω)e jωtdω = |
1 |
∞ |
Sg (ω) |
|
e jωtdω . (3.2.20) |
|
|
|
|
||||||||
опт.н |
|
2π −∞∫ |
опт |
|
2π −∫∞ Sg (ω)+ Sn (ω) |
|
||||
|
|
|
|
|
||||||
К сожалению, вычисленная по формуле (3.2.20) весовая функция физически нереализуема, поскольку согласно правой части выражения (3.2.19) является четной, то есть не равной нулю при t<0. Следовательно, нереализуема и передаточная функция (3.2.19).
Строгое решение уравнения Винера−Хопфа (3.2.17) слишком громоздко, поэтому ограничимся нестрогим выводом, воспользовавшись результатами работ Г. Боде и К. Шеннона.
Представим спектральную плотность суммы стохастических процессов на входе системы Sg+n (ω) в виде произведения двух сопряжённых
функций (в большинстве случаев такое возможно)
181
Sg+n (ω)= ψ(jω)ψ(− jω). |
(3.2.21) |
Представление спектральной плотности в виде такого произведения называется факторизацией.
Тогда, в соответствии с формулой (3.2.6) можно считать, что воздействие на входе системы является результатом прохождения белого шума
сединичной интенсивностью через некоторое фиктивное линейное звено
счастотной передаточной функцией ψ(jω).
Доказано, что физически реализуемо лишь то звено, у которого полюсы частотной передаточной функции ψ(jω) расположены в верхней части комплексной плоскости. В дальнейшем считаем, что для функции
ψ(jω) это условие выполняется. Тогда у сопряженной функции ψ(− jω) полюсы расположены в нижней части комплексной плоскости.
Частотная передаточная функция оптимальной системы, полученная в виде (3.2.18) или (3.2.19), физически нереализуема, поэтому обозначим её Wопт.н (jω) («н» − нереализуемая).
С учетом вышеизложенного структурная схема системы, учитывающая передаточную функцию ψ(jω), приведена на рис. 3.5.
|
Wн(jω) |
||||
Белый |
|
|
g(t)+n(t) |
|
y(t) |
|
ψ(jω |
Wопт.н(jω) |
|||
шум |
|
|
Sg+n(ω) |
|
|
|
|
|
|
|
|
Рис. 3.5. Оптимальная нереализуемая система
Общая передаточная функция Wн(jω) также, очевидно, нереализуема и имеет вид с учетом (3.2.18), (3.2.19) и (3.2.21)
182
W |
(jω)= ψ(jω)W |
(jω)= ψ(jω) |
Sg (ω) |
|
= |
Sg (ω) |
|
|
Sg+n (ω) |
ψ(− jω) |
|||||||
н |
опт.н |
|
|
|||||
|
|
|
|
|||||
Функция веса, соответствующая этой передаточной функции, тоже, конечно, нереализуема и равна
|
1 ∞ |
jωt |
|
1 ∞ |
Sg (ω) |
jωt |
|
|
|||
wн (t)= |
|
−∞∫Wн (jω)e |
|
dω = |
|
−∞∫ |
|
e |
|
dω . |
(3.2.22) |
2π |
|
2π |
ψ(− jω) |
|
|||||||
Реализуемая весовая функция wр (t) совпадает с нереализуемой функцией (3.2.22) при t>0 и тождественно равна нулю при t<0
wр |
(t)= wн (t) |
при t > 0, |
|
0 |
при t < 0. |
Фурье-изображение весовой функции wр (t)есть частотная передаточная функция реализуемой системы, и выражается формулой
∞ |
∞ |
∞ |
|
1 |
∞ |
S |
g |
(ω) |
|
|
Wр (jω)= ∫wр (t)e− jωtdt = ∫wн (t)e− jωtdt = ∫e− jωt |
|
∫ |
|
|
|
e jωtdω dt . |
||||
|
|
|
|
|
||||||
−∞ |
0 |
0 |
|
2π −∞ ψ(− jω) |
|
|||||
Найденная передаточная функция состоит из последовательно соединённых звена с передаточной функцией ψ(jω) и оптимальной реализуемой системы, передаточная функция которой равна
W |
(jω)= |
Wр (jω) |
= |
1 |
∞ e− jωt |
|
1 |
∞ |
Sg (ω) |
e jωtdω dt . (3.2.23) |
|
|
|
|
|
∫ |
|
|
|||||
опт.р |
ψ(jω) |
|
|
∫ |
|
|
|
|
|||
|
|
|
ψ(jω) 0 |
|
2π −∞ ψ(− jω) |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
183 |
Применение обратного преобразования Фурье к последнему выражению дает функцию веса физически реализуемой системы
|
1 |
∞ |
|
|
wопт.р (t)= |
∫Wопт.р (jω)e jωtdω. |
(3.2.24) |
||
|
||||
|
2π −∞ |
|
||
Полученная весовая функция (3.2.24) является решением уравнения Винера−Хопфа (3.2.17) и удовлетворяет условиям физической реализуемости.
Вычисления по формуле (3.2.23) не настолько сложные, как это может показаться на первый взгляд. Рассуждаем следующим образом. Если бы
Sg |
было Фурье-изображением физически реализуемой весовой |
ψ(− jω) |
функции w′(t), то операция в квадратной скобке получает оригинал w′(t) из изображения, а в следующей операции − интегрирования от 0 до ∞ нижний предел можно поменять на −∞, поскольку w′(t)=0 при t<0, и мы опять получаем изображение из оригинала. Таким образом, обе операции компенсировали бы друг друга, и оптимальная передаточная функция
соответствовала бы выражению (3.2. 18). В действительности |
Sg |
не |
|
ψ(− jω) |
|
||
Фурье-изображение физически реализуемой весовой функции. Но оказывается, что для широкого круга задач можно разложить это выражение на простые дроби
Sg |
|
Sg |
|
Sg |
|
|
||
|
= |
|
|
+ |
|
|
, |
(3.2.25) |
ψ(− jω) |
ψ(− jω) |
ψ(− jω) |
||||||
|
|
|
+ |
|
|
− |
|
|
184
где первое слагаемое в правой части представляет собой функцию с полюсами в верхней комплексной полуплоскости, а второе – с полюсами в нижней полуплоскости1.
Первое слагаемое с полюсами в верхней полуплоскости является Фурье-изображением физически реализуемой функции w′(t), равной нулю при t<0. Можно показать, что второе слагаемое, имеющее полюсы только в нижней полуплоскости, является Фурье-изображением такой весовой функции w′′(t), которая равна нулю при t>0. Но тогда операция в квадратной скобке переводит второе слагаемое в w′′(t), а следующая операция – интегрирование от 0 до ∞ − в нуль, так как w′′(t)=0 при t>0. Поэтому после подстановки суммы (3.2.25) в формулу (3.2.23) первое слага-
емое после двух преобразований даст исходное выражение |
|
Sg (ω) |
, а |
|
|
|
|
||
ψ(− jω) |
||||
|
|
|
+ |
|
второе слагаемое даст нуль.
Окончательно выражение (3.2.23) для передаточной функции оптимальной реализуемой системы будет иметь вид
W |
(jω)= |
1 |
Sg |
. |
(3.2.26) |
|
|
|
|
||||
ψ(jω) |
|
|||||
опт.р |
|
ψ(− jω) |
|
|||
|
|
|
|
|
+ |
|
Сравнив полученное выражение с формулой (3.2.18) с учетом выражения (3.2.21)
W (jω)= |
1 |
|
S |
|
|
|
|
|
|
g |
|
, |
(3.2.27) |
||
|
|
|
|||||
опт.н |
ψ(jω) |
ψ(− jω) |
|
|
|||
|
|
|
|
|
|
|
|
1 Такая операция называется расщеплением.
185