像matlab中的最小二乘的内核

Question

好的，所以我有四个100x100的2d矩阵（让我们称之为A，B，C和D）我希望解决线性方程组（Y = Xt）让我们说3x3的内核（在代码中）下面你可以看到iA，...，iB是3x3矩阵。

总而言之，我希望以优化的方式执行以下操作：

for i=2:size(A,1)-1
   for j=2:size(A,2)-1
      iA=A(i-1:i+1,j-1:j+1);
      iB=B(i-1:i+1,j-1:j+1);
      iC=C(i-1:i+1,j-1:j+1);
      iD=D(i-1:i+1,j-1:j+1);

      X=[iA(:),iB(:),ones(size(iA(:)));ones(size(iA(:))),iA(:),iB(:)]
      Y=[iC(:);iD(:)]

      t(i,j,:)=X\Y
   end
end

有没有更好的方法呢？类似内核的方法可能吗？