?? test_joint_diag.m
字號:
% Calling the joint approximate diagonalization function.m=5 % dimensionn=3 % number of matricesseuil = 1.0e-12; % precision on joint diag compteur=0;while 1 ; compteur=compteur+1;% drawing a `random' unitary matrixU= randn(m)+i*randn(m) ; [U,to_waste]=eig(U+U'); % Drawing a random set of commuting matricesA=zeros(m,m*n);for imat=1:n cols = 1+(imat-1)*m:imat*m; A(:,cols) = U*diag(randn(m,1)+i*randn(m,1))*U';end;% Perturbation of the joint structure ?% A = A + 0.001*randn(m,m*n);%% Do it[ V , DD ] = joint_diag(A,seuil);%% should be permutation matrix abs(V'*U)end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -