近代光信息处理第5章广义傅里叶变换及其光学实现.ppt
2023/6/26,1,第五章,广义傅里叶变换及其光学实现,2023/6/26,光学信息处理,2,第五章 广义傅里叶变换及其光学实现,5.1 引言5.2 广义傅里叶变换的定义及性质5.3 广义傅里叶变换的本征函数5.4 用透镜系统实现广义傅里叶变换 的基本光学单元5.5 基本光学单元的组合5.6 用自聚焦效应光波导实现广义傅里叶变换5.7 维格纳变换,2023/6/26,光学信息处理,3,5.1 引 言,二维傅里叶变换(u,v)=Fo=-o(x,y)exp-i2(ux+vy)dxdy 可以用光学系统近似实现 在本章中将研究当物体到透镜的距离d1及输出图像到透镜的距离d2不等于透镜的焦距f 时透镜或透镜系统对输入图像的变换,2023/6/26,光学信息处理,4,5.1 引 言,研究表明,d1和d2 满足一定的条件时,输出平面上将出现 o 的广义傅里叶变换:(2)又称为分数阶傅里叶变换(fractional Fourier transform),当=/2时,分数阶傅里叶变换显然变为常规傅里叶变换,2023/6/26,光学信息处理,5,数学家的贡献,早在1937年,Condon提出了广义傅里叶变换的初步概念到1980年,Namias 完整地提出了广义傅里叶变换的数学定义、性质,讨论了变换的本征函数,并用于处理谐振子的薛定谔方程、格林函数问题、在均匀磁场中的自由电子的能级、在含时间变量的均匀磁场中自由电子薛定谔方程的求解等1987年,McBride 和 Kerr 进一步研究了广义傅里叶变换,把变换看作是充分光滑的函数构成的向量空间(Frechet 空间)中的算子,在此框架内建立了广义傅里叶变换更为严谨、完整的理论系统,这两篇文章至今仍是广义傅里叶变换的理论基础,2023/6/26,光学信息处理,6,物理学家的贡献,直到90年代,光学科学家和工程师开始关注广义傅里叶变换与光学的关系,与三十年前常规傅里叶变换与光学的结合产生了傅里叶光学的情况非常相似1993年,Ozaktas 和Mendlovic 提出用平方折射率光波导(GRIN)来实现广义傅里叶变换;Lohmann,Bernardo等则用透镜系统成功地实现了这一变换;Lohmann还设计了阶数连续可变的广义光学傅里叶变换系统;Bernardo等认为应正确地称这一变换为广义傅里叶变换,而不是分数阶傅里叶变换,因为阶数既可以是整数、分数,还可以是复数,2023/6/26,光学信息处理,7,广义傅里叶变换与其他变换关系,Lohmann,Mendlovic阐明了广义傅里叶变换与维格纳变换的关系,指出可以用维格纳空间中的旋转来一般地定义广义傅里叶变换,这一定义与光波在梯度折射率介质中的传播的定义是等价的。Mendlovic等进一步讨论用广义傅里叶变换来表征信号的新方法,以及分数阶光学相关;Dorsch,Bernardo等分别提出了用光学系统实现任意阶傅里叶变换的方案;,2023/6/26,光学信息处理,8,广义傅里叶变换与其他变换关系,Ozaktas 等研究了广义傅里叶变换与小波变换的关系,他们认为广义傅里叶变换可以表为小波变换,小波函数具有h(x)=exp(ix2)的形式然而该函数是分布在(-,)上的振荡函数,并不具备小波的特点易证h(x)的傅里叶变换H(u)=exp(iu),而H(0)0,不符合小波变换的相容性条件因而我们认为广义傅里叶变换只是形式上与小波变换相似 Mendlovic等对变换的形式稍加改换,定义了广义余弦变换,该变换适用于非相干光,在数字成像、非相干光信息处理方面都有潜在的应用众所周知,夫琅和费衍射可以实现常规的傅里叶变换,Pellat-Finet则探讨了菲涅耳衍射与广义傅里叶变换的关系,2023/6/26,光学信息处理,9,傅里叶变换在科学技术的许多领域中有广泛的应用,因此我们可以预料广义傅里叶变换的应用领域将更为宽广目前,它已成为数学、量子力学中重要的应用工具 本章将研究广义傅里叶变换的数学定义、性质及实现广义傅里叶变换的光学系统,并讨论与广义傅里叶变换有密切关系的维格纳变换,2023/6/26,光学信息处理,10,5.2 广义傅里叶变换的定义及性质,5.2.1 广义傅里叶变换的定义 仅讨论一维函数的广义傅里叶变换,有关的定义和性质可以直接推广到二维的情况 函数g()的广义傅里叶变换定义为(1)通常称它为g()的广义傅里叶谱,记为G(x),2023/6/26,光学信息处理,11,以-代替上式中的,得到(2)称为广义傅里叶变换的阶可证明F-是F 的逆变换,即:F-F g()=g(x)(4),2023/6/26,光学信息处理,12,广义傅里叶变换的主值 区间为(-,)。当 超出主值区间时,相应的变换可以化成在该区间内的变换。下面将证明这一点因此F-(0)实质上只是负阶数的广义傅里叶变换。广义傅里叶变换的一个性质,在于当=/2 以及=-/2 时化成常规的傅里叶变换及逆变换:(5)(6)注意这里傅氏变换的表达式与其他各章有所不同,2023/6/26,光学信息处理,13,当=0 时没有意义,因而Fo也必须另行定义.由于 0时,sin,tan,所以有(7)其中用到极限意义下的 函数的定义:(8)从而可用上述极限过程来定义:F og()=g(x)用类似的方法还可定义:F g()=g(-x)以上两式表明:0 阶广义傅里叶变换给出输入图像本身,阶广义傅里叶变换则给出它的倒像.,2023/6/26,光学信息处理,14,虽然广义傅里叶变换仍然是线性变换,即:FAg()+Bh()=AFg()+BFh()(13)式中A,B为常数,但由于变换公式中出现二次相因子,所以它的性质和常规的傅里叶变换有了很大的差别例如,它不再满足缩放规律,2023/6/26,光学信息处理,15,5.2.2 基本性质和运算法则,g()的广义傅里叶变换谱记为 G(x),并用 g G 表示变换对(1)位移(shift)g(+)exp isin(x+cos/2)G(x+cos)当=/2 时即化为傅里叶变换的位移公式(2)宗量乘积(multification)设 D d/dx为微分算符,m 0.F m g()=(x cos+i sin D)m G(x),2023/6/26,光学信息处理,16,F m g()=(x cos+i sin D)m G(x)例如设 m=2,有(xcos+i sinD)2=x2 cos2+x cos i sin D+i sinD xcos-sin2D2=x2 cos2+i x cos sin D+i sincos+i x sincos D-sin2D2=cos(x2cos+i sin)+i x sin2D-sin2D2(3)微分(differentiation)F Dm g()=(ix sin+cos D)m G(x),2023/6/26,光学信息处理,17,(4)宗量微分混合积(mixed product)F(Dm)g()=-(sin-i x2 cos)sin+xcos 2 D+i sin cos D2 m G(x)(5)指数(exponential)F ei b g()=exp-ibcos(x-bsin/2)G(x-sin)(6)可加性(additivity)F F g()=F+g()(28)对称性:F F g=F F g=F+g逆变换:F F-g=F-F g=F o g=g,2023/6/26,光学信息处理,18,(7)周期性(periodicity)由于在广义傅里叶变换的定义中出现tan 及sin,所以变换关于 具有周期性,周期为2,这样就有以下结果:F 2n g()=g(x)F(2n+1)g()=g(-x)F 2n+g()=F g()这样当(-,时的变换F 均可化为主值区间内的变换设=p/2,阶广义傅里叶变换还可表为F(p)g,p的定义域为(-2,2,2023/6/26,光学信息处理,19,图5.1 广义傅里叶变换的周期性,例如:F(1)F(1)g()=F(2)g()=F g()=g(-x)Fg()即常规的傅里叶变换得到的输入图像的傅里叶谱,可用2f系统实现;而F(2)g()则表示两次傅里叶变换,得到输入图像的倒像,可用4f 系统实现,它们都是广义傅里叶变换的特例,其中=/2 或p=l表示常规傅里叶变换,即F/2 或 F1,F 1 常简写为F=-/2 或 p=-l 则表示常规的傅里叶逆变换,即 F-/2 或F-1,2023/6/26,光学信息处理,20,5.2.3 广义傅里叶变换群,变换算符Fo 具有如下性质,对于任意的F,有 F Fo=Fo F=F 因此可称为单位算符或恒等元 对于 F 存在,满足F F-=F-F=Fo即F-是F 的逆算符或逆元 对于任意的实数、,有F F=F F=F+F+依然是广义傅里叶变换算符因此变换算符对于乘法是闭合的 结合律:F(F F)=(FF)F=F+因而所有的广义傅里叶变换算符对于(28)(可加性)所定义的乘法构成群,可称为广义傅里叶变换群,2023/6/26,光学信息处理,21,5.3 广义傅里叶变换的本征函数,广义傅里叶变换算符F 的本征函数为 n(x)=Hn(x)exp(-x2/2)本征值为 exp(-i n)其中Hn(x)为n 阶厄米多项式,exp(-x2/2)为高斯函数,所以n(x)常称为高斯-厄米型函数(GH函数)n(x)构成区间(-,)内的完备正交函数组,因此任何平方可积的函数g(x)都可以用它展开:g(x)=-n n(x)=-n Hn(x)exp(-x2/2)其中系数n 可用厄米函数的正交性得到:(11),2023/6/26,光学信息处理,22,Fg()=-nFn()=-nexp(-in)n(x)=-n exp(-in)Hn(x)exp(-x2/2)上式又称广义傅里叶变换的级数表达式在节5.6中将讨论它在渐变折射率介质光波导中的应用,我们知道,高斯厄米函数正是量子力学一维谐振子的本征函数,图5.2给出前几阶GH函数的图形,2023/6/26,光学信息处理,23,5.4 用透镜系统实现广义傅里叶变换的基本光学单元,5.4.1 第一类基本光学单元 由广义傅里叶变换的“可加性”可知:连续执行N个阶数为 n(n=1,2,N)的变换的结果,相当于执行阶数为 1+2+N 的一次变换,亦即 F=F n(式中=n),2023/6/26,光学信息处理,24,5.4 用透镜系统实现广义傅里叶变换的基本光学单元,5.4.1 第一类基本光学单元例:=/2的常规傅里叶变换,既可由一个焦距为f 的透镜来实现,也可由两个相同规格的透镜构成的透镜组来实现,它们的焦距为(3)间距为2d,(4),2023/6/26,光学信息处理,25,图5.3 用两个透镜实现傅里叶变换,2023/6/26,光学信息处理,26,几何光学的计算还可证明,N 个焦距为(5)的透镜按图5.3的方式串联起来,间距参数(6)则该系统的合成焦距,且前焦面位于第一个透镜前d 处,后焦面位于第 N 个透镜后d 处。透镜系统能否实现傅里叶变换?必须首先证明:当单色光波通过一个透镜单元,即经过两次距离为 d 的菲涅耳衍射,并经过一次透镜相位变换,其效应相当于=/2N 阶广义傅里叶变换,才能通过变换的可加性得到该系统实现傅里叶变换的普遍结论,2023/6/26,光学信息处理,27,第一类基本光学单元,广义傅氏变换与常规傅氏变换的主要差别:不仅包含线性相位项,还包括二次相位项广义傅氏变换重新记为式中=1/d1+1/d2-1/f(8)其中的二次相位因子是由两次菲涅耳衍射及透镜相位变换的合成效果,(7),2023/6/26,光学信息处理,28,若d1=d2=f,则即二次项消失,变成常规的傅里叶变换,这是我们已熟知的傅里叶变换的光学实现方法也就是说,薄透镜单元在特别的输入距离,输出距离的配置下产生了=/2 的常规的傅里叶变换效应,(10),2023/6/26,光学信息处理,29,设想d1和d2 不等于f,看看在这种情况下有没有可能产生广义傅里叶变换,其阶数/2。将(7)式与节5.2(1)式比较,发现必须满足条件 d1=d2=d,但d 不一定等于 f。,(7),5.2(1),2023/6/26,光学信息处理,30,再仿照(4)及(3)式,设(12)(13)代入=1/d1+1/d2-1/f,得到,(3),(4),把上面的结果代入(7)式得到(15)式引入归一化坐标(无量纲)(16)其中(17)(15)式变成(18)式(18)式中常数(19),2023/6/26,光学信息处理,32,一维广义傅立叶变换,二维广义傅立叶变换,2023/6/26,光学信息处理,33,正透镜实现广义傅里叶变换,将(18)式与节5.2(1)式相比较,发现除积分号前的常数因子外,它就是二维 阶广义傅里叶变换,即(20)式中Co 和C 有所不同,在条件(12)及(13)式成立时,薄透镜在单色光的照射下,将透镜前面d 处的输入图像o 变成它的广义傅里叶谱,形成在透镜后d 处,,正透镜,2023/6/26,光学信息处理,34,称为族参数,d 称为间距参数 光学广义傅里叶变换表达式与变换的数学定义式的最大差别在于光学系统中存在族参数。很明显,只有族参数相同的光学广义傅里叶变换才能组成群,不同族参数的变换不具备可加性。族参数仅取决于及 f,然而在 确定后,透镜的焦距就确定了,这对光学单元按可加性组合带来许多限制。,2023/6/26,光学信息处理,35,在用透镜系统实现广义傅里叶变换时,我们仍用公式(15),即定义(21)其中积分号前的系数(22)与数学定义中的归一化系数并不相同,但由于我们只能探测光强的分布,因此这一差别并不带来实质性的影响,2023/6/26,光学信息处理,36,广义傅里叶的阶数另一种定义,在光学中常常使用令 p=2/,则(21)式为(24)(21)及(24)式都是经常运用的 显然,F/2 或F(1)即常规的傅里叶变换,F-/2 或F(-1)即傅里叶逆变换,分别按习惯的记法,记为F 及F 1 常规傅里叶变换仅能用正透镜实现,2023/6/26,光学信息处理,37,负透镜实现广义傅里叶变换,在 式中,当 f 0,则 0,则得到负阶数的广义傅里叶变换,它可以用图5.5所示的负透镜单元实现。由,得 d 0,,阶数 由下式决定:cos 1-d/f(25)这里d 和f 都是负数,o和 1 分别是光学中的虚物和虚像.,2023/6/26,光学信息处理,38,5.4.2 第二类基本光学单元,实现广义傅里叶变换的第二类基本光学单元如图所示,两个规格相同的正透镜焦距为f,间距为d.,设在两个透镜之间光波的传播遵循菲涅耳衍射的规律,则在1面上的场,按(25)及(13)式设 及 d,并规定d 的方向从0 指向1 向右为正,则有 将上式与(15)式相比较,发现只要设(27)(28)则有 1(u,v)=C F(o(x,y)(30)亦即第二类基本单元也能实现广义傅里叶变换,注意两种基本光学单元的族参数的定义不相同,2023/6/26,光学信息处理,40,将图5.6中两个正透镜均改为负透镜,并设 f 和 均为负值,此 仍为正值,而间距d 0,说明 o 在 1 左面,如图所示,它能实现负阶数的广义傅里叶变换由此可知,用两个透镜构成的系统也能实现广义傅里叶变换,2023/6/26,光学信息处理,41,5.5 基本光学单元的组合,在讨论用透镜或透镜组实现广义傅里叶变换时,引入了族参数第一类单元:(1)第二类单元:(2)显然,同一类型的广义光学傅里叶算符,仅当族参数相等时才有可加性:F(p2)F(p1)=F(p2+p1)(3)换言之,族参数相同的广义光学傅氏算符属于同一群(3)式暗示同属一群的光学广义傅氏算符对应的光学单元具有互相组合成复杂系统的性能,2023/6/26,光学信息处理,42,讨论:第一类光学单元的组合,当 取某一常数时,p 或 可取两个值:1=(p1=p=2/)(4)2=-(p2=2-p)(5)由此导出两个不同的d:d1f(1-cos)f 1-cos(p/2)(6)d2f(1+cos)f 1+cos(p/2)(7),2023/6/26,光学信息处理,43,举例:系统由两焦距均为f 的透镜单元构成,,第一个光学单元:d11=d12=d1=f(1-cos)=f 1-cos(p/2)(8)第二个光学单元:d21=d22=d2=f(1+cos)=f 1+cos(p/2)(9)p1=p,1=,p2=2 p,1=-(10)两个透镜的问距为 d12+d21=d1+d2=2f(11)它们的共同族参数为(12),2023/6/26,光学信息处理,44,设在0平面输入图像o(x,y),则在单色光的照射下,第一个光学单元对 o 进行阶数为p1的广义傅里叶变换,1=F(p1)o;第二个光学单元再对1进行阶数为p2 的广义傅里叶变换,2=F(p2)1,最后在输出平面2上得到输出图像2(x,y).,2023/6/26,光学信息处理,45,1=F(p1)o,2=F(p2)1,2(x,y)=F(p2)F(p1)o(x,y)=F(p2+p1)o(x,y)=F(2)o(x,y)=o(-x,-y)(13)表示最终得到输入图像的倒像当 p=1 时得到 4f 系统,其区别在于,在上述系统中,谱平面1上呈现输入信号的广义傅里叶谱而在4f 系统中,谱平面上呈现输入信号的傅里叶谱,2023/6/26,光学信息处理,46,包含负透镜的广义傅里叶变换系统,(14)上式成立的条件为 1=3=-,p1=p3=2 p(/2)2=-,p2=-p,系统由3个第一类单元构成,Ll 和L3 为同样规格的正透镜,焦距为f,而L2为负透镜,焦距为-f,它们具有共同的族参数为,2023/6/26,光学信息处理,47,对于中间的负透镜,我们有 d21=d22=(-f)(1-cos)=-f 1-cos(p/2)d21,d22 为负值,表示第二个单元的输入平面位于透镜右方1 处,而输出平面位于透镜左方2 处,关于透镜对称分布,2023/6/26,光学信息处理,48,对于两个正透镜,则有 di1=di2=f(1+cos)=f1+cos(p/2)(i=1,2,3)第一个光学单元的输出平面为1,恰为第二个单元的输入平面;而第二个单元的输出平面为2,又是第三个单元的输入平面 由于具有共同的族参数,因此变换具有可加性:F(p3)F(p2)F(p1)o=F(p3+p2+p1)o=F(4-3p)o 当p=2/3 时有F(4-3p)o=F(2)o=o(-x,-y)同样得到倒像上例表明,正、负透镜可以适当组合,其条件是它们有共同的族参数,且正、负透镜焦距的绝对值相等,2023/6/26,光学信息处理,49,事实上,属于同一群的广义傅里叶算符只要求族参 相同,而 f 和 均可不同例如设(21)两个单元对应的广义傅里叶算符分别为F1 和 F2,它们仍属于同一群,具有可加性:F1F2=F1+2(22)所以两个光学单元之可以串接,在条件(21)满足时 i 有两个解,分别记为 1=1,1=-1 2=2,2=-2(23),2023/6/26,光学信息处理,50,与之相应,有两组可能的d:d1=f1(1-cos 1),d1”=f1(1+cos 1)d2=f2(1-cos 2),d2”=f2(1+cos 2)(24)由上式还可得到:didi”=fi2(1cos2i)=fi2sin2i=(i=1,2)(25)亦即属于同一族的N个光学单元必须满足的条件 d1d1”d2d2”dNdN”=(26),2023/6/26,光学信息处理,51,举例:设用两个第一类光学单元构成系统,使第一个透镜的输出平面(即系统的“谱平面”)的复振幅分布为输入图形 o 的F/6o,而系统的输出平面为o 的常规傅里叶变换,即,2023/6/26,光学信息处理,52,两个不同规格透镜组合实现常规傅氏变换,几何光学的计算表明该系统的输出平面确是系统的焦平面若要求输出图像为o(-x,-y)(倒像),即要求:1+2=(29)则有 sin2=sin(-1)=sin1(30)从而有 f1=f2(31),2023/6/26,光学信息处理,53,亦即两个第一类单元串联,并要求输出图像为输入图像的倒像时,两个透镜焦距相等,但 和d 仍有两种不同的组合:1=,1=-2=-,2=d1=f(1-cos),d1”=f(1+cos)d2=f(1+cos),d2”=f(1-cos)注意在两种情况下,透镜的间隔均为2f用几何光学易验证在任一情形下,输出平面与输入平面共扼,且放大率为-l(倒像)由光线的回溯性定律立即可以把这一结论推广到第二种情况,2023/6/26,光学信息处理,54,从理论上讲,N 个相同规格的第一类单元串联起来,如每个单元对应的变换为F(1/N),则整个系统对应的变换为常规傅里叶变换。由此可见,单个薄透镜并不是能实现常规傅里叶变换的唯一光学模型。一般来讲,=/2 的透镜组合均能实现常规傅里叶变换。然而随着N 的增大,光能损失也会变大,且杂散光的效应会变得越来越严重。,2023/6/26,光学信息处理,55,5.6 用自聚焦效应光波导实现广义傅氏变换,5.6.1 渐变折射率介质波导的自聚焦效应 圆柱形渐变折射率波导的折射率形式为n2(r)=n121-(n2/n1)g(r)(1)柱坐标系rz的z 轴与波导的对称轴一致,又称光轴通常只考虑平方折射率介质(称为渐变折射率介质)的情形,此时有 n2(r)=n121-(n2/n1)r2,r a 其中n1,n2是两个参数,a 是波导的半径,2023/6/26,光学信息处理,56,5.6.1 渐变折射率介质波导的自聚焦效应,图5.11给出n(r)随r的变化,当ra时n 趋于常数,即波导的包层,2023/6/26,光学信息处理,57,考虑 r a 的近轴情况,n(r)=n11-(n2/n1)r21/2 n11-(n2/2n1)r2=n1-n2 r2/2 设单色光波沿z 轴方向穿过相距为 z 的两个平面z zo 及 z zo+z,在 z zo 上光波分布为o,则在z zo+z 上光波复振幅的分布为 1(x,y)=o(x,y)exp(i2n z/)o(x,y)exp(i2n1 z/1)exp(-in2 z(x2+y2)/)式中 为光波在真空中的波长,并假定z 很小,2023/6/26,光学信息处理,58,1(x,y)o(x,y)exp(i2n1 z/1)exp(-in2 z(x2+y2)/)将此公式与薄透镜的相位效应公式:tl=expik(-f)expik(x2+y2)/2f/f(4)相比较,发现它们都具有关于坐标的二次相位因子,而它恰恰是球面波的近轴表达式,所以一个厚度为 z 的很薄的介质波导,相当于一个薄透镜,其等效焦距为 l/n2 z 这样,渐变介质折射率波导应具有对平行入射光的聚焦效应,2023/6/26,光学信息处理,59,渐变折射率自聚焦效应的证明,设A和B是介质中的两点,它们之间的光程(5)根据费马原理,上式的变分为 0,且(6)设 r 是变分时光线的位移,则有 n=r n dl=dr lo,式中 lo 是光线路径在给定点的单位切线向量,2023/6/26,光学信息处理,60,代入(6)式,运用分部积分法,得到(9)由于r在两个端点A,B处为0,上式第三项消失又因为 dr=lo dl(10)所以有(11)因此有(12),2023/6/26,光学信息处理,61,将n(r)的表达式代入上式右边考虑到近轴近似,d/dl d/dz,上式化作(dn/dz)(dr/dz)+nd2r/dz2=-n2 r(13)由于 n 不随z 的变化而变化,上式左边第一项为0用n1近似代替左边第二项中的n,得到微分方程 d2r/dz2+(n2/n1)r=0(14)上式的解为(15)式中ro=r(0),ro=r(0)分别是 r 和dr/dz 在端口z=0的值,n(r)=n1-n2 r2/2(r a),2023/6/26,光学信息处理,62,设一束平行于z 轴的光照射z=0 端面,则 r00,(16)上式表明光线在渐变折射率波导中的传播路径具有周期性,如图5.13所示 设一个周期的长度为2L,则有而焦距为这就是自聚焦透镜的原理,光线在渐变折射率波导中的传播,2023/6/26,光学信息处理,63,在5.5节中曾讲过,常规的傅里叶变换效应可以由N个第一类光学基本单元串联而成的系统来实现随着N的变大,每个单元透镜的焦距 也随之变大因此每个薄透镜也就更接近于一个平板这样一系列接近密接的透镜组与上述自聚焦光波导薄层系列非常相似 因此我们可以设想一段长度为pL 的渐变折射率光波导能够实现 p 阶广义傅里叶变换,2023/6/26,光学信息处理,64,5.6.2光波在渐变折射率介质波导中的传播,为简单起见,只讨论一维情况,其结论可以直接推广到二维情况仅考虑TE波 设光波沿z 轴方向传播,电场在y方向振动,用E(x,z)表示,则它应满足亥姆霍兹方程(2/x2+kx2)E(x,z)=0(19)式中kx 是x方向的波矢量的值,kx2=k2-kz2(20)令ko为真空中的 k,kz=称为传播常数,n 为折射率,则有 k2=n2 ko2(21)代入(19)式,得到2E/x2+(n2 ko2-2)E=0(22),2023/6/26,光学信息处理,65,在二次渐变折射率介质中,n(x)n1 n2 x2/2(23)这种情况相当于量子力学中无限深势阱中的一维谐振子将上式代入(22)式得到2/2+(-2)=0(24)式中=2E,=x,=(n12 ko2-2)/2,=(n1 n2)1/2 ko(25)从(24)式解出式中Hn为n 阶厄米多项式由此可见,在渐变折射率波导中传播的本征模就是GH函数,恰恰是广义傅里叶变换的本征态,2023/6/26,光学信息处理,66,在波导中的场为式中 称光斑尺寸,(30)式中 为真空中的波长,2023/6/26,光学信息处理,67,5.6.3 用渐变介质波导实现广义傅里叶变换,引入归一化坐标(31)并记En(x)为n(x),由于n 是(-,)区间的本征函数,所以任何平方可积的函数g(x)都可以用n(x)展开:(32)系数,2023/6/26,光学信息处理,68,设在z=0处放置由g(x)表达的一维图像,则在 z=pL 处的场为(35)上面曾得到结论:长度为L(p=1)的光波导能实现光波的聚焦我们将讨论在p1时,长度为pL 的光波导对光波的变换关系 注意到焦距L的表达式(18),就有,2023/6/26,光学信息处理,69,代入(35)式,得到(37)将广义傅里叶变换对本征函数的变换关系式(节5.2(9)式)代入(37),得到(39)(39)式表明:一段长度为pL 的渐变折射率介质波导,也能实现广义傅里叶变换在z=0 处放置图像 g,并用单色光波照射,在 z=pL 处即可得到它的广义傅里叶变换谱 F(p)g。,2023/6/26,光学信息处理,70,显然,长度为p1L和p2L 的两段光波导串接起来,必然能实现 p1+p2 阶广义傅里叶变换,而变换算符的可加性及可易性是由光波导长度的可加性及可易性自然满足的。这样,我们又找到一种实现广义傅里叶变换的光学器件。串联的透镜组只能得到P/Q阶广义傅里叶变换,P和Q都是正整数,而由于技术上的原因,Q无法做到太大,所以我们只能得到阶数的步进的变化。由于光波导中阶数正比于波导的长度,因此能够实现阶数的连续的变化。光波导的缺点:空间带宽积比较小,也就是说它能处理的信息量比较小。,2023/6/26,光学信息处理,71,5.7 维格纳变换,5.7.1 维格纳变换的定义函数g(x)的维格纳(Wigner)变换定义为(1)如果把 x,v 分别理解为信号g 的空间变量和空间频率变量,则维格纳变换就是信号函数g 的空间和频谱特征的综合表现或同时反映,因此在光学信号及其他信号的处理中有广泛的应用 维格纳变换的缺点:它并不是线性变换,所以将它应用到线性系统中去有一定的困难,2023/6/26,光学信息处理,72,设G(v)是g(x)的傅里叶变换,则有(2)(3)以(2)及(3)式代入(1)式得到(4)其中的积分(5)代入(4)式,经整理后得到维格纳变换在频域中的表达式(6),2023/6/26,光学信息处理,73,容易证明逆变换表达式(7)(8)维格纳变换对于v 和x 的积分分别得到信号在空域和频域中的功率密度(9)(10)因而信号的总功率Eo可以由下面的积分来表达:(11),2023/6/26,光学信息处理,74,下面讨论当信号变化时其维格纳谱的相应变化首先把(1)式改写成更加容易处理的形式:(13)式中为光波的波长,f 是维格纳变换的一个特征量,以后我们会看到它相当于广义傅里叶变换中 f;和 x 都具有长度量纲,分别是频域和空域中的长度坐标;/f 则是空间频率变量。,2023/6/26,光学信息处理,75,5.7.2 基本性质和运算法则,(1)反演(inversion)当g(x)g(-x)时,相应的维格纳变换成为(16),可见信号在 x 空间的反演引起维格纳变换在(x,)空间的反演(参见图5.14):W(A)W(A)(17),2023/6/26,光学信息处理,76,上述变换可以用矩阵来表示,设(x,)和(x,)分别是变化前后的维格纳空问坐标,(18),则反演的结果可用下式表示:(19),2023/6/26,光学信息处理,77,(2)缩放(scaling),设 g(x)g(x/M)式中M为缩放倍率,相应的维格纳谱的变化为 W(x,)W(x/M,M)(21)或(22)即信号坐标 x 的放大将引起维格纳空间中 x 坐标按同样倍率的放大,同时 坐标按同样比例缩小,2023/6/26,光学信息处理,78,(3)薄透镜的相位变换,g(x)g(x)exp(-i x2 Q/f)设透镜的焦距为:f=f/Q 相应的维格纳变换为(25)变换表达式为(26)在维格纳空间中表现为坐标沿平行轴的移动Q x.,2023/6/26,光学信息处理,79,(4)在自由空间中的传播菲涅耳衍射,信号在自由空间中传播的距离为z 时,相应傅里叶谱G 如下变化 G()G()exp(-i z2)(27)这样,我们就有 G(+/2)G(+/2)exp-iz(+/2)2(28)G*(-/2)G*(-/2)expiz(-/2)2(29)代入(13)式,得到(30),2023/6/26,光学信息处理,80,(4)在自由空间中的传播菲涅耳衍射,相应的矩阵变换式为(32)在维格纳空间中相当于坐标沿平行于 x 轴的平移.,2023/6/26,光学信息处理,81,(5)夫琅和费衍射,将它代入(12)式,得到令R=1 得到(35)(36),2023/6/26,光学信息处理,82,在维格纳表象中,夫琅和费衍射表示的常规傅里叶变换表示绕原点旋转/2,见图5.17,2023/6/26,光学信息处理,83,5.7.3 广义傅里叶变换的维格纳变换定义,我们自然要问,当/2 时的旋转能否表示广义傅里叶变换?可以证明:在维格纳空间中,广义傅里叶变换相当于坐标系旋转 角。从另一角度来看,由于 阶广义傅里叶变换可以由第一类光学系统实现,它是两次菲涅耳衍射,以及在其间的透镜相位变换的结果。,2023/6/26,光学信息处理,84,矩阵变换法证明,设两次光波传播对应的距离为 z,透镜的焦距为 f,f=f/Q,z=R f 按广义傅里叶变换的参数定义,令R=tg(/2)=(1-cos)/sin(44)Q=sin(45)则信号的广义傅里叶变换对应的维格纳变换可以由矩阵变换式得到(46),2023/6/26,光学信息处理,85,5.7.4 维格纳变换在光学系统设计中的应用,因为信号(即光学图像)的变化可用矩阵表示或用维格纳表象中的简单变换表示,因而可以用维格纳变换进行光学系统的总体设计 例如,要设计一个变焦距的光学系统,亦即要求变换矩阵的参数Q可变透镜的焦距是固定的,因而这一要求只能通过透镜间距离的调节来达到,2023/6/26,光学信息处理,86,方案一:,前后矩阵表示常规傅里叶变换,中间的矩阵是光波在自由空间中的传播,传播距离为z则式中 表示生成倒像,而则表示三个矩阵的积相当于一个透镜,它的合成焦距 可以通过自由传播的距离 z 来调节在这个例子中,参数 f 相当于单个透镜的焦距,即 f=f1=f2,2023/6/26,光学信息处理,87,用第一类或第二类单元构造ZOOM系统第一类第二类,