Книги по разным темам Pages:     | 1 |   ...   | 11 | 12 | 13 | 14 | 15 |   ...   | 18 |

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

Одноканальная аппаратура требует применение регулируемого устройства задержки, являющегося весьма сложным в технической реализации. Другой недостаток одноканальной аппаратуры заключается в том, что с ее помощью осуществляется последовательная оценка значений корреляционной функции (1),.., ( ). Это существенно увеличивает время анализа и, n R R x X следовательно, при помощи одноканальной аппаратуры невозможно осуществить оперативную оценку корреляционной функции.

Многоканальная аппаратура в принципе позволяет осуществить оперативный анализ, так как в ней оценки значений корреляционной функции получаются одновременно. Однако такая аппаратура требует применения в п раз больше устройств задержки, чем одноканальная (n-число каналов). Реализация большого числа устройств задержки является трудной технической задачей Кроме отмеченного, многоканальной аппаратуре присущ тот существенный недостаток, что оценка корреляционной функции получается в виде совокупности ее значений, что затрудняет ее использование для проведения корреляционного анализа объекта исследования (управления).

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

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

4.2 Аппроксимативная оценка корреляционной функции Основная идея аппроксимативного способа оценки корреляционных функций заключается в выборе модели M (, 0, 1,..., ) нормированной N корреляционной функции ( ) этого или иного критерия. Другими словами, X должна быть решена задача аппроксимации функции ( ) функцией X M = M (, 0, 1,..., ). В зависимости от того, что выбирается в качестве N критерия аппроксимации, получаются различные способы оценки корреляционной функции.

Наиболее распространенным на практике является способ, основанный на обеспечении минимума величины квадратической погрешности.

= ( ) - ( )]2 d. (4.3) M X [ Величина погрешности является функцией параметров модели 0, 1,..., N. Если величина имеет один экстремум в заданном диапазоне изменения параметров 0, 1,..., N и этот экстремум является минимум, то значения параметров модели могут быть определены из системы уравнений = m (4.4) (m = 0,1,..., N) Обеспечение одного экстремума (минимума) достигается выбором соответствующего вида модели и диапазона изменения ее параметров.

Подставляя в систему (4.4) величину из уравнения (4.3), получим систему уравнений в следующем виде:

M ( ) M ( ) ( ) d - M ( ) ( )d = X X. (4.5) X X 0 (m = 0,1,.., N) Введем обозначение M ( ) = hM (, 0, 1,.., ) = hm ( ) ; (4.6) N X M ( ) M ( ) d = m (0, 1,..., ) = m. (4.7) N X Тогда система уравнений (4.4) примет вид hm ( ) ( )d -m = X. (4.8) (m = 0,1,.., N) o o M[X (t) X (t - )] Поставив в выражение (4.7) ( ) =, после ряда X o M[X (t)] преобразований получим o o o M[[X (t){ ( ) X (t - )d - m X (t)}] = m h. (4.9) (m = 0,1,.., N) Соотношение (4.8) может рассматривать как алгоритм для оценки аппроксимативным способом корреляционной функции. Такая аппаратура будет состоять из N+1 идентичных каналов. На рисунке 4.2 в качестве примера показана блок-схема m-го канала. Этот канал включает в свой состав фильтр с регулируемыми параметрами 0, 1,..., N с импульсной характеристикой hm ( ) (4.6), линейный преобразователь с коэффициентом преобразования m (4.7), функционально зависящим от параметров фильтра, вычитающее устройство (ВУ), множительное устройство (МУ), блок усреднения (БУ) и нуль - индикатор (НИ).

Процесс оценки параметров модели корреляционной функции осуществляется так. При подаче на вход аппаратуры анализируемого o процесса X (t) на выходах блоков усреднения каналов будут действовать сигналы, являющиеся оценками левых частей соответствующих уравнений системы (4.8). Путем вариации параметров фильтров добиваются нулевых значений выходных сигналов блоков усреднения всех каналов. Полученные значения параметров принимаются в качестве оценок параметров корреляционной функции исследуемого процесса.

Различают два типа модели корреляционной функции. Один из типов моделей представляет ряд по той или иной системе функций. Ко второму типу относят модели, представляющие функции заданного вида.

o X (t) 0, 1...

N ВУ m БУ НИ МУ Рисунок 34 - Блок - схема устройства для оценки корреляционной функции по аппроксимативному способу Последний тип характерен тем, что модель корреляционной функции является нелинейной функцией параметров 0, 1,..., N. Аппаратура, осуществляющая оценку параметров такой модели, имеет взаимосвязанные каналы, что при большом числе параметров приводит к существенному увеличению длительности анализа. Поэтому с целью повышения оперативности оценки параметров модели, представляющих собой функции заданного вида, целесообразно брать модели с небольшим числом параметров, не превышающем двух, трех.

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

M 1( ) = exp- ;

M 2 ( ) = exp- cos 1 ;

M 3 ( ) = exp- (cos 1 + sin 1 ) ;

M 4 ( ) = exp- (cos 1 - sin 1 ).

Как конкретно проводится синтез аппаратуры, осуществляющей аппроксимацию корреляционной функции функцией заданного вида, рассмотрим на примере.

Пусть на основании каких-то соображений в качестве модели корреляционной функции выбрана следующая:

M 1( ) = exp-. (4.10) Эта модель содержит неизвестный параметр 0, который должен быть определен на основании критерия (4.3).

Для данной модели из соотношений (4.5) и (4.6) находим h0 ( ) = - exp-, 0 = -.

С учетом определенных величин система уравнений (4.9) для данного случая примет вид o o o M[[X (t){- exp- X (t - )d + X (t)}] = или o o o M[[X (t){4 0 exp- X (t - )d - X (t)}] = 0. (4.11) На рисунке 35 представлена блок - схема аппаратуры для оценки параметра модели 0 синтезированная на основании (4.11).

Эта аппаратура состоит из двух одинаковых фильтров нижних частот с передаточными функциями с регулируемыми постоянными времени (1 + Tp) T =, вычитающего устройства (ВУ), множительного устройства (МУ), (1 + Tp) блока усреднения (БУ) и нуль - индикатора (НИ). Заметим, что в данной аппаратуре отсутствует линейный преобразователь с регулируемым коэффициентом преобразования, что обусловлено видом модели. Процесс оценки вида параметра 0, может быть легко автоматизирован, если применить устройство, осуществляющее регулировку постоянных времени фильтров на входе в зависимости от величины выходного сигнала блока усреднения.

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

Модели, представляющие собой ряды по той или иной системе функций, применяются тогда, когда имеется много априорных сведений о корреляционных свойствах анализируемого процесса. Эта модель имеет вид M ( ) = ( ). (4.12) K K K =2 (1 + Tp) (1 + Tp) ВУ БУ НИ МУ Рисунок 35 - Блок - схема аппаратуры для оценки параметра модели (4.11) корреляционной функции процесса.

Характерной ее особенностью является линейная зависимость от параметров 0,1,...,N.

Наиболее просто осуществлять оценку параметров модели (4.12) можно, если в качестве функции ( ) выбрать ортогональные, K удовлетворяющие условию 0, k m (4.13) K m ( ) ( )d =, k = m.

m В этом случае выражения (4.6), (4.7) примут вид hm ( ) = ( ), (4.14) m m = mm, (4.15) а система уравнений (4.9) преобразуется к виду o o o M[[X (t){ ( ) X (t - )d - mm X (t)}] = m. (4.16) (m = 0,1,.., N) Каждое из уравнений этой системы является линейным уравнением с одним неизвестным. В результате решения получим o o m = M[X (t) ( ) X (t - )d ]. (4.17) m m x o Если анализируемый процесс X (t) предварительно пронормирован по дисперсии то формула (4.17) примет вид o o m = M[X (t) ( ) X (t - )d ], (4.18) H m H m o o где X (t) = X (t).

H x Блок - схема аппаратуры, реализующей алгоритм (4.18), показана на рисунке 36. Аппаратура состоит из нормализатора (Н), фильтра с импульсной переходной характеристикой ( ), множительного устройства (МУ), блока m усреднения (БУ) и масштабирующего преобразователя с коэффициентом передачи. Заметим, что если из состава аппаратуры исключить m нормализатор, то будет осуществляться оценка параметров модели ненормализованной корреляционной функции RX ( ). На рисунке 4.показана лишь один из каналов. Всего таких каналов должно быть N+1.

Основным элементом каждого канала, определяющим сложность аппаратуры в целом, является фильтр с импульсной характеристикой ( ).

m o X (t) m У У m Рисунок 36 - Блок - схема аппаратуры, реализующей алгоритм (4.18) Для того, чтобы реализовать фильтр с такой импульсной характеристикой, должно быть выполнено условие ( )d < для любого m m. Это условие необходимо иметь в виду при выборе системы ортогональных функций. Кроме того, для простоты технической реализации аппаратуры желательно выбирать такую систему ортогональных функций, которая приводила бы к канонической структуре фильтров, когда в составе аппаратуры используется лишь один фильтр с импульсной характеристикой ( ). Остальные же фильтры оказываются составными частями указанного N фильтра. Для получения канонической структуры фильтров должно быть выполнено условие t (t) = (t - )m ( )d, (4.19) m+1 m где m ( ) - импульсная характеристика элементарного звена фильтра.

Если выбранная система ортогональных функций удовлетворяет условию (4.19), то блок - схема аппаратуры для оценки параметров модели будет иметь вид, показанный на рисунке 4.5. Не поясняя состава блоков аппаратуры, который понятен из рисунка, скажем лишь, что она включает в состав один фильтр с импульсной характеристикой ( ) и N одинаковых фильтров с импульсными характеристиками ( ). Все эти последовательно включенные фильтры образуют фильтр с импульсной характеристикой ( ).

N Широкое применение для оперативного анализа корреляционных функций получили модели, представляющие собой ряды по системе ортогональных функций Лагерра ( ( ) = L( ) ).

m m q q Lm ( ) = exp (4.20) (-1) Cm q!( )2.

q=Для этих функций условие (4.13) имеет вид 0, k m Lm ( )Lk ( )d = 1, (4.21), k = m ЕЕЕЕ (t) (t) (t) 0 1 N ЕЕЕЕ МУ МУ МУ ЕЕЕЕ БУ БУ БУ ЕЕЕЕ 1 1 ЕЕЕЕ 0 1 N 1 N Рисунок 37 - Блок - схема аппаратуры для оценки параметров модели корреляционной функции (канонического вида) а параметры модели, как следует из (4.18) определяются по формуле o o m = M[X (t) ( ) X (t - )d ] = Rx ( )Lm ( )d. (4.22) m L 0 Применение модели для аппроксимативной оценки корреляционных функций позволяет просто реализовать аппаратуру, во - первых, благодаря тому, что выполняется условие (4.13), что приводит к канонической структуре фильтров, а во - вторых, в следствии того, что элементарная ячейка такого фильтра представляет собой простейшее звено с передаточной функцией первого порядка p.

p + Погрешность аппроксимации корреляционной функции при использовании функции Лагерра, как это следует из выражений (4.3), (4.20) и (4.21), равна N = Rx ( )d -. (4.23) m m=Отсюда видно, что при прочих равных условиях величина погрешности аппроксимации зависит от величины параметра функции Лагерра.

Принимая во внимание формулы (4.20) и (4.22), можно показать, что при =0 и при = величина =, т.е. существует оптимальное значение, при котором минимальна.

Найдем это оптимальное значение. Рассматривая как функцию, получим следующее условие минимума величины :

=. (4.24) > Из уравнения (4.23) находим N N m 1 = -. (4.25) m m m=0 m=Из выражений (4.22) и (4.20) имеем m m 1 m + = - m-1 + m + m+1. (4.26) 2 2 m Подставив из уравнения (4.26) в уравнение (4.25), получим N + = - +. (4.27) N N + Отсюда видно, что оптимальное значение параметра должно находиться из уравнения = 0. (4.28) N +При этом должно быть выбрано такое решение уравнения, при котором > 0. Последнее условие с учетом (4.26) - (4.28), примет вид (N + 1) > (N + 2). (4.29) N N +В техническом аспекте задача выбора оптимального значения параметра функции Лагерра сводится к оптимальной подстройке постоянной времени фильтра Лагерра. Эта подстройка должна осуществляться в соответствии с алгоритмами (4.28), (4.29).

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

4.3 Оценка моментов корреляционной функции Корреляционная функция всесторонне характеризует корреляционные свойства анализируемого процесса. Однако во многих случаях достаточно ограничиться частичными сведениями о корреляционных свойствах сигнала.

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

Моментом корреляционной функции порядка К называется величина k = Rx ( )d. (4.30) k Применяются также и нормированные моменты k = ( )d. (4.31) н.k x Отличие последней характеристики от предыдущей состоит лишь в том, что она определяется через нормированную корреляционную функцию.

2 При этом, так как Rx ( ) = ( ), то =.

x x k x н.k Как видно из формулы (4.30), определение величины требует знания k корреляционной функции анализируемого процесса. В то же время эти характеристики применяются именно для того, чтобы описать процесс без знания корреляционной функции. Отсюда возникает необходимость поиска, с целью повышения оперативности анализа, таких способов оценки моментных характеристик корреляционных функций, которые не требовали бы предварительной оценки корреляционной функции.

Pages:     | 1 |   ...   | 11 | 12 | 13 | 14 | 15 |   ...   | 18 |    Книги по разным темам