FTF

基于物体三维形貌测量的傅里叶变换测量轮廓术及相关精度的计算。

clear all
clc
N=100;
M=2; %滤波器阶数
del=0.001;
a1=1.558;
a2=-0.81;
v=randn（1N）;
u（1）=v（1）;
for n=2:N-1
    u（2）=a1*u（1）+v（2）;
    u（n+1）=a1*u（n）+a2*u（n-1）+v（n+1）; 
    d（2）=0.2*u（2）+0.7*u（1）;
    d（n+1）=0.2*u（n+1）+0.7*u（n）;
end
%初始化
am=zeros（NM）;
bm=zeros（NM）;
wm=zeros（NM）;
gm=zeros（NM）;
k=0;
r（1:）=ones（1M+1）;
r（2:）=ones（1M+1）;
epsf=del;
epsb=del;
km=zeros（1M）;  
for n=3:N  
    uf=[u（n-1） u（n-2）];
    ub=[u（n） u（n-1）];
    ef（nn-1）=u（n）-uf*am（n-1:）‘;
    ef（nn）=r（n-1M）*ef（nn-1）;%ef是标量
    epsf0=epsf;
    epsf=epsf+ef（nn-1）;
    am（n:）=am（n-1:）+ef（nn-1）*gm（n-1:）;
    r（nM+1）=epsf0*r（n-1M）/epsf;
    km1=[0 gm（n-1:）]+ef（nn）/epsf*[1 -am（n:）];
   for m=1:M
       km（m）=km1（m）;
   end
   k=km1（end）;
    eb（2n-1）=ub*bm（n-1:）‘;
    eb（n+1n-1）=u（n+1-M）-ub*bm（n-1:）‘;
    r（nM）=1/（1-k*eb（nn-1））*r（nM+1）;
    eb（nn）=r（nM）*eb（nn-1）; 
    epsb=epsb+eb（nn）*eb（nn-1）;
    gm（n:）=（km+k*bm（n-1:））*r（nM）/r（nM+1）;
    bm（n:）=bm（n-1:）+gm（n:）*eb（nn-1）;
    e（nn-1）=d（n）-ub*wm（n-1:）‘;
    wm（n:）=wm（n-1:）+gm（n:）*e（nn-1）;  
end
plot（wm）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      17447  2008-06-04 20:11  FTF\1.jpg

     文件       2163  2008-05-28 20:43  FTF\FTF.fig

     文件       1135  2008-05-28 20:30  FTF\FTF.m

     文件       1089  2013-11-11 16:32  FTF\FTF2.asv

     文件       1089  2013-11-11 16:35  FTF\FTF2.m

     目录          0  2013-11-11 16:35  FTF

----------- ---------  ---------- -----  ----

                23141                    7