Question

关于imagesc作图 MATLAB

N=8;T = inv(dctmtx(N));%8×8的DCT变换矩阵X = T(:,1)*T(:,1)';A = ones(8,8)*0.125;figure(1);colormap('gray');imagesc(A);%显示图像axis square%画图区域是方形axis off %不显示轴线和标号figure(2);colormap('gray');imagesc(X);%显示图像axis square%画图区域是方形axis off %不显示轴线和标号A = 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250X = 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250为甚么两者绘图不一样？





Answer 1

在命令行输入下列指令，你会发现有趣的事情：

>> A(1,1) == X(1,1)
ans =
     0
>> A(1,1) == 0.1250
ans =
     1
>> X(1,1) == 0.1250
ans =
     0

说明看起来X里都是0.1250，但实际上并不是。

这是因为 matlab 显示结果只显示有限位的小数，看起来是0.1250，实际可能是0.125000000000000000000000000000000000091，即便你用format long，后面也是看不到的。所以实际上X里的都是0.1250的近似值，而且不一样，所以图像有差异。

可以调用roundn，只保留到小数点后5位，这样就看起来一样。

X = roundn(T(:,1)*T(:,1)',-5);