
книги из ГПНТБ / Френкс, Л. Теория сигналов
.pdfФункции Эрмита
Для Т — (—оо, оо) и w (f) = е~р полиномы
Ф„ (t) = (2« н1 У я ) - 1/2 Нп(0, п = 0, |
1,2 ,... |
(3.42) |
образуют ортонормальную систему. Нп — это |
полиномы |
Эрмита, |
которые задаются в виде |
|
|
Я Л 0 Н - 1 ) " е ^ ( е - Я |
|
(3.43) |
или рекуррентной формулой |
|
|
Нп (t) = 2 Ш п . 1 (/) — 2 ( л - 1) Я„_а (0- |
(3.44) |
|
Функции Эрмита |
|
|
яМ*)=(2"л! / Й Г 1/2 e~t!l2Hn(t) |
(3.45) |
ортонормальны с единичным весом на (— оо, оо); они могут также быть получены применением процедуры Грама — Шмидта к последователь ности {lne~i!!/2; п = 0, 1, 2, ...}.
Функции Уолша
Для Т — [0, 1] и w (0 = 1 можно построить полную ортонормальную систему функции типа «прямоугольных волн». Для определения
таких функций удобно использовать два индекса. Функции <p„A) (О, называемые функциями Уолша [9—11], определяются следующим об разом:
Ф0(0 = 1; 0 < * < 1 ;
1; 0 < / <
Фх ( 0 :
- 1 :
1; |
< / < 1, |
фУ> (0 =
, + < , < *
ф < « ф ,
Фi2)(t) =
(24—1) |
ф™’ (20; о < * < - ± - |
|
(3.46) |
||
Фт+1 ОФ |
||
|
(-1 )* + 1ф!») (2 /-1 ); |
|
|
6i |
где т = 1, 2, 3, ... и k = \ , 2, 3, .... 2т~ 1.
Эта система весьма важна для практики, поскольку функции ку сочно-постоянны и принимают только два значения (+1 и —1). Подоб ные сигналы легко могут быть получены с помощью двоичных логи-
¥ 2 ¥
Рис. 3.5. Функции Уолша, пронумерованные в соответствии с количеством перемен знака на интервале 0 < ^< 1 .
ческих схем. Перемножение таких сигналов для получения, скажем, скалярного произведения производится очень просто, с помощью ключа, изменяющего полярность и включаемого в нужные моменты
времени. Заметим, что функции ц>т(т~ !)) {t) соответствуют обычным прямоугольным волнам. Мы можем перейти к соответствующей одноиндексной системе, упорядочив функции (3.46) так, чтобы п-я функ ция п раз пересекала нулевой уровень на интервале 0 < t < 1 (т. е. п раз меняла знак). Это достигается при следующих обозначениях:
®о (0 = Фо |
(0. %(0 = Фх (0. |
(0 = фт( (0, |
(3.47) |
причем п = 2т~ 1 + |
k — 1 (рис. 3.5). |
|
|
3.4. ОПЕРАТОР РАЗЛОЖЕНИЯ СИГНАЛА В АППАРАТУРНОЙ РЕАЛИЗАЦИИ
Чтобы завершить процесс отыскания приближенного численного представления для произвольного сигнала конечной энергии, выберем некоторую систему базисных функций {<рг} (и взаимный базис {0*}), и установим эквивалентность
|
П |
|
х ~ |
2 а гф;> |
(3.48а) |
где |
г = 1 |
|
|
|
|
а,- - hi (х) = |
(х>0j); *' = 1.2,..., 11. |
(3.486) |
Компонента!
Рис. 3.6. Обобщенный анализатор спектра.
По аналогии с обычным преобразованием Фурье будем рассматри вать (3.48) как пару преобразований. В (3.48 а) сигнал представлен определенной линейной комбинацией {<р;}, а формула (3.48 б) дает линейное правило для отыскания разложения сигнала на выбранные компоненты. Последнюю операцию можно рассматривать как обрат ную по отношению к разложению, и именно ее мы хотим реализовать с помощью соответствующих физических устройств. Используя метод реализации линейных функционалов, показанный на рис. 2.7, получаем устройство, которое является обобщением лабораторного прибора — анализатора спектра (рис. 3.6). Эту схему можно выполнить иначе, заменив ячейки «перемножитель — интегратор» на ячейки «квантова тель — фильтр», что соответствует рис. 2.8.
Во многих случаях, генерируя в приборе «подходящие» колебания, не удается точно сформировать нужную систему функций {0г}. На пример, в некоторых системах обработки оптических сигналов опера-
63
ция умножения выполняется путем освещения через материал с раз ной степенью прозрачности.
В таких условиях желательно, чтобы функции, на которые умно жается сигнал, были неотрицательными, в частности может использо ваться система вещественных экспоненциальных функций. В качестве другого примера приведем случай, когда «эталонные» функции должны быть порогового типа («да» — «нет»), чтобы умножение могло произ водиться с помощью соответствующего ключа. Может случиться, что у проектировщика системы обработки нет никакого выбора и нужно использовать заданные функции. Это приводит к следующей задаче.
Имеется |
набор устройств, реализующих |
линейные функционалы |
{fUl, fUt, |
•••> fuml Как наилучшим образом |
использовать эти устрой |
ства, чтобы получить произвольный функционал /д? Для ответа на этот вопрос воспользуемся обсуждавшимся в § 2.6 соответствием между сигналами и линейными функционалами.
Задача состоит в том, чтобы найти наилучшее (в смысле минимума расстояния в сопряженном пространстве) приближение искомого функ
ционала путем ортогонального проектирования /е £ |
[L2 (Т)]* на под |
||
пространство из [L2 |
(Т)]*, натянутое на |
fu„ fu„ ..., fum}. |
|
Приближение /е |
линейной комбинацией |
{/„.} |
реализуется си |
стемой параллельных цепей, рассмотренных в § 2.6; на выходе каждой из цепей имеется звено, регулирующее коэффициент усиления, как показано на рис. 3.7. Коэффициенты усиления устанавливаются так,
|
ТП |
|
|
чтобы минимизировать ||/0 — |
ah |
||; они определяются через ска- |
|
|
k=\ |
пространстве, согласно |
(2.61). |
лярные произведения в сопряженном |
|||
Применяя ортогональное проектирование, имеем |
|
||
|
ТП |
|
|
|
А=1 |
й |
(3.49а) |
где |
|
||
|
|
|
|
ah = |
(fe,frft) = |
(vft,0). |
(3.496) |
При этом {vft; k = 1, 2, ..., m) есть взаимный базис в подпространстве Nm, натянутом на {ufe; k = 1,2, ..., m}. Нетрудно показать, что такая аппроксимация для /0 в точности соответствует минимуму расстояния до 0, получаемому при проектировании на Nm■Учитывая соответствие
норм, имеем для любого сигнала х единичной энергии из L2 |
(Т): |
||
е - 2 atUk |
/е — 2 |
ahfu = e=>- |
|
k=i |
k=i |
|
|
Ы * )— 2 |
%/«fe(x) |
;С е. |
(3.50) |
*=i |
й |
|
|
Можно сказать также, что 100 г есть наибольшая ошибка отображения в процентах по отношению к |[ х ||. Отсюда прямо вытекает схема реали зации представления сигнала с помощью имеющихся конкретных устройств. Мы хотим представить сигнал х точкой в пространствеМ п, натянутом на {<рг; г = 1, 2, ..., п}\ можем же мы вычислять только
64
скалярные произведения (х, иг), где (и,; i = 1, 2, |
т} есть базис |
ВN т,’
Всоответствии с (3.49) устройство, в котором реализуется точная
верхняя грань ошибки определения компонент, должно давать
т |
|
|
|
|
« г= 2 |
а, |
й |
(х, и*); *' = 1 ,2 ,...,« , |
(3.51) |
£=1 |
|
|
|
где
aih= (Vft, 6;).
Рис. 3.7. Аппроксимация в сопряженном пространстве; реализация fe (х) = = (х, 0) с помощью заданной системы функционалов.
Анализатор формы |
сигнала должен |
содержать, таким образом, |
|||
звенья, реализующие заданные линейные |
функционалы, |
и блок взве |
|||
шенного суммирования, |
описываемый |
матрицей |
А |
с элементами |
|
{aik} из (3.51). Такая схема показана на рис. 3.8. |
|
|
|||
Рис. 3.9 условно иллюстрирует приближение желаемого пред |
|||||
ставления. Вначале находится ортогональная |
проекция х на N m, |
||||
представимая /п-мерным вектором р = {рь р2, |
..., |
рт }, |
а затем полу |
ченная в Nm точка ортогонально проектируется на М п. Получается представление х в М п, аппроксимирующее желаемое представление х,
которое получилось бы при прямом ортогональном проектировании х
на М п. Приближенное представление от = |
{а1( а 2....... ап} получается |
||
при умножении матрицы А на вектор р — представление |
х в Nm, по |
||
отношению к базису (v;; / |
== 1, 2, ..., т}: |
|
|
«1 |
а11 а 12 • • • а 1т |
“ Р Г |
|
а 2 = |
а 21 а 22 ■■■а 2т |
р2 |
(3.52) |
(. |
|
|
|
___ |
|
|
|
й |
_O nl ^ п 2 - ^ n m _ |
|
|
1 з |
|
65 |
|
3 З а л . S27 |
|
|
Разумеется, если Мп есть подмножество из Nm, обе проекции сов падают, и полученный анализатор в точности соответствует желаемому. Если М п и N m совпадают, то матрица А неособенная, и предыдущие соотношения просто указывают на замену базиса в Мп (упражнение
2.16).
Пример 3.3. Проиллюстрируем изложенное на примере, с которым часто имеют дело в теории сигналов. Пусть х — сигнал, ограниченный по длительности в интервале [О, Т) и имеется переключатель, с помощью которого снимаются^значения|сигнала в m отдельных моментов вре мени на указанном интервидение. 3.10). Предположим также, что мы
Схема весового суммирования
a =(v. в.)
г* <'
Рис. 3.8. Модификация заданного анализатора для получения желаемого разложения сигнала.
хотим измерить коэффициенты разложения сигнала с периодом Т в ряд Фурье. Следовательно,
П
~. 2nit
x (t)~ S |
ai Q’ T , n ^ .m — 1, |
(3.53) |
. |
Я |
|
‘= - T
есть искомое разложение. Базисные функции имеют вид
. 2nit |
. 2nit |
|
Фг (0 = е ■Г |
; 0, (/) = -!_ е* г . |
(3.54) |
Полагая, что ползунок переключателя касается в каждый данный мо мент только одного полюса, и что время между касаниями равно нулю, примем в качестве базисной ортогональную систему прямоугольных функций
= |
(k— 1) |
; k = l, 2 , ... ,m, |
66
где |
|
|
|
|
u1(i) = . |
1 |
для О ^ t < |
т |
(3.55) |
|
О |
в других случаях. |
|
|
Отсюда следует, что взаимный по отношению к {ufe} базис |
есть |
|||
М 0 = -!|гМ 0 ; k = \,2 ...... т. |
(3.56) |
|||
Приближенное преобразование |
выборочных |
значений |3ft = |
(х, uft) |
в коэффициенты Фурье а г = (х, 0*) осуществляется с помощью мат рицы А, соответствующей (3.51):
Рис. 3.9. Иллюстрация приближенного пред ставления сигнала с помощью заданного анализатора.
|
,п |
t == 0, ± |
1, ± |
2, ..., ± |
ti |
|
(3.57) |
а г — S |
— , |
|
|||||
где |
k=i |
|
|
|
2 |
|
|
|
|
|
. 2 я it |
|
|
|
|
|
|
Г |
_ |
|
|
|
|
|
«4л = (ул Л ) = ~ j M O e 1 Т dt = |
|
|
||||
|
|
о |
|
|
|
|
|
|
^ |
. 2 nit |
|
2«1{к— т ) |
|
|
|
= — |
Г е |
' тdt — ^— |
шип |
е ' |
" |
. |
(3.58) |
Т |
J |
Т |
|
|
|
|
( k - l) T
Поскольку комплексные весовые коэффициенты физически не реализуемы, мы вместо них определяем действительную и мнимую
части отдельно: at = pt + jq^ Этому соответствует представление х (t) (в случае вещественных сигналов) тригонометрическим рядом:
^/2 |
9тг// п12 |
Quit |
(3.59) |
x(t) ~ Р 0+ 22PiCOs— ---- 2 |
2^ sin —- . |
||
i=i |
т i= i |
т |
|
Тем не менее удобнее рассуждать, пользуясь комплексной формой ряда Фурье, поскольку вывод соотношения (3.58) в этом случае зна
3* |
67 |
чительно проще. Измерение коэффициентов ряда (3.59) производится
спомощью схем весового суммирования, соответствующих матрицам В
иG размером (я/2) X т (рис. 3.10). Эти матрицы получены из дейст вительной и мнимой частей матрицы А:
тт
Р | = 2 |
|
bik$k> Я |
2 |
Ci hf o, |
k=l |
|
|
|
|
|
|
т |
j |
т |
Ро = |
|
2 aoh Ph = "IT 2 |
||
|
ft=l |
T |
k=l |
Отсчетные
Рис. 3.10. Анализатор гармоник.
где
bih= |
1 |
sin (ni/m) |
2ni |
k‘- T ) |
|
|
COS |
|
|||
|
|
ni/m |
|
|
|
|
1 |
sin (ni/m) |
2ni ^ k — |
||
Cib ■— |
' |
(3.60) |
|||
|
~ S i n |
||||
|
Г |
ni/m |
|
m |
|
i — 1, 2, ..., я/2; |
6 = 1, 2, ..., m. |
Упражнение 3.6. Реализуя гармонический анализ, как описано в примере 3.3, мы, естественно, интересуемся точностью, с которой отдельные коэффициенты Фурье могут быть найдены при идеальных переключателе, интеграторе и весовом сумматоре. Рассмотрев все сигналы единичной энергии из Z.2 (Т), показать, что верхняя грань ошибки в определении г-го коэффициента Фурье составляет
I at— at I2 |
1— |
sin (ni/m)\ 2“I |
|
nilm |
|||
|
|
68
У к а з а н и е :
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
|
m |
|
|
|
|
|
I К ; — |
a.i |
| = |
(х, |
Of)— |
2 |
(vft, 6;) (x, iift) |
<l|x||- |
0i— 2 |
(°i> vft) “ft |
|
||||||||||||
|
|
|
|
|
|
|
|
|
k = |
i |
|
|
|
|
|
|
II |
|
ft= 1 |
|
|
|
|
|
Предположить, что для t = |
0, ± |
1, ± |
2...... ± |
я/2 нужно обеспечить: |
|
|||||||||||||||||
|а г- — а; | < |
0,01/"l/Y |
и |
|| х || |
< |
1; |
каково |
должно |
быть яг? Как изменятся |
|||||||||||||||
результаты, если при том же количестве контактов и скорости вращения |
пере |
||||||||||||||||||||||
ключателя |
время его нахождения |
на каждом контакте уменьшить вдвое? |
|||||||||||||||||||||
|
|
|
СПИСОК ЛИТЕРАТУРЫ |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
1. |
P a p o u l i s |
A. |
The Fourier |
integral |
and |
its |
applications. |
M cGraw-Hill, |
|||||||||||||||
|
1962. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2. |
R i e s z |
F. |
and S z . - N a g y |
B. Functional |
analysis. |
Frederick |
Ungar, |
1955. |
|||||||||||||||
3. |
A x и e з e p |
H. |
|
И. и |
Г л а з м а н |
И. М. Теория линейных операторов в |
|||||||||||||||||
4. |
гильбертовом |
пространстве. М., Гостехиздат, |
1950. |
|
|
|
физики. |
||||||||||||||||
К у р а н т |
Р. |
и |
Г и л ь б е р т |
Д. |
Методы |
математической |
|||||||||||||||||
|
М .—Л., |
ГИТТЛ, |
1951. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
5. |
А х и е з е р |
Н. |
И. Лекции |
по |
теории |
аппроксимации. |
М., |
«Наука», |
|||||||||||||||
|
1965. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6. |
L е е |
Y. |
W. Statistical |
theory |
of |
communication. |
John |
Wiley and |
Sons, |
||||||||||||||
|
1960. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7. |
К a u t z |
W. |
H. Transient synthesis in the |
time domain, — «Trans. |
IRE», |
||||||||||||||||||
|
1954, v. CT-1, |
№ 3, |
p. 29—39. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
8. |
Y o u n g |
T. |
Y. |
a n d |
H u g g i n s |
W. |
H. Complementary signals and |
||||||||||||||||
9. |
orthogonalized exponentials.— «Trans. IRE», |
1962, v. CT-9, № 4, p. 362—370. |
|||||||||||||||||||||
W a 1 s h |
J. |
Z. A closed set of normal orthogonal |
functions. — «American |
||||||||||||||||||||
|
Jour., |
Math.», |
1923, v. 45, p. 5—24. |
|
|
|
|
|
|
|
|
|
|
||||||||||
10. |
H a m m о n d |
J. |
L. and |
J о h s о n |
R. |
S. Review of orthogonal |
square- |
||||||||||||||||
|
wave |
functions |
and |
their |
application |
to |
|
linear |
networks. — «J. Franklin |
||||||||||||||
11. |
Institute», 1962, v. 273, p. 211—225. |
concept |
of |
frequency |
and |
some |
applica |
||||||||||||||||
H a r m u t h |
H. |
F. A |
generalized |
||||||||||||||||||||
|
tions. — «Trans. IEEE», |
1968, v. IT-14, № 3, p. |
375—382. |
|
|
|
|
4
ИНТЕГРАЛЬНЫЕ ПРЕДСТАВЛЕНИЯ СИГНАЛОВ
4.1. НЕПРЕРЫВНЫЕ ПРЕДСТАВЛЕНИЯ
Представления сигналов, рассмотренные в предыдущих главах, называются дискретными потому, что каждый сигнал представляется счетной последовательностью величин. Основное преимущество ди скретных представлений в том, что от них просто перейти к приближен ным конечным представлениям, нужным для численных расчетов и при физических измерениях сигналов. Но в теоретических исследованиях, таких, например, как обсуждаемые в гл. 5 вопросы линейных преоб разований, желательно иметь точное представление сигнала, с кото рым можно оперировать, не заботясь о допускаемой ошибке. Для ши-
69
роких классов сигналов, таких, как L2 (— оо, оо), точные дискретные представления получаются достаточно громоздкими и с трудом подда ются простым физическим истолкованиям из-за того, что применяемые полные ортонормальные системы сложны по|своей'природе; в большей или меньшей мере все они «надуманы». Подходящее, казалось бы, раз ложение по смещенным во времени функциям {ф; (t) = ф ( t — гТ); i = 0 ,± 1, ± 2, ...} не является полным в L2 (Т) при любой исходной функции ф (0- Однако для практических сигналов ошибку приближе ния можно произвольно уменьшить, если выбрать Т достаточно малым. Мы приходим, естественно, к такому представлению сигнала, при ко тором параметр iT распределен непрерывно вдоль оси времени, а не в изолированных, равномерно расставленных точках.
Обобщение этой идеи непрерывного представления приводит к раз личным интегральным преобразованиям, формирующим соответствую щие интегральные представления сигнала [1]. Термин «непрерывный» не обозначает в данном случае обычного понятия непрерывности ото бражения, а указывает на то, что множество параметров, характери зующих используемый в представлении базис, плотно.
Данная глава посвящена некоторым общим вопросам интегральных представлений, рассматриваемых на основе аналогии или обобщения дискретных представлений. Стремясь сделать такую аналогию полной, мы сталкиваемся с рядом математических трудностей. Но несмотря на эти трудности, такой подход дает глубокое понимание сущности проблемы на физическом уровне. В части математической строгости следует полагаться на специальные исследования применительно к кон кретным интегральным преобразованиям [2]. Преобразования Фурье и Гильберта, широко используемые в теории сигналов, можно легко интерпретировать с рассматриваемой точки зрения. Действительно, вводные рассуждения по интегралам Фурье обычно базируются на ря дах Фурье. Мысль о том, что многие интегральные преобразования можно просто исследовать путем непосредственного обобщения ди скретных преобразований, оказывается весьма продуктивной и при водит к очень интересным результатам, которые обычно вполне приме нимы в практических приложениях.
4.2. БАЗИСНЫЕ И СОПРЯЖЕННЫЕ БАЗИСНЫЕ ЯДРА*’
Непрерывный аналог рассматриваемого в гл. 3 конечномерного представления получается, если заменить в (3.2) дискретный индекс i — номер базисной функции — континуальной переменной s ( S. Множество S обычно представляет собой некоторый интервал дейст вительной оси. Базис теперь выражается функцией ф (t, s), зависящей от двух переменных, и формула (3.2) принимает вид
x(t) = ^u(s)<p(t,s)ds;t£T, |
(4.1а) |
s |
|
*> Сопряженные ядра называют также обратными или взаимными. —
Прим. ред.
70