Updating a matrix based on conditions in other matrixes

Question

Ulrik William Nash 2017-7-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/347994-updating-a-matrix-based-on-conditions-in-other-matrixes

编辑： Ulrik William Nash 2017-7-10

Updating_a_matrix_based_on_conditions_in_other_matrixes.txt

I have coded the attached loop to update matrix Lr but I would ideally like to use code that can achieve what the shown loop is doing in one sweep (i.e. without needing to loop); the actual matrixes I have can be extremely large.

The code is also shown below, but the layout is off:

k = 1;
pos_neg = [1;1;0];
E = [round(rand(1,3));round(rand(1,3));round(rand(1,3))];
Lr = zeros(numel(pos_neg(:,1)),numel(E(:,1)));
L = 0.5;
for i = 1:numel(E(:,1)) 
    for ii = 1:numel(E(:,1)) 
        if  sum(E(i,:)) <= k
            if  E(i,ii) == 0
                Lr(i,ii) = L;
            end
            if  E(i,ii) == 1
                Lr(i,ii) = -L;
            end    
        end
        if  sum(E(i,:)) >= k
            if  E(i,ii) == 1
                Lr(i,ii) = L;
            end
            if  E(i,ii) == 0
                Lr(i,ii) = -L;
            end    
        end
    end
end

2 个评论
显示无隐藏无

Ngonidzashe Nedziwe 2017-7-10

hi i would like to help but im confused about what u exactly mean or are trying to achieve when u say in one sweep

Ulrik William Nash 2017-7-10

Thank you, Ngonidzashe, what I mean is without looping.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2017-7-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/347994-updating-a-matrix-based-on-conditions-in-other-matrixes#answer_273505

在 MATLAB Online 中打开

If I work through the logic correctly,

Lr = L*(1 - 2 * bsxfun( @eq, sum(E,2) <= k, E));

Or if you have R2016b or later,

Lr = L*(1 - 2 * ((sum(E,2) <= k) == E));

4 个评论
显示 2更早的评论隐藏 2更早的评论

Ulrik William Nash 2017-7-10

Thank you Walter, no problem. In the mean time I will look at your approach. Maybe I can figure it out.

Ulrik William Nash 2017-7-10

编辑：Ulrik William Nash 2017-7-10

Hi Walter,

I think your mistake was very small. As far as I can tell, all that is required is the following change:

From: Lr = L*(1 - 2 * bsxfun( @eq, sum(E,2) <= k, E))

To:

Lr = L*(2 * bsxfun( @eq, sum(E,2) >= k, E)-1)

请先登录，再进行评论。

Answer 2

Jan 2017-7-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/347994-updating-a-matrix-based-on-conditions-in-other-matrixes#answer_273503

编辑：Jan 2017-7-10

在 MATLAB Online 中打开

Start with cleaning the code by avoiding all repeated calculations inside the loop:

k = 1;
pos_neg = [1;1;0];
E  = [round(rand(1,3));round(rand(1,3));round(rand(1,3))];
Lr = zeros(numel(pos_neg(:,1)),numel(E(:,1)));
L = 0.5;
sE = size(E, 1);
for i = 1:sE
    sumEi = sum(E(i,:));
    for ii = 1:sE
        if sumEi <= k
            if  E(i,ii) == 0
                Lr(i,ii) = L;
            end
            if  E(i,ii) == 1
                Lr(i,ii) = -L;
            end    
        end
        if sumEi >= k
            if  E(i,ii) == 1
                Lr(i,ii) = L;
            end
            if  E(i,ii) == 0
                Lr(i,ii) = -L;
            end    
        end
    end
end

Now you see that the "if sumEi <= k" test does not depend on the "for ii" loop. Then move it outside this loop:

sE = size(E, 1);
for i = 1:sE
    sumEi = sum(E(i,:));
    if sumEi <= k
      for ii = 1:sE
        if  E(i,ii) == 0
           Lr(i,ii) = L;
        elseif E(i,ii) == 1
           Lr(i,ii) = -L;
        end    
      end
    elseif sumEi >= k
      for ii = 1:sE
        if  E(i,ii) == 1
           Lr(i,ii) = L;
        elseif E(i,ii) == 0
           Lr(i,ii) = -L;
        end    
      end
    end
end

Now see, shat you want to achieve actually (this should be contained as a clear and exhaustive comment in the code): Lr is set to L or -L. So start with set all elements to L and swap the sign for the specific elements:

Lr = repmat(L size(pos_neg, 1), size(E,1));

size(E, 1) is more efficient than numel(E(:,1)), because the latter creates a temporary vector only to measure its length.

sumEi = sum(E, 1);
index = ((sumEi <= k) & (E == 1)) | ...
         (sumEi >= k) & (E == 0));    % >= R2016b

Note that sumEi == k belongs to both cases. I assume one of them should be a > or < .

If the values of E are either 1 or 0, this can be expressed in one logical comparison only. I cannot test this currently, but it is something like:#

Lr    = repmat(L, size(pos_neg, 1), size(E,1));
index = xor(sum(E, 1) <= k, E == 1);  % Or E==0? Or < instead of <= ? Try this by your own
Lr(index) = -L;