ベクトルの内積の分散の単位ベクトルのすべての分布の最小値はいくつですか?


10

nx1,,xnkn>kmaxijVar(xiTxj)E[xiTxj]=0

私はいくつかの分布を試しましたが、ほとんどすべての分散はです。例えば、それぞれの各座標た分布の両方独立かつ一様から選択される及び分布の各は、分散持つ次元の単位球上の独立した均一ベクトルです。1/kxi{1/k,1/k}xik1/k

されてすべてのディストリビューションの中で最小の分散を?1/k


どれだけきつい限界に興味がありますか?つまり、n> 100kでのみ機能する1 / 100kの下限は興味深いものでしょうか?
daniello 2014年

@ daniello、n> ckの1 / ckの下限、つまりcは定数です。これを証明するには?
2014年

質問で理解できないこと:最初に、単位ベクトルの分布と言いますが、すべての分布で単位ベクトルを生成しようとしたわけではありません...すべての、? E [ | x i | ] = 1xiE[|xi|]=1
daniello 2014年

@deniello、私はすべてのベクトルを「ユニット」にするつもりでした。申し訳ありませんが、「ガウス」ベクトルで正規化を行うのを忘れてしまいました。正規化後は、均一なベクトルと同じになります。この誤りを指摘していただきありがとうございます。
peng

回答:


8

同等だが簡単に見える問題の定式化を提示し、(n / k − 1)/(n −1)の下限を示します。 また、量子情報の未解決問題との関連も示します。 [リビジョン3の編集:以前のリビジョンでは、複雑なケースの類似の質問にはSIC-POVMに関する未解決の問題が含まれているため、以下に示す下限が達成されるケースの正確な特徴付けは難しいと主張しました量子情報。ただし、SIC-POVMへのこの接続は正しくありませんでした。詳細については、以下の「量子情報のSIC-POVMへの誤った接続」のセクションを参照してください。]

同等の定式化

最初に、danielloの回答ですでに指摘したように、Var(x i T x j)= E [(x i T x j2 ] − E [ x i T x j ] 2 = E [(x i T x j2 ]。したがって、残りの答えでは、分散を忘れて、代わりにmax ij E [(x i T x j2 ]を最小化します。

次に、目標がmax ij E [(x i T x j2 ] を最小化することであると決定したら、E [ x i T x j ] = 0 の制約を無視できます。これは、単位ベクトルx 1、…、x nの場合、目的関数max ij E [(の値を変更せずにE [ x i T x j ] = 0 を満たすために、確率1/2でそれぞれを個別に否定するx i T x j2 ]。

さらに、目的関数をmax ij E [(x i T x j2 ]から(1 /(nn -1)))∑ ij E [(x i T x j2 ]に変更する最適値は変更されません。平均が最大で最大であるため、後者は最大で前者です。ただし、E((x i T x j2 ] の値は、(ij)(ijn個のベクトルx 1、…、x nをランダムに並べ替えることで等しくなります。

したがって、任意のnおよびkについて、問題の問題の最適値は(1 /(nn -1)))∑ ij E [(x i T x j2 ] の最小値に等しく、ここでx 1、…、x nは、variables kの単位ベクトルを値として取る確率変数です。

ただし、期待の線形性により、この目的関数は期待値E [(1 /(nn −1)))∑ ijx i T x j2 ]に等しくなります。最小値は最大でも平均であるため、確率分布を考慮する必要はありません。つまり、上記の問題の最適値は、次の最適値と等しくなります。

単位ベクトルを選択し、xが1、...、X N ∈ℝのkは最小にするが(1 /(NN -1)))Σ IJX I T X J2

下限

この同等の定式化を使用して、最適値が少なくとも(n / k − 1)/(n −1)であることを証明します。

1≤ため私はN、聞かせてX I = X I、X I Tは単位ベクトルに対応するランク1プロジェクタであるxはIを。次に、(x i T x j2 = Tr(X i X j)が成り立つ。

してみましょうYI X I。次に、∑ ij Tr(X i X j)= ∑ ij Tr(X i X j-n = Tr(Y 2-nであることが成り立つ。

コーシー=シュワルツの不等式は暗示そのTrの(Y 2)≥(TR Y2 / K = N 2 / K、したがってΣ IJ Trを(X I X J)= Trと(Y 2) - NN 2 / knnn −1)で割ると、目的の値は少なくとも(n / k − 1)/(n −1)になります。

特に、n = k +1の場合、danielloの答えは最適値の2倍以内です。

この下限はいつ達成できますか?

この下限(n / k − 1)/(n −1)を達成することは、Y =(n / kIを作成することと同じです。正確な特性が得られる時期はわかりませんが、次の十分な条件があります。

  • n = k +1の場合、原点を中心とする通常のkシンプレックスを形成するk +1単位ベクトルを考慮し、danielloの回答の2 /(kk +1))から最適な1 / kに改善することで、これを達成できます。
  • 場合、Nの倍数であるK、それは明らかℝ用の正規直交基底固定することにより達成可能であるKをとする基底ベクトルの各々を割り当てるN / KV 1、...、V N
  • 最後の箇条書きよりも一般的には、kを選択してn = n 1n = n 2の両方で達成できる場合、同じkn = n 1 + n 2でも達成できます。特に、もし達成され、N = K + BBが満たす整数である≥ Bを ≥0。

詳細は確認していませんが、球形の2設計はこの下限を達成するソリューションを提供するようです。

量子情報のSIC-POVMへの誤った接続

以前の改訂では、次のように述べています。

これに完全に答えるのは難しい質問だと思います。その理由は、私たちが代わりに複素ベクトル空間ℂを考慮するならばということであるkは、この質問は、量子情報で開いている問題に関連しています。

しかし、この関係は正しくありませんでした。その理由を説明します。

より正確には、次の問題を検討してください。

単位ベクトルを選択し、X 1、···、X N ∈ℂのkは(1 /(最小限に抑えるためのnnは -1)))Σ IJ | x i * x j |

上記の下限は、この複雑なバージョンにも等しく当てはまります。複合バージョンでn = k 2の場合を考えます。次に、下限は1 /(k +1)に等しくなります。

これまでのところ、それは正しかった。

一連のK 2単位ベクトルX 1、···、X K 2 ∈ℂのKと呼ばれる下限達成SIC-POVM次元でkは

この部分は間違っていました。A SIC-POVMでの集合であるK 2単位ベクトルX 1、···、X N ∈ℂのK |用 x i * x j | 2 = 1 /(k +1)すべてのijに対して。ここでの要件は、すべてのペアのために保持しなければならないことに注意Ijのではなく、すべてのペアを超えるだけの平均IJ。「等価定式化」セクションでは、最大値の最小化と平均値の最小化の間の等価性を示しましたが、これはx 1、…、x nは、そこに単位ベクトルをとる確率変数でした。ここでx 1、…、x nは単なる単位ベクトルであるため、同じトリックを使用することはできません。


5

「いいえ、分散はより小さくなる可能性があります」という質問への迅速で汚い答えは、標準的な基準とし、次のランダムプロセスを考慮します。異なる整数i、jのペアを選択、およびを設定します。他のベクトル()の場合、に1対1で割り当てます。次に、ごとに、をそのままにするか、確率で()フリップします。v1,v2,,vk{1,2,,k+1}xi=xj=v1xtt{i,j}v2,,vkt{1,,k+1}xixi12

参照することは容易である -いずれかと直交している、またはそれらが確率で同一の/逆の方向を指すそれぞれ。、X A 、X B 1E[xaxb]=0xaxb12

一方、ます。我々はそれを持っている場合にのみ確率で発生。それ以外の場合。したがって、我々はすべてのためにそれを持っている、: X AX B 2 = 1 { B } = { I J }Var[xaxb]=E[(xaxb)2](xaxb)2=1{a,b}={i,j}1(k+12)B V R [ X AX B ] = E [ X AX B 2 ] = 1(xaxb)2=0ab

Var[xaxb]=E[(xaxb)2]=1(k+12)

私の直感では、これはかなり悪い(小さい)ものですが、証明はありません。さらに興味深いのは、この構成がn >> kで故障するように見えることと、独立して(おそらく異なる分布から)選択する必要がある場合です。xi

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.