k=[P1(:) P2(:) P3(:) P4(:)]\(I(:)-A(:));

得到的k是有4個值的向量，對應於k1 k2 k3 k4

實際上，I,A,P1,P2,P3,P4是大小相等的矩陣

計算I=A+k1*P1+k2*P2+k3*P3+k4*P4的時候，是對應位置的元素得到對應位置的結果

所以，完全可以用I(:)、A(:)、P1(:)、.........將他們都線化為列向量

（這裡每個向量長度都是3003）

那麼原來的方程可以可以看作是一個很大的方程組

k1*P1(1)+k2*P2(1)+k3*P3(1)+k4*P4(1)=I(1)-A(1)

k1*P1(2)+k2*P2(2)+k3*P3(2)+k4*P4(2)=I(2)-A(2)

k1*P1(3)+k2*P2(3)+k3*P3(3)+k4*P4(3)=I(3)-A(3)

.........

k1*P1(3003)+k2*P2(3003)+k3*P3(3003)+k4*P4(3003)=I(3003)-A(3003)

這3003個方程的最優解

這是個線性方程組用矩陣表達是

Ak=b

其中 A=[P1(:) P2(:) P3(:) P4(:)] b=I(:)-A(:)

利用matlab的左除算符和容易的到最小二乘最優解 k=A\b;