作业帮matlab怎么将点云数据用最小二乘方法拟合出平面
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/07/17 16:09:27
matlab怎么将点云数据用最小二乘方法拟合出平面
数据已经导进去了
本人小白,
数据已经导进去了
本人小白,
%点X,Y,Z到平面Ax+By+Cz+D=0的距离为
%d(ABCD,XYZ)=|AX+BY+CZ+D|/sqrt(A^2+B^2+C^2)
%ABCD四个变量只有三个是互相独立的
%设A=cos(a);B=sin(a)*cos(b);C=sin(a)*sin(b)
%那么A^2+B^2+C^2=1,距离公式化简为
%d(abc,XYZ)=|cos(a)*X+sin(a)*cos(b)*Y+sin(a)*sin(b)*Z+c|
%现在有已知点序列X,Y,Z,求参数 a b c
%先构造一个函数fun(p) 输入参数为p,其中p(1)=a,p(2)=b,p(3)=c
%使用 lsqnonlin求得p,使得sum((fun(p))^2)最小
fun=@(p) cos(p(1))*X+sin(p(1))*cos(p(2))*Y+sin(p(1))*sin(p(2))*Z+p(3);
p = lsqnonlin(fun,[0 0 0]);
A=cos(p(1));B=sin(p(1))*cos(p(2));C=sin(p(1))*sin(p(2));D=p(3);
再问： ��л����~~~��������ܼ�һ�������ʾ���������ƽ���ͼ����
再答： ��ʵ���������ش�Ҳ�ǶԵģ�ֻ�Dz�������һЩ������� WM_THU�Ļش���D=0ʱ���������� tianxiawulang�Ļش��� C=0ʱ��Ҳ�������� һ������£����ȷ��ƽ�治�Ǻ�z��ƽ�еĻ��� ��ϣ���ͼ����tianxiawulang�Ļ���һ�� fun=@(p,x,y) p(1)*x+p(2)*y+p(3); para=lsqcurvefit(@(p,data) fun(p,data(:,1),data(:,2)),[1 1 1],[X Y],Z); %���Ϊ Z=para(1)*X+para(2)*Y+para(3) plot3(X,Y,Z,'.');hold on; h=ezsurf(@(x,y)fun(para,x,y),[min(X) max(X) min(Y) max(Y)]); hold off; set(h,'edgecolor','none','facealpha',0.3,'FaceColor','r');
%d(ABCD,XYZ)=|AX+BY+CZ+D|/sqrt(A^2+B^2+C^2)
%ABCD四个变量只有三个是互相独立的
%设A=cos(a);B=sin(a)*cos(b);C=sin(a)*sin(b)
%那么A^2+B^2+C^2=1,距离公式化简为
%d(abc,XYZ)=|cos(a)*X+sin(a)*cos(b)*Y+sin(a)*sin(b)*Z+c|
%现在有已知点序列X,Y,Z,求参数 a b c
%先构造一个函数fun(p) 输入参数为p,其中p(1)=a,p(2)=b,p(3)=c
%使用 lsqnonlin求得p,使得sum((fun(p))^2)最小
fun=@(p) cos(p(1))*X+sin(p(1))*cos(p(2))*Y+sin(p(1))*sin(p(2))*Z+p(3);
p = lsqnonlin(fun,[0 0 0]);
A=cos(p(1));B=sin(p(1))*cos(p(2));C=sin(p(1))*sin(p(2));D=p(3);
再问： ��л����~~~��������ܼ�һ�������ʾ���������ƽ���ͼ����
再答： ��ʵ���������ش�Ҳ�ǶԵģ�ֻ�Dz�������һЩ������� WM_THU�Ļش���D=0ʱ���������� tianxiawulang�Ļش��� C=0ʱ��Ҳ�������� һ������£����ȷ��ƽ�治�Ǻ�z��ƽ�еĻ��� ��ϣ���ͼ����tianxiawulang�Ļ���һ�� fun=@(p,x,y) p(1)*x+p(2)*y+p(3); para=lsqcurvefit(@(p,data) fun(p,data(:,1),data(:,2)),[1 1 1],[X Y],Z); %���Ϊ Z=para(1)*X+para(2)*Y+para(3) plot3(X,Y,Z,'.');hold on; h=ezsurf(@(x,y)fun(para,x,y),[min(X) max(X) min(Y) max(Y)]); hold off; set(h,'edgecolor','none','facealpha',0.3,'FaceColor','r');