计算出球面的球度误差,貌似cftool不能进行三维的曲面拟合,小弟弄了好长时间一直没有进展,求大神帮助啊 ,不胜感谢了!
我把球面拟合出来了,也画出来了,但不知道什么是球度误差,你可以用数据和拟合出的球面方程自己算。
clear;
clc;
[x,y,z]=sphere(5);
data=10*unique([x(:)-0.1,y(:)+0.2,z(:)],'rows');
data=data+rand(size(data,1),size(data,2))%我随便生成的,你用你的数据替代
f=@(p,data)(data(:,1)-p(1)).^2+(data(:,2)-p(2)).^2+(data(:,3)-p(3)).^2-p(4)^2;
p=nlinfit(data,zeros(size(data,1),1),f,[0 0 0 1]')%拟合的参数
hold on
plot3(data(:,1),data(:,2),data(:,3),'o')
[X,Y,Z]=meshgrid(linspace(-14,14));
V=(X-p(1)).^2+(Y-p(2)).^2+(Z-p(3)).^2-p(4)^2;
isosurface(X,Y,Z,V,0);
alpha .5;camlight;axis equal;grid on;view(3);
title(sprintf('(x-%f)^2+(y-%f)^2+(z-%f)^2=%f',p(1),p(2),p(3),p(4)^2))
球度误差就是模拟出一个球,球外距球面最远的点与球内距球面最远的点的距离的差值。
数据上传不上去 麻烦你有youxiang吗 我把数据和一些资料发你看看 谢谢你啦