沃爾什轉換

沃爾什轉換的轉換式為

${\displaystyle F[m]=\sum _{n=0}^{N-1}f[n]W[m,n]}$

其中${\displaystyle W[m,n]}$是沃爾什轉換矩陣的第（m,n）個元素。 舉例來說，一個8點沃爾什轉換的轉換矩陣如下：

${\displaystyle W_{8}={\begin{pmatrix}1&1&1&1&1&1&1&1\\1&1&1&1&-1&-1&-1&-1\\1&1&-1&-1&-1&-1&1&1\\1&1&-1&-1&1&1&-1&-1\\1&-1&-1&1&1&-1&-1&1\\1&-1&-1&1&-1&1&1&-1\\1&-1&1&-1&-1&1&-1&1\\1&-1&1&-1&1&-1&1&-1\\\end{pmatrix}}}$

後面會解釋沃爾什轉換矩陣是如何產生，而沃爾什轉換的反轉換式為

${\displaystyle f[m]={\frac {1}{N}}\sum _{n=0}^{N-1}F[n]W[m,n]}$

注意到正轉換式與反轉換式只差了一個常數，這是由於沃爾什轉換矩陣的反矩陣就是自己的轉置矩陣乘上一個常數的緣故。

${\displaystyle 2^{k}}$點的沃爾什矩陣可以用下面的遞迴方式產生：

起始值${\displaystyle k=1}$（2點沃爾什轉換矩陣）

${\displaystyle W_{2}={\begin{pmatrix}1&1\\1&-1\\\end{pmatrix}}}$

假設我們已經有一個${\displaystyle 2^{k}}$點的沃爾什轉換矩陣${\displaystyle W_{2^{k}}}$ 則我們可以藉由下面的方法來產生${\displaystyle 2^{k+1}}$點的沃爾什轉換矩陣${\displaystyle W_{2^{k+1}}}$

Step 1 定義${\displaystyle V_{2^{k+1}}={\begin{pmatrix}W_{2^{k}}&W_{2^{k}}\\W_{2^{k}}&-W_{2^{k}}\\\end{pmatrix}}}$

Step 2 根據變號次數把${\displaystyle V_{2^{k+1}}}$的列（row）重新排列成為${\displaystyle W_{2^{k+1}}}$

以下舉一個使用4點沃爾什轉換矩陣產生8點沃爾什轉換矩陣的例子：

${\displaystyle W_{4}={\begin{pmatrix}1&1&1&1\\1&-1&1&-1\\1&1&-1&-1\\1&-1&-1&1\\\end{pmatrix}}}$

${\displaystyle V_{8}={\begin{pmatrix}W_{4}&W_{4}\\W_{4}&-W_{4}\\\end{pmatrix}}={\begin{pmatrix}1&1&1&1&1&1&1&1\\1&-1&1&-1&1&-1&1&-1\\1&1&-1&-1&1&1&-1&-1\\1&-1&-1&1&1&-1&-1&1\\1&1&1&1&-1&-1&-1&-1\\1&-1&1&-1&-1&1&-1&1\\1&1&-1&-1&-1&-1&1&1\\1&-1&-1&1&-1&1&1&-1\\\end{pmatrix}}}$

接著對${\displaystyle V_{8}}$的列做排序即可得上面的${\displaystyle W_{8}}$。

在不同的應用上，我們較常使用的沃爾什矩陣的列的排列順序也不同，以下以一個表來區分：

序數順序（沃爾什順序）	雙積順序（培力順序）	自然順序（哈德碼得順序）	${\displaystyle W[m,n]}$
Sequency Ordering(Walsh Ordering)	Dyadic Ordering(Paley Ordering)	Natural Ordering（Hadamard Ordering）
較常用於信號處理	較常用於控制工程	較常用於數學
第0列	第0列	第0列	${\displaystyle {\begin{pmatrix}1&1&1&1&1&1&1&1\\\end{pmatrix}}}$
第1列	第1列	第4列	${\displaystyle {\begin{pmatrix}1&1&1&1&-1&-1&-1&-1\\\end{pmatrix}}}$
第2列	第3列	第6列	${\displaystyle {\begin{pmatrix}1&1&-1&-1&-1&-1&1&1\\\end{pmatrix}}}$
第3列	第2列	第2列	${\displaystyle {\begin{pmatrix}1&1&-1&-1&1&1&-1&-1\\\end{pmatrix}}}$
第4列	第6列	第3列	${\displaystyle {\begin{pmatrix}1&-1&-1&1&1&-1&-1&1\\\end{pmatrix}}}$
第5列	第7列	第7列	${\displaystyle {\begin{pmatrix}1&-1&-1&1&-1&1&1&-1\\\end{pmatrix}}}$
第6列	第5列	第5列	${\displaystyle {\begin{pmatrix}1&-1&1&-1&-1&1&-1&1\\\end{pmatrix}}}$
第7列	第4列	第1列	${\displaystyle {\begin{pmatrix}1&-1&1&-1&1&-1&1&-1\\\end{pmatrix}}}$

若使用二進位來表示各種順序的列的編號，則雙積順序的二進位編號是序數順序的格雷碼編碼，而自然順序的二進位編號是雙積順序的位元反轉。

（1） 正交性質

沃爾什轉換矩陣的每個列是互相正交的，即如果${\displaystyle m_{0}\neq m_{1}}$則${\displaystyle \sum _{n=0}^{N-1}W[m_{0},n]W[m_{1},n]=0}$

（2） 零交（zero-crossing）性質

${\displaystyle 2^{k}}$點沃爾什轉換矩陣的每個列的變號次數都不相同，分別為變號0次到變號${\displaystyle 2^{k}-1}$，這個性質是沃爾什轉換可以用來做頻譜分析的原因之一，不同的變號次數相當於不同的頻率。

（3） 奇偶性質

沃爾什轉換矩陣（沃爾什順序）中，編號為偶數的列是偶對稱，編號為奇數的列是奇對稱。（有第0列）

（4） 線性性質

若${\displaystyle f[n]\Rightarrow F[m]}$，${\displaystyle g[n]\Rightarrow G[m]}$，（${\displaystyle \Rightarrow }$表沃爾什轉換）則有${\displaystyle af[n]+bg[n]\Rightarrow aF[m]+bg[m]}$。

（5） 加法性質

${\displaystyle W[m,n]W[l,n]=W[m\oplus l,n]}$，${\displaystyle \oplus }$表示邏輯互斥或（exclusive or）

（6） 平移性質

若${\displaystyle f[n]\Rightarrow F[m]}$則${\displaystyle f[n\oplus k]\Rightarrow W[k,m]F[m]}$

（7） 調變性質

若${\displaystyle f[n]\Rightarrow F[m]}$則${\displaystyle W[k,n]f[n]\Rightarrow F[k\oplus m]}$

（8） 巴斯瓦定理（Parseval's Theorem）

若${\displaystyle f[n]\Rightarrow F[m]}$則${\displaystyle \sum _{n=0}^{N-1}|f[n]|^{2}={\frac {1}{N}}\sum _{n=0}^{N-1}|F[m]|^{2}}$

（9） 摺積性質

若${\displaystyle f[n]\Rightarrow F[m]}$，${\displaystyle g[n]\Rightarrow G[m]}$，則${\displaystyle \sum _{l=0}^{N-1}f[l]g[((n\oplus l))_{N}]=\sum _{l=0}^{N-1}g[l]f[((n\oplus l))_{N}]\Rightarrow F[m]G[m]}$，在這裡${\displaystyle \sum _{l=0}^{N-1}f[l]g[((n\oplus l))_{N}]=\sum _{l=0}^{N-1}g[l]f[((n\oplus l))_{N}]}$代表邏輯摺積（logical convolution）。

由於一個${\displaystyle 2^{k}}$點沃爾什轉換矩陣可以由${\displaystyle 2^{k-1}}$點的沃爾什轉換矩陣堆疊後做變號與排序產生，因此一個${\displaystyle 2^{k}}$沃爾什轉換可以由做兩次${\displaystyle 2^{k-1}}$的沃爾什轉換及一些加減法和排序產生，可以得到一個類似快速傅立葉變換的蝴蝶結構。

1. 運算數值皆為實數不存在複數
2. 不需要應用到乘法運算
3. 頻譜分析( spectrum analysis)
4. 正向轉換跟反向轉換結構相似
   • Forward :${\displaystyle F[m]=\textstyle \sum _{n=0}^{N-1}\displaystyle f[n]W[m,n]}$
   • Inverse : ${\displaystyle f[m]=1/N\textstyle \sum _{n=0}^{N-1}\displaystyle F[n]W[m,n]}$
5. 跟DFT有一樣的以下性質
   1. 正交性質
   2. 零交（zero-crossing）性質
   3. 線性運算性質
   4. 調變性質
   5. 巴斯瓦定理（Parseval's Theorem)

沃爾什轉換適合做頻譜分析，但未必適合做摺積

因此不適合分析線性非時變系統(LTI system)

證明:

• Linear convolution (standard form of convolution)

${\displaystyle \textstyle \sum _{n=0}^{N-1}\displaystyle f[n]g[n-l]}$

• Circular convolution

${\displaystyle \textstyle \sum _{n=0}^{N-1}\displaystyle f[n]g[(n-l)_{N}]}$

但是在Walsh Transform中Convolution Property 代表著"logic convolution"

• Logic convolution

${\displaystyle h[n]=f[n]*g[n]=\textstyle \sum _{n=0}^{N-1}\displaystyle f[n]g[(n\oplus l]=\textstyle \sum _{n=0}^{N-1}\displaystyle f[n\oplus l]g[l]}$

* 代表 "Logic convolution"

${\displaystyle \oplus }$ 代表 "logic addition" (similar to XOR) ,for example 3${\displaystyle \oplus }$7=4

${\displaystyle {\begin{array}{lcr}3&011\\\oplus 7&111\end{array}}}$

• XOR(1,0)=XOR(0,1)=1

• XOR(0,0)=XOR(1,1)=0

以訊號處理的角度,Circular convolution 可以滿足 Linear convolution ,即可分析LTI系統

假使Logic convolution 與 Circular convolution 結果不一致則無法分析 LTI系統

舉例來說:

• 當N=8
• Circular convolution

H[2]=f[0]g[2]+f[1]g[1] + f[2]g[0] + f[3]g[7] + f[4]g[6] + f[5]g[5] + f[6]g[4]+ f[7]g[3]

Logical convolution

h[2] = f[0]g[2] + f[1]g[3] + f[2]g[0] + f[3]g[1] + f[4]g[6] + f[5]g[7] + f[6]g[4]+ f[7]g[5]

由上面可知，兩者結果不一致

CDMA 領域

• 主要應用在多工，其中CDMA為主要應用
  若使用N點沃爾什轉換，則可以對N個通道做多工
  而且使用沃爾什轉換的好處是不需要同步
  其他正交轉換則需要同步
  
  舉例：CDMA使用沃爾什轉換做多工的方法
  假設現在有兩組資料要傳，分別是[1,0,1],[1,1,0]
  並且使用8點沃爾什轉換${\displaystyle W_{8}}$的第一行與第二行來當作通道一與通道二的正交基底
  1.將0變為-1
  [1,0,1]→[1,-1,1]
  [1,1,0]→[1,1,-1]
  2.調變
  對於第一組資料拿通道一來調變
  第一組資料為[1,-1,1],通道一為[1,1,1,1,1,1,1,1]
  →[1*通道一, -1*通道一, 1*通道一]
  →[1,1,1,1,1,1,1,1, -1,-1,-1,-1,-1,-1,-1,-1, 1,1,1,1,1,1,1,1]
  
  對於第二組資料拿通道二來調變
  第二組資料為[1,1,-1],通道二為[1,1,1,1,-1,-1,-1,-1]
  →[1*通道二, 1*通道二, -1*通道二]
  →[1,1,1,1,-1,-1,-1,-1, 1,1,1,1,-1,-1,-1,-1, -1,-1,-1,-1,1,1,1,1]
  3.相加
  [1,1,1,1,1,1,1,1, -1,-1,-1,-1,-1,-1,-1,-1, 1,1,1,1,1,1,1,1]+[1,1,1,1,-1,-1,-1,-1, 1,1,1,1,-1,-1,-1,-1, -1,-1,-1,-1,1,1,1,1]
  →[2,2,2,2,0,0,0,0,0,0,0,0,-2,-2,-2,-2,0,0,0,0,2,2,2,2]
  
  

• 解調
  1.如果使用N點沃爾什轉換，則把接收的訊號每隔N點拆開來
  [2,2,2,2,0,0,0,0,0,0,0,0,-2,-2,-2,-2,0,0,0,0,2,2,2,2]
  →[2,2,2,2,0,0,0,0], [0,0,0,0,-2,-2,-2,-2], [0,0,0,0,2,2,2,2]
  2.將每段訊號與通道做內積
  若大於0,則解調為1
  若小於0,則解調為0
  對於通道一
  [2,2,2,2,0,0,0,0]·[1,1,1,1,1,1,1,1]=8 → 1
  [0,0,0,0,-2,-2,-2,-2]·[1,1,1,1,1,1,1,1]=-8 → 0
  [0,0,0,0,2,2,2,2]·[1,1,1,1,1,1,1,1]=8 → 1
  通道一接收的資料為[1,0,1]
  對於通道二
  [2,2,2,2,0,0,0,0]·[1,1,1,1,-1,-1,-1,-1]=8 → 1
  [0,0,0,0,-2,-2,-2,-2]·[1,1,1,1,-1,-1,-1,-1]=8 → 1
  [0,0,0,0,2,2,2,2]·[1,1,1,1,-1,-1,-1,-1]=-8 → 0
  通道二接收的資料為[1,1,0]
  

• Bandwidth reduction
• High resolution
• Information coding
• Feature extraction
• ECG(心電圖) signal (in medical signal processing) analysis
• Hadamard spectrometer
• Avoiding quantization error

• 阿达马矩阵
• Hadamard變換
• Hadamard矩陣

• Jian-Jiun Ding, Advanced Digital Signal Processing class note,the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2008.
• H. F. Harmuth,“Transmission of information by orthogonal functions,”1970.
• Moon-Hu. Lee,“A new reverse Jacket transform and its fast algorithm,”IEEE Trans. Circuits Syst.-II, vol. 47, pp.39-46, 2000.
• K.G.Beauchamp, "Walsh Functions and Their Applications," Academic Press,1975.
• H. F. Harmuth, "Transmission of Information by Orthogonal Functions," Springer, 1969.