function y = matrix_3D(x,y)
% 三维矩阵按页相乘(not @times点乘),页数相同
     	[row_x,col_x,p] = size(x);
		[row_y,col_y,p] = size(y);% col_x == row_y
		x1 = reshape(x,row_x,col_x,1,p);
		y1 = reshape(y,row_y,1,col_y,p);
		y2 = permute(y1,[2 1 3 4]); % (1,row_y,p,col_y) y需要和x相乘的对应列拆成一列一列的,不相等的维数放在两边分别对应对方的1
		y = reshape(sum(bsxfun(@times,x1,y2),2),row_x,col_y,p);
end