我正在尝试获取一个计算二分状态(两个密度矩阵的张量)的施密特分解的函数,以并排打印施密特数(X)和基向量(e_i,e_j)的总和。 for 循环。例如:(x'*'e_1'\otimes' e_1+ ... + x'*'e_i'\otimes' e_j)。注释掉的 fprintf 语句是我已经尝试过的东西。
代码:
function [u, v] = Schmidt(m,n,v)
m = 3 %size of basis vector
n = 4 %size of basis vector
e_m = eye(m)
e_n = eye(n)
%can have a = input('enter an array>');
v = [1
2
3
4
5
6
7
8
9
10
11
12]
for i = 1:m
for j = 1:n
e_i = e_m(:,i)
e_j = e_n(:,j)
K = kron(e_i, e_j)
x = ctranspose(K)*v
W(i,j) = x
%fprintf('%d %f %f %f\n',j,e_i,e_j,x);
%fprintf(1,'Header 1\tHeader 2\n');
%fprintf(1,'%f\t%f\n','e_i e_j');
%fprintf('the basis are %d\n',e_i, e_j)
end
end
end