白光散斑法自动测量微小位移的理论分析
- 格式:pdf
- 大小:229.41 KB
- 文档页数:4
白光散斑法自动测量微小位移的理论分析
陈 鹏 林铁生
(北方交通大学理学院
,北京
100044) 张世良
(北京联合大学信息学院
,北京
100101)
摘 要 对白光散班法自动测量物体微小位移的理论依据进行了详细的分析
,通过两次傅氏变换
导出了包含位移信息的信号峰值的完全表达式
,阐述了自动化测量的物理过程
.
关键词 白光散斑
CCD成像 傅氏变换
分类号
TN27
TheoreticalAnalysisforAutomaticMeasurements
ofMinuteDisplacementbyUsingMethodof
WhiteLightSpeckle
ChenPeng
LinTiesheng
(
CollegeofNaturalScience,NorthernJiaotongUniversity,Beijing100044)
ZhangShiliang
(
InformationCollegeofBeijingUnionUniversity,Beijing100101)
Abstract
Thispaperpresentsadetailedtheoreticalanalysisforautomaticmeasure2
mentsofminutedisplacementbyusingmethodofwhitelightspeckle.Acompleteex2
pressioninvolvingdisplacementinformationisderivedbyusingFouriertransformatikon
twice.Physicalprocessofautomaticmeasurementsisinterpreted.
Keywords
whitelightspeckle
CCDimaging
Fouriertransformation
1 实验简介
在运用白光散斑数字化技术对物体面内位移进行测量的实验中
,通过
CCD采集系统对样
品位移前后分别拍摄其数字散斑图
,如图
1所示
.再使用计算机和图象卡将这两张数字散斑图
图
1 白光散斑数字化测量系统进行重叠
,将重叠后的散斑图分割成许多子图(
64×
64pixels)
.将这些图作为分析对象
,对每一幅子图进行
快速傅氏变换(
FFT)
,得到数字杨氏条纹图
,对杨氏条
纹图再进行
FFT,发现在中心最大峰值的两侧产生对
称的信号峰值
,这两个对称峰值相对于中心最大峰值
的位移是分析子图对应样品上该点的位移值[1,2].
本文收到日期
19992
012
05 陈鹏男
1968年生工程师
emailbfxb@center.njtu.edu.cn1999年
10月
第
23卷第
5期 北 方 交 通 大 学 学 报
JOURNALOFNORTHERNJIAOTONGUNIVERSITY Oct.1999
Vol.23No.52 理论分析[3]
如图
2所示
,对子图进行分析的物理过程是基于两次
FFT,考虑到子图象的尺寸比较小,
图
2 对重叠子图进行
FFT的分析过程
可认为样品在子图象范围内各点的位移大小相等
,设其为
d=u
i+v
j(1)
设样品位移的子图为
h(
x,y)
,位移后的子图为
h(
x-
u,y-
v)
,重叠后的子图为
f(
x,
y)
,H(ω
x,ω
y)、
F(ω
x,ω
y)分别为
h(
x,y)、
f(
x,y)的傅氏变换谱
,Δ
为子图的面积
,ω
x、ω
y分
别为
x、
y的空间频率
,
则
H(ω
x,ω
y)
=κ
Δh(
x,y)exp[-i2π
(ω
xx+ω
yy)
]d
xd
y(2)
F(ω
x,ω
y)
=κ
Δ[h(
x,y)
+h(
x-u,y-v)
]exp[-i2π
(ω
xx+ω
yy)
]d
xd
y(3)
由傅氏变换位移定理知
κ
Δh(
x-u,y-v)exp[-i2π
(ω
xx+ω
yy)
]d
xd
y=
κ
Δh(
x-y)exp[-i2π
(ω
xx+ω
yy)
]d
xd
y・
exp[-i2π
(ω
xu+ω
yv)
],
则
F(ω
x,ω
y)
=H(ω
x,ω
y)・
{1
+exp[-i2π
(ω
xu+ω
yv)
]}(4)
又因为
1
+exp[-i2π
(ω
xu+ω
yv)
]=1
+cos[2π
(ω
xu+ω
yv)
]-isin[2π
(ω
xu+ω
yv)
]=
1
+cos2
[π
(ω
xu+ω
yv)
]-sin2
[π
(ω
xu+ω
yv)
]-
i2sin[π
(ω
xu+ω
yv)
]・
cos[π
(ω
xu+ω
yv)
]=
2cos2
[π
(ω
xu+ω
yv)
]-i2sin[π
(ω
xu+ω
yv)
]・
cos[π
(ω
xu+ω
yv)
]=
2cos[π
(ω
xu+ω
yv)
]exp[-iπ
(ω
x+uω
yv)
],
所以可得
F(ω
x,ω
y)
=H(ω
x,ω
y)・
2cos[π
(ω
xu+ω
yv)
]exp[-iπ
(ω
xu+ω
yv)
](5)
又因为
cosα
=
2
π1+2
3cos2α
-2
15cos4α
+…
+(-1)n+
12
(2
n-1)(2
n+1)・
cos2
nα
+…
,611北 方 交 通 大 学 学 报 第
23卷cosα
≈
2
π1+2
3cos2α
=2
π1+2
3(2cos2
α
-1)=
2
π1+4
3cos2
α
-2
3=2
3π(
1+4cos2
α
)1
所以可得
F(ω
x,ω
y)≈
H(ω
x,ω
y)・4
3π{1+4cos2
[π
(ω
xu+ω
yv)
]}・
exp[-iπ
(ω
xu+ω
yv)
](6)
令
A
s(ω
x,ω
y)、
A
h(ω
x,ω
y)分别为
F(ω
x,ω
y)、
H(ω
x,ω
y)的幅值
,
则
A
h(ω
x,ω
y)
=|H(ω
x,ω
y)
|(7)
A
s(ω
x,ω
y)
=|F(ω
x,ω
y)
|≈
H(ω
x,ω
y)・
{1
+4cos2
[π
(ω
xu+ω
yv)
]}(8)
令
d=
u
i+
v
j, ω
=ω
x
i+ω
y
j,则式(8)变为
A
s(ω
x,ω
y)
=A
h(ω
x,ω
y)・4
3π[1+4cos2
(π
d・ω
)
](9)
当
d・ω
=
n+1
2时
,为暗条纹
,当
d・ω
=
n时
,为亮条纹
,n=0
,±
1
,±
2
,…1
设
d・ω
1=
n,
d・ω
2=
n+1
,
则
d・(ω
1-ω
2)=1
,
d・Δ
ω
′
=11
令Δ
ω
′在
d上的投影为Δ
ω
,则
|
d|・
|Δ
ω
|=11
因为
d=
u
i+
v
j是确定的
,所以
|Δ
ω
|为定值
,即条纹的间隔1设ω
1为一定
,则Δ
ω
′对应
的一系列点都是亮点
,它们构成亮条纹1即知
,条纹是垂直
d且与
|
d|成反比的1故式(9)表
示的物理含义正是数字杨氏振幅条纹1
引入可见度
V=(
I
max-I
min)
/(
I
max+I
min)
,
式中
I
max、
I
min分别为光强的最大值、最小值1
由式(9)可知
,对于
A
s(ω
x,ω
y)来说
,V=2
31故杨氏振幅条纹的普遍表达式为
A
s(ω
x,ω
y)
=A
h(ω
x,ω
y)・
[1
+V・
cos(
2π
d・ω
](10)
式(10)对应的光强图即为数字杨氏条纹图1
现在对上述数字杨氏条纹再作
FFT分析1
令
G
As(ζ
,η
)、
G
Ah(ζ
,η
)分别为
A
s(ω
x,ω
y)、
A
h(ω
x,ω
y)的傅氏频谱
,
则有
G
Ah(ζ
,η
)
=κ
ΩA
h(ω
x,ω
y)exp
[-i2π
(ω
xζ
+ω
yη
)
]dω
xdω
y,
G
As(ζ
,η
)=κ
Ω4
πA
h(ω
x,ω
y)
{1
+Vcos[2π
(
uω
x+vω
y)
]}・
exp
[-i2π
(ω
xζ
+ω
yη
)
]
dω
xdω
y=
4
πG
Ah(ζ
,η
)
+4
πκ
ΩA
h(ω
x,ω
y)・
Vcos[2π
(
uω
x+vω
y)
]・711第
5期 陈鹏等
:白光散斑法自动测量微小位移的理论分析