How to add values to already existing ones in a matrix using for loop?

Question

Vitaly 2013-4-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73016-how-to-add-values-to-already-existing-ones-in-a-matrix-using-for-loop

Hi, my function looks like this:

nEN=[1 2 4 5];
KG=[];
for i1=1:4
KG((2*(nEN(i1))-1),(2*(nEN(i1))-1))=Ke((2*i1-1),(2*i1-1));
KG(2*(nEN(i1)),2*(nEN(i1)))=Ke((2*i1),(2*i1));
end

Basically it assigns values from matrix 'Ke' (8x8 with constants) to matrix 'KG' which is initially filled with zeros. Numbers of rows and columns are given by the array nEN, which varies in another loop which is not shown here. The thing is sometimes nEN can have the same values as those obtained in the previous iteration, and so my loop rewrites the values which have been already present in 'KG', but I want them to be added together instead. How can I do that?? Thanks.

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

请先登录，再进行评论。

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

Answer 1

Vitaly 2013-4-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73016-how-to-add-values-to-already-existing-ones-in-a-matrix-using-for-loop#answer_83001

在 MATLAB Online 中打开

Nevermind, I solved the problem by accumulating values in KG with:

KG((2*(nEN(i1))-1),(2*(nEN(i2))-1)) = KG((2*(nEN(i1))-1),(2*(nEN(i2))-1)) + Ke(((2*i1)-1),((2*i2)-1));

and so on..

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

请先登录，再进行评论。

Answer 2

Cedric 2013-4-21

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73016-how-to-add-values-to-already-existing-ones-in-a-matrix-using-for-loop#answer_83003

编辑：Cedric 2013-4-21

在 MATLAB Online 中打开

Like this:

 KG = zeros(size(Ke)) ;           % Prealloc (are they same size?)
                                  % and initialize with 0's.
 for ...                          
    nEN = [1 2 4 5] ;             % Defined in the external loop.
    for i1 = 1 : 4
        bi_KG = 2 * nEN(i1) ;     % Base index for KG.
        bi_Ke = 2 * i1 ;          % Base index for Ke.
        KG(bi_KG-1,bi_KG-1) = KG(bi_KG-1,bi_KG-1) + Ke(bi_Ke-1,bi_Ke-1) ;
        KG(bi_KG,bi_KG)     = KG(bi_KG,bi_KG) + Ke(bi_Ke,bi_Ke);
    end
 end