CN107703626A - 一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法 - Google Patents
一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法 Download PDFInfo
- Publication number
- CN107703626A CN107703626A CN201710979945.6A CN201710979945A CN107703626A CN 107703626 A CN107703626 A CN 107703626A CN 201710979945 A CN201710979945 A CN 201710979945A CN 107703626 A CN107703626 A CN 107703626A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- lens
- mfrac
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 239000000463 material Substances 0.000 title claims abstract description 61
- 238000013461 design Methods 0.000 title claims abstract description 30
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000005452 bending Methods 0.000 claims description 10
- 230000003287 optical effect Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 7
- 238000012545 processing Methods 0.000 abstract description 2
- 230000005684 electric field Effects 0.000 description 22
- 238000004088 simulation Methods 0.000 description 8
- 238000001514 detection method Methods 0.000 description 6
- 238000009738 saturating Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000011664 signaling Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical group C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Eyeglasses (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明是一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，属于变换材料和光学分数阶傅里叶变换信号处理领域。本装置主要基于传统GRIN透镜扩大输入面，保持输出面不变或缩小，通过保角变换等方法获得其与传统GRIN透镜的变换关系，再利用变换材料的方法获得本装置的折射率的共形变换装置；本方法的步骤为：1在传统GRIN透镜的基础上保证透镜输出面和透镜厚度不变，弯曲透镜且角度在0度和360度之间；2保持透镜输出面不变，拉伸透镜至各点各向同性；3利用变换材料的方法计算本装置的材料参数。本装置在输入低调频率chirp信号的情况下具有低测量误差和高分辨率；本装置为各向同性材料，有利于工程制备。
Description
技术领域
本发明是一种基于变换材料(Transformation Material)的光学分数阶傅里叶变换 (Fractional Fourier Optics Transformation)的装置和设计方法，属于变换材料以及光学分数阶 傅里叶变换信号处理领域。
背景技术
变换材料的雏形在2003年已经出现，当目前学术界公认是在2006年由J.B.Pendry和 U.Leonhardt分别独立提出的变换光学理论。Pendry和Leonhardt分别独立提出的这一理论，可 以从已知的参数通过“变换”而得来的参数通常用来表示材料性质，所以这种通过变换方法获 得其参数的材料，被称为“变换材料”。在电磁波、光学领域，这一方法被称为“变换光学”。
早在1937年，Condon在他的论文中就提出了分数阶傅里叶变换的初步概念。自此之后， 分数阶傅里叶变换在各个领域得到了广泛的应用。1980年分数阶傅里叶变换首次作为一种数 学手段应用于光学，之后分数阶傅里叶变换在光学中一直没有得到重视，直到1993年Ozaktas 和Mendlovic首先提出了将分数阶傅里叶变换应用于将光束通过二次折射率介质中。分数阶傅 里叶变换成为一种解释光学现象的新方法。在光学中比较常见的chirp信号主要是牛顿环，对 牛顿环的调频率检测在工程中有很多应用，比如检测球形期间表面光滑度、薄膜材料厚度等。 因此精确、方便、快速地测定其厚度十分必要，本发明主要是利用光学分数阶傅里叶变换透 镜快速精确定测量信号的调频率。
利用二次折射率透镜检测chirp信号调频率的文献目前还没有检索到，本发明对传统二次 折射率进行了改进，降低了低调频率chirp信号的测量误差，提高了装置的性能。
发明内容
本发明的目的是为了改善传统的二次折射率(Graded Index，GRIN)透镜在输入调频率 较低的chirp信号情况下测量误差大和输入多个相同幅度chirp信号检测精度低的问题，提出了 一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法。
本发明的核心思想是：以介质的特性作为手段完成对光场的控制，具体主要是基于光波 在介质中的传播的性质，当预先知道输入chirp信号在分数阶傅里叶变换装置中的处理效果， 即明确其光场分布，则可以利用变换材料的方法确定装置的介质参数，从而达到预期的光场 分布；
一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，包括分数阶傅里叶变换 透镜的装置和装置的设计方法，其中，基于变换材料的分数阶傅里叶变换透镜的装置，简称 本装置，本装置的设计方法，简称本方法；
其中，本装置主要包括输入面、基于变换材料设计的二次折射率透镜和输出面；
本装置中基于变换材料设计的二次折射率透镜，简称透镜；
输入面是输入信号进入二次折射率透镜的入射面，由于本装置是基于变换材料的原则， 虽然本装置的输入面与原始的GRIN透镜输入面相比，形状发生了拉伸与改变，但在数学意义 与物理意义上，本装置的输入面的点与原始的GRIN透镜输入面为一一对应的关系；
本装置中透镜与传统的GRIN透镜相比，通过变换材料的方法，改变传统的GRIN透镜中 原有的轴向与径向梯度折射率的分布，使透镜的折射率呈类似于环状分布，主要表现为从输 入面到输出面，透镜折射率依次呈环状增大，并在此过程中透镜的环状折射率分布所形成的 环形半径逐渐减小；透镜折射率的改变导致了光在透镜的传播路径的改变；本装置在满足折 射率分布的情况下，其形状可以不受本装置示出形状的限制；在信号输入到本装置的输入面， 进入透镜后，在透镜中是一个分数阶傅里叶变换的过程，最终在输出面形成一个阶次为1的输 入信号的分数阶傅里叶变换，即为普通的傅里叶变换，在透镜内部则形成了阶次不同的输入 信号分数阶傅里叶变换；本装置是在传统的GRIN透镜的基础上扩大输入面，保持输出面不变 或缩小的共形变换，本装置可以通过复变函数的保角变换等方法获得本装置和传统二次折射 率透镜的变换关系，再利用变换材料的方法获得本装置的折射率；
输出面：即频谱面，输入光学信号经过二次折射率透镜之后，出射波为输入波的一个傅 里叶变换信号；
本方法，通过以下步骤实现：
步骤1：在二维平面下，以传统的GRIN透镜为基础，保证透镜输出面和透镜的厚度不变， 弯曲透镜；
其中，透镜的高度所指方向z方向，且步骤1假设z方向是无限延伸的；步骤1中的透镜弯 曲角度为γ，γ弯曲程度可为0度到360度；
步骤1中透镜在二维平面主伸长表示为公式(1)：
其中，h为透镜的厚度，θ和r为步骤1透镜弯曲后透镜中的点的极坐标，分别为 传统的GRIN透镜在步骤1弯曲后在极坐标下θ和r的主伸长；
步骤2：透镜在步骤1的基础上，沿透镜中心到坐标r的方向拉伸，拉伸至透镜中的各点 各向同性；
经过步骤2后透镜在二维平面的主伸长表示为公式(2)：
其中，二维平面是r和θ组成的二维平面；θ′和r′为步骤2拉伸后透镜中的点的极坐标， 分别为步骤1中的透镜在步骤2拉伸后在二维平面的主伸长，表示θ′对θ求偏导， 表示r′对r求偏导，；
其中，步骤2要拉伸至各点各向同性，主伸长要满足公式(3)：
其中，λ为分别经过步骤1弯曲和步骤2拉伸后透镜中各点在各个方向上的主伸长；
由于步骤1弯曲和步骤2拉伸过程中始终保持透镜输出面不变，则GRIN透镜和经过步骤1 和步骤2后的透镜的变换关系满足公式(4)和(5)：
其中，R1为输出面所在圆弧的半径，f(r)为步骤1弯曲和步骤2拉伸后透镜中的各点与传 统的GRIN透镜中各点的函数关系，exp[]为指数函数，f-1(r')为f(r)函数的逆函数，ln为以e 为底的对数；表示以e为底，对取对数；公式(4)也体现了本装置中输入面与原始 的GRIN透镜输入面上的点是一一对应的对应关系；
步骤3：通过变换材料的方法，基于步骤2计算出变换后新的透镜的介电常数和磁导率；
其中，变换材料的方法中材料介电常数和磁导率变化前后关系为公式(6)和公式(7)：
其中，ε′和μ′为变换后材料的介电常数和磁导率，ε0和μ0为变换前材料的介电常数和 磁导率，λ1,λ2,λ3为三维坐标系下的主伸长；
由于步骤1中假设z方向上无限延伸，即假设在z方向没有任何变形，则在z方向上的主伸 长λz＝1，通过公式(6)和(7)，可以得出经过步骤1弯曲和步骤2拉伸后透镜的介电常数、 磁导率和折射率的关系，即公式(8)：
其中，变形主要指弯曲和拉伸，也可以指压缩公式(8)中的ε0、μ0、n0为变换前透 镜的介电常数、磁导率和折射率，εr′εθ′εz′为r，θ和z方向的介电常数，μr′,μθ′,μz′为r，θ和 z方向的磁导率，nr′,n′θ为r，θ方向的折射率；
至此，经过步骤1到步骤3，完成了本方法。
有益效果
本发明一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，与传统GRIN傅里 叶透镜装置相比，具有如下有益效果：
1.在输入低调频率chirp信号的情况下具有较低测量误差；
2.当输入为多个同幅值不同调频率的chirp信号尤其是调频率较低的chirp信号，较传统的 GRIN透镜，能够测出调频率较接近的chirp信号，具有更高的分辨率；
3.本发明采用了共形映射的方法，保证了材料的各向同性，有利于工程上的制备。
附图说明
图1为本发明“一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法”在实施例 1中的本装置组成的示意图；
图2为本发明“一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法”在实施例 2中的本装置设计方法图；
图3为本发明“一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法”在实施例 3中的本装置应用于检测低调频率chirp信号与传统GRIN透镜的仿真结果对比图；
图4为本发明“一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法”在实施例4 中的本装置应用于检测多个同幅chirp信号与传统GRIN透镜仿真结果对比图；
图示说明：
图2中(a)原始的GRIN透镜(b)经过步骤1变换后的透镜(c)经过步骤2变换后的 透镜
图3中(a)GRIN透镜场图(b)本装置仿真场图(c)GRIN透镜横轴在y＝0截线的电场图(d)本装置横轴在y＝0截线的电场图；
图4中(a)输入面点场图(b)GRIN透镜场图(c)本装置仿真场图(d)GRIN透镜横轴在y＝0截线的电场图(e)本装置横轴在y＝0截线的电场图。
具体实施方式
下面根据附图及实施例对本发明进行详细说明，但本发明的具体实施形式并不局限于此。
实施例1
本实施例阐述了将本发明“一种基于变换材料的光学分数阶傅里叶变换的装置和设计方 法”中本装置的组成和具体实施。
从图1中所示，图1是本装置的组成图。
从图1中可以看出，CD是输入面，ABCD是本装置在传统GRIN透镜保持输出面不变，弯 曲角度为45度的保形映射条件下经变换材料方法得到的透镜的折射率分布，AB是输出面。
具体到本实施例中，光波从CD入射，在输入面的光学信号与传统的二次折射率透镜相比， 当chirp信号输入到本装置中，由于本装置的输入面的点与原始的GRIN透镜输入面为一一对 应，所以与同一chirp信号输入到传统的GRIN透镜相比，表面上该信号进入本装置的入射面与 信号进入传统的GRIN透镜的入射面信号是一样的，实际上，信号在输入面空间被压缩，信号 实际进入到透镜的调频率增大，相当于一个更高频率的信号进入到透镜中，这样改善了传统 的GRIN透镜在输入调频率较低的chirp信号情况下测量误差大的问题；信号从CD进入二次折 射率透镜，最终在AB面输出，在AB面输出的光学信号是CD面输出的光学信号的一个傅里叶 变换。
实施例2
本实施例阐述了将本发明“一种基于变换材料的光学分数阶傅里叶变换的装置和设计方 法”中本方法的流程和具体实施。
从图2为本发明“一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法”中本方 法的流程和具体实施。
从图2中可以看出，本方法包含如下步骤：
步骤A：保持透镜厚度和输出面不变，弯曲透镜；
具体到本实施例中，如图2(a)所示x，y坐标系二维下笛卡尔坐标系，传统的GRIN透镜 的高度为h，厚度为L，输入面CD距离y轴的距离为R2，输出面AB距离y轴的距离为R1，将图(a)中透镜弯曲到如图(b)所示极坐标系下，变换后的透镜输出面A′B′为h，厚度为L， 输入面CD距离圆心O的距离为R2，输出面AB距离圆心O的距离为R1，透镜的弯曲角度为γ；
其中，步骤A中透镜在二维平面主伸长表示为公式(9)：
步骤B：透镜输出面不变，沿坐标r方向拉伸透镜至各点各向同性；
具体到本实施例中，如图(c)所示极坐标系，变换后的透镜输出面A″B″保持不变，仍 为h，厚度为L′，输入面CD距离圆心O的距离为R2′，输出面AB距离圆心O的距离为R1，透镜的弯曲角度为γ；
其中，经过步骤A和步骤B中变换后透镜中各点的主伸长为公式(10)：
其中，图(a)GRIN透镜与图(c)透镜中的各点关系如下：
步骤C：经过步骤B得出变换后的透镜与原始透镜的关系，根据变换材料的方法，可以得 出新的透镜的材料参数为：εr′＝εθ′＝ε0，μr′＝μθ′＝μ0，εz′＝(1/λ)2ε0，μz′＝(1/λ)2μ0， nr′＝n′θ＝(1/λ)n0；
至此，从步骤A到步骤C，完成了本实施例中一种基于变换材料的光学分数阶傅里叶变换 装置设计的方法。
实施例3
本实施例阐述了将本发明“一种基于变换材料的光学分数阶傅里叶变换的装本置和设计 方法”中本装置应用于检测低调频率chirp信号与传统GRIN透镜的仿真结果对比图。本实施例 利用COMSOL Multiphysics进行计算机仿真，本实施例输入信号为cos(5y2)，传统的GRIN透 镜和本装置的电场分布。
本施例中的参数设置为：输出信号的波长0.09cm，按照实施例2中的参数为h＝3.5m， L＝3.0377m，GRIN透镜的折射率为n(y)＝1.5[1-(1/2)η2y2]且η＝0.5171，L＝π/(2η)， γ＝1radian，R1＝3.5m，R2＝8.3368m，此时，本装置输入面与传统GRIN透镜输入面的压缩比为2.3819。
假设特殊的chirp输入信号为公式(12)：
其中exp()为指数函数，f(x)为关于x的函数，i为复数，m为chirp信号的调频率，m0为 一次项系数；
则调频率与阶次的关系为公式(14)：
其中，s2＝λξ/n1，λ为波长，ξ，n1为GRIN透镜的参数，p为chirp信号的阶次；在传 统GRIN透镜中阶次与透镜厚度的关系表示为公式(15)：
其中，l为透镜内部某一平面距离输入面的距离，L为透镜的厚度；
由施例2，可以看出本装置与传统GRIN透镜输入面的压缩率为公式(16)：
则输入公式(12)的chirp信号，本装置与传统GRIN透镜之间输入信号调频率的关系为公 式(17)：
其中为本装置中输入公式(12)信号实际的调频率；
则本装置中调频率与阶次的关系为公式(18)：
在传统的GRIN透镜输入面z方向输入信号cos(5y2)，仿真结果如图3(a)(c)所示：其中 图(a)为输入该信号电场绝对值的场图，X轴、Y轴为透镜所在二维笛卡尔坐标系下的坐标 轴，可以看出GRIN透镜的矩形区域中，入射波在透镜中弯曲，并且在信号所对应的阶次聚焦， 即为透镜区域电场的最大值；图(c)X轴为图(a)中的X轴，Y轴为的电场的绝对值，图中 曲线表示的是在y＝0时所对应的截线的电场的绝对值分布，从图(c)可以测出电场绝对值最 大值所对应的X轴的位置，即为公式(15)中的l，则通过公式(14)和(15)可以计算出输 入信号的阶次和调频率如表2所示；
在本装置z方向输入信号cos(5y2)，仿真结果如图3(b)(d)所示：其中图(b)为输入该信号在本装置中的场图，X轴、Y轴为透镜所在二维笛卡尔坐标系下的坐标轴，可以看出在本装置的扇形区域内，入射波在本装置中弯曲，并且在其所对应的阶次聚焦，即为透镜区域电场的最大值，通过聚焦位置可以看到本装置与传统的GRIN透镜相比聚焦位置在透镜的中部；图(d)X轴为图(b)中的X轴，Y轴为的电场的绝对值，图中曲线表示的是在y＝0时所 对应的截线的电场的绝对值分布，从图(d)可以测出电场绝对值最大值所对应的X轴的位置， 即为公式(15)中的l，则通过公式(18)可以计算出输入信号的阶次和调频率如表1所示；
表1 GRIN透镜和本装置在输入低调频率chirp信号仿真结果对比
从表1中，可以看出本装置测量输入信号的调频率相对误差未0.91％，比传统的GRIN透镜 检测的调频率的相对误差22.31％要低很多，所以本装置在检测低调频率信号的调频率较传统 的GRIN透镜更精确。
实施例4
本实施例按照实施例3中所述的参数设置，利用利用COMSOL Multiphysics进行计算机仿 真。
本施例的输入信号为f(y)＝cos(8y2)+cos(6y2)，如图4(a)所示，其中X轴为输入平面 的长度，Y轴表示电场的绝对值值，图中曲线表示输入平面每一点所对应的电场的绝对值； 当图4(a)中信号输入到传统GRIN透镜和本装置中，其仿真结果如图4(c)(d)所示，其中X轴和Y轴分别为两个透镜所在区域的笛卡尔坐标系的坐标轴，矩形和扇形分别表示GRIN透镜和本装置的形状，颜色表示该区域的场强的绝对值；从图(c)中可以看出，当输入f(y)到GRIN透镜中，该信号聚焦为一个点，图4(d)中可以看出本装置在输入同一信号后，聚 焦为两个点；由于输入的信号是关于X轴对称的信号，因此再两个透镜中分别做一条Y＝0的平行于横轴的截线，截线的电场图如图4(e)所示，其中，图(e)中X轴为两个透镜所在 区域中的X轴，Y轴为的电场的绝对值，图(e)中曲线表示的是GRIN透镜中平行于横轴的Y＝0的截线上每一点的电场的绝对值分布，从图(e)可以测出电场绝对值最大值所对应的X 轴的位置，从图(e)可以看到电场值有两个包络，分别测出两个包络中电场绝对值最大值所 对应的X轴的位置，通过公式(18)可以计算出输入信号的阶次和调频率所示；两个透镜的 仿真结果如表2所示：
表2输入信号为f(y)＝cos(8y2)+cos(6y2)时GRIN透镜和本装置仿真结果对比
从表2可以看出本装置可以测量出输入信号的两个调频率，而传统的GRIN透镜只能测量 出一个调频率，在分辨同幅的低调频率chirp信号上，本装置具有更大的优势。
以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明 的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之 内。
Claims (4)
1.一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，其特征在于：核心思想是：以介质的特性作为手段完成对光场的控制，具体主要是基于光波在介质中的传播的性质，当预先知道输入chirp信号在分数阶傅里叶变换装置中的处理效果，即明确其光场分布，则可以利用变换材料的方法确定装置的介质参数，从而达到预期的光场分布；
一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，包括分数阶傅里叶变换透镜的装置和装置的设计方法，其中，基于变换材料的分数阶傅里叶变换透镜的装置，简称本装置，本装置的设计方法，简称本方法；
其中，本装置主要包括输入面、基于变换材料设计的二次折射率透镜和输出面；
本装置中基于变换材料设计的二次折射率透镜，简称透镜；
输入面是输入信号进入二次折射率透镜的入射面，由于本装置是基于变换材料的原则，虽然本装置的输入面与原始的GRIN透镜输入面相比，形状发生了拉伸与改变，但在数学意义与物理意义上，本装置的输入面的点与原始的GRIN透镜输入面为一一对应的关系；
本装置中透镜与传统的GRIN透镜相比，通过变换材料的方法，改变传统的GRIN透镜中原有的轴向与径向梯度折射率的分布，使透镜的折射率呈类似于环状分布，主要表现为从输入面到输出面，透镜折射率依次呈环状增大，并在此过程中透镜的环状折射率分布所形成的环形半径逐渐减小；透镜折射率的改变导致了光在透镜的传播路径的改变；本装置在满足折射率分布的情况下，其形状可以不受本装置示出形状的限制；在信号输入到本装置的输入面，进入透镜后，在透镜中是一个分数阶傅里叶变换的过程，最终在输出面形成一个阶次为1的输入信号的分数阶傅里叶变换，即为普通的傅里叶变换，在透镜内部则形成了阶次不同的输入信号分数阶傅里叶变换；本装置是在传统的二次折射率透镜的基础上扩大输入面，保持输出面不变或缩小的共形变换，本装置可以通过复变函数的保角变换等方法推导出本装置和传统二次折射率透镜的变换关系，再利用变换材料的方法获得本装置的折射率；
输出面：即频谱面，输入光学信号经过二次折射率透镜之后，出射波为输入波的一个傅里叶变换信号。
2.根据权利要求1所述的一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，其特征在于：本方法，通过以下步骤实现：
步骤1：在二维平面下，以传统的GRIN透镜为基础，保证透镜输出面和透镜的厚度不变，弯曲透镜；
其中，透镜的高度所指方向z方向，且步骤1假设z方向是无限延伸的；步骤1中的透镜弯曲角度为γ，γ弯曲程度可为0度到360度；
步骤1中透镜在二维平面主伸长表示为公式(1)：
<mrow>
<msubsup>
<mi>&lambda;</mi>
<mi>&theta;</mi>
<mi>I</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mi>&gamma;</mi>
</mrow>
<mi>h</mi>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>r</mi>
<mi>I</mi>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，h为透镜的厚度，θ和r为步骤1透镜弯曲后透镜中的点的极坐标，分别为传统的GRIN透镜在步骤1弯曲后在极坐标下θ和r的主伸长；
步骤2：透镜在步骤1的基础上，沿透镜中心到坐标r的方向拉伸，拉伸至透镜中的各点各向同性；
步骤3：通过变换材料的方法，基于步骤2计算出拉伸后透镜的介电常数和磁导率；
至此，经过步骤1到步骤3，完成了本方法。
3.根据权利要求2所述的一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，其特征在于：经过步骤2后透镜在二维平面的主伸长表示为公式(2)：
<mrow>
<msubsup>
<mi>&lambda;</mi>
<mi>&theta;</mi>
<mrow>
<mi>I</mi>
<mi>I</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<msup>
<mi>d&theta;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<mi>r</mi>
<mi>d</mi>
<mi>&theta;</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mi>r</mi>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>r</mi>
<mi>I</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>dr</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，二维平面是r和θ组成的二维平面；θ′和r′为步骤2拉伸后透镜中的点的极坐标，分别为步骤1中的透镜在步骤2拉伸后在二维平面的主伸长，表示θ′对θ求偏导，表示r′对r求偏导，；
其中，步骤2要拉伸至各点各向同性，主伸长要满足公式(3)：
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>dr</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mi>r</mi>
</mrow>
<mi>h</mi>
</mfrac>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，λ为分别经过步骤1弯曲和步骤2拉伸后透镜中各点在各个方向上的主伸长；
由于步骤1弯曲和步骤2拉伸过程中始终保持透镜输出面不变，则GRIN透镜和经过步骤1和步骤2后的透镜的变换关系满足公式(4)和(5)：
<mrow>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mfrac>
<mi>&gamma;</mi>
<mi>h</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>r</mi>
<mo>=</mo>
<msup>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>h</mi>
<mi>&gamma;</mi>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，R1为输出面所在圆弧的半径，f(r)为步骤1弯曲和步骤2拉伸后透镜中的各点与传统的GRIN透镜中各点的函数关系，exp[]为指数函数，f-1(r')为f(r)函数的逆函数，ln为以e为底的对数；表示以e为底，对取对数；公式(4)也体现了本装置中输入面与原始的GRIN透镜输入面上的点是一一对应的对应关系。
4.根据权利要求2所述的一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法，其特征在于：步骤3中，变换材料的方法中材料介电常数和磁导率变化前后关系为公式(6)和公式(7)：
<mrow>
<msup>
<mi>&epsiv;</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msub>
<mi>&epsiv;</mi>
<mn>0</mn>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&mu;</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，ε′和μ′为变换后材料的介电常数和磁导率，ε0和μ0为变换前材料的介电常数和磁导率，λ1,λ2,λ3为三维坐标系下的主伸长；
由于步骤1中假设z方向上无限延伸，即假设在z方向没有任何变形，则在z方向上的主伸长λz＝1，通过公式(6)和(7)，可以得出经过步骤1弯曲和步骤2拉伸后透镜的介电常数、磁导率和折射率的关系，即公式(8)：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&epsiv;</mi>
<mi>r</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&epsiv;</mi>
<mi>&theta;</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>&epsiv;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&mu;</mi>
<mi>r</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&mu;</mi>
<mi>&theta;</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&epsiv;</mi>
<mi>z</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>&epsiv;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&mu;</mi>
<mi>z</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>n</mi>
<mi>r</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>n</mi>
<mi>&theta;</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>n</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，变形主要指弯曲和拉伸，也可以指压缩公式(8)中的ε0、μ0、n0为变换前透镜的介电常数、磁导率和折射率，εr′εθ′εz′为r，θ和z方向的介电常数，μ′r,μ′θ,μ′z为r，θ和z方向的磁导率，n′r,n′θ为r，θ方向的折射率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710979945.6A CN107703626B (zh) | 2017-10-19 | 2017-10-19 | 一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710979945.6A CN107703626B (zh) | 2017-10-19 | 2017-10-19 | 一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107703626A true CN107703626A (zh) | 2018-02-16 |
| CN107703626B CN107703626B (zh) | 2020-04-17 |
Family
ID=61182896
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710979945.6A Expired - Fee Related CN107703626B (zh) | 2017-10-19 | 2017-10-19 | 一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107703626B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109101727A (zh) * | 2018-08-13 | 2018-12-28 | 厦门大学 | 一种基于变换光学的共形天线设计方法 |
| CN109781153A (zh) * | 2018-12-17 | 2019-05-21 | 北京理工大学 | 物理参数估计方法、装置和电子设备 |
| CN111207780A (zh) * | 2019-12-26 | 2020-05-29 | 内蒙古科技大学 | 一种基于cft的牛顿环参数估计方法 |
| CN111628264A (zh) * | 2020-06-11 | 2020-09-04 | 太原理工大学 | 一种用普通电介质构建近零折射率材料的方法 |
| CN113433690A (zh) * | 2021-06-25 | 2021-09-24 | 北京理工大学 | 基于变换材料生成涡旋光束的光学器件装置和设计方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3425770A (en) * | 1965-12-01 | 1969-02-04 | Technical Operations Inc | Superimposed photostorage and separation |
| CN1866077A (zh) * | 2005-05-20 | 2006-11-22 | 清华大学 | 傅里叶变换光学系统及体全息存储傅里叶变换光学系统 |
| CN101223747A (zh) * | 2005-07-13 | 2008-07-16 | 文卡达·古鲁普拉赛德 | 使用均匀采样测谱法的依赖于距离的波谱 |
-
2017
- 2017-10-19 CN CN201710979945.6A patent/CN107703626B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3425770A (en) * | 1965-12-01 | 1969-02-04 | Technical Operations Inc | Superimposed photostorage and separation |
| CN1866077A (zh) * | 2005-05-20 | 2006-11-22 | 清华大学 | 傅里叶变换光学系统及体全息存储傅里叶变换光学系统 |
| CN101223747A (zh) * | 2005-07-13 | 2008-07-16 | 文卡达·古鲁普拉赛德 | 使用均匀采样测谱法的依赖于距离的波谱 |
Non-Patent Citations (1)
| Title |
|---|
| VINCENT CROZATIER.ETL: "High-Resolution Radio Frequency Spectral Analysis", 《IEEE JOURNAL OF QUANTUM ELECTRONICS》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109101727A (zh) * | 2018-08-13 | 2018-12-28 | 厦门大学 | 一种基于变换光学的共形天线设计方法 |
| CN109101727B (zh) * | 2018-08-13 | 2020-06-26 | 厦门大学 | 一种基于变换光学的共形天线设计方法 |
| CN109781153A (zh) * | 2018-12-17 | 2019-05-21 | 北京理工大学 | 物理参数估计方法、装置和电子设备 |
| CN109781153B (zh) * | 2018-12-17 | 2020-08-11 | 北京理工大学 | 物理参数估计方法、装置和电子设备 |
| US11156450B2 (en) | 2018-12-17 | 2021-10-26 | Beijing Institute Of Technology | Method, device and electronic apparatus for estimating physical parameter by discrete chirp fourier transform |
| CN111207780A (zh) * | 2019-12-26 | 2020-05-29 | 内蒙古科技大学 | 一种基于cft的牛顿环参数估计方法 |
| CN111628264A (zh) * | 2020-06-11 | 2020-09-04 | 太原理工大学 | 一种用普通电介质构建近零折射率材料的方法 |
| CN113433690A (zh) * | 2021-06-25 | 2021-09-24 | 北京理工大学 | 基于变换材料生成涡旋光束的光学器件装置和设计方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107703626B (zh) | 2020-04-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107703626A (zh) | 一种基于变换材料的光学分数阶傅里叶变换的装置和设计方法 | |
| Su et al. | On tilt and curvature dependent errors and the calibration of coherence scanning interferometry | |
| Hirayama et al. | Rigorous electromagnetic analysis of diffractive cylindrical lenses | |
| Mikuła et al. | Imaging with extended focal depth by means of lenses with radial and angular modulation | |
| Zhou et al. | Improved method for rapid shape recovery of large specular surfaces based on phase measuring deflectometry | |
| Mejía-Barbosa et al. | Object surface for applying a modified Hartmann test to measure corneal topography | |
| CN111999042B (zh) | 检测光学系统任意波长透射波前的方法 | |
| Yan et al. | Experimental study on subaperture testing with iterative triangulation algorithm | |
| Worku et al. | Decomposition of a field with smooth wavefront into a set of Gaussian beams with non-zero curvatures | |
| Xu et al. | Fabrication of continuous phase plates with small structures based on recursive frequency filtered ion beam figuring | |
| Yao et al. | Three-dimensional shape measurement with an arbitrarily arranged projection moiré system | |
| Wang et al. | Calibration of geometrical aberration in transmitted wavefront testing of refractive optics with deflectometry | |
| Shikder et al. | Measurement of the fractional topological charge of an optical vortex beam through interference fringe dislocation | |
| Pant et al. | Weighted spline based integration for reconstruction of freeform wavefront | |
| Gutiérrez et al. | Chromatic aberration in metalenses | |
| Feng et al. | Optimal strategy for fabrication of large aperture aspheric surfaces | |
| Zhang et al. | Speckle pattern shifting deflectometry based on digital image correlation | |
| Cai et al. | Analysis of wavefront errors introduced by encoding computer-generated holograms | |
| CN106895784B (zh) | 一种基于液晶材料的计算全息图及其制作方法 | |
| CN103777263B (zh) | 亚波长矩形单周期光栅结构的制作方法及偶数分束器 | |
| Gonzalez-Utrera et al. | Modeling, fabrication, and metrology of 3D printed Alvarez lenses prototypes | |
| CN109884788B (zh) | 消球差连续变焦双胶合液芯柱透镜系统及应用 | |
| KR101166961B1 (ko) | 3차원 빔측정을 통한 광학계의 파면수차 복원방법 | |
| Liao et al. | Rapid fabrication technique for nanometer-precision aspherical surfaces | |
| Ji et al. | Analysis of the misalignment effect and the characterization method for imprinting continuous phase plates |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200417 Termination date: 20201019 |