function outprod.m

Question

NOR AZIERA 2022-12-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1878217-function-outprod-m

编辑： Matt J 2022-12-14

Cannot run this coding. Please help me

function [X]=outprod(V,D)
    n=3;
    A=randn(n); %random square matrix
    A=A'*A;
    [V,D]=eig(A);%find eigenvalue and eigenvector
    X=zeros(n); %define the matrix
    for i=1:n
        X=X+D(i,i)*outprod(V(:,i),V(:,i));
        display(X);
        display(A');  
    end
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Jan 2022-12-14

移动：Matt J 2022-12-14

You forgot to mention, what the problem is. The error message should mention this already.

This function overwrites the inputs V and D and calls itself recursively until the recursion limit is exhausted. The readers cannot guess reliably, what the intention of the function is. Therefore it is fragile to suggest a modification of the code.

Please explain, what you want to calculate. In future questions attach a copy of the complete error message.

请先登录，再进行评论。

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

Answer 1

Matt J 2022-12-14

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1878217-function-outprod-m#answer_1127842

编辑：Matt J 2022-12-14

在 MATLAB Online 中打开

Perhaps this is what you want?

    n=3;
    A=randn(n); %random square matrix
    A=A'*A;
    [V,D]=eig(A);%find eigenvalue and eigenvector
    
    X=outprod(V,D);
        display(X);
X = 3×3
    0.2455    0.1739    0.6999
    0.1739    0.8567    1.4576
    0.6999    1.4576    3.4055
        display(A');    
    0.2455    0.1739    0.6999
    0.1739    0.8567    1.4576
    0.6999    1.4576    3.4055
        
function [X]=outprod(V,D)
   n=length(D);
    X=0; %define the matrix
    for i=1:n
        X=X+D(i,i)*V(:,i)*V(:,i)';
    end
end