利用MATLAB计算三维坐标序列距离误差程序
2014-05-07 15:09
686 查看
利用MATLAB计算三维坐标序列距离误差程序
1.三维坐标储存在文件中,格式如下:
各坐标间的距离真值是一定值,计算相邻距离的标准差。
2.MATLAB程序如下:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 名称:caculateAccuracy.m
% 功能:读取三维世界坐标,计算精度
% 作者:LYC
% 单位:中科院苏州医工所
% 日期:2014.5.5
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 清理各变量
clc, clear
%% 设定值
dataFile = 'v20131114-2.txt'; % 设定文件名称
errorFile = 'v20131114-2_error.txt';
realValue = 1.00; % 设定真值
%% 读取文件
fid = fopen(dataFile, 'r'); % 打开文件
k = 0; % 循环计数变量
curr = [0, 0, 0]; % 临时变量
a = []; % 数据变量
while ~feof(fid) % 判断是否到了文件末尾
curr = fscanf(fid, '(%f, %f, %f)\n', 3);
if ~isempty(curr)
k = k+1;
a = [a, curr];
end
end
fclose(fid);
%% 计算各量之间的欧氏距离
dist = []; % 储存欧氏距离
for i = 1:k-1
curr2 = (a(1,i)-a(1,i+1))^2 + (a(2,i)-a(2,i+1))^2 + (a(3,i)-a(3,i+1))^2; % 计算当前距离的平方
curr3 = sqrt(curr2); %计算当前距离
dist = [dist, curr3];
end
%% 各距离减去距离真值
distError = dist - realValue; % 储存距离误差
%% 利用贝塞尔公式计算标准差
distError2 = distError.^2; % 计算误差平方
sumError2 = sum(distError2);
std = sqrt(sumError2/(k-1))
%% 将结果写入文件
fid2 = fopen(errorFile, 'w'); % 打开文件
fprintf(fid2, 'standard deviation = %f\r\n', std);
fprintf(fid2, 'distance errors are:\r\n');
for i = 1:k-1
fprintf(fid2, '%f\r\n', distError(i)); %输入当前数据与一空格
end
fclose(fid2);
disp('误差写入文件成功!');
3.误差分析结果储存在文件中,如下所示
1.三维坐标储存在文件中,格式如下:
各坐标间的距离真值是一定值,计算相邻距离的标准差。
2.MATLAB程序如下:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 名称:caculateAccuracy.m
% 功能:读取三维世界坐标,计算精度
% 作者:LYC
% 单位:中科院苏州医工所
% 日期:2014.5.5
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 清理各变量
clc, clear
%% 设定值
dataFile = 'v20131114-2.txt'; % 设定文件名称
errorFile = 'v20131114-2_error.txt';
realValue = 1.00; % 设定真值
%% 读取文件
fid = fopen(dataFile, 'r'); % 打开文件
k = 0; % 循环计数变量
curr = [0, 0, 0]; % 临时变量
a = []; % 数据变量
while ~feof(fid) % 判断是否到了文件末尾
curr = fscanf(fid, '(%f, %f, %f)\n', 3);
if ~isempty(curr)
k = k+1;
a = [a, curr];
end
end
fclose(fid);
%% 计算各量之间的欧氏距离
dist = []; % 储存欧氏距离
for i = 1:k-1
curr2 = (a(1,i)-a(1,i+1))^2 + (a(2,i)-a(2,i+1))^2 + (a(3,i)-a(3,i+1))^2; % 计算当前距离的平方
curr3 = sqrt(curr2); %计算当前距离
dist = [dist, curr3];
end
%% 各距离减去距离真值
distError = dist - realValue; % 储存距离误差
%% 利用贝塞尔公式计算标准差
distError2 = distError.^2; % 计算误差平方
sumError2 = sum(distError2);
std = sqrt(sumError2/(k-1))
%% 将结果写入文件
fid2 = fopen(errorFile, 'w'); % 打开文件
fprintf(fid2, 'standard deviation = %f\r\n', std);
fprintf(fid2, 'distance errors are:\r\n');
for i = 1:k-1
fprintf(fid2, '%f\r\n', distError(i)); %输入当前数据与一空格
end
fclose(fid2);
disp('误差写入文件成功!');
3.误差分析结果储存在文件中,如下所示
相关文章推荐
- 利用MATLAB计算三维坐标序列距离误差程序
- 求两三维坐标点间的欧几里得距离的程序(MATLAB语言)
- 利用php 计算两点地理经纬坐标之间的距离
- 各地图gps坐标系统比较与转换以及经纬度距离计算MATLAB脚本
- 实现 计算两个三维坐标点的 距离
- 如何利用MATLAB并行计算缩短程序运行时间
- 利用java、js或mysql计算高德地图中两坐标之间的距离
- 定义一个“点”(Point)类用来表示三维空间中的点(有三个坐标)。要求如下: (1)可以生成具有特定坐标的点对象。 (2)提供可以设置三个坐标的方法。 (3)提供可以计算该“点”距原点距离平方的方法。 (4)编写主类程序验证。
- oracle函数计算两个坐标点距离
- 利用Modeller构建未知序列三维结构(利用cryo-EM map)的详细流程
- 用c语言编写:输入两点坐标(X1,Y1),(X2,Y2),计算并输出两点间的距离。
- 三维场景计算任意两点的空间距离
- PHP根据两点间的坐标计算距离
- 大规模三维空间纤维距离计算程序[CUDA校园程序设计大赛]
- matlab 利用while循环计算平均值和方差
- 一个C语言写的三维坐标的小程序
- lightoj 1240 计算机图论(计算三维点到线段的距离)
- 摄像机标定用于机械臂抓举等(利用标定将图像上的点映射为三维坐标)
- Matlab投影仿真及三维曲面重构实现及演示程序
- 输入两点坐标(X1,Y1),(X2,Y2),计算并输出两点间的距离。