CC BY
44
ЭЛЕКТРОННЫЕ И ЭЛЕКТРОМАГНИТНЫЕ
УСТРОЙСТВА
УДК 621.396:681.323
С. И. Зиатдинов
ИНТЕРПОЛИРОВАНИЕ ФУНКЦИИ НА ОСНОВЕ ДИСКРЕТНОГО ПРЕОБРАЗОВАНИЯ ФУРЬЕ
Исследуются линейные искажения при интерполяции случайного процесса на основе отсчетов его спектральной плотности, полученных дискретным преобразованием Фурье. Анализируются ошибки интерполяции для различных спектрально-корреляционных характеристик случайного процесса. Показано, что в точках взятия отсчетов ошибки интерполяции равны нулю и принимают максимальные значения в середине периода дискретизации.
Ключевые слова: спектр, дискретизация, восстановление, ошибки.
При цифровой обработке информации непрерывная функция х(^) представляется последовательностью ее отсчетов х[п], взятых через период дискретизации Т, при этом п=0, 1, 2, ...
На практике для получения спектральной плотности исследуемой функции х(^) широко используется дискретное преобразование Фурье (ДПФ), позволяющее по пачке из N отсчетов функции х(^) получить N отсчетов спектральной плотности [1]:
N -1
5[к]=£ х[п]е~ 1Шпк , (1)
п=0
где П = 2п/N1, к = 0...(-1).
В то же время существует обратное ДПФ, которое по полученным отсчетам спектральной плотности (1) однозначно определяет исходную импульсную последовательность [1]:
N -1
Nk=»
Рассмотрим задачу восстановления исходной функции х(^) в любой точке временного интервала ^ = 0...( - 1)Т на основе отсчетов спектральной плотности (1).
Очевидно, что в точках I = пТ интерполирующая функция совпадает с исходной функцией х(^), а в точках I ^ пТ возникают ошибки интерполяции, оценка которых и составляет цель настоящей статьи.
Для произвольного момента времени I в пределах временного интервала ^ = 0...( - 1)Т соотношение (2) в общем виде становится комплексным и записывается следующим образом:
1 N-1
/ (*) = N I 5[к]е^ . (3)
к=0
х [п] = 1 I 5 [к]е^Тпк . (2)
44 С. И. Зиатдинов
Положим ^ = £Т + АТ, где £ = 0... (N - 2) — номер временного отсчета в пределах интервала ^ = 0...( - 1)Т; АТ = 0...Т — точка интерполирования в пределах периода дискретизации. Тогда выражение (3) принимает вид
i N-1 ,
/ (') = 1X »[k j('T+AT>■ (4)
N k=0
После подстановки соотношения (1) в формулу (4) получим
AT" T
1 N -1N-1 jQTk \l-n+
y (0 = N XX*[n]e v J . (5)
N k=0 n=0
В дальнейшем для практических случаев рассмотрим только вещественную составляющую выражения (5):
1 N-1N-1 ( AT Л
y (t) = 77 H *[n] cos QTk U - n + ~t 1. (6)
N k=0 n=0 v 1 J
Поскольку интерполирующая функция (6) при AT ^ 0 отличается от исходной функции x(t), то возникающие ошибки интерполяции оценим коэффициентом линейных искажений [2-6]
Kл.и W = V 1 - R12 (т) ,
где R12 (т) — коэффициент взаимной корреляции интерполирующей функции y(t) и исходной функции x(t).
Определим коэффициент линейных искажений для более простого случая при т = 0 . Пусть исходная функция x(t) представляет стационарный случайный процесс с нулевым математическим ожиданием. Тогда коэффициент взаимной корреляции R12 (0) может быть найден из следующего выражения:
(0) = = ^ I I В[( - -)Т + АТ] еОГ/ - п + Т-),
а у а х N а у а хк=0 Г Т )
где а у, ах — среднеквадратические значения функций у(() и х(^); В(т) — корреляционная
функция исходного процесса х(^).
Для нахождения коэффициента взаимной корреляции (0) необходимо знать средне-
квадратическое значение а у интерполирующей функции, которое определяется соотношением
11/2
а v = —
y N
N-1N-1N-1N-1 ( AT Л ( AT 1
X X X X В [n - m]cos QTk I £ - n +-J cos QTp U - m +-
k=0 n=0 p=0 m=0 v T J v T J
Пусть корреляционная функция процесса х(^) описывается соотношением
В (т) = а2 ехр (-|т|А/),
где А/ — параметр, определяющий ширину спектральной плотности случайного процесса х(^); данной корреляционной функции соответствует пологая медленно спадающая спектральная плотность.
Интерполирование функции на основе дискретного преобразования Фурье
45
Результаты расчетов коэффициента линейных искажений Кли (0) для различных значений параметров £ и АТ / Т при значении произведения А/Т = 0,005 и числе отсчетов функции х(^) N = 17 представлены в табл. 1.
_Таблица 1
AT / T 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
Kл.и ,% i = 0 0 5,99 11,47 15,32 16,59 15,24 12,07 8,42 5,49 3,28
L = 8 0 3,74 6,23 7,80 8,46 8,49 8,14 7,38 5,94 3,65
5 N 0 3,28 5,49 8,42 12,07 15,24 16,59 15,32 11,47 5,99
Как следует из анализа полученных данных, в точках взятых отсчетов функции х(^) коэффициент линейных искажений равен нулю. При этом минимальные значения коэффициента Кли (0) имеют место в середине интервала времени (N -1) Т (£ = 8) и в точках, примыкающих к моментам взятия отсчетов функции х(^).
В табл. 2 приведены результаты расчетов коэффициента линейных искажений Кли (0) для корреляционной функции процесса х(^) вида
В (т) = а2 exp (-т2А/2),
соответствующей резко падающей спектральной плотности. Вычисления произведены при прежних исходных данных.
_Таблица 2
AT / T 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
Kл.и ,% i = 0 0 1,50 3,00 4,02 4,28 3,75 2,65 1,37 0,32 0,18
L = 8 0 0,56 0,93 0,99 0,73 0,26 0,28 0,65 0,74 0,50
5 II 0 0,18 0,32 1,37 2,65 3,75 4,28 4,02 3,00 1,50
Сопоставляя результаты, представленные в табл. 1 и 2, можно отметить, что для резко падающей спектральной плотности функции x(t) ошибки интерполяции в 4—7 раз меньше, чем в случае медленно падающей спектральной плотности. В целом же характер поведения коэффициента линейных искажений Кли (0) для различных значений параметров £ и AT / T остается прежним.
список литературы
1. ГолдБ., Рейдер Ч. Цифровая обработка сигналов. М.: Сов. радио, 1973.
2. Зиатдинов С. И., Жуков А. Д. Искажение сигнала в узкополосных фильтрах // Изв. вузов. Приборостроение. 2006. Т. 49, № 12. С. 44—47.
3. Зиатдинов С. И. Линейные искажения сигнала фильтром Баттерворта // Там же. 2007. Т. 50, № 1. С. 35—39.
4. Зиатдинов С. И. Линейные искажения сигнала экстраполяторами // Там же. 2007. Т. 50, № 5. С. 57—60.
5. Зиатдинов С. И. Линейные искажения сигнала интерполятором // Там же. 2007. Т. 50, № 10. С. 50—53.
6. Зиатдинов С. И., Гирина Н. В. Анализ ошибок при тригонометрической интерполяции // Там же. 2008. Т. 51, № 5. С. 42—45.
Сведения об авторе
Сергей Ильич Зиатдинов — д-р техн. наук, профессор; Санкт-Петербургский государственный
университет аэрокосмического приборостроения, кафедра информационно-сетевых технологий; E-mail: [email protected]
Рекомендована кафедрой Поступила в редакцию
информационно-сетевых технологий 29.11.07 г.
