단일 원소 교란 시 행렬식·영구량 최소 변화 확률

초록

본 논문은 무작위로 선택된 행렬 원소 중 하나를 교란했을 때 행렬식(또는 영구량)의 변화가 최소가 되는 경우의 확률을 구한다. 원소는 비영(0이 아닌 값)일 확률이 고정된 연속형 혹은 이산형 확률변수로 모델링한다. 연속형과 이산형 경우에 따라 최소 변화 조건이 달라지며, 이를 코팩터와 영구량의 구조적 성질을 이용해 분석한다.

상세 분석

논문은 먼저 n×n 행렬 A의 원소 a_{ij}를 독립적인 확률변수로 두고, 각 원소가 0이 될 확률을 p, 비0이 될 확률을 1‑p 로 설정한다. 여기서 p는 고정된 실수이며, 연속형 모델에서는 a_{ij}가 연속 확률밀도함수 f(x)·(1‑p) 로, 이산형 모델에서는 a_{ij}가 유한 집합 {0, c_1,…,c_k} 중 하나를 선택하는 확률질량함수로 정의된다.

행렬식 det(A) 를 전개하면 det(A)=∑{σ∈S_n}sgn(σ)∏{i=1}^n a_{iσ(i)} 로 표현된다. 특정 원소 a_{pq} 를 ε만큼 변동시켰을 때 det(A) 의 변화량 Δ는 a_{pq} 가 포함된 모든 전개 항에만 영향을 미친다. 즉 Δ = ε·C_{pq}, 여기서 C_{pq}는 (p,q) 위치의 코팩터이다. 따라서 det(A) 가 최소 변화(즉 Δ=0) 를 보이려면 C_{pq}=0 이어야 한다.

연속형 경우에는 a_{pq} 가 연속적으로 변할 수 있기 때문에, C_{pq}=0 일 확률을 직접 계산한다. C_{pq}=0 은 (n‑1)×(n‑1) 소행렬의 행렬식이 0 인 경우와 동치이며, 이는 해당 소행렬이 특이행렬일 확률과 같다. 독립성 가정 하에, 각 소행렬의 원소가 0이 될 확률 p 로 동일하므로, 소행렬이 특이일 확률은 p와 (1‑p)의 조합으로 표현되는 다항식 형태가 된다. 논문은 이를 전이 행렬과 그래프 이론을 이용해 정확히 구하고, 특히 p→0 (희소 행렬) 일 때 특이 확률이 O(p) 로 선형적으로 감소함을 보인다.

이산형 경우에는 a_{pq} 가 0 혹은 고정된 비0값 c_k 로만 취할 수 있다. 여기서 Δ=0 은 두 경우로 나뉜다. 첫째, a_{pq}=0 이면서 C_{pq}가 어떤 값이든 상관없다(이미 교란 전후 차이가 0). 둘째, a_{pq}=c_k 이고 C_{pq}=0 인 경우이다. 따라서 전체 최소 변화 확률은 P(a_{pq}=0) + P(a_{pq}=c_k)·P(C_{pq}=0) 로 분해된다. C_{pq}=0 의 확률은 연속형과 동일하게 소행렬이 영구량(permanent) 혹은 행렬식이 0 인 경우와 연결되지만, 이산형에서는 영구량이 0 일 확률이 행렬식이 0 일 확률보다 더 복잡하게 나타난다. 특히 영구량은 부호가 없으므로, 영구량이 0 일 확률은 소행렬에 적어도 하나의 행이 전부 0이 되는 경우와 동치이며, 이는 조합론적 계산을 통해 p^{n‑1} 형태로 간단히 표현된다.

논문은 또한 영구량에 대한 최소 변화 문제를 다룬다. 영구량 per(A) 를 전개하면 부호가 없으므로, a_{pq} 의 교란이 영구량에 미치는 영향은 Δ = ε·P_{pq}, 여기서 P_{pq}는 (p,q) 위치의 영구 코팩터이다. 연속형에서는 P_{pq}=0 일 확률이 거의 없으므로, 최소 변화는 a_{pq}=0 인 경우에만 발생한다. 반면 이산형에서는 P_{pq}=0 일 확률이 p^{n‑1} 로 계산되며, 이는 행이 전부 0이 되는 경우와 일치한다.

핵심 통찰은 “코팩터가 0”이라는 조건이 연속형과 이산형에서 서로 다른 확률 구조를 만든다는 점이다. 연속형에서는 코팩터가 연속적인 함수이므로 0이 될 확률이 0에 수렴하지만, 이산형에서는 0이 되는 경우가 명시적으로 존재한다. 따라서 최소 변화 확률을 구하는 절차는 연속형에서는 특이성 확률을 직접 적분하거나 전이 행렬을 이용한 재귀식으로 계산하고, 이산형에서는 조합론적 사건의 확률을 곱셈 법칙으로 구하는 두 단계로 나뉜다.

마지막으로 논문은 이러한 결과를 이용해 무작위 행렬의 안정성 분석, 신호 처리에서의 잡음 민감도 평가, 그리고 양자 컴퓨팅에서의 오류 전파 모델링 등에 응용 가능함을 제시한다.