The Code Does not Compute the other elements of Matrix

Question

Gözde Üstün 2020-7-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/564185-the-code-does-not-compute-the-other-elements-of-matrix

评论： Gözde Üstün 2020-7-13

Hello,

I am trying to compute value of probabilitiy distribution The formula is ok, code is working I have no error but some values are coming as 0. However I should have a value except zero

Here is my code:

function P = Prob(A_ax,B_by,rho)
d=size(A_ax,1);
m=size(A_ax,4);
bcount=size(A_ax,3);
P = zeros(d,d,bcount,m);
for x=1:bcount %m->bcount
    for y=1:bcount %m->bcount
        for a=1:m %d
            for b=1:m %d
                 for dim=1:d
                    for dim2=1:d
                       %P(dim,dim2,x,y) = real( trace(kron(A_ax(dim,dim2,x,a),B_by(dim,dim2,y,b))*rho));
                          P(a,b,x,y) = real( trace(kron(A_ax(:,:,x,a),B_by(:,:,y,b))*rho));
                     end
                 end
            end
        end
    end
end

As you can see the other values are 0. Where am I doing wrong?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

KSSV 2020-7-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/564185-the-code-does-not-compute-the-other-elements-of-matrix#answer_465089

编辑：KSSV 2020-7-13

在 MATLAB Online 中打开

function P = Prob(A_ax,B_by,rho)
d=size(A_ax,1);
m=size(A_ax,4);
bcount=size(A_ax,3);
P = zeros(d,d,bcount,bcount);   %% Changed here 
for x=1:bcount %m->bcount
    for y=1:bcount %m->bcount
        for a=1:d
            for b=1:d                                                                          
                          P(a,b,x,y) = real( trace(kron(A_ax(:,:,x,a),B_by(:,:,y,b))*rho));                                   
            end
        end
    end
end