0

作为我们作业的一部分，我们需要将 imageA 映射到 iamge B
我正在使用此代码进行逆映射
A = zeros(8,8);
A = [imageB(1,1) imageB(2,1) 0 0 1 0 (-1)*(imageB(1,1))imageA(1,1) (-1)(imageB(2,1))imageA(1,1);
0 0 imageB(1,1) imageB(2,1) 0 1 (-1)(imageB(1,1))imageA(2,1) (-1)(imageB(2,1))imageA(2,1);
imageB(1,2) imageB(2,2) 0 0 1 0 (-1)(imageB(1,2))imageA(1,2) (-1)(imageB(2,2))imageA(1,2);
0 0 imageB(1,2) imageB(2,2) 0 1 (-1)(imageB(1,2))imageA(2,2) (-1)(imageB(2,2))imageA(2,2);
imageB(1,3) imageB(2,3) 0 0 1 0 (-1)(imageB(1,3))imageA(1,3) (-1)(imageB(2,3))imageA(1,3);
0 0 imageB(1,3) imageB(2,3) 0 1 (-1)(imageB(1,3))imageA(2,3) (-1)(imageB(2,3))imageA(2,3);
imageB(1,4) imageB(2,4) 0 0 1 0 (-1)(imageB(1,4))imageA(1,4) (-1)(imageB(2,4))imageA(1,4);
0 0 imageB(1,4) imageB(2,4) 0 1 (-1)(imageB(1,4))imageA(2,4) (-1)(imageB(2,4))*imageA(2,4)];
imageA(1,4); imageA(2,4)];
B = reshape(imageA,[],1);
D = pinv(A)*B;


然而，我的朋友使用此代码：
D = inv(A'*A)*A'*B;


我得到与他相同的值，除了我得到 0.0000 而他得到 -0.0000 所有其余的正数或负数都是相同的。

为什么是这样 ？正确答案是什么？
你建议我可以做些什么来轻松检查我的结果？

4

1 回答 1

2

这些差异似乎源于数值的不稳定性。就机器精度而言，它们有多大（eps）？你为什么不使用backslash？

>> D = A \ B;
于 2012-12-06T15:04:59.210 回答