基于曲率信息曲线重构程序??
大家好!我根据浙江大学 刘刚 文章(见附件),编写了一个基于离散曲率信息重构曲线的程序,请高手指点一下程序中的问题:clc
clear
format long
%k=;
%k=;
k=;
N=length(k);
Z=1;
KK=(N-1)*Z;
C=zeros(1,KK);
D=1;
ds=11.45/4;
ddx=zeros(1,KK);
ddy=zeros(1,KK);
nn=zeros(1,Z);
r=zeros(1,N);
AX=zeros(1,KK);
AY=zeros(1,KK);
AX(1)=0;
AY(1)=3;
for m=1:(N-1)
for n=1:Z
D=D+1;
nn(n)=(k(m+1)-k(m))*(n-1)/Z+k(m);
C(D)=nn(n);
if nn(n)==0
r(n)=1000000;
else
r(n)=abs(1/nn(n));
end
if nn(n)>0
alfa=ds/(2*3.1415928*r(n))*360;
else
alfa=ds/(2*3.1415928*r(n))*360+90;
end
ddx(D-1)=2*r(n)*sind(alfa/2)*cosd(alfa/2);
ddy(D-1)=2*r(n)*sind(alfa/2)*sind(alfa/2);
Tran=;
Trot=;
TTT=*Tran*Trot;
AX(D)=TTT(1);
AY(D)=TTT(2);
end
end
plot(AX,AY,'r-'); %重构曲线
hold on
%拟合对象
AAX=;
AAY=;
plot(AAX,AAY,'b-');
{:{40}:}请大家指点,基于离散曲率信息,重构曲线目前有很多研究,附件根据前人研究成果编写程序,没有成功,请高手予以指点一、二。 {:{33}:}没有人关注啊 Ref:
5.建议提问的网友分清 编程问题 和 专业问题 http://forum.vibunion.com/thread-36746-1-1.html
6.提问的智慧!!!!(发帖前请认真阅读) http://forum.vibunion.com/thread-21991-1-1.html
页:
[1]