:Բ

ʽ:
[a,m,er]=CylindFit(x,y,z)
x:ṩxݾ
y:ṩyݾ
z:ṩzݾ

a:صԲ淽̲,ɿ˵
m:صĸ
er:صľ

˵:

1ϵԲ淽Ϊ
$$
(x-a_1z-a_2)^2 + (y-a_3z-a_4)^2 = a_5^2
$$
ĸΪm=[m(1) , m(2) , m(3) ]  a(1)=m(1)/m(3) , a(3)=m(2)/m(3)Բ뾶Ϊ a(5)̿Ϊxƽ-a(2),yƽ-a(4),zƽ0Ȼĸmó淽.

2ӵ1֪,ĸе m(3)Ϊ0,ȷǷΪ0,ɷֱxyzһ,CylindFit(x,y,z)CylindFit(x,z,y)CylindFit(y,z,x)Ȼѡ񷵻ֵerСĽΪսɡ

ο:ι.淽̵һע[J].ߵѧо,2012,15(2):25-26

:

s=0:0.1:6;
s=s';
x=10*cos(s);
y=10*sin(s);
z=x*y;
//ĴԵ,ʵʹӦֱӵⲿ

[a,m,er]=CylindFit(x,y,z)//سִõ½,Ժԭʼʵ x^2+y^2=10^2 ϵ÷ǳ
a =
[ 0.00000000000000   -1.2212453270E-15    0.00000000000000   -9.2981178312E-16    10.0000000000000 ]
m =
[ 0.00000000000000    0.00000000000000    1.00000000000000 ]
er =
[ 4.0358827391E-12 ]