陈家璧 光学信息技术原理及应用习题解答章

1、第七章 习题解答1. 某种光盘的记录范围为内径80mm,外径180mm的环形区域,记录轨道的间距为2um.假设各轨道记录位的线密度均相同,记录微斑的尺寸为0.6um,试估算其单面记录容量. (注: 内、外径均指直径)解： 记录轨道数为 单面记录容量按位计算为 bits = 17 Gb. 按字节数计算的存储容量为 2.1GB.2. 证明布拉格条件式（7-1）等效于（7-17）式中位相失配d = 0的情形， 因而（7-18）式描述了体光栅读出不满足布拉格条件时的位相失配。证明:将体光栅读出满足布拉格条件时的照明光波长(介质内) 和入射角 (照明光束与峰值条纹面间夹角)分别记为l0和0, 则根据布拉

2、格条件式（7-1）有: 2Lsin0= l 0 其中L为峰值条纹面间距.对于任意波长la (空气中) 和入射角r (介质内), 由(7-17)式, 位相失配 d 定义为: 其中n0为介质的平均折射率, K= 2p/L为光栅矢量K的大小，f为光栅矢量倾斜角，其值为 ，qr 为再现光束与系统光轴夹角 (参见图7-9)当 d = 0 时,有 即: l为介质中的波长. 由于角度恰为照明光与峰值条纹面的夹角, 以上结果亦即布拉格条件2L sin = l.当读出光偏离布拉格角o和布拉格波长lo的偏移量分别为D和Dl时，有 利用布拉格条件式(7-17), 以及D和Dl很小时的近似关系 cosD1 和 sin

3、DD, 立即可得: d =DqKsin(f-q0) - DlK2/4pn0 即(7-18)式原题得证。3. 用波长为532nm的激光在KNSBN晶体中记录非倾斜透射光栅,参考光与物光的夹角为30o(空气中).欲用波长为633nm的探针光实时监测光栅记录过程中衍射效率的变化,计算探针光的入射角.(假设在此二波长晶体折射率均为2.27) 解： 532nm为空气中激光波长记作la1, 在晶体外的入射角为a1，其在晶体中波长为l1, 入射角为1；633nm为空气中激光波长记作la2, 在晶体外的入射角为a2，其在晶体中波长为l2，入射角为a2。本题中涉及非倾斜光栅, 光束入射角即为与光栅峰值条纹面的夹

4、角, 按题意, 则波长532nm和633nm激光应分别满足布拉格条件： 晶体中： 2Lsin1=l1 2Lsin2=l2 (1)由折射定律，换算成空气中角度和波长为： 空气中： 2Lsina1=la1 2Lsina2=la2 (2)由（2）式得： a2 = arcsin (la2 ´ sin15°/la1 )= arcsin (633 ´ sin15° /532 ) = 17.936 °故探针光的入射角应为17.936°。4. 为了与实验测量的选择角相比较,需要有体光栅在空气中的选择角的表达式. 试对小调制度近似(n<<1

5、),导出一个计算非倾斜透射光栅空气中的选择角的表达式 (所有角度均应为空气中可测量的值).解：注意我们将对应着h-x 曲线的主瓣全宽度定义为选择角, 体光栅晶体中选择角表达式为： (1)n <<1时，对非倾斜透射光栅，有： (2)设空气中参考光入射角为qro, 选择角为DQo. 由折射定律有 sin(DQo+qro) = nsin(DQ+qr)(3)展开为 ： sinDQocosqro + cosDQosinqro= n (sinDQcosqr + cosDQsinqr) (4)因为DQo和DQ很小，有如下近似：cosDQocosDQ/21, sinDQoDQo, sinDQDQ.

6、 因此(4)式可化简为： DQocosqro+ sinqro= n(DQcosqr+ sinqr)由折射定律，有sinqro = nsinqr，可得： DQo = nDQcosqr/ cosqro = nla cosqr / (nd ú sinqrú cosqro) = 2la(n2-sin2 qro)1/2 / (d sin2qro)此式可作为空气中选择角的表达式。当sin2 qro<< n2时,还可进一步简化为:以最常用的铌酸锂晶体为例, n=2.2-2.3, 当qro < 45°时, 用估算空气中的选择角, 误差只有5%左右. 5. 铌酸锂

7、晶体折射率 n =2.28, 厚度d = 3mm, 全息时间常数之比E/W = 4, 饱和折射率调制度Dnmax=5´10-5, 用l = 532nm的激光在晶体中记录纯角度复用的全息图, 物光角度取为s=30°, 参考光角度范围r = 20°-40°. 若要求等衍射效率记录且目标衍射效率设定为10-5, 试分析影响存储容量的主要因素.为了提高存储容量, 应当在哪些方面予以改进?解：本题仅涉及纯角度复用技术, 且无页面容量的数据, 故主要讨论每个空间区域内复用存储的数据页面数即角度复用度. 影响存储容量的主要因素有：（1） 有限的角度选择性及实际选择角增

8、宽（2） 光学系统对存储容量的限制（3） 噪声对存储容量的限制本题可近似为准对称的透射式光路，以平均的参考光角度r=s=30° 并利用第4题的结果, 可估算出平均选择角为:0.0535°由光学系统决定的参考光入射角范围Q=20°, 故允许的角度复用度为： Ma » Q/ DQ =20/0.0535 = 374 由于系统存在噪声, 要求有一定的目标衍射效率, 而记录材料具有有限的动态范围, 因而根据(7-77)式,角度复用度限制为： Ma = = 4´3.14´5´10-5´3 / (532´10-6

9、80;cos30°´10-5/2) =1293由以上结果可见，本题中存储容量主要受到有限的角度选择性和光学系统有限的孔径角的限制。要提高存储容量，需首先增大晶体厚度, 适当增大写入光的夹角和光学系统的孔径角, 采用较短的激光波长. 6. 用作组页器的空间光调制器为(24´36)mm2的矩形液晶器件, 含有480´640个正方形像元. 用焦距为15mm的傅里叶变换透镜和633nm激光记录傅里叶变换全息图, 问允许的参考光斑最小尺寸为多少?解：空间光调制器的复振幅透过率可描述为二维rect函数阵列. 当记录傅立叶变换全息图光路中包括空间光调制器时，4f系统中

10、频谱面上的频谱是与像元间距有关的sinc函数。记录物体基本信息的必要条件是参考光光斑至少包含物体的零频和正负基频频谱。正方形像元间距为下式所求，二者取其小： 24/480 = 0.05mm 36/640 = 0.056mm由于黑白相邻的两个像元构成物面上最小可分辨的周期, 故方波条纹的最高频率为fmax = 1/(2´0.05) = 10 mm-1.在焦距为f的傅里叶变换透镜的频谱面上, 物体+1级频谱分量的间距为： DH = 2lf fmax = 2´633´10-6´15´10= 0.19 mm所以，允许的参考光斑最小尺寸是直径为0.19

11、mm的圆斑。7. 用一种高级语言编写计算机程序, 计算在光折变晶体中角度复用存储100幅全息图的曝光时间序列. 计算用到的参数为: 饱和衍射效率0.5, 目标衍射效率10-5, 写入时间常数60s, 擦除时间常数700s. 透射式准对称光路.解: 对于透射式准对称光路，由(7-24)和(7-26)式, 光栅衍射效率表达式为： h = sin2n = 相应的折射率调制度为 Dn = , 则饱和折射率调制度Dnsat与饱和衍射效率hsat的关系为: Dnsat =而相应于目标衍射效率h0 的折射率调制度Dn0为: Dn0 = 由 D n0 = Dnsat(1-) 可得出 tw = (1)tw即为不

12、会经受擦除效应的最后一个全息图(N=M)的曝光时间. 由(7-72)式, 对于足够大的N, 曝光时序的近似递推公式为: tN = (2)或由(7-70)式, 更精确的递推公式为:tN = W (3)按(3)式建立起计算模型, 从(1)式出发, 可从后往前逆推出每个全息图的曝光时间. 应当指出的是: a) 当计算出的tN+1足够大, (3)式右边的对数函数可能会不收敛,导致不能得出有意义的tN. 此时有意义的曝光时间个数(M-N)即为在题设条件下可以达到的角度复用度.b) 若不指定全息图的个数, 按本题所给条件可以计算出角度复用度远大于100. 亦即本题中的100个全息图是大量全息图序列中的最后

13、100个, 因此可以使用近似递推公式(2)进行计算. 也就是说, 当目标衍射效率足够小时, 用(2)式计算也是相当精确的.习 题81利用4f系统做阿贝波特实验，设物函数t（x1，y1）为一无限大正交光栅 其中a1、a2分别为x、y方向上缝的宽度，b1、b2则是相应的缝间隔。频谱面上得到如图8-53（a）所示的频谱。分别用图8-53（b）（c）（d）所示的三种滤波器进行滤波，求输出面上的光强分布（图中阴影区表示不透明屏）。. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

14、. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . （a） （b） （c） （d）图8.53（题8.1 图）解答：根据傅里叶变换原理和性质，频谱函数为 T ( fx , fy ) = t ( x1 , y1 ) = · * · 将函数展开得T ( fx , f

15、y ) = * （1） 用滤波器（b）时，其透过率函数可写为 1 fx = + 1/ b1 fy = 0 F ( fx , fy ) = 0 fx ¹ 1/ b1 fy = 任何值 滤波后的光振幅函数为 T·F = 输出平面光振幅函数为 t（x3，y3）= -1 T·F = = 输出强度分布为 I（x3，y3）= = - C其中C是一个常数，输出平面上得到的是频率增加一倍的余弦光栅。（2）用滤波器（c）时，其透过率函数可写为 1 fx ，fy ¹ 0 F ( fx , fy ) = 0 fx = fy = 0 滤波后的光振幅函数为 T·F =

16、* 输出平面光振幅函数为 t（x3，y3）= -1 T·F = * - × * - 输出强度分布为 I（x3，y3）= | t（x3，y3）|2有两种可能的结果，见课本中图89和图810。（3）用滤波器（d）时，输出平面将得到余弦光栅结构的强度分布，方向与滤波狭缝方向垂直，周期为b，它与物光栅周期b1、b2的关系为 82 采用图8-53（b）所示滤波器对光栅频谱进行滤波，可以改变光栅的空间频率，若光栅线密度为100线/mm，滤波器仅允许 + 2级频谱透过，求输出面上干板记录到的光栅的线密度。解答：根据对81题的分析，当滤波器仅允许+ 2级频谱通过时，输出平面上的光振幅应表达

17、为 t（x3）= -1 = 其振幅分布为一周期函数，空间频率是基频的2倍。而干板记录到的是强度分布： I = = - C其中C是一个常数。答：干板上记录到的光栅频率是基频的4倍，即400线/mm。83 在4f系统中，输入物是一个无限大的矩形光栅，设光栅常数d = 4，线宽a =1，最大透过率为1，如不考虑透镜有限尺寸的影响， （a）写出傅里叶平面P2上的频谱分布表达式；（b）写出输出平面复振幅和光强分布表达式；（c）在频谱面上作高通滤波，挡住零频分量，写出输出平面复振幅和光强分布表达 式；（d）若将一个位相滤波器 exp（j） x2，y2 x0，y0 H（x2，y2）= 0 其它放在P2平面的

18、原点上，写出输出平面复振幅和光强分布表达式，并用图形表示。解答：将81题结果代入，其中b1 = d = 4，a1 = a = 1，除去与y分量有关的项，可得（a）P2平面上的频谱分布为：（b）输出平面： 复振幅 t（x3）= -1 T（fx）若不考虑透镜尺寸的影响，它应该是原物的几何像，即 t (x3) = *光强分布 I (x3) = | t (x3)| 2 = *(c)挡住零频分量，输出平面情况与81题（3）相同，即 t (x3) = *- I = | t (x3) | 2 由于a = d / 4，所以强度将出现对比度反转，像光栅常数仍为d = 4，线宽为a= 3，见下图 t（x3） I（

19、x3） 0 x3 ··· ··· ··· ··· 0 x3（d）将一个p位相滤波器置于零频上。滤波器可表达为 exp（j p） 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）= = 输出平面上的复振幅为 t (x3) = -1T（f x）·H（f x） = -* - 84 图8-54所示的滤波器函数可表示为： 1 fx0 H（ff ，fy）= 0

20、fx0 -1 fx0此滤波器称为希尔伯特滤波器。 证明希尔伯特滤波能够将弱位相物体的位相变化转变为光强的变化。 y L1 fy L2 x fx x' y' 图8.54（题8.4 图）解答：位相物可表达为t0（x1，y1）= A·exp j f（x1，y1） 对于弱位相物有f < 1弧度，上式近似为（忽略A）t0（x1，y1） 1+ j f（x1，y1）滤波平面得到T（fx，fy）= t0（x1，y1）=d（fx，fy）+ jF（fx，fy）其中 F（fx，fy）= f（x1，y1）。 经希尔伯特滤波器，频谱面后的光 分布为T（fx，fy）= T（fx，fy）&#

21、183;H（ff ，fy） jF（fx，fy） fx > 0 = 0 fx = 0 - jF（fx，fy） fx < 0 像平面光场复振幅为 （以下无把握） t（x3，y3）= -1T（fx，fy） jf（-x3，-y3） x3 > 0 = 0 x3= 0 - jf（-x3，-y3） x3< 0 光强分布为 I = t· t -f 2（-x3，-y3） x3 > 0 = 0 x3= 0 f 2（-x3，-y3） x3< 0（此结论和于美文书上的答案不一样，建议取消此题）85 如图8-55所示，在激光束经透镜会聚的焦点上，放置针孔滤波器，可以提供一个

22、比较均匀的照明光场，试说明其原理。 针孔 L 激光器图8.55（题8.5 图）86 光栅的复振幅透过率为 t（x）= cos 2f0 x 把它放在4f 系统输入平面P1上，在频谱面P2上的某个一级谱位置放一块/ 2位相板，求像面的强度分布。解答：将复振幅透过率函数变换为 t（x）= cos 2f0 x = 1+cos 2f0 x / 2 其频谱为 T（fx）= t（x） d（fx）+ cos 2f0 x = d（fx）+ d（fx- f0）+ d（fx+ f0） 其中第一项为零级谱，后两项以次为+1级和-1级谱。设将/ 2位相板放在+1级谱上，其透过率表达为 H（fx）= exp（j） 则频谱

23、面P2后的光振幅变为 T= T·H = d（fx）+ d（fx- f0）·exp（j）+ d（fx+ f0） = d（fx）- d（fx- f0）+ d（fx+ f0）像平面光场复振幅为 t（x）= -1 T = - exp（j2f0x3）+ exp（-j2f0x3） = - j sin（2f0x3） 像平面强度分布为I = ç t（x）ç2 = t（x）· t（x） =1- j sin（2f0x3）1+ j sin（2f0x3） =+ sin2（2f0x3） 像平面得到的仍是一周期函数，其周期缩小1倍，振幅减小4倍，本底也有所变化，并且出现图

24、形的横向位移，位移量为1/2周期。87 在用一维正弦光栅实现两个图象相加或相减的相干处理系统中，设图象A、B置于输入平面P1原点两侧，其振幅透过率分别为：tA（x1- l，y1）和 tB（x1+ l，y1）；P2平面上光栅的空间频率为f0，它与l的关系为：f0 = l /f，其中和f 分别表示入射光的波长和透镜的焦距；又设坐标原点处于光栅周期的1/4处，光栅的振幅透过率表示为： 试从数学上证明： 1）在输出平面的原点位置得到图象A、B的相减运算； 2）当光栅原点与坐标原点重合时，在输出平面得到它们的相加运算。88 如何实现图形O1和O2的卷积运算？画出光路图并写出相应的数学表达式。 解答：第一

25、步，制作O1的傅里叶变换全息图，光路如下： (x1，y1) L (x2，y2) O1 R -b H f f全息图H的透过率为tH = | m1 |2 + R02 + R0m1（fx，fy）exp_ -j 2fx b +R0m1（fx，fy）exp_ j 2fx b 其中m1= O1，R0为平面参考波的振幅，b为参考点源的横向位移量。第二步，进行卷积运算。在4f系统中，将O2置于输入平面（x1，y1），全息图置于频谱平面（x2，y2），如图 x1,y1 x2,y2 L1 H L2 O1H O2 O2(x1,y1) O2几何像 O1 O2 f1 f1 ' f2 f2 ' x3,y3

26、 频谱面后的光场为 UH'= O2·tH= m2·| m1 |2 + R02 + R0m1（fx，fy）exp_ -j 2fx b +R0m1（fx，fy）exp_ j 2fx b输出平面光场为 O2 -1 tH = R02O2 + O1 HO1 O2 + R0O1(x3-b，y3)O2 + R0O1(-x3-b，-y3)O2式中第三项即为O1 和O2的卷积运算，位置在x3 = b处。89 在4f系统中用复合光栅滤波器实现图象的一维微分 g / x ，若输入图象g在x方向的宽度为l，光栅频率应如何选取？解答：设复合光栅的空间频率为f0和f0+d，则滤波的结果使像平面上得到两套物的三重像，两个正一级像的位相差等于，它们离零级像的角间距q1、q2分别由下式确定 sin q1 = f0l，