最小二乘法拟合球面源代码（matlab）

转自：https://blog.csdn.net/zhuoyueljl/article/details/53675460

原

最小二乘法拟合球面源代码（matlab）

2016年12月15日 20:13:56

阅读数：3200

1、源代码

[html]  view plain  copy

1. <span style="font-size:14px;">[x,y,z]=sphere(5);  
2. x = [0,0,1,0,0,0.7].';   %自己的数据  
3. y = [0,1,0,0,-1,0.7].';  %自己的数据  
4. z = [1,0,0,-1,0,0].';    %自己的数据  
5. % 拟合  
6. x = x(:); y = y(:); z = z(:);  
7. data=unique([x(:)-0.1,y(:)+0.2,z(:)],'rows');  
8. f=@(p,data)(data(:,1)-p(1)).^2+(data(:,2)-p(2)).^2+(data(:,3)-p(3)).^2-p(4)^2;  
9. p=nlinfit(data,zeros(size(data,1),1),f,[0 0 0 1]')%拟合的参数  
10. hold on  
11. plot3(data(:,1),data(:,2),data(:,3),'o')  
12. [X,Y,Z]=meshgrid(linspace(-14,14));  
13. V=(X-p(1)).^2+(Y-p(2)).^2+(Z-p(3)).^2-p(4)^2;  
14. isosurface(X,Y,Z,V,0);  
15. alpha .5;camlight;axis equal;grid on;view(3);  
16. title(sprintf('(x-%f)^2+(y-%f)^2+(z-%f)^2=%f',p(1),p(2),p(3),p(4)^2))</span>  

在matlab的输入窗口输入源代码后按Enter键如下图所示

2、实验结果