УДК 517.518.85:517.443
Г.С. Ханян
ТЕОРЕМА ОТСЧЕТОВ ДЛЯ СИГНАЛА КОНЕЧНОЙ ДЛИТЕЛЬНОСТИ С НЕОБЯЗАТЕЛЬНО НУЛЕВЫМ ИНДЕКСОМ ЧАСТОТНОЙ ПОЛОСЫ
Рассматривается теорема отсчетов во временной области применительно к сигналам ограниченной длительности в полосе частот с верхней и нижней (необязательно нулевой) границами среза. Ставится цель найти более широкие достаточные условия справедливости теоремы отсчетов для финитных процессов и, тем самым, расширить сферу приложений теоремы в цифровой обработке сигналов. Проводится доказательство теоремы для базового сигнала - фрагмента гармонического колебания, и на основе принципа суперпозиции делается обобщение на полигармонический процесс с частотами гармоник, определяемыми в процессе исследования.
Теорема отсчетов; частота дискретизации; индекс полосы частот; длина цифровой реализации; гармонический сигнал; интерполяционная формула; дельта-функция.
G.S. Khanyan
SAMPLING THEOREM FOR FINITE DURATION SIGNAL WITH A NON-OBLIGATORY ZERO INDEX OF THE FREQUENCY BAND
The sampling theorem is considered in time domain with respect to signals of limited duration in a frequency band with upper and lower (non-obligatory zero) cutoff boundaries. The aim is to find wider sufficient conditions of sampling theorem validity for finite processes and thereby to broaden application sphere of this theorem in digital signal processing. A proof of the theorem is conducted for the base signal - a time fragment of a harmonic oscillation, - and on the superposition principle basis a generalization is made to polyharmonic process with harmonic frequencies determined in the course of the study.
Sampling theorem; sampling frequency; index of frequency band; length of digital realization; harmonic signal; interpolation formula; delta-function.
Введение. Одним из фундаментальных результатов, установленных в теории
информации, является теорема отсчетов Уиттекера-Котельникова-Шеннона,
в которой утверждается о возможности и способе полного восстановления аналогового сигнала s (t) по его цифровой реализации s (tn), полученной с частотой дискретизации F, если спектр сигнала ограничен полосой частот 0 < f < F/2.
В данной, классической формулировке [1] теоремы сигнал восстанавливается преобразованием
/ N sin nF(t - tn ) ,1Л
s(t) = Y s(tn )---- ' , ; ; tn = n/ F . (1)
nt~L nF(t - tn )
Обобщение теоремы на случай полосовой фильтрации выглядит так:
s(t) = Y s(t ) Sin П(G + 1) F (t - tn ) - Sin ПGF (t - tn ) , (2)
n=~ n nF (t - tn )
и рассматривается в [2], где для индекса полосы частот G указаны границы
Fbw = ([G +1/2] + | G - [G +1/2] |) F/2, FUp = ([G +1/2] +1- | G - [G +1/2] |) F/2(3)
диапазона частот протекания процесса
Flow < f < Fup , (4)
в котором преобразование (2) тождественное:
s'(t) = s(t) . (5)
Формулировка (2) теоремы с целочисленным индексом G рассматривалась, например, в [3]. В этом случае из (3) следует, что границами диапазона являются Flow = GF/2, Fup = (G+1)F/2, и индекс G = 0, 1, 2, ... нумерует в частотной области последовательно расположенные полосы (4), в которых выполняется утверждение теоремы (5). Классической формулировке (1) соответствует при этом индекс G = 0 с присоединением точки f = 0 - частоты постоянной составляющей процесса.
Постановка задачи. Теорема отсчетов приобретает практическую ценность, когда преобразование (2) имеет лишь конечное число слагаемых. Тогда оно может, в математическом обеспечении той или иной измерительно-вычислительной системы, служить инструментом интерполяции ограниченного числа данных, получаемых при проведении различного рода испытаний и тестов - электронных устройств, приборов, аппаратуры, деталей и узлов машин, метеонаблюдений и т.д.
Ясно, что при таком ограничении временной области в ядре преобразования (2) наряду с параметрами G и F должен фигурировать третий параметр - число N интерполируемых отсчетов s (tn) сигнала s (t), называемое также длиной цифровой реализации сигнала. При нулевом индексе G вид такого преобразования известен:
/(,) = у s(t) sin * F (t - tn) . (6)
tí n N sin *F (t -1,)/N
Так, в работах [4], [5] утверждение (5) в формулировке (6) теоремы доказывается для периодической функции s (t) = s (t+T), но только с нечетным числом отсчетов N = FT, укладывающихся на протяжении периода T.
В настоящей работе ставится цель найти более широкие достаточные условия справедливости теоремы отсчетов для финитных процессов - с необязательно нулевым индексом полосы частот. Тем самым будет расширена сфера приложений теоремы в цифровой обработке сигналов, например, в задаче интерполяции поля давления при малом (несильно возмущающем поток) и технически жестко обусловленном числе датчиков на входе в авиационный двигатель.
Метод исследования. При переходе от (1) к (6) ядро преобразования - функция sinc(x)=(sinx)/x заменяется ее «неполным циклическим» аналогом
. . sin х 1 2к - N +1 ,п,
SinCN (х) =--------=---- у cos---------х, (7)
N N sin х / N N k=í N
получаемым суммированием N косинусоид и стремящимся к sinc(x) при N^^.
Руководствуясь этой идеей, модификацию преобразования (2) для сигнала ограниченной длительности с произвольным индексом G представим в виде
s'(t) = !>n)sin*(G(t-^>; tn = tí ♦ n/F. (8)
,=o N sin *F (t - tn)/ N
Формулировка (8) теоремы отличается от (1), (2) и (6) еще и тем, что дискретизация в общем случае проводится со смещением: начальный момент времени t0 необязательно равен нулю или кратен интервалу дискретизации 1/F. На практике такая ситуация может иметь место, например, при аналого-цифровом преобразовании нескольких сигналов с применением мультиплексирования.
В качестве базового объекта преобразования примем гармоническое колебание с амплитудой а, частотойf начальной фазой ф и длительностью T:
s(t) = a cos (2п ft + ф); а > 0, f > 0, - п < ф < п; 10 < t < 10 + T . (9)
Как только будут установлены параметры сигнала (9) и преобразования (8), отвечающие утверждению (5) теоремы, более сложный тип сигналов можно будет порождать суперпозицией тригонометрических функций с такими параметрами.
Доказательство теоремы. Дискретизация сигнала (9) с частотой ¥ и подстановка цифровой реализации 5 (П в правую часть преобразования (8) дает, с использованием формулы выражения разности синусов через удвоенное произведение синуса полуразности на косинус полусуммы аргументов:
N 1 £ _ £ £ _ £
5'(£) = п£ со8(2п/£п + ф)8шеN/2 п¥------- со8п(2О +1)¥---------- • (10)
п=0 2 2
Вычислить сумму (10) можно, воспользовавшись разложением функции (7) с последующим представлением косинусов в экспоненциальной форме и применением формулы суммирования геометрической прогрессии:
1 _ N 1 а а ф р'2 и/ N
1 _ а ’ а
4 (11)
N дкш, а = р2 пиК / N
I чп =
и =±1; к = 0,±1,±2,...; I = 0,±1,±2,...,
где 8 - символ Кронекера, и - знак, К - заданное, I - неопределенное целое число в дискретной дельта-функции. Но тогда мы вынуждены сделать предположение, что при ненулевом индексе О число отсчетов N - четное (число слагаемых N/2 -целое). В случае О = 0 число отсчетов N может быть любым натуральным числом, поскольку оно, в соответствии с (6), совпадет с числом слагаемых разложения (7).
Подставив правую часть (7) в (10) с заменой в этой формуле N на N/2, имеем:
_ N _1 N /2_1
5'(£) = _П_ £ £ р’(2 п/гп+ф) £ £ ри(2к_N/2+1)ж¥(£_£„)/N £ р'£(О+1/2)^(£_£„) (12)
4N п=0 ’=±1 к =0 и =±1 ¿=±1
Поменяв местами индексы суммирования п и к, расположив их правее индексов суммирования’, и, £ слагаемых формулы Эйлера и введя обозначение, записываем
К = и/Г _ иЬМО/2 _ ULN/4 + N/4 _ 1/2, (13)
записываем преобразование (12) в виде, в котором слева (в комплексной форме) фигурирует исходный сигнал 5(0 вместо его цифровой реализации 5(£п):
П N /2_1 ( 1 N _1 Л
^'(£) = а £ рЧ(2п/£+ф) ££ £ р'2и(г_г°)(к_К)/N _ £ р'2пи(К_к)п/N (14)
4 у =±1 и =±1£=±1 к =0 V N п=0 у
Ясно теперь, что дальнейшее упрощение возможно, если крайняя справа сумма подпадет под второй случай формулы (11), т.е. будет дискретной дельтафункцией. Но тогда К- к должно быть целым числом, или, что то же самое, целым должен быть параметр К. Приняв это предположение, представляем (14) в виде
П N /2_1
5 V)=П £ р’(2п/н" £ £ Ё р'2 и “ _'0),к _К ”" гК-.,1». (15)
4 ’=±1 и=±1£=±1 к=0
В силу определяющего свойства дельта-функции в сумме по к остается всего одно слагаемое, и преобразование сводится лишь к сумме по индексам ’, и, ¿:
51 ) = П £££ р'‘ 2’^ 2'1Г“^«“=К-IN£N; I = 0, ± 1, ± 2....С16)
4 ’=±1 и =±1£=±1
Издержкой столь радикального упрощения является то, что результат преобразования 5 (£) становится многозначной функцией - с ветвями, нумеруемыми отмеченным в имени функции индексом I. Здесь и далее через АВ обозначена логическая дельта-функция от булевой переменной В, равная единице при истинности этой переменной и нулю при ее ложности (символ Кронекера представляет собой частный случай этой функции: 8 к , п = А к = п). В формуле (16) в роли В выступает неравенство, указывающее множество значений индекса к в (15).
Подставив в (16) К из (13), получаем, наконец, результат преобразования:
5 7 (£ ^ N >2,0>0 = "Т £ А (О№ +1)/2<| /Г _IUW|<(ОW+N _1)/2 со§(2п( / _ )£ + (ф + 2 ЖиРг0)) , (17)
2 и =±1
означающий, что исходное колебание расщепляется на две компоненты с частотами/ ± № и фазами ф + 2пШ г0.
Проанализируем теперь условия целочисленности параметра К. Представим безразмерную частоту сигнала как сумму целого числа т и правильной дроби ц:
/Г = т + ц; т = 0,1,2,...; 0 < ц < 1. (18)
Подставим (18) в (13) и напишем решение полученного уравнения:
ц = {_т + ’иК + ’¿N0/2 + ’Ш/4 _ Щ/4 + ’и/2}. (19)
Представим N0 в виде суммы целой (четной или нечетной) и дробной частей: N0 = 2Я + г + р; Я = [ N0/2], г = [NG]mod2, р = {N0}. (20)
Подстановка (20) в (19) и отбрасывание заведомо целочисленных термов из фигурных скобок дает, путем обобщения всех частных случаев комбинаций ’, и, ¿:
= 1 _ г + ’¿р +(1 _ ’¿)г sgn р (21)
2 '
Дробная часть ц безразмерной частоты /Г, как заданный параметр сигнала, не должен зависеть от комбинаций индексов суммирования ’ и £ в (21), что, очевидно, возможно только при р=0, т.е. когда N0 — целое число. Таким образом, ц - целое или полуцелое число в зависимости от нечетности или четности N0:
ц = (1 _ (N0 )mod 2)/2 = {(N0 _ 1)/2}; ^/2} = 0; {N0} = 0; 0 > 0. (22)
Это и есть условия целочисленности параметра К - достаточные условия представления результата преобразования сигнала в виде (17).
Рассмотрим теперь случай 0 = 0, руководствуясь той же методикой анализа.
В этом случае вместо (10) мы имеем дело с преобразованием а N _1 N _1 2 к _ N + 1
5 (£) = ^ ££-(27/гп + ф) сов------------- ----------------п ¥ (г _ ¡п), (23)
N п=0 к =0 N
которое в комплексной форме выглядит следующим образом:
г, N _1 (л N _1 А N — Л
5'(г) = П£р'1(2п/г+ф) ££р'2пир(г_г0)(к_К)/мI _£р'2пи(К_к)п/N I. к = и/Т + ___1. (24)
4 ’=±1 и =±1к=0 VN п=0 у 2
Требование целочисленности К приводит вместо (19)-(22) к соотношениям ц = {_т + ’иК _ и/2 + ’и/2} = {’и/2 _ и/2} = {^ _1)/2}; 0 = 0, (25) а преобразование - к многозначной функции, !-я ветвь которой задается формулой
SI (г)| 0=0 = ^ £ A|/Г_UN|<(N_1)/2 СОв(2п(/ _ иР) г + (ф + 2и0». (26)
2 и=±1
Параметр ц и на этот раз оказывается целым или полуцелым числом, что позволяет обобщить условия (22) и (25) на все случаи 0>0:
{N0} = 0, ^/2}0 = 0, ц = {(N0 + N _ 1)/2}. (27)
Далее, поскольку ц фиксировано, то смежные разрешенные безразмерные частоты отличаются на единицу, что позволяет обобщить формулы (17) и (26):
(г) = П £ А,и со8(2^(/ _ и)г + (ф+2яШг0));
2 и =±1 . ( )
О ,и = А(0№ _1)/2<| /Г _Ж |<(0№+ N _1)/2
Ясно, что за восстановление (5) сигнала может отвечать только нулевая ветвь функции (28) при О 0,±1=1, поскольку частоты и фазы компонент и=±1 при I=0 совпадают, и сумма компонент становится равной исходному сигналу (9). Однако
этого недостаточно, необходимо еще, чтобы Djj =0 для всех Ы0 и J = ±1, иначе при любом способе суперпозиции ветвей произойдет пересечение области определения функции Djj с областью определения функции D 0jJ, т.е. будет иметь место явление искажения сигнала из-за наложения частот.
Анализ наложения ветвей путем решения уравнения D0J =DJj = 1 (фактически, системы неравенств) показал, что на нулевую ветвь при нецелом G накладываются 2 из 4 компонент ветвей J=±[-G], что определяет для теоремы (8) диапазон частот
— 11 < f < F —— (29)
2 _|J up 2T
с параметрами F¡ow и Fup из (3). Границы в (29) отличаются от границ в (3) для теоремы (2) не более чем на половину спектрального разрешения 1/T каждая. В пределе T ^ ^ неравенство (29) приобретает вид (3), соответствующий сигналу бесконечной протяженности.
В силу принципа суперпозиции (линейности преобразования (8) относительно отсчетов сигнала) теорема верна для полигармонических процессов Q
s(t) = Z am COS(2n fmt + Фт ^ fm = (m + ^)/T , am ^ 0, “ П < Ф m < П (30)
m=P
с безразмерными частотами гармоник fmT - целыми или полуцелыми в зависимости от параметра ц из (27).
Целочисленные границы суммирования P и Q в терминах частот среза (Р+ц)/Т и (Q+^)/T совпадают с границами неравенства (29) и задаются формулами: P = [(NG - N +1)/2] + N [{G} +1/2](1 - {G}) + [N/2]
Q = [(NG + N -1)/2] + N [{G} - 1/2]{G}.
Число гармоник в процессе (30) или же число разрешенных частот гармонического сигнала (9) составляет
M = Q - P +1 = N | {G} —1/21 +{N/2}. (32)
Оно максимально (M=N-[N/2]) при целом индексе G и минимально (M=0) при по-луцелом. Последнее означает невыполнение утверждения (5) теоремы (8), так же,
как и (2). В простейшем случае (N=1) из (27), (30)-(32) следует, что можно по единственному отсчету восстановить лишь функцию s(t)=a0cosq 0, т.е. константу.
Заключение. Сформулирована и доказана при определенных достаточных условиях теорема отсчетов для финитных полигармонических процессов с выводом аналитических формул для границ дискретного спектра частот гармоник в зависимости от индекса частотной полосы и числа интерполируемых отсчетов.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Котельников В.А. О пропускной способности «эфира» и проволоки в электросвязи // Материалы к I Всесоюзному съезду по вопросам технической реконструкции дела связи и развития слаботочной промышленности. - М.: Всесоюзн. энергетич. комитет, 1933.
2. Ханян Г.С. Обобщение теоремы отсчетов на случай нецелого индекса полосы частот // Международная научно-техническая конференция к 100-летию со дня рождения В.А. Котельникова: Москва, 21-23 октября 2008 г. Тезисы докладов. - М.: Изд. дом МЭИ, 2008. - C. 35-37.
3. Петров В.В., Усков А.С. Информационная теория синтеза оптимальных систем контроля и управления (Непрерывные системы). - M.: Энергия, 1975. - 232 с.
4. Linden D.A. A Discussion of Sampling Theorems. Proc. of the IRE, July 1959. - P. 1219-1226.
5. Brown J.L., Jr. An RKHS Analysis of Sampling Theorems for Harmonic-Limited Signals // IEEE Trans. of Acoustics, Speech, and Signal Processing. - April 1985. - Vol. ASPP-33, № 2. - Р. 437-440.
Статью рекомендовал к опубликованию д.т.н., профессор С.Г. Крутчинский.
Ханян Гамлет Сократович - Центральный институт авиационного моторостроения им. П.И. Баранова; e-mail: khanian@mail.ru; 111116, г. Москва, ул. Авиамоторная, 2; старший научный сотрудник; к.т.н.
Khanyan Gamlet Sokratovich - Central Institute of Aviation Motors; e-mail: khanian@mail.ru;
2, Aviamotornaya street, Moscow, 111116, Russia; principal researcher; cand. of eng. sc.
УДК 621.372.54
В.П. Тепин, А.В. Тепин
СИНТЕЗ ЭКВИВАЛЕНТНЫХ ПЕРЕСТРАИВАЕМЫХ АНАЛОГОВЫХ И ЦИФРОВЫХ ФИЛЬТРОВ
Рассматриваются произвольный аналоговый фильтр с управляемой частотой настройки при неизменной форме частотной характеристики, а также эквивалентный ему рекурсивный цифровой фильтр, управляемый двумя способами: прямым — изменением частоты дискретизации, а также косвенным - изменением весовых коэффициентов. Предлагаются методы синтеза таких фильтров, основанные на частотном масштабировании и билинейном z-преобразовании аргумента передаточной функции базового аналогового или цифрового фильтра. Эти методы гарантируют эквивалентность не только частотных характеристик, но также и законов перестройки фильтров рассматриваемого класса.
Перестраиваемый фильтр; синтез передаточной функции; прямой и косвенный способы перестройки; частотное масштабирование; закон перестройки.
V.P. Tepin, A.V. Tepin SYNTHESIS OF EQUIVALENT TUNABLE ANALOG AND DIGITAL FILTERS
In this paper, we consider a tunable analog filter and equivalent recursive digital filter. Their adjustment frequency is controlled provided that frequency response remains invariable. Two modes of digital filter tuning are discussed: direct one — by varying a sampling frequency, and indirect one - by varying the filter weighting coefficients. We propose several methods for either filter synthesis which use frequency scaling and bilinear z-transform of basic analog or digital filter transfer function’s argument. These methods ensure both the frequency responses and control laws equivalence for the filters considered.
Tunable analog filter; tunable digital filter; transfer function synthesis; direct and indirect control modes; frequency scaling; control law.
Введение. Перестраиваемые аналоговые и цифровые фильтры широко используются в современных системах автоматического управления, а также в информационно-измерительных и телекоммуникационных системах, обеспечивая адаптивную селекцию и обработку, формирование, поиск и автоматическое сопровождение сигналов, слежение за их частотой, подавление нестационарных помех, а также решение ряда других задач. Среди них наиболее распространены фильтры, относящиеся к классу «частотно-масштабируемых», управление частотой настройки которых достигается без искажения формы частотных характеристик, путем изменения масштаба аргумента передаточной функции [1]. Известно, что при этом масштабируются не только частотные, но и временные характеристики [2]. Такие фильтры применяются также для реализации управляемых генераторов периодических и непериодических колебаний, в частности - параметрических генераторов динамического хаоса, а также систем защищенной передачи информации на их основе [3, 4].