(x,y)坐标的RMSE

时间:2013-05-10 20:15:09

标签: matlab dataset coordinates cartesian

我在笛卡尔平面上有一个坐标(X,Y)值的数据集,我想找到这些点的RMSE为1坐标(x1,y1)。基本上,(X,Y)在(x1,y1)左右,我想找到他们的RMSE。

任何人都可以帮助我,因为我不确定我做得对吗:

我这样做:

    Err = Err + sqrt[(X(i) - x1)^2 + (Y(i) - y1)^2]     - - - - Previous error + current error (Distance between points)

    RMSE = sqrt[(1/no_Of_Points)*Err^2]

这是对的吗?

我在MATLAB中这样做,所以如果需要,请随意编写代码。

这是我的代码:

    RMSEright = 0;
    countright = 0;

    for i = 1:1:size(VarName1,1)
        [x,y] = pol2cart(VarName2(i,1), VarName1(i,1));
        if x > 0
            RMSEright = RMSEright + (((featureright(1,1) - x)^2)+((featureright(1,2) - y)^2))^0.5;
            countright = countright + 1;
        end
    end

    RMSEright = ((1/countright)*RMSEright)^0.5

谢谢!

1 个答案:

答案 0 :(得分:0)

我在您的代码中看到您有兴趣测量中心x_centerx(1) > x_center(1)右边的点的RMSE。

您可以对此向量化,并避免使用循环。让z = z = [X Y];,然后让z = z(z(:,1) > z_center(1), :);

然后,计算RMSE sqrt(mean(vecnorm(z-z_center, 2, 2).^2))

这是示例代码:

% Center point
z_center = [1 1];
plot(z_center(1), z_center(2),'r*'); hold on;

% Data Points
X = [1.2256; 1.2931; 0.1213; 1.0387; 0.5158; 1.5895; 1.1485; 1.5914; 0.6550; 0.8026];
Y = [0.9002; 0.2724; 1.0856; 0.0871; 1.2780; 0.8944; 0.9888; 0.3219; 1.2895; 1.0577];
z = [X Y];
plot(z(:,1), z(:,2),'*');

% Select "right" points
z = z(z(:,1) > z_center(1), :);

% Compute RMSE
sqrt(mean(vecnorm(z-z_center, 2, 2).^2))
plot(z(:,1), z(:,2),'oc');

legend('center', 'points', 'right points');

返回

  

ans = 0.6710

和情节

rmse