3D matrix with equation
5 views (last 30 days)
Show older comments
I have 2 matrices, Y and P, both of size m by n. I need to implement the following equation:
I'm having an issue with the index changing between j and l.
My initial try was:
eq = 2*Y.*Y - Y.*P - Y.*P
However, this is wrong of course since the multiplication would be with the same indices. I was thinking I would have to expand the matrices into a third dimension to overcome this problem with an extra index and then run some similar code for the equation, but I don't really understand how to perform this expansion and how the code for the equation would look like afterwards.
If you know another solution than to use a 3D matrix, this would of course also be much appreciated.
If it helps, I need to sum over index i afterwards, so this will reduce the dimensions again afterwards.
To clarrify: i runs from 1 to m and j & l run from 1 to n.
0 Comments
Accepted Answer
Matt J
on 2 Jun 2021
Edited: Matt J
on 2 Jun 2021
If it helps, I need to sum over index i afterwards, so this will reduce the dimensions again afterwards.
Don't multiply-then-sum. Do them in one combined step using matrix multiplications (and transpositions, where appropriate).
eq = 2*Y.'*Y - P.'*Y - Y.'*P %entries are eq(j,l)
More Answers (2)
J Chen
on 2 Jun 2021
does mean ? You can calculate the result using three loops on i, j, and l.
0 Comments
David Hill
on 2 Jun 2021
Should be something like this if I understood you correctly.
[m,n]=size(Y);
N=zeros(n);
for j=1:n
for l=1:n
N(j,l)=sum(2*Y(:,j)*Y(:,l)'-Y(:,l)*P(:,j)'-Y(:,j)*P(:,l)','all');
end
end
0 Comments
See Also
Categories
Find more on Matrix Indexing in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!