白光散斑法自动测量微小位移的理论分析

  • 格式: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

+…

,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(

󰁩

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π

)

]dω

xdω

y,

G

As(ζ

)=κ

Ω4

πA

h(ω

x,ω

y)

{1

+Vcos[2π

(

x+vω

y)

]}・

exp

[-i2π

)

]

xdω

y=

4

πG

Ah(ζ

)

+4

πκ

ΩA

h(ω

x,ω

y)・

Vcos[2π

(

x+vω

y)

]・711第

5期 陈鹏等

:白光散斑法自动测量微小位移的理论分析