统计误差计算 Matlab

Question

我有两个值向量，我想对它们进行统计比较。为简单起见，假设 A = [2 3 5 10 15] 和 B = [2.5 3.1 4.8 10 18]。我想计算标准偏差、均方根误差 (RMSE)、平均值，并方便地呈现，可能是直方图。你能帮我怎么做，以便我理解吗？我知道问题可能很简单，但我是新手。非常感谢！

编辑： 这就是我想要实现 RMSE 的方式。

dt = 1;
for k=1:numel(A)
err(k)=sqrt(sum(A(1,1:k)-B(1,1:k))^2/k);
t(k) = dt*k; 
end

但是它给了我比我预期更大的价值，因为例如3 和 3.1 仅在 0.1 上有所不同。 这就是我如何计算每个周期的参考值与该周期中相应估计值之间的误差。 你能告诉我，我做得对还是错？

abs_err = A-B;

Answer 1

您循环遍历向量的方式不是逐个元素，而是通过增加向量长度，也就是说，您在每次迭代时都比较以下内容：

      A(1,1:k)    B(1,1:k)
      --------    --------
k=1   [2]         [2.5]
 =2   [2 3]       [2.5 3.1]
 =3   [2 3 5]     [2.5 3.1 4.8]
....

您绝不会比较仅 2 和 2.1！

假设A 和B 是相同长度的向量（并且都是列向量或行向量），那么您需要函数std(A-B)、mean(A-B)，如果您查看matlab 交换，您会发现用户贡献的rmse(A-B)，但您也可以将 RMSE 计算为sqrt(mean((A-B).^2))。至于显示直方图，试试hist(A-B)。

在你的情况下：

dt = 1;
for k=1:numel(A)
    stdab(k) = std(A(1,1:k)-B(1,1:k));
    meanab(k) = mean(A(1,1:k)-B(1,1:k));
    err(k)=sqrt(mean((A(1,1:k)-B(1,1:k)).^2));
    t(k) = dt*k; 
end

如果您想计算每个向量对差异A(1,1:k)-B(1,1:k) 的直方图，您还可以在循环中包含hist(A(1,1:k)-B(1,1:k))。