Файл: Френкс, Л. Теория сигналов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 116

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Отсюда следует, что вектор х — х ортогонален ко всем векторам в М п.

Для того чтобы показать, что ||х — х || — минимальная

норма из всех

|| х — х | | , рассмотрим произвольный вектор х 6 М п:

 

1х х р =! (х

х) (х х) ||2 = (х х, х —х) —

—(х—X, X —х) —(х —X, X —х) + (х—X, X —х).

Поскольку х — х £ М п,

средние слагаемые пропадают, и мы имеем

||х - х ||2 = ||х —i f + «x —i f .

(3.11)

Ясно, что минимум достигается при х = х. Назовем х ортогональной проекцией х на Мп ', а ц =^х — х — погрешностью приближения х век­ тором х. Точность приближения численно характеризуется нормой тр Положив в (3.11) х — 0,'получим

h f = ||x |M )ip .

 

(3-12)

Теперь ясно, что отношения эквивалентности,

соответствующие

разбиению L2 (Т) согласно (3.8), могут быть записаны следующим об­

разом:

 

 

х~у=*>(х, Ф{) — (у, <р4); i ■—1, 2,...,

п,

(3.13)

что соответствует приведенному во введении примеру 1.7, иллюстри­ рующему отношения эквивалентности и разбиения. Из рассмотрения этого примера следует, что (3.13) может быть записано также иначе:

х ~ у = ^ х —увМ ,

*

(3-14)

где М — {х; (х, х) = О для всех х £ Мп}.

Легко показать,

что М

в (3.14) есть линейное подпространство. М называется ортогональным дополнением М п, потому что любой вектор из L2 (Т) может быть един­ ственным образом представлен суммой вектора из Мп и вектора из Д4, причем эти векторы ортогональны, т. е. для любого х имеет место:

x = x-fz; х£7Ип, z £ М, (х, z) = 0.

(3.15)

Итак, естественно рассматривать L2 (Г) как прямую сумму подпрост­

ранств М п и М, т. е. L2 (Т) = М п +

М. Заметим,

что единственным

общим элементом М п и М является нулевой вектор.

Все эти понятия

графически поясняются на рис. 3.1.

*

 

Упражнение 3.2.^Прямые суммы и проекции. Можно сказать, что линейное

пространство L есть сумма подпространств

М и N, т. е. L = М + N, в том

случае, когда L является пространством, натянутым на наибольшее количество линейно независимых векторов, содержащихся в М и N совместно. Тогда любое г £ L может быть представлено как г = х + у, где х £ М, а у £ N. Показать,

что разложение г = х + у единственно тогда и только тогда,

когда М f) jV =

= {0}, т. е. если общим для М и N вектором является только нулевой вектор.

В этом случае говорят, что L есть прямая сумма М и N, т. е.

L =

М 0

N. За­

метьте, что М и N не обязательно ортогональны.

Если L = М ф

N и z = х +

+ у — любой вектор из L, то х £ М называется

проекцией

z на

М

вдоль N

51


(рис. 3.2). Пусть {ср;; i = I,

2, ... >

является базисом

в

L,

и пусть {в*;

i =

=

1, 2...} есть взаимный базис, и,

наконец,

пусть М натянуто на

базис

{<р;;

i =

1, 2, ..., п) и

N натянуто на базис

i = п +

1,

п +

2, ...}.

Показать,

что проекция z на М вдоль N может быть представлена в виде

 

 

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

Х =

2 (z, 90<Pi.

 

 

 

 

 

 

 

1, 2, ..., n}

 

i—1

 

взаимным базисом в М.

причем {0j; i =

не обязательно

является

Понятие неортогональной проекции будет полезно при рассмотрении некоторых вопросов обнаружения сигналов (см. гл. 10).

Пример

3.1. Пусть

М 3 есть подпространство из

L2

(Т), Т =

= [0, оо], натянутое на

действительные экспоненциальные функции

{е~*, e~2t, е-3 /}. Пусть

требуется найти

наилучшее

приближение

в М з для прямоугольного импульса: х (t) =

1 для 0 ^

t ^

х/2; л: (() =

= 0 для t >

V2. Можно записать

,,

 

 

 

со

 

 

 

 

(фь>

Фт) = ^е-(*+"г) {dt= U(k + m)\ k, m= 1,

2,

3.

 

о

 

 

 

 

Рис. 3.1. Иллюстрация ортогонального проектирования x eL 2(T) на конечномерное подпространство.

Следовательно, в (3.4) матрица G имеет вид:

- 1 1 1 ~

 

 

 

G =

72

—240

180

—240

900

—720

 

180

—720

600

Согласно (3.5) в'заимный баЗис есть:

 

 

01(/) = 72е-<—240e~w+

180e~3^.

 

02 (t) =•—240е~' + 900е-«—720е-3',

 

03 (0 = 180е-*—720е-2^ + 600е-3<

 

i/г

( « 1 = —0,912,

«г = (х, 6i)= ^

ос2 = 3,57,

 

0

1 « з= — 1.36,

52


следовательно.

 

 

х (t) = 0,912е_/ -f 3,57e_2/— 1,36e-3'

и

« = {—0,912; 3,57; — 1,36}

есть представление x (t) в R3. Погрешность аппроксимации составляет

IIТ) IP = II х IP—1| X IP =>Iti |j = 0,289.

Приближение

иллюстрирует­

 

 

 

ся на рис. 3.3.

 

 

полу­

 

 

 

 

Пример

3.2. Для

 

 

 

чения

взаимного

базиса не

 

 

 

всегда нужно находить обрат­

 

 

 

ную матрицу, что очень труд­

 

 

 

но при больших п.

Развивая

 

 

 

предыдущий пример, мы про­

 

 

 

демонстрируем

способ

оты­

 

 

 

скания

взаимного базиса для

 

 

 

Мп,

натянутого

на

(e~kt;

 

 

 

k =

1,2,

...,

п }

при

Т —

 

 

 

= [0, оо3и произвольного п. Нам нужны такие 0т , что

 

 

 

 

(Ф*. 0rn) =

S е-« 0m (t) dt = Qm(s)\s=k= 8hm,

(3.16)

 

 

 

 

 

 

о

 

 

 

где 0 т (s) есть преобразование Лапласа от 0m (t)

Поскольку 0m (t)

в общем случае есть линейная комбинация

{cpft; &=

1, 2,

..., п}, при

любом m ©m

(s) есть рациональная функция s с полюсами в точках

s =

— 1, —2,

..., —п. С другой стороны, из (3.16) следует, что ©m(s)

обращается в нуль в точках s = 1, 2, 3, ...,

п, исключая точку s = m.

Рис. 3.3.

Аппроксимация

прямоугольного импульса

х с

помощью

*> Для функций

времени,

определенных на Т — [0,

оо),

преобразование

 

 

ОО

 

 

 

Лапласа имеет

вид:

х (s) = \

х (t) e~~st dt, где s — комплексная переменная.

 

 

о

 

 

через их нули

Это преобразование применяется для представления сигналов

и полюсы (в плоскости s) и для определения самих сигналов с помощью вычетов.

Мы предполагаем,

что читатель знаком с этими методами. Для справок может

быть использована

гл. 9 в [1].

•*

53


Поэтому ©m (s) может быть записана следующим образом:

П

em(s)=

Кт П (S-A)

(3 .17)

^ -------

(■s—m) П (s + k) k= i

Согласно (3.16) 0 m (m) — 1 и мы можем определить также константу

Кт-

са, п)= п

(j—k)=>Km

( - 1 ) " C ( - m , п)

(3.18)

k=\.

 

С (т, п)

 

Ьф}

 

 

 

Наконец, найдем обратное преобразование Лапласа обычным спо­ собом, с помощью вычетов

 

 

П

 

 

 

0т (0 = 2

Уml фДО.

(3.19)

 

 

l=

1

 

где

YmZ (S 0 ©m (®) |s= —/

С (—т, п) С (—I, п)

 

(1-\-т)С (т, п)С(1, п)

 

 

 

 

 

Упражнение 3.3. Применяя (3.19),

найти 0; (t) примера

3.1.

3.2.ПОЛНЫЕ ОРТОНОРМАЛЬНЫЕ СИСТЕМЫ

Впредыдущем параграфе мы применили теорему проектирования

для нахождения наилучшего представления произвольного сигнала в конечномерном подпространстве М п, но мы обошли вопрос о том, как выбрать подходящее подпространство. Ясно, что для любого заданного подпространства более обширное подпространство из L2 (Т) не пред­ ставляется вполне адекватно: это следует из существования векторов, ортогональных к М п. Задача нахождения оптимального подпростран­ ства имеет смысл только по отношению к некоторому ограниченному (обычно компактному) исходному подпространству из L2 (Т). Такого рода задача рассматривается в гл. 6. В настоящем параграфе мы рассмотрим более простой вопрос о сходимости представления, полагая, что число измерений подпространства можно произвольно увеличивать, и приведем ряд примеров, демонстрирующих часто применяемые ме­ тоды нахождения систем базисных функций, безотносительно к опти­ мальности или неоптимальности натянутых на них подпространств.

Во многих курсах функционального анализа [2, 3] доказывается, что L2(T) есть полное сепарабельное пространство. Из этого свойства следует, что ортогональным проектированием можно получить сколь

угодно близкую аппроксимацию для любого х £ L2

(Г), если выбрать п

достаточно большим

 

=2( х , фг)срг.

(3.20)

/= 1

 

54


В (3.20) функции фг взяты из бесконечного счетного множества ортонормальных функций, удовлетворяющих некоторым условиям.

Заметим вначале, что из (3.12) следует при х = хп

I! X---Х„ II2 = ЦX Ц3--- 2 |(Х,ф*)|2.

г= 1

Поэтому, при любом п

 

2 |( х ,ф г)|2< ||х ||2.

(3.21)

г=I

 

Соотношение (3.21) известно как неравенство Бесселя, оно показывает, что сумма квадратов коэффициентов разложения (х, ф;) ограничена для

любого х ( L2

(Т). Из

неравенства

(3.21) видно также,

что {х„}

в (3.20) есть последовательность Коши, так как

для любого е > 0

при достаточно большом п0имеем

 

 

 

I I

—хт ||2 =

2

|(х, фг)|3< 82; п > т > п 0

(3.22)

 

 

 

г'=т+1

 

 

 

и, поскольку

L2

(Т) — полное пространство,

последовательность

{хп} сходится к некоторой точке в L2 (Т)

. Последовательность сходится

к х, если {фг;

i =

1, 2,

...}

есть полная ортонормальная система.

Ортонормальная система называется полной, если не существует до­ полнительных, отличных от нуля ортогональных векторов, которые можно было бы прибавить к этой системе. Благодаря сепарабельности / L2 (Г) полная ортонормальная система является счетной. Очевидно, полная ортонормальная система является для L2 (Т) аналогом базиса

в конечномерном пространстве. Произвольная бесконечная ортонор­ мальная система не обязательно полна; например, использованная в (1.34) для разложения во временной ряд система функций

У 2W sin 2nW [t— {H2W)\

i 0, + 1, + 2,...

2nW [t~(ii2W)}

:

является ортонормальной, но не полной в L2 (— оо, оо), поскольку функции с полосой больше W не принадлежат подпространству, на­ тянутому на эту систему. Заметим, что не следует смешивать понятия полноты ортонормальной системы и полноты метрического простран­ ства.

Для полной ортонормальной системы неравенство Бесселя (3.21) переходит в равенство

оо

 

2 |(х, фг)|2 = ||х ||2 для любого х £ L2(T),

(3.23)

/= 1

 

которое известно так же, как условие полноты. На основании сказан­ ного выше можно утверждать, что для любого е > 0 и любого х£ L2 (Т) имеется такое п0, что

П

(х,ф;)фг < е при п > я 0.

(3.24)

/ = 1