본문으로 이동

아다마르 변환

한울위키, 우리 모두의 백과사전.
imported>InternetArchiveBot님의 2025년 11월 3일 (월) 08:08 판 (1 개의 출처 구조, 0 개의 링크를 깨진 것으로 표시) #IABot (v2.0.9.5)
(차이) ← 이전 판 | 최신판 (차이) | 다음 판 → (차이)
파일:1010 0110 Walsh spectrum (single row).svg
행렬 곱불 함수아다마르 행렬월시 스펙트럼:[1]
(1, 0, 1, 0, 0, 1, 1, 0) × H(8) = (4, 2, 0, −2, 0, 2, 0, 2)
파일:1010 0110 Walsh spectrum (fast WHT).svg
(1, 0, 1, 0, 0, 1, 1, 0)의 월시 스펙트럼을 계산하는 더 빠른 방법인 고속 월시-아다마르 변환
원래 함수는 월시 스펙트럼을 통해 산술 다항식으로 표현될 수 있다.

아다마르 변환(영어: Hadamard transform, 월시-아다마르 변환(Walsh–Hadamard transform), 아다마르-라데마허-월시 변환(Hadamard–Rademacher–Walsh transform), 월시 변환(Walsh transform), 또는 월시-푸리에 변환(Walsh–Fourier transform)이라고도 함)은 일반화된 푸리에 변환의 한 예시이다. 이는 2m개의 실수 (또는 복소수 또는 초복소수)에 대해 직교, 대칭, 대합적, 선형 연산을 수행한다 (아다마르 행렬 자체는 순수하게 실수이다).

아다마르 변환은 크기 2의 이산 푸리에 변환(DFT)으로 구성된 것으로 간주될 수 있으며, 실제로는 2 × 2 × ⋯ × 2 × 2 크기의 다차원 DFT와 동일하다.[2] 이는 임의의 입력 벡터를 월시 함수의 중첩으로 분해한다.

이 변환은 프랑스수학자자크 아다마르(프랑스어: [adamaʁ]), 독일계 미국인 수학자 한스 라데마허, 그리고 미국인 수학자 조세프 월시의 이름을 따서 명명되었다.

정의

아다마르 변환 Hm은 2m × 2m 행렬인 아다마르 행렬(정규화 인자로 스케일링됨)로, 2m개의 실수 xn을 2m개의 실수 Xk로 변환한다. 아다마르 변환은 두 가지 방식으로 정의될 수 있다: 재귀적으로 정의하거나, 인덱스 n과 k의 이진 (기수-2) 표현을 사용하여 정의한다.

재귀적으로, 우리는 1 × 1 아다마르 변환 H0단위 H0 = 1로 정의하고, m > 0에 대해 Hm을 다음과 같이 정의한다: Hm=12(Hm1Hm1Hm1Hm1) 여기서 1/2는 때때로 생략되는 정규화 요소이다.

m > 1의 경우, Hm은 다음과 같이 정의될 수도 있다: Hm=H1Hm1 여기서 크로네커 곱을 나타낸다. 따라서, 이 정규화 인자를 제외하면, 아다마르 행렬은 완전히 1과 −1로 구성된다.

동등하게, 아다마르 행렬은 (k, n)-번째 항목을 다음과 같이 정의할 수 있다: k=i=0m1ki2i=km12m1+km22m2++k12+k0n=i=0m1ni2i=nm12m1+nm22m2++n12+n0

여기서 kj와 nj는 각각 k와 n의 비트 요소(0 또는 1)이다. 왼쪽 상단 모서리의 요소에 대해 k=n=0으로 정의한다. 이 경우, 다음이 성립한다: (Hm)k,n=12m/2(1)jkjnj

이는 입력과 출력을 각각 nj와 kj로 인덱싱된 다차원 배열로 간주할 때, 정확히 2×2××2×2의 다차원 DFT이며, 유니터리가 되도록 정규화되었다.

아다마르 행렬의 몇 가지 예시는 다음과 같다. H0=+(1)H1=12(1111)H2=12(1111111111111111)H3=123/2(1111111111111111111111111111111111111111111111111111111111111111)(Hn)i,j=12n/2(1)ij 여기서 ij는 i와 j의 이진 표현의 비트 단위 스칼라곱이다. 예를 들어, n2인 경우, (Hn)3,2=(1)32=(1)(1,1)(1,0)=(1)1+0=(1)1=1이 되며, 이는 위 결과와 일치한다(전체 상수 무시). 행렬의 첫 번째 행, 첫 번째 열 요소는 (Hn)0,0으로 표시된다.

H1은 정확히 크기 2의 DFT이다. 이는 또한 Z/(2)의 두 요소 덧셈 그룹에 대한 푸리에 변환으로 간주될 수 있다.

아다마르 행렬의 행은 월시 함수이다.

월시-아다마르 변환의 장점

실수

행렬 H의 위 정의에 따르면, 여기서 H = H[m,n]으로 둔다. H[m,n]=(1111)

월시 변환에서는 행렬에 1과 -1만 나타난다. 1과 -1은 실수이므로 복소수 계산을 수행할 필요가 없다.

곱셈이 필요 없음

DFT는 무리수 곱셈이 필요한 반면, 아다마르 변환은 그렇지 않다. 부호 반전만으로 충분하므로 유리수 곱셈도 필요하지 않다.

일부 속성은 DFT와 유사

월시 변환 행렬에서 각 행과 각 열은 월시 함수이며, 부호 변화의 수가 연속적으로 증가한다. 즉, 첫 번째 행과 첫 번째 열에는 부호 변화가 없다(모든 항목이 1과 같다). 두 번째 행에는 하나의 부호 변화가 있고, 세 번째 행에는 두 개의 부호 변화가 있으며, 이런 식으로 마지막 행까지 이어진다. H[m,n]=(1111111111111111111111111111111111111111111111111111111111111111)이를 이산 푸리에 변환 ej2πmn/N과 비교하면, 각 행 i에는 i1개의 영점 교차가 포함된다.

이산 푸리에 변환에서 m이 0과 같을 때(첫 번째 행에 해당), 결과도 1이다. 이후 행에 대해서는 행렬의 특성인 신호 주파수가 첫 번째 원본 행렬에서는 낮게 시작하여 다음 행으로 갈수록 증가하다가 마지막 행에서 최고조에 달하는 것을 관찰할 수 있다.

푸리에 변환과의 관계

아다마르 변환은 사실 2 × 2 × ⋯ × 2 × 2 크기의 다차원 DFT와 동일하다.[2]

또 다른 접근 방식은 아다마르 변환을 불 대수군 (/2)n에 대한 푸리에 변환으로 보는 것이다.[3] [4] 유한 (아벨) 군에 대한 푸리에 변환을 사용하여, 함수 f:(/2)n의 푸리에 변환은 다음과 같이 정의된 함수 f^이다. f^(χ)=a(/2)nf(a)χ¯(a) 여기서 χ(/2)n지표이다. 각 지표는 어떤 r(/2)n에 대해 χr(a)=(1)ar의 형태를 가지며, 여기서 곱셈은 비트 문자열에 대한 불 스칼라곱이다. 따라서 f^의 입력을 r(/2)n(폰트랴긴 쌍대성)와 동일시하고, f^:(/2)n를 다음과 같이 정의할 수 있다: f^(r)=a(/2)nf(a)(1)ra

이는 ff^의 입력을 불리언 문자열로 간주할 때 f의 아다마르 변환이다.

아다마르 변환이 2n개의 복소수 벡터 v에 왼쪽에서 아다마르 행렬 Hn을 곱하는 위의 공식화 측면에서, 이 동등성은 fv의 요소 인덱스에 해당하는 비트 문자열을 입력으로 받고, fv의 해당 요소를 출력하도록 취함으로써 확인된다.

이는 벡터 v에 적용될 때 순환군 /2n의 지표를 사용하는 일반적인 이산 푸리에 변환과 비교된다.

계산 복잡도

고전적인 영역에서, 아다마르 변환은 고속 아다마르 변환 알고리즘을 사용하여 nlogn 연산(n=2m)으로 계산할 수 있다.

양자 영역에서, 아다마르 변환은 양자 논리 게이트로서 병렬화될 수 있으므로 O(1) 시간에 계산할 수 있다.

양자 컴퓨팅 응용

아다마르 변환은 양자 컴퓨팅에서 광범위하게 사용된다. 2 × 2 아다마르 변환 H1은 아다마르 게이트로 알려진 양자 논리 게이트이며, n-큐비트 레지스터의 각 큐비트에 아다마르 게이트를 병렬로 적용하는 것은 아다마르 변환 Hn과 동일하다.

아다마르 게이트

양자 컴퓨팅에서 아다마르 게이트는 하나의 큐비트 회전으로, 큐비트-기저 상태 |0|1를 계산 기저 상태 |0|1의 동일한 가중치를 갖는 두 개의 양자 중첩 상태로 매핑한다. 일반적으로 위상은 다음과 같이 선택된다: H=|0+|120|+|0|121|

디랙 표기법에서. 이는 변환행렬 H1=12(1111) 에 해당한다. 이는 |0,|1 기저, 즉 계산 기저에서도 알려져 있다. 상태 |0+|12|0|12는 각각 |+|로 알려져 있으며, 함께 양자 컴퓨팅에서 극좌표 기저를 구성한다.

아다마르 게이트 연산

H(|0)=12|0+12|1=:|+H(|1)=12|012|1=:|H(|+)=H(12|0+12|1)=12(|0+|1)+12(|0|1)=|0H(|)=H(12|012|1)=12(|0+|1)12(|0|1)=|1

0 또는 1 큐비트에 아다마르 게이트를 한 번 적용하면 (처음 두 연산에서 볼 수 있듯이) 관측될 경우 0 또는 1이 동일한 확률로 나올 양자 상태가 생성된다. 이는 표준 확률적 계산 모델에서 공정한 동전을 던지는 것과 정확히 같다. 그러나 아다마르 게이트를 연속으로 두 번 적용하면 (마지막 두 연산에서 효과적으로 수행되는 것처럼), 최종 상태는 항상 초기 상태와 동일하다.

양자 알고리즘의 아다마르 변환

양자 아다마르 변환은 아다마르 변환의 텐서곱 구조 때문에 각 큐비트에 아다마르 게이트를 개별적으로 적용하는 것과 같다. 이 간단한 결과는 양자 아다마르 변환이 log2N 연산을 필요로 하며, 고전적인 경우 Nlog2N 연산과 비교된다는 것을 의미한다.

n-큐비트 시스템의 경우, 각 n 큐비트(각각 |0으로 초기화됨)에 작용하는 아다마르 게이트N=2n 형태일 때 균일한 양자 중첩 상태를 준비하는 데 사용될 수 있다. 이 경우 n 큐비트를 사용하면 결합된 아다마르 게이트 Hnn개의 아다마르 게이트의 텐서곱으로 표현된다: Hn=HHHn times

결과적인 균일 양자 중첩 상태는 다음과 같다: Hn|0n=12nj=02n1|j 이는 모든 N=2n에 대해 아다마르 게이트를 사용하여 균일 양자 상태를 준비하는 것을 일반화한다.[5]

이 균일 양자 상태의 측정|0|N1 사이의 무작위 상태를 초래한다.

많은 양자 알고리즘은 아다마르 변환을 초기 단계로 사용한다. 앞서 설명했듯이, 이는 |0으로 초기화된 n 큐비트를 |0,|1 기저의 모든 2n 직교 상태의 동일한 가중치를 갖는 중첩으로 매핑하기 때문이다. 예를 들어, 이는 도이치–조사 알고리즘, 사이먼의 알고리즘, 번스타인-바지라니 알고리즘, 그리고 그로버 알고리즘에서 사용된다. 쇼어 알고리즘은 초기 아다마르 변환과 양자 푸리에 변환을 모두 사용하는데, 둘 다 유한군에 대한 푸리에 변환의 일종이다. 첫 번째는 (/2)n에 대한 것이고 두 번째는 /2n에 대한 것이다.

일반적인 경우 N2n일 때 균일 양자 중첩 상태를 준비하는 것은 비자명하며 더 많은 작업이 필요하다. 균일 중첩 상태를 준비하기 위한 효율적이고 결정론적인 접근 방식 |Ψ=1Nj=0N1|j 는 모든 N에 대해 게이트 복잡도와 회로 깊이가 O(log2N)에 불과하며 최근에 발표되었다.[6] 이 접근 방식은 오직 n=log2N 큐비트만 필요하다. 중요한 점은 이 접근 방식에서 균일 중첩 상태 |Ψ를 생성하는 데 보조 큐비트나 다중 제어 양자 게이트가 필요하지 않다는 것이다.

분자 계통학(진화생물학) 응용

아다마르 변환은 분자 데이터를 사용하여 계통수를 추정하는 데 사용될 수 있다.[7][8][9] 계통학은 유기체 간의 관계를 이해하는 데 초점을 맞춘 진화생물학의 하위 분야이다. DNA 다중서열정렬에서 얻은 위치 패턴 빈도의 벡터(또는 행렬)에 적용된 아다마르 변환은 나무 위상에 대한 정보를 담는 또 다른 벡터를 생성하는 데 사용될 수 있다. 계통 발생학적 아다마르 변환의 가역적 특성은 나무 위상 벡터로부터 위치 가능성을 계산할 수 있게 하여, 아다마르 변환을 최대 가능도 추정에 사용할 수 있도록 한다. 그러나 후자의 응용은 위치 가능성을 계산하는 다른 훨씬 더 효율적인 방법이 있기 때문에 위치 패턴 벡터에서 나무 벡터로의 변환보다 덜 유용하다.[10][11] 그러나 계통 발생학적 아다마르 변환의 가역적 특성은 수학적 계통 발생학을 위한 우아한 도구를 제공한다.[12][13]

계통 발생학적 아다마르 변환의 메커니즘은 나무 T에 대한 위상 및 분지 길이 정보를 제공하는 벡터 γ(T)사이트 패턴 벡터 또는 행렬 s(T)를 사용하여 계산하는 것을 포함한다.

γ(T)=H1(ln(Hs(T))) 여기서 H는 적절한 크기의 아다마르 행렬이다. 이 방정식은 해석을 단순화하기 위해 세 개의 방정식 시리즈로 다시 작성될 수 있다:

r=Hs(T)ρ=lnrγ(T)=H1ρ

이 방정식의 가역적 특성은 다음과 같이 예상되는 사이트 패턴 벡터(또는 행렬)를 계산할 수 있게 한다:

s(T)=H1(exp(Hγ(T)))

DNA에 대한 Cavender–Farris–Neyman (CFN) 두 상태 치환 모델을 사용하여 뉴클레오타이드를 이진 문자(퓨린 A와 G는 R로, 피리미딘 C와 T는 Y로 인코딩)로 인코딩할 수 있다. 이렇게 하면 다중 서열 정렬을 위 사이트 패턴 벡터 s(T)로 인코딩할 수 있으며, 이는 다음 예시에서와 같이 나무 벡터 γ(T)로 변환될 수 있다:

특정 나무에 대한 아다마르 변환을 보여주는 예시 (작업 예시에 대한 값은 Waddell et al. 1997[14]에서 인용)
인덱스 이진 패턴 정렬 패턴 γ(T) ρ=H1γ(T) r=exp(ρ) s(T)=H1ρ
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

이 표에 제시된 예시는 간략화된 세 가지 방정식 체계를 사용하며, Newick 형식으로 ((A,B),(C,D))로 작성될 수 있는 4개 분류군 나무에 대한 것이다. 사이트 패턴은 ABCD 순서로 작성된다. 이 특정 나무는 두 개의 긴 말단 가지(사이트당 0.2개의 트랜스버전 치환), 두 개의 짧은 말단 가지(사이트당 0.025개의 트랜스버전 치환), 그리고 짧은 내부 가지(사이트당 0.025개의 트랜스버전 치환)를 가지고 있다. 따라서 Newick 형식으로는 ((A:0.025,B:0.2):0.025,(C:0.025,D:0.2))로 작성될 것이다. 이 나무는 데이터가 최대 단순성 기준으로 분석될 경우 긴 가지 끌림을 나타낼 것이다(분석된 서열이 관찰된 사이트 패턴 빈도가 s(T)=H1ρ 열에 표시된 예상 빈도에 근접할 만큼 충분히 길다고 가정). 긴 가지 끌림은 인덱스 6의 사이트 패턴(나무 ((A,C),(B,D))를 지지하는)의 예상 수가 참 나무(인덱스 4)를 지지하는 예상 사이트 패턴의 수를 초과한다는 사실을 반영한다. 분명히, 계통 발생학적 아다마르 변환의 가역적 특성은 나무 벡터 γ(T)가 정확한 나무에 해당한다는 것을 의미한다. 따라서 변환 후의 단순성 분석은 통계적으로 일치하며,[15] 올바른 모델(이 경우 CFN 모델)을 사용하는 표준 최대 우도 분석과 동일하다.

0 인덱스의 사이트 패턴은 (뉴클레오타이드를 퓨린 또는 피리미딘으로 인코딩한 후) 변경되지 않은 사이트에 해당한다. 별표가 있는 인덱스(3, 5, 6)는 "parsimony-informative"하며, 나머지 인덱스는 단일 분류군이 다른 세 분류군과 다른 사이트 패턴을 나타낸다(따라서 표준 최대 우도 계통수에서 말단 가지 길이와 동등하다).

RY 데이터 없이 뉴클레오타이드 데이터를 사용하려면(궁극적으로 0과 1로 재코딩), 사이트 패턴을 행렬로 인코딩할 수 있다. 4개 분류군 나무를 고려하면 총 256개의 사이트 패턴(4개의 뉴클레오타이드의 4제곱)이 존재한다. 그러나 기무라 3매개변수(또는 K81) 모델의 대칭성을 통해 DNA에 대한 256개의 가능한 사이트 패턴을 64개 패턴으로 줄일 수 있으며, 이로 인해 4개 분류군 나무에 대한 뉴클레오타이드 데이터를 8 × 8 행렬로 인코딩할 수 있다.[16] 이는 위에서 전이(RY) 사이트 패턴에 사용된 8개 요소 벡터와 유사하다. 이는 클라인 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) 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)의 예시[16]를 사용하여 설명할 수 있으며, 이는 네 가지 영장류 헤모글로빈 유사유전자의 다중 서열 정렬을 기반으로 한다:

인코딩된 서열 정렬 예시 (Hendy et al. 1994[16]에서 발췌) (값은 9879개 사이트 중 개수임)
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열이 전이 차이에 해당하기 때문이며, 이는 거의 모든 유전체 영역 비교에서 트랜스버전 차이보다 더 빠르게 축적된다(그리고 이 예시에 사용된 헤모글로빈 유사유전자에서는 확실히 더 빠르게 축적된다[17]). AAGG 사이트 패턴을 고려하면 클라인 군 비트 쌍의 두 번째 요소에 대해 이진 패턴 0000이 되고, 첫 번째 요소에 대해 0011이 된다. 이 경우 첫 번째 요소를 기반으로 한 이진 패턴은 인덱스 3(테이블에서 단일 별표로 표시)에 해당한다. GGAA, CCTT, TTCC 사이트 패턴도 정확히 같은 방식으로 인코딩된다. AACT 사이트 패턴은 두 번째 요소에 대해 이진 패턴 0011, 첫 번째 요소에 대해 0001로 인코딩된다. 이는 첫 번째 요소에 대해 인덱스 1, 두 번째 요소에 대해 인덱스 3을 생성한다. 두 번째 클라인 군 비트 쌍을 기반으로 한 인덱스는 8을 곱하여 열 인덱스(이 경우 24열)를 산출한다. AACT 사이트 패턴의 개수가 포함될 셀은 두 개의 별표로 표시되어 있다. 그러나 예시에서 숫자가 없는 것은 서열 정렬에 AACT 사이트 패턴이 없음을 나타낸다(마찬가지로 CCAG, GGTC, TTGA 사이트 패턴도 같은 방식으로 인코딩되며, 이들도 없음).

기타 응용

아다마르 변환은 데이터 암호화, JPEG XRMPEG-4 AVC와 같은 많은 신호 처리데이터 압축 알고리즘에서도 사용된다. 영상 압축 응용 프로그램에서는 일반적으로 변환된 차이의 절대값 합계 형태로 사용된다. 또한 양자 컴퓨팅에서 상당수의 알고리즘의 중요한 부분이다. 아다마르 변환은 핵자기 공명, 질량 분석법결정학과 같은 실험 기술에도 적용된다. 또한 일부 버전의 지역 민감 해싱에서도 의사 난수 행렬 회전을 얻기 위해 사용된다.

같이 보기

외부 링크

각주

  1. Compare Figure 1 in Townsend, W.J.; Thornton, M.A. (2001). 〈Walsh spectrum computations using Cayley graphs〉. 《Proceedings of the 44th IEEE 2001 Midwest Symposium on Circuits and Systems (MWSCAS 2001)》. MWSCAS-01. IEEE. 110–113쪽. doi:10.1109/mwscas.2001.986127. ISBN 0-7803-7150-X. 
  2. Kunz, H.O. (1979). 《On the Equivalence Between One-Dimensional Discrete Walsh–Hadamard and Multidimensional Discrete Fourier Transforms》. 《IEEE Transactions on Computers》 28. 267–8쪽. doi:10.1109/TC.1979.1675334. S2CID 206621901. 
  3. “Fourier Analysis of Boolean Maps– A Tutorial –, pp. 12–13” (PDF). 2021년 5월 1일에 원본 문서 (PDF)에서 보존된 문서. 2025년 7월 17일에 확인함. 
  4. Lecture 5: Basic quantum algorithms, Rajat Mittal, pp. 4–5
  5. Nielsen, Michael A.; Chuang, Isaac (2010). 《Quantum Computation and Quantum Information》. Cambridge: 케임브리지 대학교 출판부. ISBN 978-1-10700-217-3. OCLC 43641333. 
  6. Alok Shukla and Prakash Vedula (2024). 《An efficient quantum algorithm for preparation of uniform quantum superposition states》. 《Quantum Information Processing》. 23:38. 38쪽. arXiv:2306.11747. Bibcode:2024QuIP...23...38S. doi:10.1007/s11128-024-04258-4. 
  7. Hendy, Michael D.; Penny, David (December 1989). 《A Framework for the Quantitative Study of Evolutionary Trees》. 《Systematic Zoology》 38. 297쪽. doi:10.2307/2992396. JSTOR 2992396. 
  8. Hendy, Michael D.; Penny, David (January 1993). 《Spectral analysis of phylogenetic data》 (영어). 《Journal of Classification》 10. 5–24쪽. doi:10.1007/BF02638451. ISSN 0176-4268. S2CID 122466038. 
  9. Székely, L. A., Erdős, P. L., Steel, M. A., & Penny, D. (1993). A Fourier inversion formula for evolutionary trees. Applied mathematics letters, 6(2), 13–16.
  10. Felsenstein, Joseph (November 1981). 《Evolutionary trees from DNA sequences: A maximum likelihood approach》 (영어). 《Journal of Molecular Evolution》 17. 368–376쪽. Bibcode:1981JMolE..17..368F. doi:10.1007/BF01734359. ISSN 0022-2844. PMID 7288891. S2CID 8024924. 
  11. Stamatakis, Alexandros (2019), Warnow, Tandy (편집), “A Review of Approaches for Optimizing Phylogenetic Likelihood Calculations” (영어), 《Bioinformatics and Phylogenetics》, Computational Biology (Cham: Springer International Publishing) 29, 1–19쪽, doi:10.1007/978-3-030-10836-6, ISBN 978-3-030-10836-6, S2CID 145834168, 2020년 10월 10일에 확인함 
  12. Chor, Benny; Hendy, Michael D.; Holland, Barbara R.; Penny, David (2000년 10월 1일). 《Multiple Maxima of Likelihood in Phylogenetic Trees: An Analytic Approach》 (영어). 《Molecular Biology and Evolution》 17. 1529–1541쪽. doi:10.1093/oxfordjournals.molbev.a026252. ISSN 1537-1719. PMID 11018159. 
  13. Matsen, Frederick A.; Steel, Mike (2007년 10월 1일). Ané, Cécile; Sullivan, Jack (편집). 《Phylogenetic Mixtures on a Single Tree Can Mimic a Tree of Another Topology》 (영어). 《Systematic Biology》 56. 767–775쪽. arXiv:0704.2260. doi:10.1080/10635150701627304. ISSN 1076-836X. PMID 17886146. 
  14. Waddell, Peter J; Steel, M.A (December 1997). 《General Time-Reversible Distances with Unequal Rates across Sites: Mixing Γ and Inverse Gaussian Distributions with Invariant Sites》 (영어). 《Molecular Phylogenetics and Evolution》 8. 398–414쪽. Bibcode:1997MolPE...8..398W. doi:10.1006/mpev.1997.0452. PMID 9417897. 
  15. Steel, M. A.; Hendy, M. D.; Penny, D. (1993년 12월 1일). 《Parsimony Can Be Consistent!》 (영어). 《Systematic Biology》 42. 581–587쪽. doi:10.1093/sysbio/42.4.581. ISSN 1063-5157. 
  16. Hendy, M. D.; Penny, D.; Steel, M. A. (1994년 4월 12일). 《A discrete Fourier analysis for evolutionary trees.》 (영어). 《Proceedings of the National Academy of Sciences》 91. 3339–3343쪽. Bibcode:1994PNAS...91.3339H. doi:10.1073/pnas.91.8.3339. ISSN 0027-8424. PMC 43572. PMID 8159749. 
  17. Miyamoto, M. M.; Koop, B. F.; Slightom, J. L.; Goodman, M.; Tennant, M. R. (1988년 10월 1일). 《Molecular systematics of higher primates: genealogical relations and classification.》 (영어). 《Proceedings of the National Academy of Sciences》 85. 7627–7631쪽. Bibcode:1988PNAS...85.7627M. doi:10.1073/pnas.85.20.7627. ISSN 0027-8424. PMC 282245. PMID 3174657.