Преобразование Адамара (также известное как преобразование Уолша – Адамара, Адамара– Преобразование Радемахера – Уолша, преобразование Уолша или преобразование Уолша – Фурье ) является примером обобщенного класса преобразований Фурье. Он выполняет ортогональную, симметричную, инволютивную, линейную операцию с 2 действительными числами (или комплексные или гиперкомплексные числа, хотя сами матрицы Адамара являются чисто действительными).
Преобразование Адамара можно рассматривать как построенное из дискретных преобразований Фурье размера 2 (ДПФ), и фактически оно эквивалентно многомерному ДПФ размером 2 × 2 × ⋯ × 2 × 2. Он разлагает произвольный входной вектор на суперпозицию функций Уолша.
Преобразование названо в честь французского математика Жака Адамара, немецко-американский математик Ганс Радемахер и американский математик Джозеф Л. Уолш.
Преобразование Адамара H m - это матрица 2 × 2, матрица Адамара (масштабированная с помощью коэффициента нормализации), которая преобразует 2 действительных числа x n в 2 действительных числа X k. Преобразование Адамара можно определить двумя способами: рекурсивно или с помощью двоичного (base -2) представления индексов n и k.
Рекурсивно, мы определяем преобразование Адамара 1 × 1 H 0 с помощью тождества H0= 1, а затем определяем H m для m>0 по:
где 1 / √2 - нормализация, которая иногда опускается.
Для m>1 мы также можем определить H m следующим образом:
, где представляет произведение Кронекера. Таким образом, за исключением этого коэффициента нормализации, матрицы Адамара полностью состоят из 1 и -1.
Аналогично, мы можем определить матрицу Адамара ее (k, n) -м элементом, записав
где k j и n j - это двоичные цифры (0 или 1) для k и n соответственно. Обратите внимание, что для элемента в верхнем левом углу мы определяем: . В этом случае имеем:
Это как раз многомерный ДПФ, нормализованное до унитарного, если входы и выходы рассматриваются как многомерные массивы, индексированные n j и k j соответственно.
Далее следуют некоторые примеры матриц Адамара.
где - побитовое скалярное произведение двоичных представлений чисел i и j. Например, если , то , согласившись с вышеизложенным (игнорируя общую константу). Обратите внимание, что первая строка, первый элемент столбца матрицы обозначается .
H1это именно ДПФ размера 2. Его также можно рассматривать как преобразование Фурье на двухэлементной аддитивной группе Z / (2).
Строки матриц Адамара - это функции Уолша.
В классической области преобразование Адамара может быть вычислено за операций () с использованием быстрого преобразования Адамара алгоритм.
В квантовой области преобразование Адамара может быть вычислено за времени, поскольку это квантовая логика вентиль, который можно распараллелить.
В квантовой обработке информации преобразование Адамара, которое в данном случае чаще называют вентилем Адамара контекст (см. квантовый вентиль ), является одним кубитом вращением, отображающим базисные состояния кубита и в два состояния суперпозиции с равным весом вычислительного базисного состояний и . Обычно фазы выбираются так, чтобы
в нотации Дирака. Это соответствует матрице преобразования
в основа, также известная как. Состояния и известны как и соответственно, и вместе они составляют in квантовые вычисления.
многие квантовые алгоритмы используйте преобразование Адамара в качестве начального шага, поскольку оно отображает m кубитов, инициализированных с помощью к суперпозиции всех двух ортогональных состояний в на основе равного веса.
Примечательно, что вычисление квантового преобразования Адамара представляет собой простое применение логического элемента Адамара к каждому кубиту индивидуально из-за структуры тензорного произведения преобразования Адамара. Этот простой результат означает, что квантовое преобразование Адамара требует log n операций по сравнению с классическим случаем n log n операций.
Одно применение вентилей Адамара к 0 или 1 кубиту приведет к созданию квантового состояния, которое, если будет наблюдаться, будет равно 0 или 1 с равной вероятностью (как видно из первых двух операций). Это похоже на подбрасывание справедливой монеты в стандартной вероятностной модели вычислений . Однако, если вентиль Адамара применяется дважды подряд (как это фактически делается в последних двух операциях), то конечное состояние всегда совпадает с исходным.
Преобразование Адамара можно использовать для оценки филогенетических деревьев на основе молекулярных данных. Филогенетика - это подполе эволюционная биология сосредоточена на понимании взаимоотношений между организмами. Преобразование Адамара, применяемое к вектору (или матрице) частот паттернов сайтов, полученному из ДНК выравнивания множественных последовательностей, можно использовать для создания другого вектора, который несет информацию о топологии дерева. Обратимый характер филогенетического преобразования Адамара также позволяет вычислять вероятности сайта из вектора топологии дерева, что позволяет использовать преобразование Адамара для оценки максимального правдоподобия филогенетических деревьев. Однако последнее приложение менее полезно, чем преобразование вектора шаблона сайта в вектор дерева, потому что есть другие способы вычисления вероятностей сайта, которые намного более эффективны. Однако обратимая природа филогенетического преобразования Адамара действительно представляет собой элегантный инструмент для математической филогенетики.
Механика филогенетического преобразования Адамара включает вычисление вектора , который предоставляет информацию о топологии и длинах ветвей для дерева с использованием вектора или матрицы шаблона сайта .
где - матрица Адамара подходящего размера. Это уравнение можно переписать в виде серии из трех уравнений, чтобы упростить его интерпретацию:
Обратимый характер этого уравнения позволяет рассчитать вектор (или матрицу) ожидаемого шаблона сайта следующим образом:
Мы можем использовать модель замещения с двумя состояниями Кавендера-Фарриса-Неймана (CFN) для ДНК, кодируя нуклеотиды в виде двоичных знаков (пурины A и G кодируются как R, а пиримидины C и T кодируются как Y). Это позволяет кодировать множественное выравнивание последовательностей как вектор, который может быть преобразован в вектор дерева, как показано в следующем примере:
Индекс | Двоичный шаблон | Шаблоны выравнивания | ||||
---|---|---|---|---|---|---|
0 | 0000 | RRRR и YYYY | -0,475 | 0 | 1 | 0,6479 |
1 | 0001 | RRRY и YYYR | 0,2 | -0,5 | 0,6065 | 0,1283 |
2 | 0010 | RRYR и YYRY | 0,025 | -0,15 | 0,8607 | 0,02 |
3 * | 0011 | RRYY и YYRR | 0,025 | -0,45 | 0,6376 | 0,0226 |
4 | 0100 | RYRR и YRYY | 0,2 | -0,45 | 0,6376 | 0,1283 |
5 * | 0101 | RYRY и YRYR | 0 | -0,85 | 0,4274 | 0,0258 |
6 * | 0110 | RYYR и YRRY | 0 | -0,5 | 0.6065 | 0,0070 |
7 | 0111 | RYYY и YRRR | 0,025 | -0.9 | 0.4066 | 0.02 |
В примере, показанном в этой таблице, используется упрощенная схема с тремя уравнениями, и это для дерева четырех таксонов, которое может быть записано как ((A, B), (C, D)); в формате newick. Шаблоны сайта записываются в порядке ABCD. Это конкретное дерево имеет две длинные концевые ветви (0,2 трансверсионных замен на сайт), две короткие концевые ветви (0,025 трансверсионных замен на сайт) и короткую внутреннюю ветвь (0,025 трансверсионных замен на сайт); таким образом, он будет записан как ((A: 0,025, B: 0,2): 0,025, (C: 0,025, D: 0,2)); в формате newick. Это дерево будет демонстрировать притяжение длинных ветвей, если данные анализируются с использованием критерия максимальной экономии (при условии, что анализируемая последовательность достаточно длинна, чтобы наблюдаемые частоты паттернов сайтов были близки к ожидаемым частотам. показано в столбце ). Привлекательность длинной ветви отражает тот факт, что ожидаемое количество шаблонов сайтов с индексом 6, поддерживающих дерево ((A, C), (B, D)); - превышает ожидаемое количество шаблонов сайтов, поддерживающих истинное дерево (индекс 4). Очевидно, обратимая природа филогенетического преобразования Адамара означает, что древовидный вектор означает, что древовидный вектор соответствует правильному дереву. Таким образом, экономичный анализ после преобразования статистически согласован, как и стандартный анализ максимального правдоподобия с использованием правильной модели (в данном случае модели CFN).
Обратите внимание, что шаблон сайта 0 соответствует сайтам, которые не изменились (после кодирования нуклеотидов как пурины или пиримидины), индексы со звездочками (3, 5 и 6) являются «информативными для экономии» и. остальные индексы представляют собой паттерны участков, в которых один таксон отличается от трех других таксонов (таким образом, они эквивалентны длине конечных ветвей в стандартном филогенетическом дереве максимального правдоподобия).
Если кто-то желает использовать нуклеотидные данные без перекодирования как R и Y (и в конечном итоге как 0 и 1), можно закодировать данные как матрицу. Если мы рассмотрим дерево из четырех таксонов, то всего будет 256 паттернов сайтов (четыре нуклеотида в четвертой степени). Однако симметрии трехпараметрической модели Кимуры (или K81) позволяют уменьшить количество шаблонов сайтов с 256 до 64 шаблонов, что позволяет кодировать нуклеотидные данные для дерева с четырьмя таксонами как 8 x 8, аналогично вектору из 8 элементов, использованному выше для шаблонов сайтов трансверсии (RY). Это достигается путем перекодирования данных с использованием четырехгруппового кодирования Клейна :
Нуклеотид 1 | Нуклеотид 2 | Нуклеотид 3 | Нуклеотид 4 |
---|---|---|---|
A (0,0) | G (1,0) | C (0,1) | T (1, 1) |
C (0,0) | T (1,0) | A (0,1) | G (1,1) |
G (0,0) | A (1,0) | T (0,1) | C (1,1) |
T (0,0) | C (1,0) | G (0,1) | A (1,1) |
Как и в случае данных RY, шаблоны сайтов индексируются относительно к базе в произвольно выбранном первом таксоне с базами в последующих таксонах, закодированных относительно этой первой базы. Таким образом, первый таксон получает битовую пару (0,0). Используя эти битовые пары, можно создать два вектора, похожие на векторы RY, а затем заполнить матрицу, используя эти векторы. Это можно проиллюстрировать на примере Hendy et al. (1994), который основан на множественном выравнивании последовательностей четырех псевдогенов гемоглобина приматов:
0 | 8 | 16 | 24 | 32 | 40 | 48 | 56 | |
---|---|---|---|---|---|---|---|---|
0 | 8988 | 9 | 10 | 12 | 24 | 90 | ||
1 | 41 | 9 | ** | |||||
2 | 45 | 13 | ||||||
3 | 54 * | 14 | 3 | |||||
4 | 94 | 20 | ||||||
5 | 1 | |||||||
6 | 2 | 2 | ||||||
7 | 356 | 1 | 1 | 75 |
Гораздо большее количество паттернов сайтов в столбце 0 отражает тот факт, что столбец 0 соответствует переходным различиям, которые накапливаются быстрее, чем различия трансверсии практически во всех сравнениях геномных областей (и определенно быстрее накапливаются в псевдогенах гемоглобина, использованных для этого рабочего примера). Если мы рассмотрим шаблон сайта AAGG, он будет иметь двоичный шаблон 0000 для второго элемента битовой пары группы Клейна и 0011 для первого элемента. в этом случае двоичный шаблон, основанный на первом элементе первого элемента, соответствует индексу 3 (так, строка 3 в столбце 0; обозначена одной звездочкой в таблице). Шаблоны сайтов GGAA, CCTT и TTCC будут закодированы точно так же. Шаблон сайта AACT будет закодирован двоичным шаблоном 0011 на основе второго элемента и 0001 на основе первого элемента; это дает индекс 1 для первого элемента и индекс 3 для второго. Индекс, основанный на второй битовой паре группы Клейна, затем умножается на 8, чтобы получить столбец 24 строки 1 (ячейка обозначена двумя звездочками, но отсутствие числа в примере указывает на то, что выравнивание последовательностей не включает шаблоны сайтов AACT (аналогично, Шаблоны сайтов CCAG, GGTC и TTGA отсутствуют).
Преобразование Адамара также используется в шифровании данных, а также во многих обработка сигналов и алгоритмы сжатия данных , такие как JPEG XR и MPEG-4 AVC. В сжатие видео приложений, он обычно используется в форме суммы абсолютных преобразованных разностей. Он также является важной частью алгоритма Гровера и алгоритма Шора в квантовых вычислениях. Преобразование Адамара также применяется в экспериментальных методах, таких как ЯМР, масс-спектрометрия и кристаллография. Оно дополнительно используется в некоторых версиях местность-сенсити ve хеширование для получения псевдослучайных поворотов матрицы.