Броуновское движение для КУА
Броуновское движение может быть создано при помощи инфинитезимальной образующей. Броуновское движение на множестве унитарных матриц может быть создано инфинитезимальными образующими из гауссова унитарного ансамбля. Стационарное распределение этого броуновского движения тогда становится идентичным распределению КУА.
In[1]:=
![Click for copyable input](assets.ru/brownian-motion-on-cue/In_56.png)
mats = RandomVariate[GaussianUnitaryMatrixDistribution[0.1, 2],
100000];
mats = Table[MatrixExp[I mat], {mat, mats}];
Сгенерируйте броуновскую траекторию с начальными точками, выбранными в КУА.
In[2]:=
![Click for copyable input](assets.ru/brownian-motion-on-cue/In_57.png)
initial = RandomVariate[CircularUnitaryMatrixDistribution[2]];
res = FoldList[#2.#1 &, initial, mats];
Рассчитайте фазы собственных значений и сравните их с плотностью вероятности собственных значений матриц из КУА.
In[3]:=
![Click for copyable input](assets.ru/brownian-motion-on-cue/In_58.png)
phases = RandomSample /@ Arg[Eigenvalues /@ res];
код на языке Wolfram Language целиком
Out[4]=
![](assets.ru/brownian-motion-on-cue/O_27.png)