УДК 519.2
ОЦЕНКА СПЕКТРА МОЩНОСТИ СТАЦИОНАРНОГО СЛУЧАЙНОГО ПРОЦЕССА СПЛАЙНОМ ПЕРВОГО ПОРЯДКА
Устинова Ирина Георгиевна,
канд. техн. наук, доцент каф. высшей математики физико-технического института ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», Россия, 634050, г. Томск, пр. Ленина, д. 30. E-mail: igu@sibmail.com
Лазарева Елена Геннадьевна,
канд. физ.-мат. наук, доцент каф. общей математики ФГБОУ ВПО «Национальный исследовательский Томский государственный университет», Россия, 634050, г. Томск, пр. Ленина, д. 36.
E-mail: lazareva-e-g73@mail.ru
Подберезина Елена Ивановна,
канд. физ.-мат. наук, доцент каф. высшей математики физико-технического института ФГБОУ ВПО «Национальный исследовательский Томский политехнический университет», Россия, 634050, г. Томск, пр. Ленина, д. 30. E-mail: pei@tpu.ru
Актуальность работы обусловлена тем, что спектр мощности, так же как и функция корреляции, является одной из важнейших характеристик второго порядка случайного процесса, так как он показывает, какого рода гармоники преобладают в данном процессе, какова структура процесса; позволяет получить оценки спектрального состава изучаемых полезных сигналов и помех. По спектрам возможен синтез (восстановление) сигнала, а также построение линейных, в том числе оптимальных фильтров, а также получение оценки погрешности линейной фильтрации.
Цель работы: Оценка спектра мощности стационарного случайного процесса в виде сплайна первого порядка при следующих схемах измерений: когда в каждый момент времени производится ровно одно измерение и когда в каждый момент времени производится случайное число измерений. Исследование статистических характеристик полученных оценок Методы исследования: Для расчетов использованы методы теории вероятностей и математической статистики.
Результаты: Получена несмещенная оценка спектра мощности в виде сплайна первого порядка в двух схемах измерения: в каждый момент времени производится ровно одно измерение; в каждый момент времени производится случайное число измерений.
Показано, что дисперсия оценок асимптотически ведет себя как —, где T - время наблюдения.
Спектр мощности, корреляционная функция, сплайн первого порядка, оценки параметров, статистические свойства оценок
Введение
Цифровой спектральный анализ является мощным современным методом обработки сигналов [1]. В его основании лежит оценка и изучение спектра мощности сигнала, принимаемого за некоторый случайный процесс [2]. Функция плотности спектра мощности определяет распределение дисперсии случайного процесса по частотам [3].
Для нахождения оценки функции плотности спектра существует следующие основные способы:
S(а) = — J^[t]cos(«t)dT затруднительно, ибо этот интеграл не всегда вычисляется в элементарных функциях, и корреляционная функция известна, как правило, на дискретном множестве переменной т [4-8].
£(а) = 5С2Р , где с - скейлинговая константа,
зависящая от амплитуды сигнала; Л - фрактальная размерность, причем эта взаимосвязь носит оценочный характер, ибо она является асимптотическим решением интегрального уравнения [13-16].
Целью данной работы является нахождение оценки спектра мощности стационарного случайного процесса в виде сплайна первого порядка. Обычно сплайн определяют следующим образом [17-19]. Функцию Б(Ь), определенную и непрерывную на отрезке [а,Ь], называют полиномиальным сплайном порядка т с узлами е{а=х0<х1<е2<...<хп=Ь}, если на каждом из отрезков [хн,х^, ]=1,п Б(Ь) есть алгебраический полином степени т.
Постановка задачи
Пусть значения процесса у(Ь) измеряются на отрезке времени [0;Т]. Моменты измерений Ь, 1=1,N известны точно. Будем полагать, что у=у(Ь )+£;, где £ - независимые случайные величины, причем, М[<^]=0, Б[х1]=аг, М[у(Ь)]=0, МЫЬМЬ )]=Ф-Ь ].
По результатам наблюдений требуется построить оценку 8(а) спектра мощности Б(а) в виде сплайна первого порядка.
Построение оценки спектра мощности в виде сплайна первого порядка
Разобьем всю ось частот а на отрезки [0;П], [0;20], [20;30],.... Рассмотрим статистику
г,]; г ф]
где для функции ф(т) справедливо условие:
<р(т) = |ф (ю)созютс!ю. (2)
Заметим, что
М[2] = -’- X КУ] - ^М(] - Ь].
ПГ /,]; г Ф ]
Можно считать, что имеет место следующее приближение:
М[2] * — |^Я[ы - V]^(и - v)dиdv.
пГ 0 0
Сделав замену переменных, получим
JJR[u - v ]^(u - v )dudv =
= Jdv J R(z)^(z)dz.
u - v = z du = dz
Поменяем порядок интегрирования, используя область интегрирования, изображенную на рис. 1.
Fig. 1. Integration domain
Это позволит от исходного двойного интеграла перейти к однократному интегралу
J R[ z]<p( z)dz J dv + J R[ z ]ф( z)dz J dv =
— T —z 0 0
= J (T + z )R[ z ]ф( z )dz + J (T — z )R[ z ]ф( z )dz =
= Г |^1 - ГГ|я[2]ф(2).
В асимптотическом случае, когда Т^х>, выра-12 „
жение ^ ^ 0 и, соответственно,
Г ( II ^ "
| ^ 1 - ^ Я[2]ф(2)йЪ * |Д[2]ф(2)йЪ.
Таким образом
М[2] * — | Д[г]ф(2)<£ = - |Я[г]ф(2)<£. Воспользовавшись выражением (2), получим
+» +х
M[Q] ^ JR[t] JФ(a) cosardadr
+» +х
Jф(a)da JR[r]cosardT.
Используя известное соотношение Вине-ра-Хинчина [2], окончательно получим
M[Q] - J^(a)S(a)da.
Для построения оценки £(ю) спектра мощности Б(ю) в виде сплайна первого порядка будем полагать, что на отрезке [(к-1)П;Ш] функция 5(ю) представима в виде
кО-ю ю-(к -1)0
^(ю) ^к-1 — + ^
Рассмотрим интеграл
I = J S(a) I A
kQ- a a - (k - 1)Q ,
+ B----——— Ida,
і 2 kQ - a Q
a - (k - 1)Q
e [(k - 1)Q; kQ],
о a g[(k - 1)Q; kQ].
Используя выражение (2), находим p(r) = J — f" 3k -1 - 3 al cosarda =
& 3 іcos(k - 1)Qr -)
Qt2 I-cos kQr J
a e[(k - 1)Q, kQ],
а £[(к -1)0, кО].
Полагаем, что при к=1 <р(т)=ф0(т), при к=2 ф(т)=^1(т),.. Подставляя в выражение (1) найденные значения <р(г), находим последовательно БоБ!,..., которые соединяем отрезками прямых, что и дает оценку Б(а) в виде сплайна первого порядка. Из построения оценок коэффициентов сплайна следует, что полученные оценки узлов сплайна являются несмещенными.
Оценка дисперсии полученных оценок
спектра мощности
Найдем теперь асимптотическую оценку дисперсии Д0] статистики 2 при Т^да. Для этого рассмотрим
q2 =-1г X X X Ху^у/ул^— t Mt —)^ t j=i ,■=i i=i k=i
Усредняя по реализациям процесса при фиксированных моментах измерений, получим
M [ ] = M [ .y (tj) y(ti) 7(t) Я4)] =
= R[tj—h ]R[t,— tk] + R[tj — t ]R[t — 4] +
+R[tj— tk ]R[h— h]Тогда
м [Q2] =
—2^2 п T
вычисление которого приводит к выражению
-{£*_1(2 А + В) + £ (А + 2 В)}.
Заметим, что приА=2 и Б=-1 I = — £к-1. Тогда Ф(а) выражения (3) будет иметь вид:
ХХХХр(?/ -p(t, - У/=1 i=1 1=1 k=1
іR[t/ - ti]R[t1 - tk] + R[t,. -1,]R[ti - tk] +)
v+R[t/ - tk ]R[ti - t1]
Далее следует рассмотреть усреднение по моментам измерений каждого слагаемого, содержащего сомножитель вида й[^-^ ] Д [ЬгУ выражения (4). Здесь возможны следующие варианты:
п2Г2 00001хф(м - и&)ф^ - v&)dudvdu&dv&J’
вычисление которого дает при T^<ю квадрат математического ожидания статистики 2, а именно
—2^ R[z]ф(г)а2^ , которое вычитается при нахождении дисперсии.
=к, то есть тогда рассматриваемое слагаемое
выражения (4) будет иметь вид
П Г 0 0 0
После замены переменных и изменения порядка интегрирования последнее выражение может быть приведено к виду
Г |R[г]ф(г)а2 | R[^]ф(Ь)Л +
z-T T —z
-J R[z ]p( z )dz J tR[t ]p(t )dt
откуда видно, что при T^<ю основной вклад в диспер
сию будет вносить слагаемое —— J R[г]ф(2)с£
Аналогичный результат получаем при равенстве индексов j и к, і и і, j и і.
/гг1гГп т
— J JR[u - u &]R[u - u &]p(u - u &)p(u - u &)dudu&.
Делая замену переменных и меняя порядок интегрирования, этот интеграл приводим к виду
_2т->2
T J|^1 — ~~jR2[z]ф2(z)dz, который при Т^да
стремится к -1— [ R2[z]ф2(z)dz.
Рассмотрев точно так же два других слагаемых выражения (4), после проведенного исследования можно утверждать, что дисперсия статистики Q
при Т^к> убывает как —.
Оценка спектра мощности в виде сплайна
первого порядка при случайном числе данных
в моментах измерений
К условиям предыдущей задачи добавим еще одно, а именно: в каждый момент времени производится случайное число измерений щ, где щ - случайные величины, распределенные по закону Пуассона с параметром Я. Аналогичная задача рассматривается в работе [20], однако в ней используется другой подход к получению сплайна, а именно подход, в котором сначала ищется оценка функции корреляции, а затем уже находится оценка спектра мощности, причем коэффициенты сплайна оцениваются все сразу.
Выясним, как в нашем случае наличие условия случайного числа измерений в каждый момент времени повлияет на оценки коэффициентов сплайна.
Для решения задачи разбиваем всю ось частот аналогично тому, как мы делали это ранее и рассмотрим статистику вида
Q = Я— X - 4). (5)
Я ПТ і,,; іФгде ф(т) задаются выражением (2). Усреднив (5) по величинам щ, получим выражение
Q = -^ X УіУрі*- - 4 )>
пТ ,• “Й, & &
которое в дальнейшем и будет использоваться для нахождения оценок коэффициентов сплайна. Все остальные выкладки будут аналогичны рассмотренным ранее в разделе «Построение оценки спектра мощности в виде сплайна первого порядка».
Имитационное моделирование полученных оценок
Все расчеты и построение графиков выполнены в системе Mathcad 2000.
В соответствии с [21] было проведено моделирование реализации стационарного случайного процесса, которое описывается разностным уравнением первого порядка
х(п) = 0,8х(п -1) + у(п), (6)
где у(п) - случайные величины, равномерно распределенные на [-1; 1]. В качестве первого элемента генерируемой реализации х(0) было взято значение 1,5. Результат моделирования реализации стационарного случайного процесса представлен на рис. 2.
Рис. 2. Модельная реализация стационарного случайного процесса
Fig. 2. Model implementation of a stationary random process
На рис. 3 построены оценки спектра мощности коррелограммным методом и в виде сплайна первого порядка.
Оценим значимость полученных оценок спектра мощности. Для изучения качества полученной модели воспользуемся критерием Фишера-Снедекора. Для этого найдем наблюдаемое значение критерия
F = Е (S1І—S )2
І=1_________
Е (Si—S іі )2
здесь m - число оцениваемых паn — m
\\5(<B)
\\ 5,(ш) * / » / S,(0) /
Рис. 3. Оценки: $(ю) спектра мощности S(m) коррелограммным методом, S(rn) - в виде сплайна первого порядка для реализации случайного процесса рис. 2
Fig. 3. Estimations: S(rn) of power spectrum S(m), S((o) as a first-order spline for implementing the random process of the Fig. 2
раметров, _п_- объем выборки, 5 - выборочное среднее, Si, 1=1, п - значения спектра для модельной реализации случайного процесса (6), который имеет вид:
Б = ■
и Su, i=1,n - оцененные значе0,41 - 0,4 ео81
ния спектра мощности в виде сплайна. Наблюдаемое значение критерия оказалось равным 23, 684, что значительно превосходит табличное значение Ра,т-1,п-т=3,6 критерия Фишера-Снедекора при уровне значимости а=0,05 и числом степеней свободы кг=т-1 и к2=п-т. Так как ¥>¥ФМ, то полученная оценка спектра мощности в виде сплайна первого порядка является значимой. Индекс корреляции
Е (Б- - Бп )2
Е (Б - Б)2
полученная оценка спектра в виде сплайна первого порядка весьма тесно связана с аналитическим выражением спектра случайного процесса (6).
Таким образом, сплайновая оценка качественно правильно отражает вид спектра мощности.
Выводы
Выполнен расчет коэффициентов сплайна первого порядка в приложении к оценке спектра мощности стационарного случайного процесса в двух ситуациях, когда в каждый момент времени производится:
Показано, что полученные оценки являются несмещенными. Установлено, что дисперсия оценок асимптотически ведет себя как -1, где Т - время наблюдения.
СПИСОК ЛИТЕРАТУРЫ
Поступила 24.03.2014 г.
UDC 519.2
EVALUATION OF THE POWER SPECTRUM OF A STATIONARY RANDOM PROCESS AS A FIRST-ORDER SPLINE
Irina G. Ustinova,
Cand. Sc., National Research Tomsk Polytechnic University, 30, Lenin Avenue,
Tomsk, 634050, Russia. E-mail: igu@sibmail.com
Elena G. Lazareva,
Cand. Sc., National Research Tomsk State University, 36, Lenin Avenue, Tomsk, 634050, Russia. E-mail: lazareva-e-g73@mail.ru
Elena I. Podberezina,
Cand. Sc., National Research Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia. E-mail: pei@tpu.ru
The urgency of the work is caused by the fact that the power spectrum as well as the correlation function, is one of the most important characteristics of the second order of stochastic process. The spectrum shows what kind of harmonics prevails in the process, its structure; allows estimating spectral composition of the studied useful signals and noises. It is possible to synthesize (recover) signal, as well as to construct linear, including optimal filters, to obtain error estimates of linear filtering by the spectra.
The main aim of the study is to evaluate the power spectrum of a stationary random process as a first-order spline under the following measurement schemes: one measurement is made every moment, a random number of measurements is made every moment. Study of statistical characteristics of the estimates.
The methods used in the study: the methods of probability theory and mathematical statistics are used for calculation.
The results: The authors obtained an unbiased estimate of the power spectrum as a first-order spline in two schemes of measurement: one measurement is made every moment, a random number of measurements is made every moment. It is shown that the variance of
estimates behaves asymptotically like T, where Tis the time of observation.
Power spectrum, correlation function, first-order spline, parameter estimations, statistical properties of estimations.
REFERENCES
Received: 24 March 2014.