Аппроксимация корреляционных функций и спектральных плотностей мощности ортогональными функциями Сонина-Лагерра Текст научной статьи по специальности «Математика»

УДК 517.587:519.216

АППРОКСИМАЦИЯ КОРРЕЛЯЦИОННЫХ ФУНКЦИЙ И СПЕКТРАЛЬНЫХ ПЛОТНОСТЕЙ МОЩНОСТИ ОРТОГОНАЛЬНЫМИ ФУНКЦИЯМИ СОНИНА-ЛАГЕРРА

С. А. Прохоров, И. М. Куликовских

Самарский государственный аэрокосмический университет им. ак. С. П. Королева,

443086, г. Самара, Московское ш., 34.

E-mails: sp@smr.ru; eirisan@rambler.ru

Исследованы аппроксимативные возможности ортогональных функций Сони-на—Лагерра при заданном параметре ортогонального базиса. Оценены параметры аппроксимирующего выражения, используемые при построении моделей корреляционной функции и спектральной плотности мощности, в соответствии с минимумом взвешенной квадратической погрешности аппроксимации.

Ключевые слова: аппроксимация, корреляционная функция, ортогональный базис Сонина—Лагерра, преобразование Фурье, случайный процесс, спектральная плотность мощности.

Применение ортогональных функций Сонина—Лагерра для аппроксимации и восстановления корреляционных функций (КФ) и спектральных плотностей мощности (СПМ) связано с их широкими аппроксимативными возможностями.

Проведём корреляционно-спектральный анализ в ортогональном базисе Сонина—Лагерра с изменяемым параметром а при значениях а = 0,1, 2.

Представим ортогональные функции Сонина—Лагерра в следующей форме:

Ьа) (т,7) = С- ± в-*, (1)

Ск+а з=0 *•

значение функций в «нуле» которых (0,7) = 1, а норма

{Г(а+1) 7

' при к = т; (2)

0, при к = т.

Представим модель КФ с параметрами 7, во, въ • • •, вт в виде

т т

Ка(т) = вк 4°° (т,7)1(т)+^ А 4“)(-Т,7)1(-т) , (3)

к=0 к=0

где а%. — дисперсия процесса, а коэффициенты разложения определяются так:

x

со

в — (а) f Kx(T)Lka) (т,7)^(“) (т) dr. (4)

ILk II о

Прохоров Сергей Антонович — зав. кафедрой информационных систем и технологий Самарского государственного аэрокосмического университета; д.т.н., проф.

Куликовских Илона Марковна — аспирант кафедры информационных систем и технологий Самарского государственного аэрокосмического университета.

Здесь ^(а) (т) = та —вес ортогональной функции.

Определив коэффициенты разложения вк, необходимо решить задачу определения параметра масштаба 7, обеспечивающего минимальную допустимую погрешность аппроксимации. Для выбранного ортогонального базиса данная задача решена в [1].

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

Д = .1

0

п 2

Kx(т) -Y^ ek^(Q° (т) Лт- (Б)

k=0

Однако коэффициенты (4) не обеспечивают выполнения основного свойства КФ:

т

ка(0) = а% ^ вкьк“)(0,7) = • (6)

к=0

Для обеспечения последнего условия вводят коэффициенты разложения Ьк, для которых

Pa(0) = Е bkLk“)(0,7) = 1, (7)

k=0

где ра(т) = Ка(т)/^2 —нормированная корреляционная функция (НКФ).

Таким образом, модель КФ (3) можно представить как модель НКФ с учётом ограничения (7) в следующем виде:

m m

Ра(т) = [Е bk 4“)(т,7)1(т) + Е bk 4“)(-т,7)1(-т) • (В)

k=0 k=0

В работе [2] приводится методика расчета коэффициентов разложения и оценка значений погрешностей аппроксимации КФ при ограничениях на её модель. Для удобства погрешность аппроксимации представляют в виде суммы Д = Ді + Д2, где Ді —погрешность аппроксимации, вызванная конечным числом членов разложения m; Д2 —погрешность аппроксимации, вызванная дополнительным ограничением, обеспечивающим выполнение основного свойства КФ.

Выражения для коэффициентов разложения bn и погрешностей аппроксимации Ді, Д2 для выбранных значений а представлены в табл. 1.

Построим ортогональную модель автокорреляционной функции (АКФ) идеального полосового шума px (т) = (sin 0,5т/0,5т) cos 5т в базисе Сонина— Лагерра с параметром а = 1 (см. рис. 1).

Учитывая выражение (3), оценим СПМ:

т

k=0

Ла)

2m

Sa(^) = JE & Re wka)(jW), (9)

где Wk )(jw) —преобразование Фурье ортогональных функций. 186

ГО

Т аблица 1

Коэффициенты разложения и погрешности аппроксимации КФ

Дъ А2

ЬП

Ьп вп +

/ т \

(1 ~ 12 рк)

у к=0 7

(т + 1)

/1

рХ(тМт - тЕв2’ 0 к=0 / т \ 2

(1 — I] вк)

' г-—л '

Д2 = “

7 (т + 1)

Ьп вп +

/ т \

2(1- — 2 рк)(п +1) 4 к=0 7________

(т + 1)(т + 2)

/1 X \

рХ(тМтМт — 72 Е

0 к=0

/ т \2

2 (1 — I] вк)

2

72 (к + 1)?

2 =

72 (т + 1)(т + 2)

Ьп рп + "

/ т \

3 ( 1 — Й (п + 1)(п + 2)

к=0

/* 4

Д1= /рХ (т )р(т Мт —

У 1„____п

в2

(т + 1)(т + 2)(т + 3)

Д

12

73 к= (к+1)(к+2)’

/ т \ 2

(1 — 2 вк)

V — л /

2 =

73 (т + 1)(т + 2)(т + 3)

а

оо

0

о

1

о

2

Рис. 1. Аппроксимация КФ функциями Сонина—Лагерра (а =1); 7 = 4, т = 250. Относительная погрешность аппроксимации 6 = 0,12

Выражения для действительной и мнимой части преобразования Фурье ортогональных функций представлены в работе [1]. Однако форма их представления имеет недостатки, связанные с необходимостью вычисления комбинаторного числа сочетаний, что не позволяет вычислять данные выражения для высоких порядков. В табл. 2 представлена альтернативная форма записи данных выражений, свободная от указанного недостатка.

Оценка СПМ, записанная по параметрам модели КФ (см. рис. 1), показана на рис. 2.

Другим способом оценки СПМ (9) является построение её аппроксимативной ортогональной модели:

£аМ = ск4°° (^,7).

к=0

(10)

С учётом чётности СПМ модель (8) должна гарантировать выполнение

Таблица 2

Преобразование Фурье ортогональных функций Сонина—Лагерра

о

— (— 1)к СОБ ф СОБ [(2к + 1)ф]

7

(—1)к+1 СОБ фБІП [(2к + 1)ф]

1

1 + ( —1) соБ[2(к + 1

7 • (к + 1)

2 /1 (—1)к сов[(2к+3)ф]

(—1)к+1 зіп [2 (к + 1) ф])

7-(к+1)(к+2) 42

7Т +

2 СОБ ф

+к+1

7 • (к + 1)

2 /(—1)к+1 віп[(2к+3)ф] tg ф^

Т- (к+1) (к+2)

2 СОБ ф

2

а

2

0

1

2

Рис. 2. Модель СПМ с характеристиками шє = 5,3, Б(о>є) = 0,6771,

Д^є = 5,6692

189

Таблица 3

Коэффициенты разложения и погрешности аппроксимации СПМ

а

Ста

ДЪ Д2

Ста — вга +

4 - £ А (-1Г ____&=0_______

т + 1

(-1)П

« 1 Ш Л1 — ^Хн М ^ - в2

п 7 к=0

д.

4 -£ А (-1)к ■=0

2—

7

т+1

СП ----- в«, +

4 - £ ■=0

вк [(к+1)шоё 2] (к+1)

Ш [(к+1)Шоё 2]

£ и

■=0

[(п + 1) шоё 2]

/1 у \

^Хн М ^ М ^ - ^2 Е

Г к=0 (к + 1)!

д

7 Ш вк [(й+1)Шоё 2]

4 - 2^-------------------

2—

■=0

(к+1)

^2 Ш [(^+1)Шоё 2]

■=0 (1+1)

■=0

Сп — вга +

Ш

2 - V 8 ^

■ =0

вк [(к+2)ё1у 2] (к+1)(к+2)

£

■ =0

[(■+2)ё1у 2]2 (■+1)(к+2)

[(п + 2) &у 2]

го

4 Ш

д1 — у ^Хн М ^ М ^ -^

в2

д

2—

8 - £

■ =0

г к=0(к +1) (к +2)

Ш вк [(к+2)ё1у 2] (■+1)(к+2)

7'

зШ

£

■ =0

[(■+2)ё1у 2] (■+1)(к+2)

Ш

0

2

1

Ш

го

1

2

1

2

Ш

2

4

Аппроксимация корреляционных функций и спектральных плотностей мощности ...

условия нормировки [2]:

ГО

г т 2

I ск—(0) = .

0 к=0

Выражения для коэффициентов разложения сп и погрешностей аппроксимации Д1, Д2 для выбранных значений а представлены в табл. 3.

Оценим КФ по параметрам модели СПМ [2]:

ШЯе Ш1т

Каху (т)— 2^^ вк Яе ■ Ке ^ Ке(^Т ) + ^ вк 1т ' 1т ^ 1т(?Т) , (11)

■ =0 ■ =0

где —преобразование Фурье ортогональных функций;

го

Яе — и і ||2 I (и)^д. тге(и) аЯе)

№Яе||2 У 0 го

АІт — ||2 I 1т (и)^А:1т(ич аІт)

Іт 1

0

Представим решение обратной задачи: построим ортогональную модель СПМ и по параметрам модели оценим КФ.

На рис. 3 представлена аппроксимация СПМ (см. рис. 2) с применением ортогональных функций Сонина—Лагерра (а — 0).

По параметрам модели СПМ оценивается КФ. На рис. 4 отражены результаты решения прямой (аппроксимация КФ и оценка СПМ) и обратной (аппроксимация СПМ и оценка КФ) задач.

Рис. 3. Аппроксимация СПМ функциями Сонина—Лагерра (а — 0); 7 — 20, т — 400. Относительная погрешность аппроксимации 6 — 0,058

Рис. 4. Наложение исходной КФ и КФ, полученной после аппроксимативного корреляционно-спектрального анализа

Совпадение моделей исходной и результирующей АКФ (см. рис. 4) свидетельствует об адекватности предложенных ортогональных моделей и правильности разработанных алгоритмов аппроксимативного корреляционно— спектрального анализа.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Прохоров, С. А. Частотные характеристики ортогональных функций Сонина—Лагерра [Текст] / С. А. Прохоров, И. М. Куликовских // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. — 2007. — № 2 (15). — C. 123-127.

2. Прикладной анализ случайных процессов [Текст] / под ред. С. А. Прохорова. Самара: СНЦ РАН, 2007. — 582 с.

Поступила в редакцию 12/V/2008; в окончательном варианте — 14/X/2008.

MSC: 93B40, 42C05, 33C45

APPROXIMATION OF CORRELATION FUNCTION AND POWER SPECTRAL DENSITY WITH SONIN-LAGUERRE ORTHOGONAL FUNCTIONS

S. A. Prokhorov, I. M. Kulikovskikh

S. P. Korolyov Samara State Aerospace University,

443086, Samara, Moskovskoe sh., 34.

E-mails: spSsmr.ru; eirisanarambler.ru

Approximate capabilities of Sonin-Laguerre orthogonal functions with predefined parameter of orthogonal basis are studied. Parameters of approximation expression are evaluated, that are used for construction of the models of correlation function and, power spectral density in compliance with the minimal weighted quadratic error of approximation.

Key words: approximation, correlation function, Sonin-Laguerre orthogonal basis, Fourier transformation, random process, power spectrum density.

Original article submitted 12/V/2008; revision submitted 14/X/2008.

Prokhorov Sergey Antonovich, Dr. Sci. (Tech), Prof. Head of Information Systems and Technologies Chair of Samara State Aerospace University.

Kulikovskikh Ilona Markovna, Postgraduate Student of Information Systems and Technologies Chair of Samara State Aerospace University.