Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. — 2008. — № 2(17). — С. 262—266
Математическое моделирование
УДК 681.5.015
ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ИСПЫТАТЕЛЬНЫХ ИМПУЛЬСОВ НА ОСНОВЕ СТОХАСТИЧЕСКИХ РАЗНОСТНЫХ УРАВНЕНИЙ
М. А. Заусаева, В. Е. Зотеев
Самарский государственный технический университет,
443100, Самара, ул. Молодогвардейская, 244.
E-mails: zausmasha@mail.ru; zoteev@pm.samgtu.ru
Рассматривается задача параметрической идентификации типовых ударных воздействий по результатам наблюдений их амплитудно-частотной характеристики. Эта задача решается на основе стохастических разностных уравнений, описывающих результаты измерений мгновенных значений амплитудного спектра.
Ключевые слова: амплитудно-частотная характеристика, линейно параметрическая дискретная модель, стохастическое разностное уравнение.
Для задания испытательного воздействия в промышленном эксперименте обычно используют записи реальных ударных процессов, полученные с помощью специальной записывающей аппаратуры в условиях нормальной эксплуатации. Для последующего их моделирования и воспроизведения на испытательных установках часто используют амплитудный спектр (амплитудно-частотную характеристику) ударного воздействия, который позволяет правильно выбрать структуру и характеристики (в частности, полосу пропускания) виброиспытательного комплекса [1].
Рассматривается задача параметрической идентификации типовых ударных воздействий по результатам наблюдений их амплитудно-частотной характеристики. Основные типы испытательных импульсов при простом и сложном колебательном ударах, их временные характеристики x(t) и амплитудные спектры а(ш) приведены в таблице [1].
Задача оценки параметров четырех различных воздействий, представленных в таблице, решается на основе стохастических разностных уравнений, описывающих результаты измерений а^ мгновенных значений амплитудного спектра реального ударного процесса.
Амплитудно-частотную характеристику ш прямоугольного импульса будем рассматривать на промежутке (0, Ip), а полусинусоидального импульса — на промежутке (0, . В этом случае амплитудные спектры могут быть описаны соответственно
функциями
2AT cos Цр . _
, , 2A . шТ . . I — 1-, при ш = т,
а(ш) =—sin — и а(ш) = s п A- п2
2
при Ш = р .
Мгновенные значения = а(Дшк) (к Є N) амплитудных характеристик, где
Зотеев Владимир Евгеньевич — доцент кафедры прикладной математики и информатики Самарского государственного технического университета; к.ф.-м.н., доцент.
Заусаева Мария Анатольевна — аспирант кафедры прикладной математики и информатики Самарского государственного технического университета.
Типичные ударные воздействия и их спектры
Воздействие ж(£) а(а>)
Прямо- угольный импульс ( 0, при 4 < 0, < А, при 0 ^ 4 ^ Т, [ 0, при 4 > Т. АТ 2 Бт ^ —Т
Пилооб- разный импульс Г 0, при 4 < 0, < А Т, при 0 ^ 4 ^ Т, [ 0, при 4 > Т.
“ А/1 2Т -1п ^Т + 81п —Т — у — Т ' —Т 2
Полусину- соидальный импульс Г 0, при 4 < 0, < А -1п Т, при 0 ^ 4 ^ Т, [ 0, при 4 < Т. Т 2 СОБ ■ЦТ
1 ^2Т2 1 п2
Сложный колеба- тельный удар Г 0, при 4 < 0, | е-“* 008 (^о4 + ф), при 4 ^ 0. — •у/—4-2(—2-а2) —2 + (—2+а2)2
Дш — шаг равномерной дискретизации, для представленных выше типовых ударных воздействий описываются соответственно дискретными функциями вида
2А Д^&Т
= дл ’ (1)
А2 - 2 й1п Дц^Т | 2 -111 Д~';' 7 ) , (2)
Л , . 1 Л , . О / 5 V /
д^2&2 V Д^т д^т
_ 2АТ 008 ; = п (3)
ак = "V 1 - А—22т2 ’ к = Д^г (3)
п2
~2 д^2&2
ак =--------------------------------------------2. (4)
Д^4&4 — 2 (^2 — а2) Д^2^2 + (^2 + а2)
2
К
С помощью несложных алгебраических преобразований построены линейно параметрические дискретные модели, связывающие в виде рекуррентных формул мгновенные значения дискретных функций (1)—(8).
При обработке экспериментальных данных формируется выборка результатов измерений ак — ак + £к, которые содержат случайную помеху £к (& = 1, 2,..., N, где N — объём выборки). В этом случае построенные линейно параметрические дискретные модели принимают вид стохастических разностных уравнений:
- для прямоугольного импульсного воздействия:
а1 — Л2 + £1, а2 = Лз + £2,
+ (1 — I) ак-+а —
ак
— 2а — I) ак )£к-2 — 2(1 —
3,4,
-1Л1 +
!) Л1£к-
(5)
1 + £к ,
где
Д^Т 2А Д^Т
------, Л2 — —-------- 81П ------,
2 ’ 2 Д^ 2 ’
А
Л3 — —— -1п Д ^Т; Да>
Л
1
для пилообразного импульсного воздействия:
= Лк+3 + 2ак £к,
к — 1, 2,
,6;
к3ак +(к — 6)3«|-б— Л1 [(к — 1)3ак-1 +(к — 5)3ак—5] +
где
+Л2 [(к — 2)3акк-2 +(к — 4)3«1—4] + Лз(к - 3)3ак-з+
+ 2(к 6) ак-6£к —6 2Л1(к 5) ак—5£к—5
—2Л2(к — 4)3ак-4£к-4 — 2Лз(к — 3)3ак—з£к—з —
—2Л2(к — 2)3ак-2£к-2 — 2Л1(к — 1)3ак—1£к—1+ +2к3ак£к, к — 7, 8, .. ., N
ДшТ\
Л1 — 2(1 + 008 ДшТ + 008 2— ) ,
ДшТ ДшТ \
Л2 — — (3 + 4 008 ДшТ + 4 008 —2----+ 4 008--2— 008 ДшТ I
ДшТ ДшТ \
Лз — 4(1 + 008 ДшТ + 008----2---+ 2 008--2— 008 ДшТ I ,
(6)
А2
Дш2к2
22 81п ДшкТ +
ДшкТ \
ДшкТ
ДшкТ
81п
Лк+з — -—2л I 1 —
- для полусинусоидального импульсного воздействия:
2
к — 1, 2,..., 6;
а1 — Л4 + £1, а2 — Л5 + £2,
ак — ак-2 — 2ак-1Л1 + [к2ак + (к — 2)2ак—2] Л2 — 2(к — 1)2ак—1Л3+ (7)
+ [1 — (к — 2)2Л^ £к-2 + 2 [(к — 1)2Лз — Л1] £к-1+
+(1 — к2Л2)£к, к — 3,4,..., N,
где
ДшТ
Л2
Дш2Т2
Лз — Л1 Л2 ,
л 2АТ 008 Дг
Л4 —-----------------2
п 1 -
Д—2Т 2
Л5 —
2АТ 008 ДшТ
п1
4Д—2Т2
для сложного колебательного удара:
а2 — Л2 + 2а1£1,
к2а|-1 — (к — 1)2а| — Дш2 (4к3 — 6к2 + 4к — 1) а|а|—
22 к- 1
22 — а2 а2
к“к—1(1 — 2к)Л1 +
+ (2к2ак-1 — 2Дш2 (4к3 — 6к2 + 4к — 1) а|ак-1 —
2Л1(1 2к)акак—1) £к-1
— (2(к — 1)2ак — 2Дш2(4к3 — 6к2 + 4к — 1)ак а|—1 —
— 2Л1(1 — 2к)акак—1) £к,
к — 2, 3, .. ., N
где
Л1
2 (ш2 — а2) , Л2
Дш2
Дш4 — 2 (ш2 — а2) Дш2 + (ш2 + а2)
2
(8)
Описанные выше линейно параметрические дискретные модели можно представить в виде обобщенной регрессионной модели:
Ь — ^Л + п, П — Рл £,
2
а
Л
1
2
2
п
2
2
п
п
где Л — вектор неизвестных коэффициентов линейно параметрической дискретной модели; є — N -мерный вектор случайной помехи в результатах наблюдений; п — N -мерный вектор эквивалентного случайного возмущения в стохастическом разностном уравнении; Ь — N -мерный вектор свободных членов системы уравнений; Р — матрица регрессоров, элементы которой содержат результаты измерений амплитудно-частотной характеристики; Рд — матрица линейного преобразования случайной помехи в стохастическом разностном уравнении эквивалентного возмущения. Элементы векторов и матриц, входящих в систему (9), несложно описать, исходя из формул (5)-(8).
В основе параметрической идентификации ударного воздействия по результатам измерений амплитудного спектра лежит среднеквадратичное оценивание коэффициентов соответствующего стохастического разностного уравнения. В этой связи наиболее эффективным является алгоритм, который включает итерационную процедуру уточнения среднеквадратичных оценок на основе преобразования разностного уравнения, описывающего результаты измерений амплитудного спектра [2]:
Рд-1Ь = Р-^Л + є.
В основе итерационной процедуры лежит минимизация функционала
Ь - Р-1. РЛ л(к) л(к)
=Ф шт,
2
2
є
где Р^) — обратная матрица, элементы которой зависят от коэффициентов Л^к). Очевидно, что данный функционал представляет собой квадратичную форму относительно искомых коэффициентов Л^-. Следовательно, он достигает своего неотрицательного минимума. При этом нетрудно показать, что минимум данного функционала достигается в точке
Л(к+1} — (рт^—к) р)—1 к) Ь, %к) — Рл(к)РТ). (10)
На первом шаге итерационной процедуры находится начальное приближение:
Л(0) — (р Тр) —1 р Т ь.
Затем по формуле (9) при Л(к) — Л(0) вычисляется следующее приближение:
Л(1) — (р Ч-1р)—1 р X>
Оно вновь подставляется в правую часть формулы (10), и находится новое приближение Л(2) и т. д.
Достаточное условие сходимости итерационной процедуры определяется на основе принципа сжимающихся отображений аналогично методу простых итераций решения систем нелинейных уравнений [3]. Проведенные численно-аналитические исследования показали хорошую сходимость метода: 2-3 итерации.
Найденные среднеквадратические оценки Л^- коэффициентов линейно параметрических дискретных моделей (5)-(8) лежат в основе вычисления оценок параметров амплитудных спектров. При определении параметров прямоугольного импульса следует воспользоваться соотношениями
- 2 Л л Л2^^
Т = ——агееов Л1, А =
Параметры пилообразного импульса можно вычислить по формулам:
1 / Лі — 2 + \/Л1 + 4Л1 + 4Л2
Т = —— агссоэ --------------------------------------
Д^ I 4
А =
Л4Д^2 агеео8(Л)
агссоэ (Л) + 2 яіп агссоэ (Л)] — 2 эш [агссоэ (Л)]'
д Л1—2+л/Л?+4 А? +4 А 2 «
где Л — -------1—4--------, параметры полусинусоидального импульсного воздей-
ствия выражаются через коэффициенты линейно параметрической дискретной модели (8) по формулам
2 Л4 [п2 — Д^2Т21
Т — ——агееов Л1, А —--------------7------,
Д^ 2пТЛ1
а параметры сложного колебательного удара вычисляются так:
2Д^
Лі + л1— Д^2 — Лі
2Д^
Лі + л2 — Д^2 + Лі
2
Проведены численно-аналитические исследования эффективности разработанного метода параметрической идентификации типовых ударных воздействий. На рисунке кривая 1 описывает истинную амплитудно-частотную характеристику, а точки 2 соответствуют смоделированным результатам наблюдений ак = = ак + єк (к = 1,2,..., 50) амплитудночастотной характеристики сложного колебательного удара при о>о = 1, а = 0,5, Д^ = 0,1, є = 10%. Кривая 3 на рисунке описывает восстановленную с помощью разработанного метода параметрической идентификации амплитудно-частотную характеристику.
Полученные результаты свидетельствуют о высокой помехоустойчивости разработанного метода, в основе которого лежит среднеквадратичное оценивание коэффициентов линейно параметрической дискретной модели в форме стохастических разностных уравнений. Такой подход к решению задачи позволяет практически устранить смещение в оценках и тем самым добиться высокой точности определения параметров ударных воздействий по экспериментальным данным.
Истинная амплитудно-частотная характеристика, смоделированные результаты наблюдений и восстановленная по ним амплитудно-частотная характеристика сложного колебательного удара
2
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Вибрации в технике [Текст]: Справочник в 6 т.—М.: Машиностроение, 1981.—Т. 5.— 496 с.
2. Зотеев, В. Е. Параметрическая идентификация линейной динамической системы на основе стохастических разностных уравнений [Текст] / В. Е. Зотеев // Матем. моделирование. — 2008. — Т. 20, № 9. — С. 120-128.
3. Волков, Е. А. Численные методы [Текст] / Е. А. Волков. —М.: Наука, 2003. —256 с.
Поступила в редакцию 29/ІХ/2008; в окончательном варианте — 09/ХІ/2008.
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. — 2008. — № 2(17). — С. 267—270
MSC: 65P40, 34C15, 37M05
DEFINITION OF TRIAL IMPULSES PARAMETERS ON THE BASIS OF STOCHASTIC DIFFERENCE EQUATIONS
M. A. Zausaeva, V.E. Zoteev
Samara State Technical University,
443100, Samara, Molodogvardeyskaya str., 244.
E-mails: zausmasha@mail.ru; zoteev@pm.samgtu.ru
Problem, of parametrical identification of sample percussion influences according to the results of their amplitude-frequent characteristic observation is studied. This problem is solved with the help of stochastic difference equations describing the results of measurements of the instantaneous values of amplitude spectrum.
Key words: amplitude-frequent characteristic, linear parametric discrete model, stochastic difference equations.
Original article submitted 29/IX/2008; revision submitted 09/XI/2008.
Zoteev Vladimir Eugenievich, Ph.D. (Phis. & Math.) Assist. Prof., Dept. of Applied Mathematics and Computer Science of Samara State Technical University.
Zausaeva Mariya Anatolievna, Postgraduate Student, Dept. of Applied Mathematics and, Computer Science of Samara State Technical University.
УДК 519.856
УСТОЙЧИВОСТЬ МЕТОДА СТОХАСТИЧЕСКОГО ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ ДИСКРЕТИЗАЦИИ
К. А. Дрозденко
Самарский государственный технический университет,
443100, Самара, ул. Молодогвардейская, 244.
E-mail: mexicanose@rambler.ru
Доказана теорема устойчивости метода стохастического динамического программирования с переменными параметрами дискретизации.
Ключевые слова: стохастическое программирование, динамическое программирование, распределение ресурсов, переменные параметры дискретизации, устойчивость метода.
Рассмотрим метод стохастического динамического программирования распределения ресурсов с переменными параметрами дискретизации. Модель распределения ресурсов в стохастических условиях имеет следующий вид:
(sj-i,q)= sup {p(Zj+i,q) [/j+1,9) + Zj+i,q]} ;
Xj ex
J i
) = 1 = C;
j=0 i=0
Дрозденко Константин Александрович — аспирант кафедры прикладной математики и информатики Самарского государственного технического университета.