CC BY
122
В.Е. Зотеев
ИТЕРАНИОННЫЙ МЕТОД СРЕДНЕКВАДРАТИЧНОГО ОЦЕНИВАНИЯ ПАРАМЕТРОВ КОРРЕЛЯЦИОННОЙ ФУНКЦИИ РАСПРЕДЕЛЕНИЯ СЛУЧАЙНОГО ПОЛЯ НЕУПРУГОЙ РЕОЛОГИЧЕСКОЙ ДЕФОРМАЦИИ
Рассматривается эффективный численный метод среднеквадратичного оценивания параметров корреляционной функции распределения случайного поля неупругой реологической деформации. В основе метода лежит линейно параметрическая дискретная модель в форме стохастического разностного уравнения, коэффициенты которой известным образом связаны с параметрами корреляционной функции. Приводятся результаты численно-аналитических исследований сходимости метода и погрешности вычисления параметров корреляционной функции. Данный метод может быть рекомендован для эффективного решения задач параметрической идентификации различных по физической природе объектов, динамические процессы в которых описываются затухающей по экспоненте синусоидой.
Рассматривается задача определения параметров корреляционной функции распределения случайного поля неупругой реологической деформации по экспериментальным данным [1]. Важность этой проблемы диктуется широким кругом методов решения краевых стохастических задач пластичности и ползучести, основанных на использовании аналитической аппроксимации корреляционной функции [2,3]. Особенностью исследуемого процесса является сильная зашумленность истинных значений корреляционной функции случайной аддитивной помехой в результатах наблюдений. Достоверное определение параметров корреляционной функции в таких условиях требует применения статистических методов обработки результатов эксперимента. Классические способы вычисления среднеквадратичных оценок параметров корреляционной функции приводят к численным итерационным методам решения нелинейной системы нормальных уравнений, эффективность которых сильно зависит от выбора начального приближения решения [4]. Как правило, применение таких способов в практике обработки результатов эксперимента вызывает серьезные затруднения.
В данной работе поставленная задача решается в формате современных информационных технологий. Принципиальным моментом для нового способа решения является разработка и применение линейно параметрической дискретной модели (ЛПДМ) корреляционной функции в форме стохастического разностного уравнения, коэффициенты которой известным образом связаны с параметрами корреляционной функции. При таком подходе задача определения параметров корреляционной функции случайного поля неупругой реологической деформации сводится к устойчивому среднеквадратичному оцениванию коэффициентов стохастического разностного уравнения [5].
Известно [1], что корреляционная функция распределения случайного поля неупругой реологической деформации описывается выражением вида
а
~ (х) = а0 ехр (- |а|х
008 (( х) + ь (( х)
(1)
где а, ( и а0 - параметры, подлежащие определению. В ходе эксперимента обычно формируются выборки небольших объемов (N = 20 + 30), а результаты наблюдений ук = ~ (хк) + е к, к = 0, N — 1, содержат случайную аддитивную составляющую ек, существенно искажающую теоретические значения корреляционной функции ~(хк). Поэтому устойчивое и достоверное оценивание параметров а, ( и а0 является важнейшей проблемой при исследовании случайных полей неупругой реологической деформации.
В данной работе рассматривается новый метод определения параметров корреляционной функции распределения случайного поля неупругой реологической деформации, в основе которого лежит итерационная процедура среднеквадратичного оценивания коэффициентов ЛПДМ в форме стохастического разностного уравнения. Разработка этого метода включает решение следующих основных задач:
—построение линейно параметрической дискретной модели, связывающей последовательные отсчеты корреляционной функции, в форме стохастического разностного уравнения;
- разработка алгоритма определения параметров корреляционной функции на основе среднеквадратичных оценок коэффициентов стохастического разностного уравнения;
—исследование эффективности итерационного метода среднеквадратичного оценивания коэффициентов ЛПДМ;
—оценка погрешности вычисления параметров корреляционной функции в процессе обработки результатов эксперимента.
Построение линейно параметрической дискретной модели корреляционной функции в форме стохастического разностного уравнения. Основой предлагаемого метода определения параметров корреляционной функции распределения случайного поля неупругой реологической деформации является линейно параметрическая дискретная модель. Следует различать линейно параметрическую дискретную модель в виде рекуррентной формулы, линейно связывающей некоторые комбинации дискретных значений исследуемой функциональной зависимости, и ЛПДМ в форме стохастического разностного уравнения, описывающего эту зависимость с учетом аддитивной случайной помехи в результатах наблюдений.
Рассмотрим построение ЛПДМ в виде рекуррентной формулы, связывающей три последо-
вательных отсчета непрерывной функции у (x) = a0 exp (— ax)
cos (fix) + a sin(ffx)
Полагая
х = тк , к = 0,1,2,* , где т - период равномерной дискретизации функции, получаем ее дискретный аналог:
yk = a0 exp(- atk)
cos {¡ftk) + sin (fftk)
(2)
Применяя к обеим частям данного выражения /-преобразование, получаем: г{~к}= а0 1 — * 1е "(<=08—3т — а3 2' 5111 3), где 1= 2еатсо^Зт, 12 = —г"2“. Отсюда
Z{yk}—z %Z{yk}—z 2A2Z{yk} = a0 —z la0e at(cosfit—afi !sinft). Возвращаясь в пространство оригиналов, используя первую теорему смещения: z —Z {pk }=yk—r, r = 0, 1,2, при условии, что yk—r = 0 при k — r < 0, получаем разностное уравнение вида
у — 1iyk—1 —l2yk—2 = dkÜ0 — 8k—1Ü0 e ~at (cos ft—af —1 sin ft), k=0, l,2,* ,
|l, при k = 0;
где Sk = < ’ ’ - символ Кронекера.
[0, при k ф 0
При k = 0 и k = 1 с учетом 14 = a0af -1e ~at sin ft отсюда соответственно получаем
у0 = a0 = 13, у = a0 e ~at cos ft + a0af 4e ~at sin ft = “Jo +14. (3)
При k > 2 приходим к линейно параметрической дискретной модели в форме разностного уравнения второго порядка
yk = ЛУ—1 + 12y?k—2 = ^ k = 2,3, * , (4)
которое рекуррентно описывает временную последовательность дискретных значений корреляционной функции (2).
Коэффициенты в модели (4) связаны с параметрами корреляционной функции соотношениями:
1 = 2e~at cos ¡ft, 12 = — e~2at, 13 = a0, 14 = a0af ~le~at sin ¡ft. (5)
При обработке эксперимента формируется выборка результатов наблюдений yk, k = 0, N — 1, где N - объем выборки. Результаты наблюдений yk содержат аддитивную случайную помеху ek: yk = yk + £k, где yk - дискретные значения, соответствующие теоретической корреляционной функции, и, следовательно, соотношениям (1) и (2). Тогда ЛПДМ в форме стохастического разностного уравнения можно представить в виде
Уо = Я3 + е о,
У\ = Я + Я4 2 Є о + е 1 ’
Ук = Я ук-1 + Я2 Ук-2 + Ц к ,
Цк = -Я2Єк-2 — Я1Єк-1 + Єк ,
к = 2,Ж -\.
В матричной форме ЛПДМ, описывающая корреляционную функцию при аддитивной случайной помехе в отсчетах наблюдений, имеет вид
ГЬ = F Я + ц,
[ц = Ре
или в развернутой форме обобщенной регрессионной модели -
Ьі = Я\ /л + Я2 /і 2 + Я3 /і 3 + Я4 /і 4 + Ц і,
Ц і = Рі\е 0 + Рі 2Є1 + рі3Є 2 + * + Рг,Ы-\£ N-2 + еЖ-\, г = \, Ж .
Здесь Я = (Я\, Я2, Я3, Я4 )Т - вектор неизвестных коэффициентов ЛПДМ; е = (е 0, е \,*
N - мерный вектор случайной помехи в результатах наблюдений; ц = (Л\, Л2,* ,ЛЖ У - N -мерный вектор эквивалентного возмущения в стохастическом разностном уравнении; Ь = (у0,у\,у2,* ,ук-1,* ,yN-\ )Т- N - мерный вектор правой части; F- матрица регрессоров размера N X 4, столбцы которой описываются формулами:
/ \Т
/,2 =(0, 0, Уо, У\, У2 , * , Ук-3,К , У N-3 )Т , /,3 =(1,0,К ,0)Т ,
(7)
(8) (9)
Є N-1 )Т -
А =
Уо
У\, У2 , * , У к-2 ’
, УN-2
/,4 =(0,1,0,* ,0)Т ; Р - матрица размера N X N в стохастическом уравнении эквивалентного возмущения- нижняя треугольная, ленточная, трех диагональная. Первый столбец матрицы Р
имеет вид: Рі\ = ются формулами:
\,-—, - Я2, о, о, 22
Т
Остальные столбцы матрицы Р, у =2^, описыва-
Р« =
0,
1,
-11,
- Я2 о,
= ЪК,
= 2, N,
= 2, N -1, = 2, N - 2,
(Ю)
« = 2, N - 3.
1 < г < У, г =
г = У +1, г = У + 2,
У + 3 < г < N,
Построенная линейно параметрическая дискретная модель в форме системы из двух стохастических разностных уравнения (7) или (8) и (9) лежит в основе предлагаемого метода определения параметров корреляционной функции распределения случайного поля неупругой реологической деформации.
Определение параметров корреляционной функции основе стохастического разностного уравнения. Исходными данными, на основе которых определяются параметры корреляционной функции, служат результаты вычислений, полученные в процессе физического эксперимента. Проводя анализ экспериментальных данных, например [1], представленных для математической обработки, следует отметить, во-первых, небольшие объемы выборок: N < 20, а, во-вторых, то, что эти выборки охватывают не более одного периода корреляционной функции:
2р
Т = -^ • Очевидно, что при этом период дискретизации т, используемый при формировании
экспериментальных данных, сравнительно мал (в относительных к периоду корреляционной функции единицах). Это отрицательно влияет на устойчивость вычислений среднеквадратичных оценок коэффициентов ЛПДМ и должно быть учтено при разработке численного метода определения параметров корреляционной функции.
Полученные в ходе эксперимента результаты вычислений ук содержат аддитивную случайную помеху ек. Будем предполагать, что случайные возмущения ек имеют нулевое математическое ожидание, дисперсию Г2, не коррелированны между собой (матрица дисперсий-
о
ковариаций: V [e ]= а2е E , где E - единичная матрица) и подчиняются нормальному закону распределения: e k е N (о, ol).
Алгоритм определения параметров корреляционной функции на основе среднеквадратичных оценок коэффициентов стохастического разностного уравнения включает следующие основные этапы.
1. Формирование на основе ЛПДМ системы линейных алгебраических уравнений. На этом этапе по приведенным выше формулам, описывающим элементы стохастического разностного
уравнения (8), вычисляются элементы матрицы F = (fy), j = 1,4 и правой части b = (b), i = 1, N , переопределенной системы линейных алгебраических уравнений относительно коэффициентов 1j, j = 1,4 .
2. Вычисление среднеквадратичных оценок коэффициентов стохастического разностного уравнения. Данный этап является определяющим в процессе формирования погрешности оценок параметров корреляционной функции. От того, насколько обоснованно и корректно будет выбран метод среднеквадратичного оценивания, существенно зависит точность вычисления коэффициентов ЛПДМ. Разработка и исследование эффективного численного итерационного метода среднеквадратичного оценивания коэффициентов стохастического разностного уравнения (8) является основной задачей, решаемой в рамках данной работы.
3. Вычисление параметров корреляционной функции распределения случайного поля неупругой реологической деформации. С учетом известных соотношений (5) по среднеквадратичным
оценкам Xi коэффициентов ЛПДМ вычисляются параметры корреляционной функции:
а =—— ln(- Я2), ß =1 arccos—^=, a0 = X3. (11)
2т Т 2л/-!
4. Оценка погрешности результатов вычислений. На заключительном этапе алгоритма определения параметров корреляционной функции оценивается погрешность полученных результатов. Так как в алгоритме вычислений применяются статистические методы обработки экспериментальных данных, то, очевидно, что погрешность носит случайный характер, и для ее описания можно использовать среднее квадратическое отклонение оценок параметров корреляционной функции. Построение доверительных интервалов для этих оценок позволит указать предельные (с заданной доверительной вероятностью) абсолютную и относительную погрешности.
Разработка и исследование итерационного метода среднеквадратичного оценивания коэффициентов стохастического разностного уравнения. Применение статистических методов при вычислении коэффициентов l j линейно параметрической дискретной модели позволяет существенно повысить точность оценивания. Однако вычисление среднеквадратичных оценок Xj по обобщенной регрессионной модели, которая описывается уравнением (8):
п2
4 N
bt — ^Xjfij+1 = 2 ^ min , приводит к значительному асимптотическому смещению в
j=1 J i=1
оценках из-за корреляции между отсчетами случайного эквивалентного возмущения hi и элементами матрицы fij, i = 1,2, [6]. Этот существенный недостаток устраняется за счет использования итерационной процедуры при вычислении среднеквадратичных оценок.
Из (10) очевидно, что всегда существует обратная матрица P—1. Тогда после преобразования обобщенной регрессионной модели (7) получаем Pxlb = P— F X + e . Вычисление среднеквадратичных оценок коэффициентов Xj на основе минимизации функционала
II —1 —1 || 2 и ц2
Jx = \PX b — PX F X\ = e ^ min позволяет устранить асимптотическое смещение, связанное
с корреляцией между случайным возмущением и результатами наблюдений, и тем самым повысить точность вычисления параметров корреляционной функции. Так как элементы матрицы P— зависят от Xj и X2, то, в общем случае, минимизация функционала Jх приводит к решению системы нелинейных относительно Xj уравнений. Рассмотрим итерационную процедуру вычисления среднеквадратичных оценок коэффициентов Xj. На первом этапе вычисляются 126
первоначальные МНК - оценки : |b - FA ||2 = Щ\2 ^ min . Затем формируется приближение матрицы Р х » Р(1 (l(0), Ai0), преобразуется обобщенная регрессионная модель (7), и вычисляются новые, уточненные оценки A(l) коэффициентов ЛПДМ. Эти оценки используются для вычисления второго приближения матрицы Р A » Р^2 (l{l), a2!)) и т.д. Процесс уточнений повторяется до тех пор, пока не будет выполнено условие останова, например: A) - 2A(i-1) + A(i 2)| < 0,0ЦA(i^. Найденные на последней итерации оценки A(i) принимаются за
истинные значения коэффициентов ЛПДМ и по ним далее вычисляются параметры корреляционной функции распределения случайного поля неупругой реологической деформации.
Проведены численно-аналитические исследования сходимости описанной итерационной процедуры. Численный эксперимент был организован следующим образом. Генерировалась выборка —к дискретных значений функции (2) с параметрами a0 = l, а = 0,5, ß = 2 . Период дискретизации t = 0,15 и объем выборки N = 20 при формировании тестового сигнала выбирались в соответствии с режимом проведения физического эксперимента, в частности:
t 2p —
T = 0,048, где T = -ß- - период корреляционной функции. В отсчеты ук, соответствующие
точным значениям функции, добавлялась аддитивная помеха £к, значения которой обеспечивали заданную величину мощности помехи 3%, 5% и 10% в относительных единицах к мощно-
1Ы
сти сигнала: ^^ 100% . Результаты вычислений параметров корреляционной функции усредня-
\\у\\
лись по l000 реализациям.
На рис. 1 и рис. 2 представлены зависимости относительной погрешности вычисления параметра а, характеризующего затухание корреляционной функции, и параметра ß - частоты корреляционной функции, от числа итераций n для различных значений мощности случайной помехи в результатах наблюдений. Кривые 1, 2 и 3 построены по результатам вычислений при мощности случайной помехи 3%, 5% и 10%. Точки на графиках при n = 0 соответствуют результатам среднеквадратичного оценивания без предварительного преобразования обобщенной регрессионной модели (7).
Р и с. 1 Зависимость относительной погрешности вычисления параметра а от числа итераций (в логарифмическом масс-штабе)
Р и с. 2 Зависимость относительной погрешности вычисления параметра Ь от числа итераций (в логарифмическом масс-штабе)
По результатам исследований можно сделать выводы, что, во-первых, применение итерационного метода среднеквадратичного оценивания позволяет существенно, в сотни раз, уменьшить погрешность вычисления параметров корреляционной функции. Во-вторых, результаты вычислений практически не изменяются после второй итерации, что свидетельствует о хорошей сходимости метода. При этом точность оценивания параметра / - частоты корреляционной функции - на два порядка выше, чем параметра а.
Проведены исследования эффективности предложенного итерационного метода вычисления среднеквадратичных оценок. Целью таких исследований являлся сравнительный анализ погрешности вычисления параметров корреляционной функции двумя различными способами: на основе среднеквадратичных оценок коэффициентов модели (8) и на основе итерационного численного метода. Для этого, как и в предыдущем случае, формировалась выборка дискретных значений тестового сигнала с параметрами а0 = 1, а = 0,5, / = 2, т = 0,15 , N = 20 , что соответствует рабочей области проведения физического эксперимента по изучению неупругой реологической деформации. В отсчеты тестового сигнала добавлялась случайная помеха, мощность которой изменялась в диапазоне от 0 до 20%.
В качестве критериев эффективности в сравнительном анализе методов среднеквадратичного оценивания использовались:
„ й \М [а ]— а|
— смещение оценок а и р относительно истинных значений (в %): Ла =---------,—,--- 100%;
\а\
— дисперсии оценок параметров корреляционной функции: 0[а ] = М [( — М [а ])2 ];
— обобщенная ошибка вычисления коэффициентов разностного уравнения - второй момент
= Б
оценки относительно истинного значения коэффициента: М
Здесь М [] - оператор математического ожидания, статистическая оценка которого в каждом случае находилась посредством усреднения М = 1000 реализаций для каждой точки эксперимента.
На рис. 3, 4 представлены зависимости смещения оценок параметров корреляционной функции (в относительных единицах) от мощности случайной помехи. Так как результаты вычислений различаются в сотни, тысячи и более раз, то для оси ординат используется логарифмический масштаб. Точки 1 соответствуют методу наименьших квадратов, примененному к модели (8) без ее преобразования, а точки 2 - итерационному методу вычислений среднеквадратичных оценок коэффициентов ЛПДМ. Очевидно, что точность вычисления при использовании итерационного метода среднеквадратичного оценивания на несколько порядков выше по сравнению с методом наименьших квадратов, примененному непосредственно к модели (8). Это достигается за счет устранения большого смещения в МНК - оценках, обусловленного корреляцией между случайной эквивалентной помехой ц и элементами матрицы Г .
1§(Аа,%) 2
1
-2
-3
1 ; ♦ ♦ • * . ♦ * 1
. ♦ ♦ і
♦ 10 ■
0 2 5 .■*■■■ 15
■ ■ ■ ■ 1
1
18(43,%)
£,%
1 ♦ ♦
\ ♦ ♦ # # * ' ♦
0 ♦ 5 ■ ■ 10 ■ ■ ■ 15 ■ ■ ■ ■ 1
2 . ■ ■ " ■ ■
■
£,%
Р и с. 3. Зависимость смещения оценки параметра р и с. 4. Зависимость смещения оценки парамет-
а от величины случайной помехи в результатах ра / от величины случайной помехи в результа-
наблюдений (в логарифмическом масштабе) ТаХ наблюдений (в логарифмическом масштабе)
128
0
Несмещенность оценок на практике не всегда гарантирует высокую точность. Большое значение имеет и дисперсия оценок, которая зависит от выбора параметров обработки экспериментальных данных. Проведены исследования зависимости дисперсий оценок коэффициентов ЛПДМ и параметров корреляционной функции от мощности случайной помехи в результатах
эксперимента. Выборочные дисперсии оценок коэффициентов ЛПДМ 52[/А] и 52
вычисля-
лись по М = 1000 реализациям для каждой точки эксперимента: 52 [А, ]= —
м ,
У (ам - А■ )2
М -1
У А г
¿—I г
Р=1
М
по формулам:
А3 = —----, 3 = 1,2. Оценки дисперсий для параметров корреляционной функции определялись
3 М
52 [а ]
да
■ А ]■
дА
в которых в соответствии с (11)
дА
[]+ 2 д3 д3
да
ООУ
дА
дА1 дА2
1 эз
2тА ’ дА
[А1,А ]+
д3
дА2
2 [а ],
(12)
(13)
тд/ - 4А2 - А
д3
и —^
эАь2
А
І2
'2 2тА^д/ - 4А - А
выборочная ковариация между оценками А1 и А вычислялась по формуле
[, А2
УМ р - А ) р - А)
_ р=1
ООУ
М -1
На рис. 5, 6 представлены зависимости (в логарифмическом масштабе) среднего квадратического отклонения оценок параметров корреляционной функции от величины случайной помехи в результатах наблюдений. Очевидно, что применение итерационного метода (точки 2) позволило на два порядка уменьшить дисперсии по сравнению со среднеквадратичным оцениванием по модели (8) без ее преобразования (точки 1).
2
2
5
1
а
Р и с. 5 Зависимость среднего квадратического Р и с. 6 Зависимость среднего квадратического
отклонения оценки параметра а от величины отклонения оценки параметра 3 от величины
случайной помехи в результатах наблюдений случайной помехи в результатах наблюдений
(в л°гарифмическ°м масштабе) (в логарифмическом масштабе)
В процессе формирования результирующей ошибки заметную роль играет погрешность среднеквадратичного оценивания коэффициентов А ЛПДМ. При фиксированной выборке результатов эксперимента эта погрешность зависит от смещения и дисперсии оценки. Причем
129
требование несмещенности на практике не всегда оправдано, так как оценка с небольшим смещением и малой дисперсией может оказаться предпочтительнее несмещенной оценки, но с большей дисперсией. Поэтому целесообразно использовать обобщенный показатель ошибки -второй момент оценки относительно истинного значения параметра. На рис. 7, 8 представлены зависимости (в логарифмическом масштабе) обобщенной ошибки вычисления коэффициентов А и А от величины случайной помехи в результатах наблюдений. Точки 1 соответствуют среднеквадратичным оценкам, полученным на основе модели (8) без ее преобразования, а точки 2 - итерационному методу. Как и в предыдущих случаях, применение итерационного метода среднеквадратичного оценивания существенно, на несколько порядков, позволяет повысить точность вычислений.
Р и с. 7 Зависимость обобщенной ошибки вы- Р и с. 8 Зависимость обобщенной ошибки вычисления коэффициента А1 от величины случай- числения коэффициента А2 от величины слу-ной помехи в результатах наблюдений (в лога- чайной помехи в результатах наблюдений (в ло-
рифмическом масштабе) гарифмическом масштабе)
Таким образом, проведенные численно-аналитические исследования помехозащищенности итерационного численного метода, в основе которого лежит линейно параметрическая дискретная модель в форме стохастического разностного уравнения, показали его высокую эффективность при вычислении параметров корреляционной функции распределения случайного поля неупругой реологической деформации.
Оценка погрешности вычисления параметров корреляционной функции в процессе обработки результатов эксперимента. Основным источником погрешности вычисления параметров корреляционной функции является случайная помеха в результатах эксперимента, поэтому для оценки погрешности целесообразно использовать статистические методы [7]. Процедуру вычисления результирующей погрешности можно разбить на следующие основные этапы. Вначале на основе статистической обработки результатов оценивания коэффициентов обобщенной регрессионной модели в форме стохастических разностных уравнений (8) и (9) находятся оценки дисперсий 52 [Ху ] для коэффициентов Аз, 3 = 1,3 . Затем с учетом формул (12) и (13) вычисляются оценки дисперсий 52[а], 52[/] и 52[а0] параметров корреляционной функции. На заключительном этапе вычисляются доверительные границы случайной погрешности, которые могут быть интерпретированы как предельные абсолютные погрешности Ла , Л/ и Ла0 (с заданной доверительной вероятностью) для оценок параметров корреляционной функции. Полученные результаты позволят также найти и предельные (с заданной доверительной вероятно. о Ла Лр ~ Ла0
стью) относительные погрешности вычисления да = —р, дВ = и да0 =-,—0 параметров
\щ |р| Ы
корреляционной функции.
При использовании итерационного метода вычисления среднеквадратичных оценок коэф-
I -1 -1 II2 II || 2
фициентов ЛПДМ минимизируется функционал Jm = \Pl b - Px FM = e ^ min. Это при-
£■ ( T 1 \—1 T 1 T
водит к оценкам l = \F Wl FJ F Wl b , где матрица WM = PXPX . Отсюда можно получить l = l + (ft W— F)-1 FT (P— )T e . Тогда M[l]=l + M (ftW — F) FT (p— J e] = l + Dl, где
Dl = M pT W —1F) FT (p— ) eJ - смещение, обусловленное корреляцией между случайным
возмущением e и элементами матрицы F . Так как при данном алгоритме оценивания смещение мало и практически равно нулю (что подтверждается результатами не только численных, но и аналитических исследований), то в первом приближении M[l]= l.
Рассмотрим ковариационную матрицу оценок l: V[l] = M -M[l) -M[l])
. Исполь-
зуя линейную аппроксимацию вектора l -I = (ftWxF)1FT [px 1J e »(ftW;-1/7 j-1 FT (^21 J e,
где F - неслучайная матрица регрессоров, элементы которой описываются теми же формулами, что и для матрицы F с учетом замены ук = ук + ek на ук, получаем:
( -1 )(-1) ] = M [( W- F j1 F t (p- JeeT P- F ( w-1F j1 ] = ( w-1F j1 M [eeT ] =
= [ftW^f)1^2 » [ftWxF^g2. Обозначив элементы матрицы C = [ftWxFj1 через ciJ- , i = 1,4 , j = 1,4, дисперсии и ковариации коэффициентов 1 можно оценить по формулам 52 [l ]= ciis2, ооу[; , 1j ]= CjS2, где s2- оценка дисперсии <72 случайного возмущения в результатах наблюдений. При обработке экспериментальных данных в соответствии с данной мето-
N
К -е)2
w w 2 к=1 дикой дисперсию случайного возмущения можно оценить по выражению s =-------------------------,
N - 4
1 N
е =— У ек , где ек - остатки, которые находятся через оценки вектора коэффициентов ЛПДМ
N к=1
по формуле е = р;1ь - p;1f; .
По известным оценкам дисперсий и ковариаций коэффициентов ЛПДМ можно найти оценки средних квадратических отклонений параметров корреляционной функции. Из формул (12)
и (13) следует s[«]= \ s[;], s[/] =------, 1 V4-l2s2[ii]-4l112cov[l1,12]+1s2 [1] и
2Л 4 2t| ;|v - 412 -1
sk]=s[13].
В качестве оценки предельной абсолютной погрешности вычисления параметров корреляционной функции (с доверительной вероятностью p) можно использовать границу доверительного интервала, которая для a, например, вычисляется по формуле: Da = tTs[a] [6]. Здесь величина tT = t (g,v) берется из таблицы распределения Стьюдента при числе степеней свободы v = N - 4 и уровне значимости g = 1 - p . При доверительной вероятности p = 0,99 и объеме выборки N > 15 принимаем tT = 2,7. Границы доверительных интервалов для остальных параметров вычисляются по аналогичным формулам.
Проведены численно-аналитические исследования, подтверждающие достоверность полученных формул для оценки предельных абсолютной и относительной погрешностей вычисления параметров корреляционной функции при обработке экспериментальных данных. На рис. 9 - 10 кривые 1 - 3 описывают границы для предельной относительной погрешности вычисления параметров a и / , построенные по десяти реализациям для каждой точки эксперимента. Кривые 1 и 3 соответствуют максимальному и минимальному значениям предельной относительной погрешности, а кривая 2 - значению, усредненному по десяти реализациям.
Точки 4 на этих графиках отображают результаты вычислений действительной относительной погрешности: по 10 значений в каждой точке численного эксперимента. Очевидно, что погрешность оценок параметров корреляционной функции практически полностью укладывается в границы для предельной относительной погрешности, вычисленные в соответствии с предложенной методикой по представленным выше формулам.
8а,%
30
20
10
; // // '*7 у/
1 3 - \ / . / Vх jr \ ; '/h / />ч / \ / / *
L > у V 1 \ * к * 4 ^ і А д А *
/ * * і А А і f і і * 1 1 . ! * 1 А А і!:іФі Ґ:,!:
8b,%
10
15
£,%
10
10
15
0 5
Р и с. 9 Зависимости оценок предельной относи- Р и с. 10 Зависимости оценок предельной относительной погрешности (кривые 1, 2 и 3) и относи- тельной погрешности (кривые 1, 2 и 3) и относительных погрешностей (точки 2) вычисления па- тельных погрешностей (точки 2) вычисления параметра а от величины случайной помехи раметра Ь от величины случайной помехи
5
Применение итерационного метода среднеквадратичного оценивания параметров корреляционной функции при исследовании ползучести материалов. Разработанный метод оценивания параметров корреляционной функции на основе итерационной процедуры среднеквадратичного оценивания коэффициентов стохастического разностного уравнения был применен в научно-техническом эксперименте при исследовании ползучести материалов. Результаты измерений, полученные в ходе физического эксперимента и лежащие в основе оценивания параметров корреляционной функции распределения случайного поля неупругой реологической деформации, представлены в таблице.
Результаты эксперимента для построения корреляционной функции распределения случайного поля неупругой реологической деформации
х (мм) 0 2 4 6 8 10 12 14 16 18
K(x) 106 195 88 51 -6 31 28 -11 -99 -129 -115
х (мм) 20 22 24 26 28 30 32 34 36
K(x) 106 -98 -52 -102 -15 -3 44 10 93 99
Применение итерационного метода среднеквадратичного оценивания к выборке результатов эксперимента объема N = 19 и с шагом дискретизации t = 2 мм позволило получить следующие оценки параметров корреляционной функции: a = 0,014, b = 0,16 и a0 = 135 10-6.
Предельные относительные погрешности вычисления этих параметров соответственно составляют 191%, 19% и 46%. На рис. 11 точки 1 соответствуют результатам эксперимента, представленным в табл. 1. Кривая 2 является графиком корреляционной функции, восстановленной по экспериментальным данным. Для сравнения на этом же рисунке представлена кривая 3, описывающая корреляционную функцию, параметры которой найдены методом наименьших квадра-
N N
тов, примененным непосредственно к нелинейной функции (2): ^(к — ук) =У ek ^ min.
k=1 к=1
При этом решение нелинейной системы нормальных уравнений в МНК находилось методом Ньютона с помощью пакета поиска решений из приложения Microsoft Excel. Однако решить
задачу таким способом оказалось возможным только на основе первоначального приближения, полученного опять же итерационным численным методом среднеквадратичного оценивания коэффициентов ЛПДМ. Пусть Н0 - нулевая гипотеза: функции, представленные на рис. 11 кривыми 2 и одинаково хорошо аппроксимируют экспериментальные данные 1. Для проверки этой гипотезы найдено отношение остаточных дисперсий Г = 1,30 , которое имеет распределение Фишера. При уровне значимости у = 0,05 и соответствующих числах степеней свободы Гкр = 2,33 > Г . Следовательно, можно считать, что по отношению к экспериментальным данным различие между этими двумя функциями мало и обусловлено только случайной помехой в результатах наблюдений.
Р и с. 11 Восстановленная по эксперименталь- Р и с. 12 Моделирование корреляционной функ-
ным данным корреляционная функция распреде- ции распределения случайного поля реологиче-
ления случайного поля реологической деформа- ской деформации при исследовании ползучести
ции при исследовании ползучести материала материала
X
Проведено численно-аналитическое моделирование результатов данного эксперимента с целью анализа статистических характеристик погрешности оценок параметров корреляционной функции. При моделировании генерировался тестовый сигнал (2) с параметрами близкими к параметрам восстановленной по экспериментальным данным корреляционной функции:
а = 0,01, Ь = 0,2, а0 = 100 • 10-6, т = 2 и N = 19. В отсчеты тестового сигнала добавлялась случайная аддитивная помеха, мощность которой составляла 50% от мощности сигнала. При этом оценки 5 среднего квадратического отклонения (с.ко.) по остаточной сумме квадратов были соизмеримы с оценкой с.к.о. 5 = 34,7 10-6, вычисленной по результатам физического эксперимента. Исследовалась 1000 реализаций смоделированных таким образом функций. Одна из этих реализаций представлена на рис. 13 кривой 1, а выборка ее значений зашумленных случайной помехой - точками 2. Графики восстановленных по этим точкам корреляционных функций представлены на рис. 13 кривыми 3 и 4. Кривая 3 построена по результатам вычисления параметров корреляционной функции с помощью итерационного численного метода среднеквадратичного оценивания коэффициентов ЛПДМ. Кривая 4 соответствует функции, параметры которой найдены с помощью численной процедуры решения нелинейной системы нормальных уравнений (с выбором первоначального приближения на основе результатов, полученных итерационным методом среднеквадратичного оценивания). Очевидно, что обе кривые практически совпадают и мало отличаются от истинной кривой 1. Относительные погрешности определения параметров соответственно составили да = 56%, 80 = 0,8% и 8а0 = 21% для первого случая и 8а = 42% , 80 = 1,0% и 8а0 = 20% при втором методе вычислений. Остаточные
суммы квадратов оказались равными 0,017110-6 и 0,0167 10-6, соответственно, что также свидетельствует о статистической эквивалентности результатов оценивания обоими методами.
С целью сравнения погрешности вычисления параметров корреляционной функции в данном физическом эксперименте с предельными оценками относительной погрешности проведено усреднение относительных погрешностей вычисления параметров по 1000 реализаций: Sa = 85%, Sß = 4,8% и Sa0 = 16%. Видно, что эти значения в 2 - 3 раза меньше предельно допустимых.
Таким образом, в работе получены следующие научные результаты:
— построена линейно параметрическая дискретная модель, описывающая последовательность отсчетов корреляционной функции распределения случайного поля неупругой реологической деформации;
— описан алгоритм определения параметров корреляционной функции на основе среднеквадратичных оценок коэффициентов линейно параметрической дискретной модели в форме стохастического разностного уравнения;
— разработан итерационный численный метод среднеквадратичного оценивания коэффициентов ЛПДМ и проведены численно-аналитические исследования его эффективности;
— разработана методика и получены формулы, описывающие предельные относительные погрешности вычисления параметров корреляционной функции;
— приведены результаты апробации итерационного метода среднеквадратичного оценивания в физическом эксперименте при исследовании ползучести материалов.
Все это позволяет говорить о том, что разработан новый эффективный метод определения параметров корреляционной функции распределения случайного поля неупругой реологической деформации, в основе которого лежит итерационная процедура среднеквадратичного оценивания коэффициентов линейно параметрической дискретной модели в форме стохастического разностного уравнения, и который может быть рекомендован для эффективного решения задач параметрической идентификации различных по физической природе объектов, динамические процессы в которых описываются затухающей по экспоненте синусоидой.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Радченко В.П., Дудкин С.А., Тимофеев М.И. Экспериментальное исследование и анализ полей неупругих микро-и макронеоднородностей сплава АД-1. // Вест. Самарского госуд. техн. ун-та. Сер. Физ.-мат. науки. 2002. Вып. 16. С. 111-117.
2. Попов Н.Н., Самарин Ю.П. Исследование полей напряжений вблизи границы стохастически неоднородной полуплоскости при ползучести. // ПМТФ. 1988. №1. С. 159-164.
3. Должковой А.А., Попов Н.Н., Радченко В.П. Решение стохастической краевой задачи установившейся ползучести для толстостенной трубы методом малого параметра // ПМТФ. 2006. Т. 47. №1. С. 161 - 171.
4. БоглаевЮ.П. Вычислительная математика и программирование. М.: Высш. шк., 1990. 544 с.
5. Зотеев В.Е. Разработка и исследование линейных дискретных моделей колебаний диссипативных систем. // Вест. Самарского госуд. техн. ун-та. Сер. Физ.-мат. науки. 1999. Вып. 7. С. 170-177.
6. ВучковИ., БояджиеваЛ., Солаков О. Прикладной линейный регрессионный анализ. М.: Финансы и статистика, 1987. 238 с.
7. Грановский В.А., Сирая Т.Н. Методы обработки экспериментальных данных при измерениях. Л.: Энергоатомиз-дат. 1990. 288 с.
Поступила 15.01.2006.
