Сворачивание по одному индексу - Matlab - Ответ 3889748

clc
clear
 
 
disp('matrix A');
A=fix(randn (4,4,4,4));
disp('matrix B');
B=fix(randn (4,4,4,4));
 
%Расчет матрицы D
for l=1:4
    for s=1:4
        for m=1:4
            sum = 0;
            for i=1:4
                for j=1:4
                    sum = A(l,s,i,j)*B(i,j,s,m)+ sum;
                end
            end
            D(l,s,m)=sum;
        end
    end
end
D
 
%Расчет единичной матрицы E
   for c1=1:4
       for c2=1:4
           for m1=1:4
               for m2=1:4
                   if ((c1==m1) && (c2==m2))
                       E(c1,c2,m1,m2) = 1;
                   else
                       E(c1,c2,m1,m2) = 0;
                   end
               end
           end
       end
   end
E
 
%Расчет матрицы F
for l1=1:4
    for l2=1:4
        for s=1:4
            for m=1:4
               sum = 0;
               for c1=1:4
                   for c2=1:4
                       sum= B(c1,c2,s,m);
                   end
               end
               F(l1,l2,s,m)=sum;
            end
        end
    end
end
F