Tối thiểu trên tất cả các phân phối của vectơ đơn vị phương sai của sản phẩm chấm của vectơ là gì?


10

Tôi cố gắng để tìm thấy một bản phân phối qua vectơ ngẫu nhiên, nói , trên cầu đơn vị chiều (nơi ) giảm thiểu tuân theo ràng buộc .nx1,Giáo dục,xnkn>ktối đaTôijVmộtr(xTôiTxj)E[xTôiTxj]= =0

Tôi đã thử một số bản phân phối và hầu hết tất cả chúng đều có phương sai 1/k . Ví dụ: cả phân phối trong đó mỗi tọa độ của mỗi xTôi đều được chọn độc lập và thống nhất từ {-1/k,1/k} và phân phối trong đó mỗi xTôi là một vectơ đồng nhất độc lập trên mặt cầu đơn vị k -chiều có phương sai 1/k .

1/k phương sai tối thiểu trong số tất cả các phân phối?


Làm thế nào chặt chẽ của một ràng buộc mà bạn quan tâm? Đó là, liệu giới hạn dưới 1 / 100k chỉ hoạt động với n> 100k có thú vị hay không?
daniello

@daniello, ý bạn là giới hạn dưới của 1 / ck cho n> ck trong đó c là hằng số? Làm thế nào để chứng minh điều này?
peng

Một cái gì đó tôi không hiểu trong câu hỏi: ban đầu bạn nói phân phối trên các vectơ đơn vị , nhưng không phải tất cả các phân phối mà bạn nói bạn đã thử tạo vectơ đơn vị ... Bạn có nghĩa là với tất cả , ? E [ | x i | ] = 1xiE[|xi|]=1
daniello

@deniello, tôi dự định biến tất cả các vectơ thành "đơn vị" .. Xin lỗi, tôi đã quên thực hiện chuẩn hóa trên vectơ "Gaussian", sau khi chuẩn hóa, nó sẽ giống như vectơ đồng nhất. Cảm ơn bạn đã chỉ ra sai lầm này.
peng

Câu trả lời:


8

Tôi sẽ trình bày một công thức tương đương nhưng có vẻ đơn giản hơn của vấn đề và hiển thị giới hạn dưới của ( n / k - 1) / ( n 1). Tôi cũng cho thấy một kết nối với một vấn đề mở trong thông tin lượng tử. [Chỉnh sửa trong phiên bản 3: Trong các phiên bản trước, tôi đã tuyên bố rằng việc mô tả chính xác các trường hợp mà giới hạn dưới hiển thị dưới đây có thể gặp khó khăn vì một câu hỏi tương tự trong trường hợp phức tạp bao gồm một vấn đề mở về SIC-POVMs thông tin lượng tử. Tuy nhiên, kết nối này với SIC-POVM không chính xác. Để biết chi tiết, hãy xem phần Kết nối không chính xác với SIC-POVM trong thông tin lượng tử mà bên dưới.]

Công thức tương đương

Đầu tiên, như đã được chỉ ra trong câu trả lời của daniello, lưu ý rằng Var ( x i T x j ) = E [( x i T x j ) 2 ] - E [ x i T x j ] 2 = E [( x i T x j ) 2 ]. Vì vậy, trong phần còn lại của câu trả lời, chúng ta quên đi phương sai và thay vào đó tối thiểu hóa tối đa ij E [( x i T x j ) 2 ].

Tiếp theo, khi chúng tôi quyết định mục tiêu của chúng tôi là để giảm thiểu tối đa ij E [( x i T x j ) 2 ], chúng ta có thể bỏ qua những hạn chế đó E [ x i T x j ] = 0. Điều này là bởi vì nếu chúng ta có ngẫu nhiên vectơ đơn vị x 1 , chụp , x n , sau đó chúng ta có thể phủ định từng độc lập với xác suất 1/2 để thỏa mãn E [ x i T x j ] = 0 mà không thay đổi giá trị của hàm mục tiêu max ij E [( x i T x j) 2 ].

Hơn nữa, thay đổi hàm mục tiêu từ max ij E [( x i T x j ) 2 ] thành (1 / ( n ( n −1))) ∑ ij E [( x i T x j ) 2 ] không thay đổi giá trị tối ưu. Cái sau nhiều nhất là cái trước vì trung bình nhiều nhất là tối đa. Tuy nhiên, chúng tôi luôn có thể làm cho các giá trị của E [( x i T x j ) 2 ] cho sự lựa chọn khác nhau của ( i , j ) ( ij ) bằng cách hoán vị n vectơ x 1 , Mạnh , x n ngẫu nhiên.

Vì vậy, đối với bất kỳ nk , giá trị tối ưu của vấn đề được đề cập bằng giá trị tối thiểu của (1 / ( n ( n −1))) ∑ ij E [( x i T x j ) 2 ] trong đó x 1 , Số, x n là các biến ngẫu nhiên lấy vectơ đơn vị tính bằng ℝ k làm giá trị.

Tuy nhiên, theo tính tuyến tính của kỳ vọng, hàm mục tiêu này bằng giá trị mong đợi E [(1 / ( n ( n −1))) ∑ ij ( x i T x j ) 2 ]. Bởi vì mức tối thiểu tối đa là trung bình, không cần phải xem xét phân phối xác suất nữa. Nghĩa là, giá trị tối ưu của vấn đề trên bằng với giá trị tối ưu của các vấn đề sau:

Chọn vectơ đơn vị x 1 , ..., x n ∈ ℝ k để giảm thiểu (1 / ( n ( n -1))) Σ ij ( x i T x j ) 2 .

Chặn dưới

Sử dụng công thức tương đương này, chúng tôi sẽ chứng minh rằng giá trị tối ưu ít nhất là ( n / k - 1) / ( n 1).

Với 1≤ in , đặt X i = x i x i T là máy chiếu hạng 1 tương ứng với vectơ đơn vị x i . Sau đó, nó giữ rằng ( x i T x j ) 2 = Tr ( X i X j ).

Đặt Y = ∑ i X i . Khi đó, nó giữ rằng ∑ ij Tr ( X i X j ) = ∑ i , j Tr ( X i X j ) - n = Tr ( Y 2 ) - n .

Bất đẳng thức Cauchy từ Schwarz ngụ ý rằng Tr ( Y 2 ) (Tr Y ) 2 / k = n 2 / k , và do đó ∑ ij Tr ( X i X j ) = Tr ( Y 2 ) - nn 2 / k - n . Bằng cách chia cho n ( n 1), chúng ta thu được giá trị mục tiêu ít nhất là ( n / k - 1) / ( n 1).

Cụ thể, khi n = k +1, câu trả lời của daniello nằm trong hệ số 2 từ giá trị tối ưu.

Khi nào giới hạn dưới này có thể đạt được?

Đạt này thấp hơn ràng buộc ( n / k - 1) / ( n -1) tương đương với làm Y = ( n / k ) tôi . Tôi không biết đặc tính chính xác khi có thể đạt được, nhưng theo các điều kiện đủ tồn tại:

  • Khi n = k +1, có thể đạt được bằng cách xem xét các vectơ đơn vị k +1 tạo thành một k -simplex thông thường ở trung tâm gốc, cải thiện từ 2 / ( k ( k +1)) trong câu trả lời của daniello thành 1 / k tối ưu 2 .
  • Khi n là bội số của k , có thể đạt được rõ ràng bằng cách sửa một cơ sở trực giao của ℝ k và gán từng vectơ cơ sở cho n / k của v 1 , Lỗi , v n .
  • Tổng quát hơn điểm đạn cuối cùng, nếu có thể đạt được với một số lựa chọn k và cả n = n 1n = n 2 , thì nó cũng có thể đạt được cho cùng kn = n 1 + n 2 . Cụ thể, có thể đạt được nếu n = a k + b trong đó ab là các số nguyên thỏa mãn ab 0 .

Mặc dù tôi chưa kiểm tra chi tiết, nhưng dường như bất kỳ thiết kế 2 hình cầu nào cũng đưa ra giải pháp đạt được giới hạn dưới này.

Kết nối không chính xác với SIC-POVM trong thông tin lượng tử

Trong các phiên bản trước, tôi đã tuyên bố:

Tôi nghi ngờ rằng trả lời này hoàn toàn là một câu hỏi khó. Lý do là nếu thay vào đó chúng ta xem xét không gian vectơ phức tạp ℂ k , câu hỏi này có liên quan đến một vấn đề mở trong thông tin lượng tử.

Nhưng mối quan hệ này là không chính xác. Tôi sẽ giải thích tại sao.

Chính xác hơn, hãy xem xét vấn đề sau:

Chọn vectơ đơn vị x 1 , ..., x n ∈ ℂ k để giảm thiểu (1 / ( n ( n -1))) Σ ij | x i * x j | 2 .

Các giới hạn dưới ở trên giữ bằng nhau trong phiên bản phức tạp này. Hãy xem xét trường hợp n = k 2 trong phiên bản phức tạp. Khi đó giới hạn dưới bằng 1 / ( k +1).

Cho đến nay, nó đã đúng.

Một tập hợp các k 2 đơn vị vectơ x 1 , ..., x k 2 ∈ ℂ k đạt các ràng buộc thấp hơn được gọi là một SIC-POVM trong chiều k ,

Phần này không chính xác. Một SIC-POVM là một tập hợp các vectơ đơn vị k 2 x 1 , Khác, x nk mà | x i * x j | 2 = 1 / ( k +1) với mọi ij . Lưu ý rằng ở đây yêu cầu phải giữ cho tất cả các cặp ij , không chỉ trung bình trên tất cả các cặp ij . Trong phần Công thức tương đương của Fem, chúng tôi đã chỉ ra sự tương đương giữa giảm thiểu tối đa và tối thiểu hóa trung bình, nhưng điều này là có thể bởi vì x 1, Mạnh, x n là các biến ngẫu nhiên lấy các vectơ đơn vị ở đó. Ở đây x 1 , chụp , x n chỉ là các vectơ đơn vị, vì vậy chúng ta không thể sử dụng cùng một mẹo.


5

v1,v2,Giáo dục,vk{1,2,Giáo dục,k+1}xTôi= =xj= =v1xtt{Tôi,j}v2,Giáo dục,vkt{1,Giáo dục,k+1}xTôi-xTôi12

E[xmộtxb]= =0xmộtxb12

Vmộtr[xmộtxb]= =E[(xmộtxb)2](xmộtxb)2= =1{một,b}= ={Tôi,j}1(k+12)(xmộtxb)2= =0mộtb

Vmộtr[xmộtxb]= =E[(xmộtxb)2]= =1(k+12)

xTôi

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.