Сжатие сигналов МТКС. Преобразования сигналов. Кодеры с Субполосные преобразованием кодеры, страница 2

C0  d∞,

при этом

tdt=0.

Непрерывное изменение параметров а и b избыточно, поэтому вводят дискретные значения: a=am0 ; b=nb0a0m ; m,nZ , a01, b0≠0.

C увеличением параметра масштаба увеличивается размер шага сдвига (при анализе на низких частотах детали не важны) . при b0=1

m,ntn

тогда по аналогии с рядом Фурье сигнал можно представить в виде набора дискретных коэффициентов разложения: ∞ dm,nstdt.

−∞

восстановленный сигнал в этом случае:

stn

mZ nZ

Основы теории кратномасштабного анализа

При КМА сигнал представляется в виде совокупности его последовательных приближений

Теория КМА базируется на теории функциональных пространств. Под КМА понимается описание пространства сигналов через вложенные не пересекающиеся подпространства:

...⊂V 2V 1V 0V1V2⊂...

V m=0,             ∪V m=L2R , mZ.

при этом если st∈V m , то   s2t∈V m1 , а также существует базисная функция

0 t∈V 0 , и ее целочисленные сдвиги образуют ортонормированный базис  V 0 ,

при этом функции m,nt=2−m/2 2−m xn, образуют ортонормированный базис V m сигнал может быть представлен посредством последовательных приближений в V m st= lim smt,

m−∞

тогда появляется возможность анализа  сигнала на различных уровнях разрешения или масштаба. При этом базисная функция удовлетворяет первому основному уравнению КМА:

t.                                                           (!)

n

масштабирующая (scaling) функция и параметры           тесно связаны между собой: hn

m1,kt. p p

При этом выражение (!) можно представить в частотной области:

,

где                           H =∑hn⋅exp−j⋅⋅n  масштабирующий множитель уравнения КМА

n

При  рекурсивном применении  масштабирующего уравнения КМА:

     .

m=1    2

Свойства коэффициентов масштабирующего уравнения для скейлинг-функции:

Свойство 1. . n

Свойство 2.                                2⋅∑hnhn2k=k.

n

Свойство 3. ∣H ∣2∣H ∣2=1

Введем понятие подпространства           как ортогональное дополнение W m пространства                    доV m   V m1

V m1=V mW m ,       W m=0, ∪W m=L2R .

пусть t    образует ортонормированный базис пространства         , тогда W 0

t∈W 0⊂V −1           и можно представить как разложение по базису V −1

tn                             (!!)

n

Выражение (!!) представляет собой масштабирующее уравнение для вейвлет-функции по аналогии со скейлинг-функциями определим семейство вейвлет-функций m,nt=2−m/2 2−m tn.

Тогда масштабирующее уравнение для вейвлет-функции в частотной области имеет вид

,         =G.

Свойства вейвлет-функции:

Свойство 1.

hngn=0 ⇒ gn=−1 hn1

n

Свойство 2. Свойство 3.

G=−H −⋅exp−j⋅

0=0, ⇒ ∑gn=0;

n

n

В этом случае сигнал можно представить в виде проекции на           (грубой V m аппроксимации) и множества деталей как проекцию на W j

st=∑cm,n⋅m,nt∑d j,k⋅j,kt

n                           j

В качестве примера вейвлет-базиса можно рассмотреть вейвлет Хаара:

                                                                   h=[ 1      1 ];                      g=[ 1 ];

2 2                  2 2

Дискретное вейвлет-преобразование (ДВП, DWT)

Как было показано выше, базисные функции основных уравнений КМА зависят от коэффициентов уравнений КМА h и g. Т.о, возможно построение  базисных функций на основе коэффициентов уравнений КМА h и g, а также возможно производить анализ дискретных сигналов без использования базисных функций, путем итеративного вычисления коэффициентов разложения, получив таким образом полностью дискретный процесс декомпозиции:

c j,kn2k ;

n

dj,k=2∑c j−1,ngn2k .

Матричное описание:      n c1,0        h0 h1 h2 h3    0    0    0    0     c 0,0 c1,1 0    0 h0 h1 h2 h3        0    0    c0,1  v j1=M jv j; c1,2          0    0    0    0 h0 h1 h2 h3 c0,2

v jv j1;

представление посредством системы банк-фильтров (субполосного кодирования)

сигнал делится на 2 ветви: НЧ (H – ветвь), которая выделяет информацию и ВЧ (G – ветвь), определяющую детали. Применяя рекурсивно к H – ветви получим:

Заменим эквивалентной структурой системы анализа

Условие полного восстановления сигналов

O =[F H exp jprE G exp jqs]S ...