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

Определение параметров испытательных импульсов на основе стохастических разностных уравнений Текст научной статьи по специальности «Математика»

CC BY
93
36
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АМПЛИТУДНО-ЧАСТОТНАЯ ХАРАКТЕРИСТИКА / ЛИНЕЙНО ПАРАМЕТРИЧЕСКАЯ ДИСКРЕТНАЯ МОДЕЛЬ / СТОХАСТИЧЕСКОЕ РАЗНОСТНОЕ УРАВНЕНИЕ / AMPLITUDE-FREQUENT CHARACTERISTIC / LINEAR PARAMETRIC DISCRETE MODEL / STOCHASTIC DIFFERENCE EQUATIONS

Аннотация научной статьи по математике, автор научной работы — Заусаева Мария Анатольевна, Зотеев Владимир Евгеньевич

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

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

Похожие темы научных работ по математике , автор научной работы — Заусаева Мария Анатольевна, Зотеев Владимир Евгеньевич

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

Definition of Trial Impulses Parameters on the Basis of Stochastic Difference Equations

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.

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

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. — 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

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

Лз — Л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.

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

УДК 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

Дрозденко Константин Александрович — аспирант кафедры прикладной математики и информатики Самарского государственного технического университета.

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