こちらの記事を読む library(BosonSampling) library(complexplus) m <- 2 # size of matrix (m x m) U <- randomUnitary(m) V <- randomUnitary(m) # 複素正方行列 # BosonSampling::randomUnitary() の中身をなぞる M <- matrix(complex(real = rnorm(m^2)…

１次元２状態量子ウォークを定める２ｘ２複素行列は行列式が±１のユニタリ行列 (行列式=1の場合) (行列式=-1の場合) #' 2x2 Unitary Matrix with determinant +/- 1 #' #' 2x2 Unitary Matrix #' @param theta, a real value. #' @param phi1, a real value.…