WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     | 1 |   ...   | 15 | 16 || 18 |

Окно Пугачева – Даниэля представляет собой прямоугольник в частотной области и определяется выражением 2 w,| w | w g(w) =, (5.73) 0,| w |> w соответствующее корреляционное окно имеет вид w sin h( ) =. (5.74) w На рисунке 23 изображены графики функций ( 5.73 ) и (5.74 ) Очевидно, что использование этого окна эквивалентно полосовой фильтрации и требует большого объема вычислительных процедур.

Оценка СПМ при этом имеет вид w T sin( ) ^ €( S w ) = R ( ) cos w d. (5.75) w Погрешности оценивания в этом случае довольно велики, особенно для процессов с АКФ, обладающих низкой колебательностью.

Оценка сверху относительной среднеквадратической погрешности оценки Пугачева - Даниэля определяется выражением m 2 m cm = + = + (5.76) cm T T S (w) или 2 m = + (5.77) 2 R( )cos(w )d.

T S(w) 3 m Для случайного сигнала с АКФ R(w) = e- | | cos w0 (5.78) и СПМ:

S(w) = + (5.79) 2 2 + (w + w0 )2 + (w + w0 ) выражение (5.77) имеет вид 2 m = + C (5.80) 4 T m где 1- 3( + x)2 1- 3( + x) C = + (5.81) [1+ ( + x)2 ]1+ ( - x) здесь = ; x =.

Из (5.80) видно, что для уменьшения дисперсии оценки (5.75) необходимо уменьшать ширину корреляционного окна =, но при этом m растет смещение. Зависимость (5.68) носит параболический характер и, следовательно, имеет единственный экстремум – минимум.

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

m Оптимальная ширина корреляционного окна (5.74) находится из условия существования экстремума следующим образом:

2 4 4 C T = 0;.... =, (5.82) m m где - оптимальное, в смысле среднеквадратического критерия, m значение ширины окна.

Погрешность оптимизированной оценки Пугачева – Даниэля имеет вид (2.1+ 0.49) C. (5.83) (T ) Окно Барлетта представляет собой треугольник во временной области 1- ;...

m h( ) = (5.84) m 0;... > m Соответствующее спектральное окно имеет вид:

sin2 ( ) m g() = (5.85) m ( )m а спектральная оценка:

m € € S() = (5.86) R( )1- )d.

cos( 0 m Графики функций (5.84) и (5.85) показаны на рисунке 43.

Как видно из (5.85) и рисунка 43, спектральная оценка имеет выраженные боковые лепестки, что крайне не желательно при спектральном анализе. Это объясняется тем, что спектр исследуемого сигнала, и особенно его слабые составляющие искажаются, так как боковые максимумы лепестков «маскируют» пики слабых гармоник обрабатываемого сигнала.

h( ) g() - m m Рисунок 44 – Корреляционное и спектральное окна Барлетта Среднеквадратическая погрешность оценки определяется выражением:

2 2 2 2 m m = + C. (5.87) R()cos( )d = + 3 T 3 T m Оптимальная ширина окна определяется соотношением:

m 8TC = и с учетом этого m 0. = (1.61+ 0.67) C2. (5.88) (T )Не смотря на простоту технической реализации, окно Барлетта довольно редко используется в современном спектральном анализе из-за явно выраженных боковых лепестков и сравнительно низких метрологических характеристиках.

Окно Хэмминга описывается во временной во временной и частотной областях уравнениями 0.54 + 0.46cos, m h( ) = (5.89) m 0, > m sin(m) sin(m - ) sin(m - ) g() = 1.08 + 0.46 + 0.46. (5.90) m (m + ) (m - ) Графики корреляционного и спектрального окон Хэмминга изображены на рисунке 44:

h( ) g( ) - m m 2m m Рисунок 45 – Корреляционное и спектральные окна Хэмминга Погрешность определяется выражением:

m 2 = 0.8 +, (5.91) 4.35 2S() R( )cos( )d T m = (0.8 + 0.5) C. (5.92) (T )Оценка, получаемая при помощи окна Тьюки, дает результаты несколько лучше, чем оценка Хэмминга. Окна Тьюки во временной и частотных областях описываются следующими соотношениями:

1+ ;

2 cos m, m h( ) = (5.93) 0; > m sin( ) sin( + ) sin( - ) m m m m m g() = + + (5.94) m 2 ( + ) 2 ( - ) m m m а оценка СПМ:

m € S() = R( )1+ cos ( )d. (5.95) cos m Графики функций (5.93) и (5.94) приведены на рисунке 45. Боковые экстремумы здесь того же порядка, что и у оценки Хэмминга, хотя погрешности оценивания меньше.

2 m = 0.75 + (5.96) 2 2 R( )cos( )d, T m = (0.95 + 0.25) C2.

(T ) h( ) g() - m m 2m m Рисунок 46 – Корреляционное и спектральное окна Тьюки Окно Парзена /4/ во временной и частотных областях определяется соотношениями:

1- 6 2 3 m + 6, m m 2 2 m < m, 1- h( ) =, m 0, > m m sin g() = m 4. (5.97) m С помощью этого окна синтезируется спектральная оценка с наименьшими боковыми максимумами (менее 0.2 % основного экстремума), хотя по своим метрологическим характеристикам она и уступает оценки Тьюки.

Спектральные оценки без боковых лепестков и более высокими метрологическими показателями по сравнению с рассмотренными могут быть получены при использовании корреляционного окна /7/ q h( ) = (5.98) Ak Qk ( ), k = где Ак – параметры окна (коэффициенты разложения);

q – порядок окна;

Qk( ) – базисные функции.

В качестве базисных функций, отвечающих требованиям:

- отсутствие «отсечки» во временной области, типичных для всех рассмотренных выше окон и неизбежно ведущей к появлению боковых максимумов;

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

В качестве таких функций в /7/ предложены следующие:

k Qk ( ) = (5.99) e- k! Qk ( ) = (5.100) e- (k = 1) Qk ( ) =. (5.101) e(k +1) Параметры Ак корреляционного окна выбираются в ходе процедуры оптимизации /7/, который может производиться на основании одного из следующих критериев:

- минимума среднеквадратической погрешности оценивания;

- минимума интегральной среднеквадратической погрешности;

- абсолютной безлепестковости получаемых оценок СПМ.

6 Методы оценки законов распределения составляющих объекта исследования Знание закона распределения составляющих объекта исследования необходимо, с одной стороны, для его информационного описания, а с другой – для обоснования и правильного выбора динамического диапазона первичных измерительных преобразователей информационноизмерительный систем и измерительно-вычислительных комплексов.

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

Основное применение получили:

- непосредственная оценка;

- аппроксимативные способы оценивания.

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

F( ) = P{X (t) } x x 0 или для простоты, F( ) = P{X } x x 0 x F( ) = f (x)dx = (x) f (x)dx, (6.1) x - здесь (x) изменяет пределы интегрирования:

1,- < x x (x) =, x 0, x > или (x) = 1(- < x ), x F( ) = M[ (t)]; {X (t)} = Z(t), x где Z(t) – выходной сигнал измерительного устройства.

Далее заменяем оператор математического ожидания на оператор усреднения:

€ € € F( ) = M[Z(t)] = M[{X (t)}] (6.2) x Соотношение (6.2) будем использовать как алгоритм для синтеза измерительного устройства, структурная схема которого приведена на рисунке 46.

€ X(t) Z(t) ФП БУ F( ) x XРисунок 47 – Структурная схема устройства для оценки функции распределения На рисунке 46 ФП – функциональный преобразователь, БУ – блок усреднения. График функции (x) изображен на рисунке 47.

(x) x XРисунок 48 – График функции преобразования ФП График, показывающий преобразование случайного входного сигнала преобразователем, представлен на рисунке 48.

X t t Рисунок 49 – График преобразования случайного входного сигнала преобразователем Оценка может быть несмещенной, ее методическая статистическая погрешность определяется величиной:

kz функциональный преобразователь ФП осуществляет = c, ст T нелинейное преобразование, поэтому спектр сигнала на его выходе расширяется по сравнению со спектром входного сигнала, а интервал корреляции – уменьшается. Поэтому kz, тогда kx kx c.

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

Функциональный преобразователь ФП представляет собой пороговое устройство с порогом X.

До сих пор мы рассматривали задачу оценивания только одного значения F(x) при x=x0, нам же необходимо оценить функцию распределения F(x) в диапазоне значений X.

Один из способов такой оценки заключается в замене F(x) совокупностью ее значения, взятые через равные интервалы x.

Для грубой оценки берут 20-30 значений F(x), для точной 200-300.

Анализ функции распределения может быть последовательным или параллельным.

При последовательном анализе аппаратура одноканальная, последовательно изменяя x0, (см. рисунок 50) оцениваем F(x) в диапазоне.

Сложность технической реализации аппаратуры определяется числом каналов C, здесь же C=1, но длительность анализа при этом максимальна:

если для оценки одного значения требуется время T, то для оценивания N значений будет необходимо затратить время Тэ=NT € F(x) x xm x xb Рисунок 50 – Оценивание функции распределения F(x) в диапазоне.

При параллельном анализе аппаратура многоканальна, все значения F(x) оцениваются одновременно(см.рисунок 51).

Сложность технической реализации здесь наибольшая: C=N, зато время минимально: Tэ = T.

XФП БУ ) F(x1) X(t) ФП БУ ) F(x2) ФП БУ ) F(x3) x N Рисунок 51 – Параллельный анализ функции распределения. Схема аппаратуры.

6.2 Непосредственный способ оценки плотности вероятности Для того, чтобы оценить плотность вероятности случайного процесса f(x) попытаемся представить ее в виде математического ожидания некоторого сигнала. И так, имеем стационарный случайный процесс X(t), нужно определить f(x) при некотором конкретном значении x=x Функция распределения: F( ) = x 0 (x) f (x)dx;

F (x)f (x)dx = M (x) f (x0)= F(x0)= = (6.3).

x0 - x0 x Оценку получим заменой в соотношении (6.3) оператора математического ожидания на оператор усреднения:

{X (t)} € f€(x0)= M (6.4) x0.

Структура устройства для определения оценки плотности вероятности показана на рисунке 51.

(x) €( X(t) f ) x xУ (x) xРисунок 52 – Структурная схема устройства для измерения плотности вероятности по алгоритму (64) Однако, в точке x0 производная =, то есть функция xпреобразования представляет собой дельта-функцию, что технически не реализуемо. Надо искать какой-то другой функциональный преобразователь ФП, но при этом возникает погрешность от смещенности.

Плотность вероятности – это такая вероятностная характеристика, оценка которой всегда принципиально смещена из-за нереализуемости алгоритма.

Представим значение плотности вероятности в точке x=x0 в виде x x F x0 + - F x0 - F(x0 + x)- F(x0 ) 2 f (x0 ) = F(x0 ) = lim = lim, (6.5) x x x 0 x F(x0 + x)- F(x0 ); lim (x0 ) = f (x0 ).

(x0 ) = (6.6) x x То есть, если оценивать (x0 ) и брать достаточно малую x, то можно получить оценку плотности распределения f€(x0 ) с достаточно высокими метрологическими характеристиками.

График функций (x0 )и (x0 + x)- (x0 ) изображены на рисунке 52.

x0 x0+ x x x0 x0+ x x Рисунок 53 - График функции преобразования ФП для синтеза технически реализуемого фильтра X В качестве оценки плотности вероятности будем брать оценку € величины : (x0 ) = f€(x0 ).

€ € € € F(x0 + x)- F(x0 ) M[1(- < x x0 + x)] M[1(- < x x0 )] € (x0 )= = - = x x x (6.7) € M[1(- < x x0 + x)-1(- < x x0 )], = x € [1(x€ f€(x0 ) = x(x0 )M x x0 + x)]. (6.8) x € Если говорить о метрологических свойствах оценки f (x0 ), то выбором x можно сделать величину погрешности от смещенности сколь угодно см малой. Статистическая методическая погрешность определяется знакомым выражением:

k C см T 6.3 Аппроксимативные способы оценки плотности вероятности Здесь задача состоит в построение модели плотности распределения.

Для этого необходимо выбрать критерий адекватности и построить модель.

При этом имеют место три этапа:

1) выбирается вид модели;

fM (x, 0, 1,..., ), N где i -параметры модели, причем, чем больше объем априорной информации о f(x), тем точнее будет вид модели;

2) выбирается критерий адекватности, который ставит в соответствие модель и саму плотность вероятности {fm (x), f (x)};

3) производится планирование и организация эксперимента по определению параметров модели 0, 1,..., обеспечивающих оптимум N критерия адекватности.

^ В настоящее время при аппроксимативном оценивании f (x) в основном используются три критерия адекватности:

1) критерий моментов;

2) критерий производных;

3) квадратический критерий.

Рассмотрим кратко особенности построения оценок f(x) по каждому из них.

6.3.1 Аппроксимативное оценивание плотности вероятности по критерию моментов Сущность критерия моментов состоит в следующем. Модель признается адекватной объекту (в данном случае - истинной плотности вероятности, если моменты сигнала, имеющего плотность распределения равную модели, равны соответствующим моментам исследуемого сигнала).

Количество моментов должно быть равно числу неизвестных параметров модели.

Возьмем начальные моменты K k kx[X ] = M[X ] = X f (x)dx - Нужно найти (N+1) моментов сигнала и модели. Начальный момент k-го порядка модели определится равенством km[X ] = xk fm(x)dx = k ( 0, 1,... ) (6.9) N - Итак, критерий моментов состоит в том, чтобы k km = kx, = 0, N (6.10) Чем больше N, тем более адекватной будет модель. Необходимо получить оценки начальных моментов (их N+1). Для этого заменим X оператор математического ожидания М оператор усреднения:

k € € km = M[X ] (6.11) и получаем случайную величину, которая, тем не менее, несет информацию о k-ом начальном моменте ak и связана с ним. Эта случайная величина представляет собой функцию от оценок параметров модели k (0, 1... ). Оценки параметров должны удовлетворять всем требованиям, N предъявляемым к статистическим оценкам (они должны быть состоятельными, несмещенными, эффективными).

Свойства параметров модели будут полностью определяться видом оператора усреднения.

Pages:     | 1 |   ...   | 15 | 16 || 18 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.