陈家璧版_光学信息技术原理及应用习题解答(8-11章)

习 题

8．1利用4f系统做阿贝—波特实验，设物函数t（x1，y1）为一无限大正交光栅 ?1xx??1yy? t(x1,y1)??rect(1)?comb(1)???rect(1)?comb(1)?

a1b1??b2a2b2??b1 其中a1、a2分别为x、y方向上缝的宽度，b1、b2则是相应的缝间隔。频谱面上得

到如图8-53（a）所示的频谱。分别用图8-53（b）（c）（d）所示的三种滤波器进行滤

波，求输出面上的光强分布（图中阴影区表示不透明屏）。

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

（a） （b） （c） （d） 图8.53（题8.1 图）

解答：根据傅里叶变换原理和性质，频谱函数为 T ( fx , fy ) = ? [ t ( x1 , y1 )]

1b1 = { 将函数展开得 T ( fx , fy ) =

a1b1? [rect(x1a1)]·? [comb(x1b1)] } ?{

1b2? [rect(y1a2)·? [comb(y1b2)]}

{sinc(a1fx)+sinc(a1b1a2b2)δ(fx-1b11)+sinc(a1b1a2b2)δ(fx+1b11b2)+???}

?

a2b2{sinc(a2fy)+sinc()δ(fy-b2)+sinc()δ(fy+)+???}

（1）用滤波器（b）时，其透过率函数可写为

1 fx = + 1/ b1 fy = 0

F ( fx , fy ) = 0 fx ? 1/ b1 fy = 任何值 滤波后的光振幅函数为 T·F =

a1b1sinc(a1b1)[δ(fx-1b1)+δ(fx+1b1)]

输出平面光振幅函数为 t’（x3，y3）= ? -1[ T·F ]

a1b1a1b12?x3b1)]+exp[-j(2?x3b1)]}

= sinc(){exp[j( =

2a1b1sinc(a1b1)?cos(2?x3b1)

输出强度分布为 I（x3，y3）=

4a1b1222sinc(2a1b1)?cos(22?x3b1)

=

2a1b12sinc2(a1b1)?cos(4?x3b1) - C

其中C是一个常数，输出平面上得到的是频率增加一倍的余弦光栅。 （2）用滤波器（c）时，其透过率函数可写为

1 fx ，fy ? 0

F ( fx , fy ) = 0 fx = fy = 0 滤波后的光振幅函数为 T·F =

a1b1a2b2{sinc(a1b1a2b2)?(fx-1b11b2)+sinc(a1b1)?(fx+1b1)+???}

?

{sinc()?(fy-)+sinc(a2b2)?(fy+1b2)+???}

输出平面光振幅函数为 t’（x3，y3）= ? -1[ T·F ]

1b1 = {

[rect(x3a1)?comb(x3b1)] -

a1b1rect(x3b1)}

× {

1b2[rect(y3a2)?comb(y3b2)] -

a2b2rect(y3b2)}

输出强度分布为

I（x3，y3）= ? t’（x3，y3）?2

有两种可能的结果，见课本中图8．9和图8．10。

（3）用滤波器（d）时，输出平面将得到余弦光栅结构的强度分布，方向与滤波狭缝方向垂直，周期为b’，它与物光栅周期b1、b2的关系为

=b’11b12+1b22

8．2 采用图8-53（b）所示滤波器对光栅频谱进行滤波，可以改变光栅的空间频率，若光栅线密度为100线/mm，滤波器仅允许 + 2级频谱透过，求输出面上干板记录到的光栅的线密度。

解答：根据对8．1题的分析，当滤波器仅允许+ 2级频谱通过时，输出平面上的光振幅应

表达为

t’（x3）= ? -1 { sinc(a1b1a1b1)[?(fx-2b1)+?(fx+2b1)]}

=

2a1b1sinc()cos4?x3b1

其振幅分布为一周期函数，空间频率是基频的2倍。而干板记录到的是强度分布： I =

4a1b122sinc22(a1b1()cos24?x3b18?x3b1

=

2a1b12sinc2a1b1)cos - C

其中C是一个常数。

答：干板上记录到的光栅频率是基频的4倍，即400线/mm。

8．3 在4f系统中，输入物是一个无限大的矩形光栅，设光栅常数d = 4，线宽a =1，最

大透过率为1，如不考虑透镜有限尺寸的影响， （a）写出傅里叶平面P2上的频谱分布表达式；

（b）写出输出平面复振幅和光强分布表达式； （c）在频谱面上作高通滤波，挡住零频分量，写出输出平面复振幅和光强分布表达 式； （d）若将一个π位相滤波器 exp（jπ） x2，y2 ≤ x0，y0

H（x2，y2）=

0 其它

放在P2平面的原点上，写出输出平面复振幅和光强分布表达式，并用图形表示。

解答：将8．1题结果代入，其中b1 = d = 4，a1 = a = 1，除去与y分量有关的项，可得

（a）P2平面上的频谱分布为：

T(fx)=14{sinc(fx)+sinc(14)?(fx-14)+sinc(14)?(fx+14)+???}

（b）输出平面：

复振幅 t（x3）= ? [T（fx）]

若不考虑透镜尺寸的影响，它应该是原物的几何像，即 t (x3) =

14[rect(x3)?comb(2

-1

x34)]

x34)]

2光强分布 I (x3) = | t (x3)| =

116[rect(x3) ?comb((c)挡住零频分量，输出平面情况与8．1题（3）相同，即 t (x3) =

14[rect(x3)?comb(x34)]-14rect(x34)

I = | t (x3) | 2

由于a = d / 4，所以强度将出现对比度反转，像光栅常数仍为d = 4，线宽为 a’= 3，见下图

t（x3） I（x3）

0 x3 ··· ··· ··· ··· 0 x3

（d）将一个?位相滤波器置于零频上。滤波器可表达为 exp（j ?） f x = f y = 0 H（f x，f y）=

1 f x，f y ? 0 只考虑一维情况，频谱变为 T’（f x）= T（f x）·H（f x）

= =

1414{sinc(fx)exp(j?)+sinc(1414)?(fx-1414)+sinc(1414)?(fx+1414)+???}

{-sinc(fx)+sinc()?(fx-)+sinc()?(fx+)+???}

输出平面上的复振幅为

t (x3) = ? [T（f x）·H（f x）] = -14rect(x3)+14[rect(x3)?comb(x34)] -

14rect(x34)

-1

8．4 图8-54所示的滤波器函数可表示为：

1 fx＞0 H（ff ，fy）= 0 fx＝0

-1 fx＜0 此滤波器称为希尔伯特滤波器。

证明希尔伯特滤波能够将弱位相物体的位相变化转变为光强的变化。

y L1 fy L2

x fx

x'

y'

图8.54（题8.4 图）

解答：位相物可表达为

t0（x1，y1）= A·exp [ j ?（x1，y1）] 对于弱位相物有? ? 1弧度，上式近似为（忽略A）

t0（x1，y1）? 1+ j ?（x1，y1）

滤波平面得到

T（fx，fy）= ? [t0（x1，y1）]

=?（fx，fy）+ j?（fx，fy）

其中 ?（fx，fy）= ? [?（x1，y1）]。 经希尔伯特滤波器，频谱面后的光

分布为T’（fx，fy）= T（fx，fy）·H（ff ，fy）

j?（fx，fy） fx ? 0

= 0 fx ? 0 - j?（fx，fy） fx ? 0

像平面光场复振幅为 （以下无把握） t’（x3，y3）= ? -1[T’（fx，fy）]

j?（-x3，-y3） x3 ? 0

= 0 x3? 0 - j?（-x3，-y3） x3? 0

光强分布为 I = t’· t’?

-? 2（-x3，-y3） x3 ? 0 = 0 x3? 0

?（-x3，-y3） x3? 0

（此结论和于美文书上的答案不一样，建议取消此题）

8．5 如图8-55所示，在激光束经透镜会聚的焦点上，放置针孔滤波器，可以提供一个比

较均匀的照明光场，试说明其原理。

针孔 L 激光器

图8.55（题8.5 图）

8．6 光栅的复振幅透过率为

t（x）= cos 2πf0 x

把它放在4f 系统输入平面P1上，在频谱面P2上的某个一级谱位置放一块λ/ 2位相板，求像面的强度分布。

解答：将复振幅透过率函数变换为

t（x）= cos 2πf0 x = [1+cos 2πf0 x ] / 2 其频谱为

T（fx）= ? [t（x）] = =

1212 2

?（fx）+ ?（fx）+

121? [cos 2πf0 x] ?（fx- f0）+

144?（fx+ f0）

其中第一项为零级谱，后两项以次为+1级和-1级谱。设将λ/ 2位相板放在+1

级谱上，其透过率表达为

H（fx）= exp（jπ） 则频谱面P2后的光振幅变为 T’= T·H

= =

1212?（fx）+ ?（fx）-

1414 ?（fx- f0）·exp（jπ）+ ?（fx- f0）+

1414?（fx+ f0）

?（fx+ f0）

像平面光场复振幅为

t’（x）= ? -1 [T’] =

12 -

1214exp（j2πf0x3）+

1214exp（-j2πf0x3）

= - j sin（2πf0x3）

像平面强度分布为

I = ? t’（x）?2 = t’（x）· t’（x）? = =

1414[1- j sin（2πf0x3）][1+ j sin（2πf0x3）] +

14 sin（2πf0x3）

2

像平面得到的仍是一周期函数，其周期缩小1倍，振幅减小4倍，本底也有所变化，

并且出现图形的横向位移，位移量为1/2周期。

8．7 在用一维正弦光栅实现两个图象相加或相减的相干处理系统中，设图象A、B置于

输入平面P1原点两侧，其振幅透过率分别为：tA（x1- l，y1）和 tB（x1+ l，y1）；P2平面上光栅的空间频率为f0，它与l的关系为：f0 = l /λf，其中λ和f 分别表示入射光的波长和透镜的焦距；又设坐标原点处于光栅周期的1/4处，光栅的振幅透过率表示为：

G(x2,y2)?1????j(2?fx?)??1?exp02?2?2??

???exp?j(2?fx?02?2? 试从数学上证明：

?)?? 1）在输出平面的原点位置得到图象A、B的相减运算；

2）当光栅原点与坐标原点重合时，在输出平面得到它们的相加运算。 8．8 如何实现图形O1和O2的卷积运算？画出光路图并写出相应的数学表达式。 解答：第一步，制作O1的傅里叶变换全息图，光路如下： (x1，y1) L (x2，y2) O1 R -b H f f 全息图H的透过率为 tH = | ?1 |2 + R02 + R0?1（fx，fy）exp_[ -j 2πfx b]

+R0?1?（fx，fy）exp_[ j 2πfx b]

其中?1= ? [O1]，R0为平面参考波的振幅，b为参考点源的横向位移量。 第二步，进行卷积运算。在4f系统中，将O2置于输入平面（x1，y1），全息图置于频谱平面（x2，y2），如图

x1,y1 x2,y2

L1 H L2 O1? O2 O2(x1,y1) O2几何像

O1? O2

f1 1 ' 2 fff ' 2x3,y3

频谱面后的光场为 UH'= ? [O2]·tH

= ?2·{| ?1 | + R0+ R0?1（fx，fy）exp_[ -j 2πfx b]

+R0?1?（fx，fy）exp_[ j 2πfx b]}

输出平面光场为

O2? ? -1[ tH]

= R02O2 + O1 ?O1? O2 + R0O1(x3-b，y3)?O2 + R0O1?(-x3-b，-y3)?O2

式中第三项即为O1 和O2的卷积运算，位置在x3 = b处。

8．9 在4f系统中用复合光栅滤波器实现图象的一维微分 ?g / ?x ，若输入图象g在x方向的宽度为l，光栅频率应如何选取？

解答：设复合光栅的空间频率为f0和f0+?，则滤波的结果使像平面上得到两套物的三重

像，两个正一级像的位相差等于π，它们离零级像的角间距?1、?2分别由下式确

定

sin ?1 = f0?， （1）

sin ?2 =（ f0 +?）? （2）

因而正一级像离零级像的线间距分别为

l1 = sin ?1·f （3） l2 = sin ?2·f （4）

其中f是透镜焦距。分析可知，得到微分结果的条件是

l1 - l2 ? l / 2 （l为物的宽度） （5）

将（1）、（2）两式代入（3）、（4）两式，再代入（5）式，得到

l1 - l2 = ? ·? f ? l / 2

2

2

? ?

l2?f （6）

因而复合光栅的空间频率差应满足（6）式关系，才能得到微分图像。

8．10 用4f系统通过匹配滤波器作特征识别，物g（x，y）的匹配滤波器为G*（fx，fy），

当物在输入平面上平移后可表示为g（x - a， y - b），求证此时输出平面上相关亮点的

位置坐标为xi = a，yi = b。

8．11 用一个单透镜系统对图象进行θ调制假彩色编码，如图8-55所示。已知调制物Om

的光栅空间频率为100线/mm，物离透镜的距离为20cm，图象的几何宽度为6 × 6cm，

试问透镜的孔径至少应多大，才能保证在频谱面上可进行成功的滤波操作。（工作波长范围为650.0—444.4nm）。

L 频谱面 O' 白光 Om d f 图8.56（题8.11图） 解答：

设：f0 = 100线/mm，d = 20cm，a×b = 6×6cm，?max= 650．0nm，?min= 444．4cm

求：透镜最小孔径 ?min

解：调制物Om的最大线度为

2l =（a2+b2） = 6√2cm

l = 3√2cm

欲在频谱面上进行成功的滤波操作，必须使所有物点的一级衍射波都能进入透 镜，最大衍射角θmax应与?max相应，即

sinθ由几何关系得到

max / f0 = max

1/2

? sinθmax = [（? / 2）- l ] / d 所以有

? = 2 [d·f0·?max + l ]

代入数据，得 ? = 110．85mm ? 111mm

答：透镜孔径至少应达到111mm，才能保证在频谱面上进行成功的滤波操作。

第九章习题解答

9-1. 用白光再现彩虹全息时，如果彩虹全息有实狭缝象，在狭缝实象处观察全息图，人眼将能观察到单色的全息象，试分析人眼在狭缝前后位置时的全息象的颜色分布情况。如彩虹全息再现的是虚狭缝，再分析人眼观察到的全息象情况。 答：

N B P M A 在图示的情况下，物的两个端点为A和B点，它们被全息记录在一条线区域上，当白光再现时，这一区域的衍射光是色散的，长波长的衍射角较大，而短波长的衍射较小，。按图示的光路结构， A点的长波长沿AM方向衍射，短波长沿AN方向衍射，B点的长波长沿BN方向衍射，短波长沿BM方向衍射。假设沿AP和BP方向衍射的波长相同，那么人眼在P处观察将看到单色象，当眼睛靠近全息图时，将看到象的上方偏蓝，而下方偏红，反之则相反。

对于虚狭缝的情况，如上图所示，P点是某一衍射波长的虚狭缝，A和B两点是两线全息图，象上的两点与它们对应，AM是线全息图A最短波长的衍射方向，BM是线全息图B的最长波长衍射方向。显然，眼睛在M点观察，将能看到A、B之间的所有象点，但它们的颜色呈光谱色分布，在图示情况下，上部是紫色，下部是红色。眼睛观察到的象的范围由眼睛离全息图的距离决定，离得越远，观察到的范围越大。

9-2. 用白光点光源再现彩虹全息时，人眼将能观察到由光谱色组成的单色象。如果用白光线光源作为再现光源，线光源的扩展方向与狭缝方向垂直，这时观察到的是消色差的黑白象，试解释其原因。 答：线光源可以看成由无数个点光源组成，每一个点光源都按光谱色排列形成一组彩色狭缝，线光源上不同点形成的狭缝的位置各不相同，它们在与狭缝垂直的方向上平移。这无数个狭缝相互迭合在一起，使人眼在该处观察时，无数个不同波长的再现象重合在一起，这也就形成了消色差的黑白象。

9-3. 在一步法彩虹全息记录光路中，物的大小为10cm，人双眼的瞳孔间距为6.5cm，透镜的孔径为20cm，对物体1：1成像，如狭缝距全息图30cm，要求人双眼能同时看见完整的象，试计算成像透镜的焦比。 解：

设透镜L对物1:1成象，由上图可见，只有人眼在POQ三角区域内才能观察到整个物体的象，显然狭缝应在使MN大于等于双眼间距处。设透镜孔径为D，焦距为f，物体AB大小为a，人眼瞳孔间距为d，狭缝距象的距离为L。由图中的几何关系可以得到

2f A D O′ B L O M d N Q P P A B Q M Q′ OO??AB?LAB?MN?aLa?d

D?a2f?aOO?

所以

Df?2a?dLDD?a

将D=20cm，a=10cm，d=6.5cm，L=30cm代入上式，得到D/f=2.2，这是一个很不切合实际的数据。实际上用一步彩虹全息是不可能获得大观察范围。

9-4. 在用横向面积分割法制作彩色彩虹全息母板的方法中，已知下列条件：三色光的中心波长分别为645.2nm、526.3nm和444.4nm；第一步记录时被记录物中心位于建在母全息图HM的坐标系的z轴，物体距HM30cm；第二步记录时参考光为平行光，入射角30?；白光再现时入射光是入射角为45?的平行光，三色再现狭缝位于z轴；设两次记录的波长均为442nm。试据以上条件，确定HM上H1、H2、和H3的位置。如果每个狭缝的光谱带宽为10nm，试确定狭缝宽度。 解：

x r SR SG SB HM l H H c O z O z x (a) (b)

图(a)和图(b)是本习题的记录和再现光路。在图(a)中设三个狭缝图象中心的入射角分别是θR、θG和θB，参考光入射角是θr，图(b)是共轭再现光路，再现光入射角是θc，图象中心三狭缝的不同色光的衍射光方向相同，设衍射角为θi。由光栅方程

sin?i?sin?c???0?sin?o?sin?r?，

得到共轭再现时的θo，

sin?o?sin?r??0??sin?i?sin?c?，

式中λ0和λ分别为记录波长和再现波长，θo的下标o分别对应于R、G和B。按题中条件及图(a)的坐标系，

?r??30?,?c?135?,?i?180?

以λ0=442.0nm和λ=444.4nm、526.3nm、645.2代入上式，得到

sin?B??0.203,?B??11.73?

sin?G??0.094,?G??5.38?

sin?R?0.016,?R?0.89?

以l=-30cm计算，得到

xB?ltan?B?6.23cm

xG?ltan?G?2.82cm

xR?ltan?R??0.47cm

由于狭缝有一定宽度而引起的带宽可对光栅方程微分得到，即 ??o??0sin?i?sin?c?2cos?o??。

由此可得到狭缝宽度

wB?0.97cm

wG?0.34cm

wR?0.22cm

9-5. 采用图示的双狭缝彩虹全息记录光路，可以在同一张底片上记录两个物体的彩虹全息图，记录的方法步骤是：第一步，用挡板挡住S2，用S1对物体O1曝光；第二步，用挡板挡住S1，用S2对物体O2曝光。然后将显影的全息图用白光照明，人眼在不同位置即可看到不同物体的再现象。

图 习题9-5 答：（1）以R的共轭光再现全息图，在原狭缝处得到实狭缝，此时在实狭缝处观察，观察

O S1 S2 （1）画出再现狭缝实象的示意图，说明再现象的特点。

（2）解释多狭缝彩虹全息图作多目标存储和假彩色编码的原理。

L R O? 到的是赝象，即凹凸与原记录物体相反的象。

（2）记录多目标时，在透镜处放置多个狭缝，分别对应于不同的物体进行全息记录，再现时用单色光，眼睛位于不同的狭缝处，将能见到不同的物体。进行假彩色编码时，将待编码的物体分别置于透镜的物面，在透镜处分别遮挡不同的狭缝，记录在同一张全息图上，再现时用白光再现，不同编码图象按不同颜色重叠在一起，结果形成假彩色图象。

9-6. 有人提出用蓝色单色激光以反射全息的方式也可以记录二维彩色照片，方法是：将三张分色片分别置于记录介质一侧，并用散射屏照明，在记录介质另一侧分别用三个不同角度的参考光入射，在同一记录介质上记录三张分色片，当白光以某一角度再现全息图时，三分色片将分别被三原色再现，呈现彩色图象。试说明其原理，并作相应的三参考光入射角设计计算。（提示：用三棱镜与记录介质用匹配液匹配的方法可以增加参考光入射角；用布喇格条件和光栅方程进行设计，参考第5章和第7章体全息部分）。

答：反射体全息有很好的波长选择性，利用布喇格条件，可以按所要求的衍射方向和衍射波长推算出两束相干记录光束的入射角，这也就是本习题提出的记录彩色二维图片的原理。体全息中记录光与再现光之间的布喇格条件关系是

cos?o?1sin?o?1?cos?r?1sin?r?1?cos?c?2sin?c?2?cos?i?2?0

???sin?i?2?0

式中角度量θ的下标o、r、c、i分别表示记录时的物光、参考光和再现时的再现光、衍射光，它们均是介质内的量，λ1和λ2分别表示记录光和满足布喇格条件的再现光波长。按题意，θc、θi和λ2均由使用条件决定，利用上式求出θo和θr，从上式可以求出

?o??i??c2?i??c2?sin?1??1?i??c?sin??2?2??1?i??c?sin??2?2?? ????。 ???r??sin?1在本习题中，θc取207.7°（在记录介质中角度，空气中为225°），θI取0°，λ2则取三原色波长，利用上式可求出在记录介质内的物光和参考光入射角。如取三原色波长为450nm、530nm和630nm，得到参考光和物光入射角如下表 波长λ2(nm) 450 参考光入射角θr 物光入射角θo 202.75° 4.95° 530 226.83° -19.13° 630 238.96° -31.26° 显然，如果在参考光入射的一边不采取其他措施，530nm和630nm的介质内入射角是无法实现的。采用45°直角三棱镜匹配的光路如下图所示。

x 散射光照明

??r z 二维照片掩膜 记录干板 三棱镜

把上表的数据折算成空气中的入射角如下表所示（表中??r为在直角三棱镜斜边的入射角），如果物光由散射角非常大的散射屏照明，表中的物光散射角没有什么意义，因照明光中总有一部分物光满足表中的入射条件。实际上大散射角的照明会降低衍射效率，采用有一定散射角的散射屏更适合物光照明，散射的中心光线因满足下表中物光入射角。作全息记录时，在全息干板前分别换三原色掩膜，对应于它们的光路参数作三次曝光。在对全息干板进行后处理后，用白光以45°角入射，就能观察到彩色二维照片。

波长λ2(nm) 450 530 参考光入射角??r 物光入射角??o

35.14° 7.53° -2.78° -29.88° 630 -21.51° -52.07° 第十章 习题解答

10.1 试比较被动三维传感和主动三维传感系统的原理、系统结构、适用范围和优缺点。（思考题）

10.2 在三角测量法中通常采用的三种坐标系统如图10.6所示。试推导三种坐标关系中，物体的距离或高度z与测量变量△x之间的关系式，即三角测量法中的测量方程。

解答：

（1）投影光轴与成像光轴平行。所构成的物三角形和像三角形是相似的直角三角形, 测量方程是。

z?bh?x

式中：b是物三角形的基线，h是像三角形的高度，△x是像三角形的基线，z是物体的距离或高度。

（2）投影光轴和成像光轴相交。θ是投影光轴与成像光轴的夹角，O是两光轴交点并作为物体高度计量的原点，I和I′是成像系统的入瞳和出瞳，线阵探测器与成像光轴垂直，与I′点的距离为f;当物距l较大时，f近似地等于成像透镜的焦距。由图中所示的几何关系可以导出。

z?l?xsin??f?cos???x

(3) 投影光轴和成像光轴相交,探测器基线与成像光轴成一倾角β，当满足Scheimpflug条件，即满足关系 tg??k?tg? 时，待测距离z和可测变量△x之间的关系式为

(OI?f)sin???xfsin??cos?sin???xz?

10.3 为什么说激光散斑对三角法测量精度具有重要影响，试解释公式(10.16)和(10.17)的物理含义，并说明如何提高激光三角法测量精度。

解答：由于物体表面的微观起伏的不确定性，在探测器上的像点的散斑分布也是不确定的，这种不确定性引起的光点中心的定位误差，因此激光散斑对三角法测量精度具有重要影响。公式(10.16)和(10.17)表明，这种不确定性与透镜的数值孔径、激光的波长和散斑的对比度有关。通过增大透镜的数值孔径，减小波长C，降低散斑的对比度可以提高激光三角法测量精度。

10.4 在位相测量剖面术中，由于变形光栅像与传统的干涉条纹图相类似，因此变形光栅像有时又被称为“干涉图”。 在干涉计量中，光波长被作为度量微观起伏的尺度，而在位相测量剖面术中与投影条纹间距有关的“等效波长”被作为度量三维宏观面形的尺度。试比较这两种方法在物理概念上和条纹处理方法上的异同性。（思考题）

10.5 由于实际得到的位相数据是一个二维的采样点阵列，所以位相展开应针对二维进行。模仿一维位相函数的位相展开过程，推导二维截断位相函数φw(i,j) 展开过程的数学表达式.

解答：可以首先沿二维数据阵列中某一列进行位相展开，然后以该列展开后的位相为基准，沿每一行进行位相展开，得到连续分布的二维位相函数。例如首先沿二维数据阵列中第一列展开，然后以该列展开后的位相为基准，沿每一行展开。

φu(i,1)?φw(i,1)?2πnini?int??φw?i,1??φw?i?1,1??/2π?0.5??ni?1n0,1?0φu?i,j??φw?i,j??2πnj?φu?i,1?ni,0?0nj?INT??φw?i,j??φw?i,j?1?2π?0.5???nj?1

10.6采用远心光路的PMP系统如图10.22所示。设图中?=300 ，?’=00，在参考平面上看到的投影正弦光栅是等周期分布的，其周斯P0=5mm，求该系统的等效波长。如果系统对条纹位相的测量精度为2?/100, 求系统的测量精度。试讨论提高系统的测量精度的方法。

解答：等效波长 ?e?P0/tg?＝8.7mm,系统的测量精度为0.087mm. 减小等效波长?e可以提高系统的测量精度。

10.7 位相测量轮廓术和傅立叶变换轮廓术是基于三角测量原理，试比较调制度测量轮廓术与上面两种方法在原理上的区别，并比较三种方法的测量精度。（思考题） 10.8

飞行时间法(TOF)是基于直接测量激光或其他光源脉冲的飞行时间来确定物

体面形的方法。图10.43是采用位相检测技术的TOF系统框图，对时间的测量可以通过对调制光波的位相测量来实现。光束经9MHz的调制器调制后投射到物接收的信号经9MHz的滤波器后与基准信号比较，然后从位相变化计算出距离的变化。假定位相的测量精度为2?/100，求系统的测量精度。如果保持位相的测量精度不变，光束的调制频率提高到90MHz, 系统的测量精度是多少。

解答：系统的测量精度为0.33m。如果光束的调制频率提高到90MHz, 系统的测量精度提高到0.03m。

第十一章习题

11.1 试证明任意两个相互统计独立的随机变量之间相关系数为零。 答：参阅《统计光学（基本概念个习题）》P21。

证明 设UV为两相互统计独立的随机变量，由于其相互独立，UV的联合概率密度函数可为边缘概率密度函数的乘积：

pUV(u,v)?pU(u)?pV(v)

因而U和V的相关函数为

?UV?????uvpUV(u,v)dudv

???uv

根据定义，U和V的协方差则为

??CUV??UV?uv?0

故相关系数?(??CUV/?U?V)为零。

11.2 若N个微小随机相幅矢量

?Nakej?k之和中每一个幅值

akN及相位?k都相互独立；

所有的ak具有相同的概率分布，数学期望与二阶矩分别为a和a?；随机位相?k均布于???,??区间内。试计算：（1）当N趋近于无穷大时这N个随机相幅矢量之和的实部和虚部的均值与方差及相关系数；（2）实部和虚部的联合概率密度函数并绘出复平面上等概率密度曲线图。如果随机位相?k均布于????????,???区间内计算结果及函数

图象有何变化？

答：参阅《统计光学》P42-45及《统计光学（基本概念个习题）》P30-34。

解：设和矢量为a（用复数表示），即

a?aej??1NN?ak?1kej?k

其实部与虚部分别为

r?Reae?j???1N1NN?ak?1Nkcos?k

i?Imae?j????ak?1ksin?k

根据N个微元矢量的相互独立性可得

r?1N1NN?ak?1Nkcos?k

i??ak?1ksin?k

对于均布与???cos?k?sin?k?????,?之中的Φk有 22?1cos?kd?k?2????/2/2?1?

????/2/2?sin?kd?k?0

所以

r?1NNa?2?2Na,i?0

??2为了要计算?2及?i，首先计算二阶矩r2及i2，

r2r?1N1NNNk??ak?1n?1NNancos?kcos?n

i2???ak?1n?1kansin?ksin?n

式中

??cos?kcos?n,cos?kcos?n??2cos?k,??4?,2?????

?/211??(1?cos2?k)d?k,???/2?2???4,当n?k???2???1,当n?k??2??sin?ksin?n,sin?ksin?n??2??sin?k,4?,2??? ???/211??(1?cos2?k)d?k,?????/2?2?0,当n?k???1?,当n?k?2因而

r2?a22?4?2(N?1)a,i22?a22

故有 ?2r?r2?r2?a22?4a2?,?2i?i?i22?a22

要计算相关函数，还需要计算实部与虚部之间相关函数，

ri?1NNNk??ak?1n?1ancos?kcos?n

其中

??cos?ksin?n,cos?ksin?n????cos?ksin?n,0,当n?k?? ??/21sin2?d?,当n?kkk????/22???0最后可得出

??ri?ri?r??0

i这说明实部与虚部之间仍是不相关的。

因为N是个很大的数，由中心极限定理知实部r与虚部i分别为高斯随机变量，即

????exp??????2???2Nr?a??????????? 22?a4a????2??2?2??????pR(r)?2?1a22?4a2?pI(i)???2?1i???exp??? 22a??a22???2??2因为r与i不相关，其联合概率密度函数仍为高斯型，可写作 pRI(r,i)2????2N???r?a???????1?exp???22228a??2aa4aa?(?)??2?22????2?????2?i? ??2??a?????图1.12为r与i的联合等概率密度曲线简图，若图中各椭圆之长短轴分别为an与bn，则有

an:bn?a:(a?228a2?)

图 pRI(r,i)的等概率密度曲线简图

11.3 若图11.3中光学成像系统的脉冲响应函数为h?x,y?，光瞳函数为P??,??，试证明全

息再现的变形前后两波面散斑场之间的相关因子的表达式（11.27b）可以转化为（11.13）式中复自相干度的表达形式。

证明：变形前后两波面散斑场之间的相关因子的表达式(11.27b)为 ???h(r?r)h?r?r?h(r?r)*0002?dr0?d2(r)dr0

该式右边分子是脉冲响应函数的自相关函数，根据自相关定理，可以表示成脉冲响应函数的傅氏变换的模平方的傅氏反变换，即

?h(r?r0)hr?r0?d2(r)dr0???P(ξ,η)exp*?????π?j(ξΔx?ηΔy)dξdη ???λz?式中Δx和Δy是d2(r)的两个分量。当Δx和Δy即d2(r)为零时，式 (11.27b) 右边分子变为分母的形式，因而有 ?h(r?r0)h*?r?r0?dr0???P(ξ,η)?dξdη

将上述两式代回相关因子的表达式中，（11.27b）即可转化为（11.13）式中复自相干度的表达形式。

μ?Δx,Δy????P(ξ,η)?exp??π?j(ξΔx?ηΔy)dξdη???λz?? （11.13）

??P(ξ,η)

dξdη11.4 在如图11.2的全息干涉记录光路中置于物面处的被测物在激振器驱动下垂直于物面

进行稳态振动，照明参考光与物面法线方向夹角可近似为不变的30°。用时间平均法记录下的全息图处理后放到图11.3的全息干涉再现光路中，再现参考光与物面之间几何关系与记录光路中完全相同。如果在再现时观察到的两条节线之间正好有五条暗纹，试问稳态振型的振幅最大值是多少？假设记录和再现的光波长为633纳米。 答：因为在时间平均法记录下的全息图重现像的振动图样中，不运动的区域即节线处显示最亮的条纹，随着条纹级次的增加，亮条纹的强度逐渐下降。如果在再现时观察到的两条节线之间正好有五条暗纹，稳态振型的振幅最大值处是第三条暗纹。因此

?πλBmax??.????

故振幅最大值为

Bmax??.?????

11.5 分离再现外差全息干涉方法中，如果被测物的空间频谱分布在?f?与f?之间，而且

两束参考光均为平面波，试问参考角?，?（参阅图11.7）需满足什么条件才能保证外差全息干涉图与噪声项完全分离？

答：分离再现外差全息干涉方法中，复原到全息记录面上的处理好的全息图振幅透过率可表示为

λ?π??.???μ

τ(rh)?τ??βR??A?(rh)???R??A?(rh)??

再现时用调制频率为?0与?0???的两束光照明全息图，则由全息图射出的光场分布由三部分组成，第一部分为再现光场中的直透项，不会产生干涉现象，第二部分为再现出的原物光场及交叉再现光场，第三部分为再现出的原物光场的共轭光场及交叉再现光场。因为第二部分和第三部分光场分布在第一部分直透光场的两侧，只要第一部分直透光场与第二部分光场不重叠，第三部分光场就不会与第二部分再现出的原物光场重叠，只需要保证第一部分直透光场与第二部分光场不重叠就可以了。如果参考角?和?的两束参考光分别记录

A?(rh)和A(rh)，两束参考光的空间频率可表示为

f??sinαλ 和 f??sinβλ 第一部分的直透项中对应于A?(rh)A?(rh)??的空间频谱分布在??f??f?与?f??f?之间，对应于

的空间频谱分布在??f??f?与?f??f?之间，其他直透项对应着频域的δ函数。

第二部分中再现出的原物光场A?(rh)和A(rh)的空间频谱，分布在?f?与f?之间。第二部分中再现出的两个交叉再现光场的空间频谱，则分布在?f??f??f?与f??f??f?之间和

?f??f??f?与f??f??f?之间。不失一般性，假设?

场和两个交叉再现光场不重叠，则要求 f???f??f??f?

欲使第二部分中再现出的原物光场和第一部分中对应于A?(rh)重叠，则要求

f????f??f? 因而需满足

α?sin-1?和A?(rh)?的直透项不

?λf?

sinα???? β?sicλ??f???

λ??

11.6 相移全息干涉技术中，已知三次相移?(t)分别取为

??，

???，

???，相应光强测量

结果为I?(ri)I?(ri)I?(ri)，试计算(11-35)式中r处的相对初位相Δ(ri)。 答：相对初位相为：Δ(ri)?tg

11.7 在图11.9的双光楔剪切散斑干涉记录光路中，光楔玻璃折射率n0为1.5，楔角?为

0.05°，物距l0为500mm，照明使用波长为633纳米的氦氖激光且入射光与观察方向

z夹角?为30°，而??x,y?及u?x,y?分别为?x,y?点沿z及x方向变形产生的位移

??I??I?I??I?。

分量。第一次拍摄两次曝光剪切散斑干涉图后，使物体绕z轴旋转180°再拍摄一次。对应同一?x,y?点在两次拍摄的两次曝光剪切散斑干涉图中的条纹级数分别为2.5和3，试问?x,y?点沿z及x方向的应变为多大？ 答： 首先用式(11.41)计算剪切量

δx??l?(n???)α?2?500??1.5-1???.????.??????.????

根据式(11.43)可以计算出x方向的面内应变为

?u?N?(??cosψ?)?N?(??cosψ?)???x?δx??sinψ?(??cosψ?)?sinψ?(??cosψ?)λ?.???????.??????????000 ???.????cos?-30sin?300????cos?-300???????cos30?

???sin?-30????cos30?0???.?????????再根据式(11.42) 计算出z方向的离面应变为

?.???.???????? ?ω?x?

?.??????.???.????????????.?????.???????

11.8 在图11.16所示的全场滤波光学系统中，透镜焦距为500毫米，使用波长为633纳米

的氦氖激光照明，放大率为1，滤波孔中心离光轴15MM，滤波孔直径有1.5、3、6MM三种。如果被测物是一个直径55MM的旋转圆盘，两次曝光之间转动角度为0.35°，记录两次曝光散斑图时的放大率也为1，而且相对口径足够大。试问用三种不同口径

滤波时得到的全场滤波条纹图上可以观察到条纹的区域直径分别为多大？ 答：该全场滤波光学系统可测量的面内位移的最大值分别为 d2??.??λfDf??.????.???????.?????????max???.???mm。

d2max???.????.????????.????.?????????????.???mm。

?? d2max????????.???mm。

因为被测物是一个直径55MM的旋转圆盘，两次曝光之间转动角度为0.35°，面内位移最大值为0.257mm、0.129mm、0.064mm的区域都是圆，其直径分别为 φ??d??????.??????.????.????????.??????.?????.??mm φ???.????????.??????.?????.??mm φ???.????????.??????.?????.??mm

本文来源：https://www.bwwdw.com/article/2wy7.html