Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Френкс, Л. Теория сигналов

.pdf
Скачиваний:
80
Добавлен:
22.10.2023
Размер:
12.79 Mб
Скачать

Функции Эрмита

Для Т — (—оо, оо) и 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

Соседние файлы в папке книги из ГПНТБ