Научная статья на тему 'Дисперсионные характеристики статических моделей стохастических объектов'

Дисперсионные характеристики статических моделей стохастических объектов Текст научной статьи по специальности «Строительство и архитектура»

CC BY
111
14
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДИСПЕРСИОННАЯ ИДЕНТИФИКАЦИЯ / СТАТИЧЕСКАЯ ПАРАМЕТРИЧЕСКАЯ МОДЕЛЬ / СТОХАСТИЧЕСКИЙ ОБЪЕКТ / НЕПАРАМЕТРИЧЕСКАЯ ОЦЕНКА РЕГРЕССИИ / DISPERSING IDENTIFICATION / STATIC PARAMETRICAL MODEL / STOCHASTIC OBJECT / A NONPARAMETRIC ESTIMATION FOR REGRESSION

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Рубан А. И.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The Dispersing Characteristics of Static Models of Stochastic Objects

The dispersing characteristics determining quality of static parametrical models of stochastic objects are constructed, and their estimations are found at the assumption, that the conditional variance of an output is constant. Parametrical models are nonlinear on the channel "input output" but they are linear concerning optimized parameters. It is shown that the nonparametric estimation of regression is close to objectively existing ideal model and according to dispersing characteristics it is possible to keep up with change of quality of suboptimal parametrical models at selection of their structure. Quality of nonparametric and parametrical models compared on several numerical tests with various level of noise.

Текст научной работы на тему «Дисперсионные характеристики статических моделей стохастических объектов»

УДК 681.5.015

Дисперсионные характеристики статических моделей стохастических объектов

А.И. Рубан*

Сибирский федеральный университет, 660041 Россия, Красноярск, Свободный, 791

Получена 15.05.2008, окончательный вариант 15.06.2008, принята к печати 10.03.2009 Построены дисперсионные характеристики, определяющие качество статических параметрических моделей стохастических объектов, и найдены их оценки при предположении, что условная дисперсия выхода постоянна. Параметрические модели, нелинейные по каналу "вход-выход", но линейные относительно оптимизируемых параметров. Показано, что непараметрическая оценка регрессии близка к объективно существующей идеальной модели и по оценкам дисперсионных характеристик можно следить за изменением качества субоптимальных параметрических моделей при подборе их структуры. На тестовых численных примерах произведено сравнение качества непараметрической и параметрической моделей при различном уровне помех.

Ключевые слова: дисперсионная идентификация, статическая параметрическая модель, стохастический объект, непараметрическая оценка регрессии.

Введение

Рассматриваем объект, который имеет один вход X и один выход У. Х и У являются непрерывными случайными величинами. Особое место среди идеализированных характеристик (моделей) стохастического объекта занимает регрессия

СЮ

М{У\х}== п(х) = J у/(у\х)3,у, (1)

которая представляет собой среднюю зависимость между выходом и входом объекта [1, 2].

Регрессия п(Х) удовлетворяет критерию минимума среднего квадратического отклонения выходов объекта и модели

П(Х) = а^шт Му,х{(У - п(Х))2}. (2)

п(х)

Минимальное значение оптимизируемого в (2) квадратичного показателя

А = Му,х {(У - п(Х ))2} (3)

названо в [1, 2] средней условной дисперсией. Там же введена средняя дисперсия регрессии

В2 = Мх{(п(Х) - ту)2} (4)

* Corresponding author E-mail address: [email protected]; [email protected]

1 © Siberian Federal University. All rights reserved

и приведена теорема о разложении дисперсии выхода объекта Пу = Му {(У — ту )2} на сумму дисперсий П\ и П2 :

Пу = Пх + П2 = Мух {(У — п(Х ))2} + Мх {(п(Х) — ту )2}. (5)

В дальнейшем для упрощения записи формул убираем индексы операторов математического ожидания М{•}, где были указаны случайные величины, по которым ведется усреднение.

На основе средней дисперсии регрессии П2 строится [1,2] безразмерная характеристика - дисперсионное отношение

Пу\х = V7 П2/Пу, (6)

которое принадлежит интервалу [0; 1]. Показано также, что дисперсионное отношение равно коэффициенту корреляции между выходом объекта и выходом модели в виде регрессии:

о М {[У — ту ][п(Х) — тп(х)]} (7)

Пу\х = Ху\п(х) =-----------------------------, (7)

ау °п(х)

где ту, тп(х) — математические ожидания, причём т^(х) = ту; ау = л/Пу, ап(х) = VП2

-- средние квадратические отклонения. Кроме того, квадрат нормированного коэффици-

М{(у—ту )(Х—тх)} лг V

ента корреляции гух = —^-ауаХ-между У и X не превосходит квадрата диспер-

сионного отношения

г\у < Пу\х ■ (8)

Дисперсионное отношение пу\х достигает предельных значений в следующих случаях.

Пух = 0, если выход объекта У не зависит от входа X. Обратное утверждение верно не всегда.

Пух = 1, если У и X связаны жесткой зависимостью в виде функции У = Э(Х). Тогда М{У\Х} = Э(Х). В этом случае другие воздействия (кроме X) не влияют на выход У и Пх = 0 (т. е. П2 = Пу).

Неравенство (8) позволило [1, 2] построить характеристику

ІУ\х

Лух2 Гух =4/1 - 4х (9)

Пу \ х V ПУ \ х

которая названа относительной степенью нелинейности по каналу вход X выход У и находится в интервале [0; 1]. По её величине можно судить о степени линейности (или нелинейности) истинных усреднённых характеристик объекта или, как говорят, "о линейности (или нелинейности) объекта". При 7у\х = 0 объект линейный, при 7у\х = 1 — нелинейный.

В работах [3, 4] на основе использования непараметрической оценки Розенблатта -— Парзена для плотностей распределения вероятностей [5, 6] построены сравнительно простые и эффективные оценки вышеуказанных дисперсионных характеристик.

Постановка задачи

При разработке моделей объектов стремятся построить хорошие — субоптимальные -параметрические нелинейные модели и количественно оценить их качество. При этом все вычисления должны быть выполнены на основе одной реализации синхронных измерений

входа и выхода объекта. Вычислительная схема количественного оценивания теоретических характеристик качества линейной (по каналу "вход-выход") параметрической модели на основе экспериментальных данных построена в [3, 4, 6]. В основе решения лежит использование состоятельной непараметрической оценки регрессии. В этой оценке важно "хорошо" настроить (для каждой конкретной реализации измерений входа и выхода объекта) так называемый коэффициент размытости. После его получения мы располагаем хорошей по своей сути (состоятельной "в статистическом смысле слова") моделью. Её в дальнейшем можно применять (и это часто с успехом делают) при решении достаточно сложных проблем, например, при синтезе эффективных алгоритмов: классификации в распознавании образов, адаптивного оптимального управления, глобальной оптимизации.

Наряду с несомненным достоинством непараметрической оценки регрессии (за счёт универсализма и состоятельности) она имеет недостатки. Первый из них состоит в том, что при расчёте её (даже для каждого фиксированного X = х значения входа) приходится непрерывно "перелопачивать" выборку (экспериментальные данные). Во-вторых, оценка представляет собой "универсальный вычислитель" , который позволяет практически для любого стохастического объекта на основе "синхронных" измерений его входа и выхода вычислять однозначную высокого качества статическую модель без привлечения априорной информации. Непараметрической модели в таком специфичном виде нельзя придать физический смысл и очень трудно в неё включить априорную информацию о ранее полученных параметрических моделях. Её трудно адаптировать к изменяющимся условиям, и она не обладает прогнозирующим свойством.

Человечество на протяжении всей своей истории с успехом применяет субоптимальные параметрические модели и не спешит заменять их другими моделями, да и острой необходимости в этом не существует. Параметрические модели удобны при адаптации к изменяющимся условиям и при прогнозировании свойств объектов. По этой причине возникает проблема количественного оценивания качества предлагаемых исследователем параметрических нелинейных моделей. В представленной работе предложен путь решения этой проблемы. Общая схема та же, что и при оценивании качества линейных (по каналу вход-выход объекта) параметрических моделей [3, 4, 6].

В предлагаемой работе построены теоретические статистические характеристики (дисперсионные характеристики) и вычислительные алгоритмы их оценивания (на основе измерений входа и выхода объекта). На численных примерах продемонстрирована правильность предложенного способа оценивания качества получаемых моделей.

Дисперсионные характеристики параметрических моделей

Построим дисперсионные характеристики нелинейных параметрических статических моделей стохастических объектов. Берём за основу широко распространённую параметрическую нелинейную (по каналу вход-выход) модель объекта с линейными параметрами

т

параметры ао, аі, ..., ат которой найдены из критерия минимума среднего квадратич-

т

ного отклонения выхода объекта У и модели ао + ^2 а 3 ^3 (X)

3=1 т

ао, а1, ..., ат = а^ шіп М {[У — а0 — а3 ^3 (X )]2}. (11)

а0,а1,...,ат

3=1

Модель (10) субоптимальная, т. е. она оптимальная (в смысле квадратичного критерия (11)) по параметрам в классе моделей выбранной структуры.

Как и для регрессии справедлива теорема о разложении дисперсии выхода объекта на две дисперсионные компоненты:

Иу = М {(у — у )2} + М {(у — ту)2 } = Иі,парам + И2,парам. (12)

Средняя дисперсия И2,парам = М{(у — ту)2} характеризует ту часть дисперсии Иу выхода объекта У,, которая вызвана влиянием входной переменной X, и это влияние описано с помощью субоптимальной модели (10).

Если объект полностью описан моделью (10) (т. е. нет дополнительных воздействий на объект и нет дополнительной стохастичности внутри объекта; в этом случае У = У), то Иі,парам =0 и И2,парам = Иу.

Если объект никак не описан моделью (10), то Иі,парам = Иу и одновременно И2, парам = °.

По аналогии с дисперсионным отношением (6), в котором использована оптимальная модель в виде регрессии п(Х) введём дисперсионное отношение, основанное на применении субоптимальной модели у (10),

Ш2, парам

Пу\х, парам — о —^{у }—' (13)

Оно также лежит в интервале [0; 1] и равно нормированному коэффициенту корреляции

Яуу между выходами объекта и модели (10) (Дуу называют часто множественным коэф-

фициентом корреляции):

М{[У - ту]у - ту]}

Пу\х, парам = Яу \у —---------оуа~------------------------------, (14)

где М{У} — ту, ау — ^Мх{[У - ту]2}.

Для дисперсионных отношений (6), (13) справедливо неравенство:

Пу\х,парам < Пу\х. (15)

Об улучшении качества субоптимальной параметрической модели У можно судить по приближению пу\х,парам к пу\X, что эквивалентно приближению ^2,парам к В2.

Введём коэффициент

, = Пу| х,парам _ /И2,парам _ М{[У — ту]2} )

?^ = пух V И2 = V М{(п(Х) — ту)2}’ ( )

который характеризует степень близости параметрической модели Y = ао+

m

+ ^2 a.j<pj(X) и регрессии n(X). Заметим, что Ху представляет собой усреднённую

j=i п

дисперсионную характеристику.

Если коэффициент Хуравен единице, то обе модели в среднем одинаково хорошо описывают объект, т. е. субоптимальная параметрическая модель (10) достигла своих предельных свойств — она лучшая.

Если коэффициент Хуменьше единицы, то параметрическая модель хуже регрессии и есть необходимость улучшать параметрическую модель за счёт замены одних базисных функций j (X) другими или добавления новых базисных функций.

Если даже Ху =1, то это не означает, что объект хорошо описан с помощью построенной модели Y. Можно лишь утверждать, что в среднем качество описания объекта на основе субоптимальной параметрической модели Y и на основе регрессии одинаковое. Если же при этом одинаковые дисперсионные отношения г/у|х,парам = Пу\х меньше единицы, то для получения лучшей модели необходимо кроме выбранных входов использовать дополнительные входные воздействия объекта и заново строить модель.

Используем также относительный коэффициент корреляции между выбранной параметрической моделью и регрессией:

^ = M{[Y - ту][n(X) - ту]} (17)

= /------~-------- /----------------- . ( )

у/М{[Y - ту]2}^ M{[n(X) - ту]2}

Модуль этого коэффициента корреляции |Жу| по своей сути эквивалентен коэффициенту Ху .

Непараметрическая оценка регрессии

Исходной информацией служит независимая выборка: (xi,yi), i = 1,n. Считаем, что условная дисперсия выхода объекта является постоянной (не зависит от входа): D{Ylx} = = Му{[Y — n(x)]2|X = x} = const. Эта гипотеза похожа на используемое в методе наименьших квадратов предположение, что измерения выхода объекта "равноточные".

Состоятельная оценка регрессии, построенная на основе использования непараметрической оценки условной плотности распределения вероятности [5, 6],

M{Ylx} = Vn(x)= ]Г Kn () yi, KN () = nK ч (18)

i=i ' ' ' ' y] K

j=1

относится к классу универсальных оценок. Здесь K(■) — ядро, h — коэффициент размытости. Ядра лучше брать усечёнными (в целях существенного сокращения объема вычислений) и достаточной степени гладкости, например параболическим:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

KЫ = ( °.75(1 - z2) N< 1

K ^=1°, 1 <|z|.

Оценка (18) представляет собой выпуклую комбинацию измерений yi выхода объекта. Веса Kn(■) в этой выпуклой комбинации определяются входом объекта. Чем ближе измерение xi к значению x, для которого рассчитывается оценка nn(x), тем значение нормированного колоколообразного ядра Kn(■) больше.

Коэффициент размытости Н вычисляется из решения экстремальной задачи, в которой квадратичный показатель соответствия выхода объекта и выхода модели на системе экспериментальных точек построен либо "методом скользящего экзамена":

П П / \

II = ЩУ - Пи(Хг ))2 =т|п, Пи(Хг) = ^ КМ( У3, (19)

г=1 3 = 1,3=*

либо прямым разбиением выборки на две части (по одной строится оценка регрессии, а по другой — минимизируемый показатель):

12 = ^ (У* - Пи2 (Хг))2 =Ш1П, ПП2 (х) = ^ Км (Х НХЛ Уз. (20)

геЫП1 3еМП2 ' '

Здесь номера первой части выборки образуют множество МП1 мощности П1 и аналогично номера второй части выборки образуют множество МП2 мощности П2. Примером множеств МП1 и МП2 служат нечётные и чётные номера измерений в упорядоченной по хг выборке:

(х*,Уг), г = 1,П.

При равномерных отсчётах по оси х, т. е. при хг+1 — хг = Д, г = 1,п — 1, целесообразно перейти от размерного коэффициента размытости Н ядра к безразмерному в = Д/Н, лежащему в интервале [0; 1] :

К„ (ХХ—Хз ) = Кг1 (Д Х^Х.) = к„ (вХ — Х

,н Д У V Д

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

Оценивание дисперсионных характеристик

Предположение, что условная дисперсия D{Y\x} = const, приводит к тому, что при расчёте математических ожиданий M{ф(X,Y)} на основе экспериментальных данных (Xi,yi), i = 1,п надо перейти к обычному оператору "среднего арифметического":

1 n

M{ф(Х, Y)}^ M{ф(Х, Y)} = - V^(xi,yi). п

i=l

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

Вычисляем оценку дисперсионного отношения (6):

Vy\х

Y,(Vn (xi) — m Y )2 i=1

n '<

Е(уі — m y )2

l n l n l n D2 = -y'(nn(xi) — mY)2, mY = -V'Уі, Dy = - У](Уі — my)2 /ji /ji /ji

i=1

i=1

i=1

Оценка дисперсионного отношения (13), основанного на применении оценки субопти-мальной модели

т т

У(х) = «0 + ^2 (х) (т. е. Уг = У(хг) = «о + ^ 0.3^3 (х*), г = 1,п), (22)

3=1 3=1

имеет вид:

'Лу\х, парам

В У

Е(л - ту)2

^---------------, (23)

Е(у - ту)2 І=1

где оптимальные параметры «о, «1,..., «т которой найдены из критерия наименьших квадратов [7]:

ао,аі,...,ат = а^ шіп -У'Ьі - «о - «З ¥>3 (хі )]2, (24)

а .Г>1 Гк™ 71. < ^ ^

хт — <х±£, _ ^.ши _ / \уг схо / <-3

ао,а1,...,ат п

г=1 3=1

т. е. параметры удовлетворяют системе линейных уравнений

1 п т 1 п

~^2 [Уг — «0 — ^ «3 ^3 (Х*)] = [Уг — У*] = 0,

П г=1 3=1 П г=1

^^^г(хг)[уг — «о — ^2 «3^3(хг)] = -Х^г(х*)[У* — Уг] = 0, I = 1,т. (25)

г=1 3=1 г=1

У(х) является оценкой параметрической модели (10), а «о, «1,..., Ат — это оценки соответствующих параметров этой модели.

Теорема (12) (о разложении дисперсии выхода объекта на две дисперсионные компоненты) справедлива и для оценок соответствующих дисперсий:

1°у = 1 1,парам + 1 2,парам. (26)

п

Оценки Iу, ]°2,парам приведены в (21), а I 1,парам = п Е [Уг — У*]2. При доказательстве

используются условия оптимальности параметров (25).

Оценка коэффициента (16)

П

І=1

Л = ІУ\х, парам = Iв 2,парам

^ Пу\х V В

2

Е [Лі - ту]2

1

Е [Пп(хі) - ту]2 І=1

(27)

характеризует степень близости выходов субоптимальной параметрической модели (22) и оптимальной непараметрической модели —- оценки регрессии М{У\х*} = Пп(х*), г = 1, п. Оценка относительного коэффициента корреляции (17)

п

Е(У* — т у )(Пп(х*) — т у)

^ = /п“‘ 1п (28)

у — ту)%/^Е(Пп(х*) — ту)2

также представляет собой коэффициент взаимосвязи между субоптимальной параметрической моделью (22) и непараметрической оценкой регрессии.

Процесс обработки информации при сравнении качества непараметрической и параметрической моделей представлен на рис. 1.

Рис. 1. Схема расчета характеристик близости непараметрической и параметрической статистических моделей стохастических объектов

Пример

В выборке (xi,yi), i = 1,100, входа и выхода объекта значения входа xi берутся последовательно из диапазона [0; 10] на равных друг от друга интервалах Д. Каждое значение выхода yi равно сумме значений выбранной сигнальной части ^(xi) и помехи:

yi = ф(х{) + а&, ф(хг) = \xi - 5| +cos xH + 3sin2xi -\x4 - 2|. (29)

Здесь a — среднеквадратичное уклонение аддитивной помехи, р — отношение "шум-сигнал" , £i — значение генератора нормально распределённой помехи с нулевым мате-

матическим ожиданием и единичной дисперсией. Величина р равна отношению интервала изменения помехи к интервалу а изменения сигнальной части выхода имитируемого объекта. При нормально распределённой помехе интервал изменения её обычно берётся равным 6a. В этом случае р = а, и тогда среднеквадратичное уклонение помехи a = ра определяется величиной отношения "шум-сигнал" р, которая варьируется от 0.01 до 2 (что

соответствует интервалу от 1 до 200%) на приводимых ниже графиках для дисперсионных характеристик.

При аддитивной равномерно распределённой в интервале [-9, в] помехе слагаемое в (29), соответствующее помехе, имеет вид в^і, Сі — значение генератора равномерно распределённой в интервале [-1; 1] помехи. Тогда р = а и в = ра.

Рис. 2. Фрагмент поведения (при 100 %-й помехе) непараметрической модели nn(x) и первого варианта параметрической модели yn(x) (слева), а также зависимость оценок дисперсионных характеристик от уровня р стохастичности объекта

Первый вариант. Структура выбранной параметрической модели имеет вид y(x) = = «о + ai\x — 5\. В ней присутствует постоянная компонента ао и базисная функция \x — 5\ с параметром ai. При наличии экспериментальных данных, представленных в левой части рис. 2 для 100 %-й нормально распределённой помехи, из квадратичного критерия рассчитаны параметры модели yn(x) = —2.973 + 0.803\x — 5\. По тем же экспериментальным данным настроена (при этом оптимальный коэффициент размытости в = 0.12) непараметрическая оценка регрессии nn(x). Она даже визуально существенно лучше оптимальной параметрической модели. Непараметрическая модель близка к сигнальной части выхода объекта ф(^.

На этом же рисунке представлено поведение оценок всех основных дисперсионных характеристик при изменении отношения "шум-сигнал" р от 1 до 200%. Оценка D2 дисперсии регрессии слабо меняется, что указывает на хорошее качество непараметрической оценки регрессии и при высоком уровне помех.

Параметрическая и непараметрическая оценки регрессии далеки друг от друга, и Ау близка к нулю. Оценка относительного коэффициента корреляции Ж у между этими моделями тоже невелика, но больше Ху. Оценка дисперсионного отношения для регрессии Пу\х существенно больше дисперсионного отношения г/у\х,парам для параметрической модели. Сравнительно большая величина оценки дисперсионного отношения г/у\х для регрессии подтверждает высокое качество непараметрической оценки регрессии.

Второй вариант. Модель y(x) = ао + ai\x — 5\ + а2 cos x имеет две базисные функции, которые взяты из набора базисных функций, присутствующих в сигнальной части выхода объекта (29). При 100 %-й нормально распределённой помехе оценка субоптимальной параметрической модели yn(x) = —3.002 + 0.911 \x — 5\ + 0.748cosx и непараметрическая оценка регрессии (при оптимальном коэффициенте размытости в = 0.11) представлены в левой части рис. 3.

Из рис. 3 видно, что параметрическая оценка регрессии находится ближе к непараметрической, чем в прошлом эксперименте. В случае 10 %-й помехи оценка коэффициента их близости поднимается с 0.1 (рис. 2) до 0.23, а величина оценки Жукоэффициента

корреляции возрастает с 0.3 до 0. 48.

Рис. 3. Фрагмент поведения (при 100 %-й помехе) непараметрической модели nn(x) и второго варианта параметрической модели yn(x) (слева), а также зависимость оценок дисперсионных характеристик от уровня р стохастичности объекта

Третий вариант. Параметрическая модель является более сложной и включает три базисные функции (специально из того же набора базисных функций, которые присутствуют в сигнальной части при имитации объекта): y(x) = ао + ai|x — 5| + a2 cosx + аз sin2x. Её конкретный численный вариант при 100 %-й нормально распределённой помехе имеет вид: УП(х) = —1.273 + 0.892|x — 5| +0.757cosx + 2.571 sin2x. Непараметрическая оценка регрессии имеет оптимальный коэффициент в = 0.11. Эти модели представлены в левой части рис. 4.

Рис. 4. Фрагмент поведения (при 100 %-й помехе) непараметрической модели nn(x) и третьего варианта параметрической модели yn(x) (слева), а также зависимость оценок дисперсионных характеристик от уровня р стохастичности объекта

Качество параметрической модели возросло. Это подтверждает и возрастание оценок для дисперсионного отношения г/у\х,парам для коэффициента близости непара-

метрической и параметрической моделей и для относительного коэффициента корреляции Жумежду моделями.

Четвёртый вариант. Параметрическая модель y(x) = ао+ai |x—5|+а2 cos x+аз sin 2x+ +a4|x — 2| по структуре совпадает с сигнальной частью выхода объекта (при его имитации). На рис. 5 при величине помехи 100% оптимальная модель yn(x) = 0.172 + 0.988|x — 5| + +0.856cosx + 2.796sin2x — 1.117|x — 2|. Она фактически накрывает сигнальную часть выхода. Отклонение вызвано только влиянием помехи. В непараметрической оценке регрессии оптимальный коэффициент размытости в = 0.10.

Как следует из рис. 5, при небольшом уровне помех параметрическая оценка бывает даже лучше непараметрической, при этом оценка Хунесколько выше 1. Оценка относительного коэффициента корреляции Ж у тоже близка к 1, но никогда не превышает её. Все дисперсионные характеристики говорят о наилучшем выборе структуры параметрической модели. С ростом уровня помех непараметрическая оценка всё-таки лучше параметрической.

Рис. 5. Фрагмент поведения (при 100 %-й помехе) непараметрической модели цп(х) и четвертого варианта параметрической модели уп(х) (слева), а также зависимость оценок дисперсионных характеристик от уровня р стохастичности объекта

Анализ всех численных результатов показывает, что оценка относительного коэффициента корреляции Ж у имеет несколько меньшую чувствительность к изменению близости

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

непараметрической и параметрической оценок регрессии по отношению к Ху .

Автор благодарен научному сотруднику М. А. Меретилову за выполненные численные расчёты.

Заключение

В многомерном случае (когда объект имеет несколько входов X и один выход У) все вышеприведенные дисперсионные характеристики сохраняются. В них лишь скалярную случайную величину X надо заменить на векторную X. Расчеты всех оценок (особенно непараметрической оценки регрессии), естественно, усложняются, но главная идея остаётся неизменной.

Вышеприведённая схема анализа моделей стохастических объектов применима и к нелинейным (по каналу "вход-выход" и по параметрам) параметрическим субоптимальным моделям У = «о + Э(х, X), оптимальные параметры ао, ~а которых удовлетворяют критерию наименьших квадратов.

Все приведённые результаты справедливы для так называемого гомоскедастичного варианта, т. е. когда условная дисперсия выхода В{У\х} является постоянной величиной. Это условие было выполнено в численных примерах при имитации объекта за счёт представления выхода объекта в виде суммы сигнальной части (некоторой функции от входа) и помехи, дисперсия которой постоянная для всех значений из рассматриваемого диапазона изменения входа. Гомоскедастичность нарушается, например, если вход объекта измеряется с аддитивной помехой (даже с постоянной дисперсией). Гетероскедастичные случаи [2] требуют специальных исследований, но общая идея сохраняется.

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

Результаты могут быть распространены и на динамические стохастические объекты.

Список литературы

[1] Райбман Н. С., Чадеев В. М. Построение моделей процессов производства. М.: Энергия, 1975. 376 с.

[2] Дисперсионная идентификация / Под ред. Н. С. Райбмана. М.: Наука, 1981. 336 с.

[3] Рубан А. И. Непараметрическая дисперсионная идентификация // Изв. АН СССР. Техническая кибернетика. 1989. № 3. С. 212.

[4] Rouban A. I. Nonparametric dispersing identification // Advances in Modeling & Analysis: Series D. Mathematical Tools; General Computer Tools. France: A.M.S.E. 1998. Vol. 1. №2. P. 43-50.

[5] Рубан А. И. Идентификация стохастических объектов на основе непараметрического подхода // Автоматика и телемеханика. 1979. №11. С. 106-118.

[6] Рубан А. И. Методы анализа данных: учеб. пособие. 2-е изд., испр. и доп. Красноярск: ИПЦ КГТУ, 2004. 319 с. (Глава 4)

[7] Себер Дж. Линейный регрессионный анализ. М.: Мир, 1980. 456 с.

The Dispersing Characteristics of Static Models of Stochastic Objects

Anatoly I. Rouban

Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041 Russia

The dispersing characteristics determining quality of static parametrical models of stochastic objects are constructed, and their estimations are found at the assumption, that the conditional variance of an output is constant. Parametrical models are nonlinear on the channel "input - output" but they are linear concerning optimized parameters. It is shown that the nonparametric estimation of regression is close to objectively existing ideal model and according to dispersing characteristics it is possible to keep up with change of quality of suboptimal parametrical models at selection of their structure. Quality of nonparametric and parametrical models compared on several numerical tests with various level of noise.

Key words: dispersing identification, static parametrical model, stochastic object, a nonparametric estimation for regression.

i Надоели баннеры? Вы всегда можете отключить рекламу.