アダマール変換

ブール関数アダマール行列そのウォルシュスペクトルである: [1] (1, 0, 1, 0, 0, 1, 1, 0) × H(8) = (4, 2, 0, −2, 0, 2, 0, 2)
高速ウォルシュ・アダマール変換は、(1, 0, 1, 0, 0, 1, 1, 0) のウォルシュスペクトルをより高速に計算する方法です。
元の関数は、ウォルシュ スペクトルを使用して算術多項式として表現できます。

アダマール変換(ウォルシュ・アダマール変換アダマール・ラーデマッハ・ウォルシュ変換ウォルシュ変換ウォルシュ・フーリエ変換とも呼ばれる)は、一般化されたフーリエ変換の一例である。2 m 個の実数(または複素数複素数とも呼ばれるが、アダマール行列自体は純実数である)に対して直交対称反転、線形演算を行う

アダマール変換はサイズ2の離散フーリエ変換(DFT)から構築されていると見なすことができ、実際にはサイズ2 × 2 × ⋯ × 2 × 2の多次元DFTと同等です。[2]任意の入力ベクトルをウォルシュ関数 の重ね合わせに分解します

この変換は、フランスの 数学者 ジャック・アダマールフランス語: [adamaʁ])、ドイツ系アメリカ人の数学者ハンス・ラーデマッハー、およびアメリカの数学者ジョセフ・L・ウォルシュにちなんで名付けられました。

意味

アダマール変換H mは、 2 m  × 2 m行列(正規化係数でスケーリングされたアダマール行列)であり、2 m 個の実数x n を2 m 個の実数X kに変換します。アダマール変換は、再帰的に定義するか、インデックスnkの2進数基数-2)表現を用いて定義する2つの方法で定義できます

再帰的に、1 × 1 アダマール変換H 0を恒等式 H 0 = 1で定義し、次にm  > 0の場合のH m を次のように定義します。 ここで、1/ √2省略されることもある正規化です。

m  > 1の場合、 H m は次のように定義することもできます。ここで、 はクロネッカー積を表します。したがって、この正規化係数を除けば、アダマール行列は 1 と −1 のみで構成されます。

同様に、アダマール行列をその( kn )番目の要素 で定義すると、

ここで、k jn jはそれぞれknのビット要素(0または1)です。左上隅の要素については、次のように定義します。この場合、次の式が成り立ちます。

入力と出力がそれぞれn jk jでインデックス付けされた多次元配列とみなされる場合、これはまさにユニタリに正規化された多次元 DFT です。

アダマール行列の例をいくつか示します。ここでは、数値 i と j の2進表現のビット単位のドットです。例えば の場合、となり、これは上記と一致します(全体の定数は無視します)。行列の最初の行、最初の列の要素は で表されることに注意してください

H 1はまさにサイズ2のDFTである。これはZ /(2) の2元加法群上のフーリエ変換とみなすこともできる。

アダマール行列の行はウォルシュ関数です。

ウォルシュ・アダマール変換の利点

本物

上記の行列Hの定義に従って、ここではH = H [ m , n ] とします。

ウォルシュ変換では、行列には​​1と-1のみが現れます。1と-1は実数なので、複素数計算を行う必要はありません。

掛け算は不要

DFTでは無理数乗算が必要ですが、アダマール変換では必要ありません。有理数乗算も不要で、符号反転だけで済みます。

いくつかの特性はDFTの特性と似ている

ウォルシュ変換行列では、各行と各列はウォルシュ関数であり、符号の変化の回数は連続的に増加します。つまり、最初の行と最初の列では符号の変化はゼロです(すべての要素が1に等しい)。2行目と3列目では符号の変化が1回、3行目では符号の変化が2回、というように続きます。これを、各行iにゼロ交差が含まれる離散フーリエ変換と比較してください

離散フーリエ変換では、m が0 に等しい場合 (最初の行に対応)、結果も 1 になります。後続の行では、信号周波数が最初の生の行列では低く始まり、最後の行まで次の行に向かって増加するという行列の特性が見られます。

フーリエ変換との関係

アダマール変換は実際には2×2×⋯×2×2の多次元DFTと同等である。[2]

別のアプローチは、アダマール変換をブール群 上のフーリエ変換と見なすものである。[3] [4]有限(アーベル)群上のフーリエ変換を用いると、関数のフーリエ変換はによって定義される関数となる。ここで指標である。各指標はに対して の形を持ち、乗算はビット列上のブール内積である。したがって、 への入力をと同一視することができポンチャギン双対性)、によって定義される。

これは、とへの入力をブール文字列と見なしたのアダマール変換です

上記の定式化において、アダマール変換は左側の複素数ベクトルをアダマール行列で乗算するものですが、 を の要素のインデックスに対応するビット文字列を入力として取りの対応する要素を出力することで、同等性がわかります

これを、複素数ベクトルに適用するときに巡回群の指標を使用する通常の離散フーリエ変換と比較してください。

計算の複雑さ

古典的領域では、高速アダマール変換アルゴリズムを使用して、アダマール変換を演算()で計算できます

量子領域では、アダマール変換は並列化できる量子論理ゲートであるため、時間内に計算できます

量子コンピューティングの応用

アダマール変換は量子コンピューティングにおいて広く用いられている。2×2アダマール変換はアダマールゲートとして知られる量子論理ゲートであり、 -量子ビットレジスタの各量子ビットに並列にアダマールゲートを適用することは、アダマール変換と同等である

アダマール門

量子計算において、アダマールゲートは1量子ビットの 回転であり、量子ビット基底状態とを、計算基底状態とに等しい重みを持つ2つの重ね合わせ状態にマッピングする。通常、位相は次のように選択される。

ディラック記法で表されます。これは、計算基底とも呼ばれる基底変換行列に対応します。状態と はそれぞれと として知られ、量子コンピューティングにおける極基底を構成します

アダマールゲート演算

0または1の量子ビットにアダマールゲートを1回適用すると、観測された場合、0または1が等しい確率で得られる量子状態が生成されます(最初の2つの操作で見られるように)。これは、標準的な確率計算モデルにおいて、公平なコインを投げるのと全く同じです。しかし、アダマールゲートを2回連続して適用した場合(最後の2つの操作で実際に行われているように)、最終状態は常に初期状態と同じになります。

量子アルゴリズムにおけるアダマール変換

量子アダマール変換の計算は、アダマール変換のテンソル積構造のため、各量子ビットに個別にアダマールゲートを適用するだけです。この単純な結果は、量子アダマール変換は、古典的な演算の場合と比較して、演算を必要とすることを意味します。


量子ビット系では量子ビット(それぞれ に初期化されている)に作用するアダマールゲートは、が の形をとるときに均一な量子重ね合わせ状態を準備するために用いられる。この場合、量子ビットの場合、複合アダマールゲートはアダマールゲートのテンソル積として表される

結果として得られる均一な量子重ね合わせ状態は次のようになる。これは、任意のに対してアダマールゲートを用いた均一な量子状態の作成を一般化したものである[5]

この均一な量子状態を測定すると、~の間のランダム状態が得られます

多くの量子アルゴリズムは、アダマール変換を初期ステップとして用います。これは、前述のように、アダマール変換は、で初期化されたn個の量子ビットを、等しい重みを持つ基底 のすべての 2 n個の直交状態の重ね合わせにマッピングするためです。例えば、これはDeutsch–Jozsa アルゴリズムSimon のアルゴリズムBernstein–Vazirani アルゴリズムGrover のアルゴリズムで使用されています。Shorのアルゴリズムは、初期アダマール変換と量子フーリエ変換 の両方を使用することに注意してください。これらはどちらも有限群 上のフーリエ変換の一種で、前者は 上、後者は 上です

である一般的な場合における均一な量子重ね合わせ状態の準備は非自明であり、より多くの作業を必要とする。 ゲート複雑度と回路深度が全てに対してのみ であるような、重ね合わせ状態を準備するための効率的かつ決定論的なアプローチが最近発表された[6] 。 このアプローチに必要なの は 量子ビットのみである。重要なのは、このアプローチでは均一な重ね合わせ状態を作成するために、補助量子ビットや複数の制御を持つ量子ゲートは必要ないということである

分子系統学(進化生物学)の応用

アダマール変換は、分子データから系統樹を推定するために用いることができる[7] [8] [9] 系統学は、生物間の関係を理解することに焦点を当てた進化生物学の分野である。DNA多重配列アライメントから得られた部位パターン頻度のベクトル(または行列)にアダマール変換を適用すると、系統樹のトポロジーに関する情報を含む別のベクトルを生成することができる。系統樹アダマール変換の可逆性により、系統樹のトポロジーベクトルから部位尤度を計算することも可能となり、アダマール変換を用いて系統樹の最大尤度推定を行うことができる。しかし、後者の応用は、部位パターンベクトルから系統樹ベクトルへの変換ほど有用ではない。なぜなら、部位尤度を計算する他の方法[10] [11]の方がはるかに効率的だからである。しかし、系統樹アダマール変換の可逆性は、数理系統学にとって優れたツールを提供する。[12] [13]

系統学的アダマール変換の仕組みには、サイト パターン ベクトルまたは行列を使用してツリーのトポロジと枝の長さに関する情報を提供するベクトルの計算が含まれます

ここで、適切なサイズのアダマール行列です。この式は、解釈を簡略化するために、3つの式を連ねた式として書き直すことができます。

この方程式の可逆性により、次のようにして予想されるサイト パターン ベクトル (または行列) を計算できます。

DNAのCavender-Farris- Neyman(CFN)二状態置換モデルでは、ヌクレオチドをバイナリ文字(プリンAとGはR、ピリミジンCとTはY)としてエンコードすることで使用できます。これにより、多重配列アライメントをサイトパターンベクトルとしてエンコードし、これをツリーベクトルに変換することが可能になります。次の例をご覧ください。

特定の木のアダマール変換を示す例(実例の値はWaddell et al. 1997 [14]から改変)
索引バイナリパターン配置パターン
00000RRRRとYYYY−0.475010.6479
10001RRRYとYYYR0.2−0.50.60650.1283
20010RRYR と YYRY0.025−0.150.86070.02
3*0011RRYYとYYRR0.025−0.450.63760.0226
40100RYRRとYRYY0.2−0.450.63760.1283
5*0101RYRYとYRYR0−0.850.42740.0258
6*0110RYYRとYRRY0−0.50.60650.0070
70111RYYYとYRRR0.025−0.90.40660.02

この表に示す例では、簡略化された 3 つの方程式スキームを使用しており、これは 4 つの分類群のツリーに対するもので、newick 形式では ((A,B),(C,D)); と記述できます。サイト パターンは ABCD の順序で記述されます。この特定のツリーには、2 つの長い末端分岐 (サイトあたり 0.2 個のトランスバージョン置換)、2 つの短い末端分岐 (サイトあたり 0.025 個のトランスバージョン置換)、および 1 つの短い内部分岐 (サイトあたり 0.025 個のトランスバージョン置換) があるため、newick 形式では ((A:0.025,B:0.2):0.025,(C:0.025,D:0.2)); と記述されます。このツリーは、データが最大節約基準を使用して分析された場合に長い分岐の魅力を示します(分析された配列が、観測されたサイト パターン頻度が列に示されている予測頻度に近くなるように十分に長いと想定)。長い枝の吸引力は、インデックス6を持つサイトパターン(ツリー((A,C),(B,D));を支持する)の期待数が、真のツリー(インデックス4)を支持するサイトパターンの期待数を超えているという事実を反映している。系統的アダマール変換の可逆性は、明らかに、ツリーベクトルが正しいツリーに対応することを意味することを意味する。したがって、変換後の節約分析は統計的に一貫しており、[15]正しいモデル(この場合はCFNモデル)を用いた標準的な最尤分析と同様に一致する。

0のサイトパターンは、ヌクレオチドをプリンまたはピリミジンとしてエンコードした後、変化していないサイトに対応することに注意してください。アスタリスク付きのインデックス(3、5、6)は「簡略化情報」であり、残りのインデックスは、1つの分類群が他の3つの分類群と異なるサイトパターンを表します(したがって、標準的な最尤系統樹における末端の枝の長さに相当します)。

ヌクレオチドデータをRとY(そして最終的には0と1)として再コード化せずに使用したい場合は、部位パターンを行列としてコード化することができます。4分類群のツリーを考えると、合計256の部位パターン(4ヌクレオチドの4乗)があります。しかし、キムラの3パラメータ(またはK81)モデルの対称性により、DNAの256の可能な部位パターンを64パターンに減らすことができ、4分類群のツリーのヌクレオチドデータを、上記の転座(RY)部位パターンに使用した8要素のベクトルと同様の方法で、8×8行列[16]としてコード化することができます。これは、クラインの4群を使用してデータを再コード化することで実現されます

系統発生アダマール変換のためのクラインの4群コーディング
ヌクレオチド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)0,1C (1,1)
T (0,0)C (1,0)G (0,1)ア(1,1)

RYデータと同様に、サイトパターンは任意に選択された最初の分類群の塩基を基準にインデックス付けされ、後続の分類群の塩基はその最初の塩基を基準にエンコードされます。したがって、最初の分類群はビットペア(0,0)を受け取ります。これらのビットペアを用いて、RYベクトルに類似した2つのベクトルを作成し、それらのベクトルを用いてマトリックスを構築することができます。これは、4つの霊長類ヘモグロビン擬遺伝子の多重配列アライメントに基づくHendy et al. (1994) [16]の例を用いて説明できます。

エンコードされた配列アラインメントの例(Hendy et al. 1994 [16]より)(値は9879サイト中のカウント)
08162432404856
08988910122490
1419**
24513
354*143
49420
51
622
73561175

列 0 のサイト パターンの数が非常に多いのは、列 0 がトランジション差異に対応し、実質的にすべてのゲノム領域の比較においてトランスバージョン差異よりも急速に蓄積する (そしてこの実例で使用したヘモグロビン擬遺伝子では確実により急速に蓄積する[17] ) という事実を反映している。サイト パターン AAGG を考えてみると、クラインのグループのビット ペアの 2 番目の要素がバイナリ パターン 0000、最初の要素が 0011 になる。この場合、最初の要素に基づくバイナリ パターンでは、最初の要素はインデックス 3 に対応する (したがって、列 0 の行 3。表では 1 つのアスタリスクで示される)。サイト パターン GGAA、CCTT、および TTCC はまったく同じ方法でエンコードされる。サイト パターン AACT は、2 番目の要素に基づいてバイナリ パターン 0011 でエンコードされ、最初の要素に基づいて 0001 でエンコードされ、これにより、最初の要素のインデックスは 1、2 番目の要素のインデックスは 3 になる。 2 番目のクラインのグループのビット ペアに基づくインデックスに 8 を掛けて列インデックスを生成します (この場合は列 24 になります)。AACT サイト パターンの数を含むセルは 2 つのアスタリスクで示されます。ただし、例に数字がないことは、シーケンス アラインメントに AACT サイト パターンが含まれていないことを示します (同様に、同じ方法でエンコードされる CCAG、GGTC、および TTGA サイト パターンも存在しません)。

その他のアプリケーション

アダマール変換は、データ暗号化だけでなく、JPEG XRMPEG-4 AVCといった多くの信号処理データ圧縮 アルゴリズムにも用いられています。ビデオ圧縮アプリケーションでは、通常、変換された差分の絶対値和の形で用いられます。また、量子コンピューティングにおける多くのアルゴリズムの重要な部分でもあります。アダマール変換は、NMR質量分析結晶構造解析といった実験技術にも応用されています。さらに、局所性に敏感なハッシュ法のいくつかのバージョンでは、擬似ランダムな行列回転を得るために用いられています。

参照

  • リッター、テリー(1996年8月)「ウォルシュ・アダマール変換:文献概説」
  • Akansu, Ali N. ; Poluri, R. (2007年7月). 「直接拡散CDMA通信におけるウォルシュ型非線形位相直交符号」(PDF) . IEEE Transactions on Signal Processing . 55 (7): 3800–6 . Bibcode :2007ITSP...55.3800A. doi :10.1109/TSP.2007.894229. S2CID  6830633.
  • Pan, Jeng-shyang 離散フラクショナル アダマール変換を使用したデータ暗号化手法 (2009 年 5 月 28 日)
  • ラコヴィッツ、パウェル博士。ウォルシュ・アダマール変換と金融収益系列のランダム性検定(2015年4月7日)
  • Beddard, Godfrey; Yorke, Briony A. (2011年1月). 「アダマール変換を用いたポンプ・プローブ分光法」(PDF) . 2014年10月18日時点のオリジナル(PDF)からアーカイブ。 2012年4月28日閲覧
  • Yorke, Briony A.; Beddard, Godfrey; Owen, Robin L.; Pearson, Arwen R. (2014年9月). 「アダマール変換を用いた時間分解結晶構造解析」. Nature Methods . 11 (11): 1131– 1134. doi :10.1038/nmeth.3139. PMC 4216935.  PMID 25282611  .

参考文献

  1. ^ Townsend , WJ; Thornton, MA (2001). "Walsh spectrum computes using Cayley graphs". Proceedings of the 44th IEEE 2001 Midwest Symposium on Circuits and Systems (MWSCAS 2001) . MWSCAS-01. Vol. 1. IEEE. pp.  110– 113. doi :10.1109/mwscas.2001.986127. ISBN 97844222223 ... 0-7803-7150-X
  2. ^ ab Kunz, HO (1979). 「1次元離散ウォルシュ・アダマール変換と多次元離散フーリエ変換の等価性について」. IEEE Transactions on Computers . 28 (3): 267–8 . doi :10.1109/TC.1979.1675334. S2CID  206621901.
  3. ^ ブール写像のフーリエ解析–チュートリアル–、pp. 12–13
  4. ^ 講義5:基本的な量子アルゴリズム、ラジャット・ミッタル、pp.4~5
  5. ^ ニールセン、マイケル・A. ;チュアン、アイザック(2010). 『量子計算と量子情報』ケンブリッジ:ケンブリッジ大学出版局. ISBN 978-1-10700-217-3. OCLC  43641333。
  6. ^ Alok ShuklaとPrakash Vedula (2024). 「均一な量子重ね合わせ状態の準備のための効率的な量子アルゴリズム」.量子情報処理. 23:38 (1): 38. arXiv : 2306.11747 . Bibcode :2024QuIP...23...38S. doi :10.1007/s11128-024-04258-4.
  7. ^ ヘンディ, マイケル・D.; ペニー, デイヴィッド (1989年12月). 「進化樹の定量的研究のための枠組み」 .系統動物学. 38 (4): 297. doi :10.2307/2992396. JSTOR  2992396.
  8. ^ Hendy, Michael D.; Penny, David (1993年1月). 「系統発生データのスペクトル解析」 . Journal of Classification . 10 (1): 5– 24. doi :10.1007/BF02638451. ISSN  0176-4268. S2CID  122466038.
  9. ^ Székely, LA, Erdős, PL, Steel, MA, & Penny, D. (1993). 進化樹のためのフーリエ反転公式.応用数学レターズ, 6 (2), 13–16.
  10. ^ Felsenstein, Joseph (1981年11月). 「DNA配列からの進化樹形図:最大尤度法」 . Journal of Molecular Evolution . 17 (6): 368– 376. Bibcode :1981JMolE..17..368F. doi :10.1007/BF01734359. ISSN  0022-2844. PMID  7288891. S2CID  8024924.
  11. ^ Stamatakis, Alexandros (2019), Warnow, Tandy (ed.), 「系統学的尤度計算の最適化のためのアプローチのレビュー」 , Bioinformatics and Phylogenetics , Computational Biology, vol. 29, Cham: Springer International Publishing, pp.  1– 19, doi :10.1007/978-3-030-10837-3_1, ISBN 978-3-030-10836-6, S2CID  145834168 , 2020年10月10日取得
  12. ^ Chor, Benny ; Hendy, Michael D. ; Holland, Barbara R. ; Penny, David (2000-10-01). 「系統樹における複数の尤度最大値:解析的アプローチ」. Molecular Biology and Evolution . 17 (10): 1529– 1541. doi : 10.1093/oxfordjournals.molbev.a026252 . ISSN  1537-1719. PMID  11018159.
  13. ^ Matsen, Frederick A.; Steel, Mike (2007-10-01). Ané, Cécile ; Sullivan, Jack (編). 「単一系統樹上の系統発生的混合は、別のトポロジーの系統樹を模倣することができる」. Systematic Biology . 56 (5): 767– 775. arXiv : 0704.2260 . doi : 10.1080/10635150701627304 . ISSN  1076-836X. PMID  17886146.
  14. ^ Waddell, Peter J; Steel, MA (1997年12月). 「サイト間で不等な速度を持つ一般的な時間可逆距離:不変サイトにおけるΓ分布と逆ガウス分布の混合」. Molecular Phylogenetics and Evolution . 8 (3): 398– 414. Bibcode :1997MolPE...8..398W. doi :10.1006/mpev.1997.0452. PMID  9417897.
  15. ^ Steel, MA; Hendy, MD; Penny, D. (1993-12-01). 「倹約は一貫性を持つ!」 . Systematic Biology . 42 (4): 581– 587. doi :10.1093/sysbio/42.4.581. ISSN  1063-5157.
  16. ^ abc Hendy, MD; Penny, D.; Steel, MA (1994-04-12). 「進化樹の離散フーリエ解析」. Proceedings of the National Academy of Sciences . 91 (8): 3339– 3343. Bibcode :1994PNAS...91.3339H. doi : 10.1073/pnas.91.8.3339 . ISSN  0027-8424. PMC 43572. PMID 8159749  . 
  17. ^ Miyamoto, MM; Koop, BF; Slightom, JL; Goodman, M.; Tennant, MR (1988-10-01). 「高等霊長類の分子系統学:系譜関係と分類」. Proceedings of the National Academy of Sciences . 85 (20): 7627– 7631. Bibcode :1988PNAS...85.7627M. doi : 10.1073/pnas.85.20.7627 . ISSN  0027-8424. PMC 282245. PMID 3174657  . 
Retrieved from "https://en.wikipedia.org/w/index.php?title=Hadamard_transform&oldid=1319241538"