3D matrix with equation

3 次查看（过去 30 天）
Wilko
Wilko 2021-6-2
评论： Wilko 2021-6-2
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.

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

采纳的回答

Matt J
Matt J 2021-6-2
编辑：Matt J 2021-6-2
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)
  1 个评论
Wilko
Wilko 2021-6-2
Thanks, this solves my issue. I did this earlier in a similar situation, but difference in notation together with some reading about 3D matrices earlier messed up my train of thought.

请先登录，再进行评论。

更多回答（2 个）

J Chen
J Chen 2021-6-2
does mean ? You can calculate the result using three loops on i, j, and l.
David Hill
David Hill 2021-6-2
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

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by