发明名称 用于改善焦线均匀度的柱面透镜列阵的设计方法
摘要 一种用于改善焦线均匀度的柱面透镜列阵的设计方法,采用了非时齐模拟退火算法,改变经典柱面透镜列阵设计上各列阵元结构的一致性,使得新型的柱面透镜列阵由宽度、焦距不等的柱面透镜列阵元组成。本发明设计的子列阵元宽度、焦距保持焦面处的子焦线长度一致,最终这些子焦线在焦面处的迭加得到一个强度分布近乎平顶的焦线,从而实现焦线的空间强度分布均匀化。
申请公布号 CN102566046A 申请公布日期 2012.07.11
申请号 CN201110342353.6 申请日期 2011.11.03
申请人 中国科学院上海光学精密机械研究所 发明人 张琥杰;周申蕾;鄔融;林尊琪
分类号 G02B27/00(2006.01)I;G02B3/06(2006.01)I 主分类号 G02B27/00(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张泽纯
主权项 一种用于改善焦线均匀度的柱面透镜列阵的设计方法,其特征在于该方法包括下列步骤:①确定入射光空间光强分布A(x,y),工作波长λ,聚焦透镜焦距F,柱面透镜列阵的宽度D,目标焦线长度L;②选定柱面透镜阵列元序列宽度为优化变量,焦线均匀度为优化目标函数;③设定柱面透镜列阵的单元数N,N的通常取值范围为5~8;④柱面透镜列阵元的初始宽度di=D/N,焦距fi=Fdi/L,设定阵列元排列方向为x方向;柱面透镜初始位相函数为: <mrow> <mi>&Phi;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mi>&Sigma;</mi> <mi>i</mi> </munder> <mi>rect</mi> <mo>[</mo> <mrow> <mo>(</mo> <mn>2</mn> <mi>x</mi> <mo>+</mo> <mi>D</mi> <mo>-</mo> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>]</mo> <mi>exp</mi> <mo>[</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mi>&pi;</mi> <mrow> <mi>&lambda;</mi> <msub> <mi>f</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>D</mi> <mo>/</mo> <mn>2</mn> <mo>-</mo> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>]</mo> </mrow>靶面上焦线光强分布为: <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mi>&lambda;</mi> <mn>2</mn> </msup> <msup> <mi>f</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <mo>&Integral;</mo> <mo>&Integral;</mo> <mi>A</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mi>&Phi;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>[</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> <mi>&lambda;F</mi> </mfrac> <mi>xu</mi> <mo>]</mo> <mi>dxdy</mi> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mrow>计算初始条件下的目标函数rms <mrow> <mi>rms</mi> <mo>=</mo> <msup> <mrow> <mo>[</mo> <mfrac> <mrow> <mo>&lang;</mo> <msup> <mrow> <mo>(</mo> <mi>I</mi> <mo>-</mo> <mover> <mi>I</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&rang;</mo> </mrow> <msup> <mover> <mi>I</mi> <mo>&OverBar;</mo> </mover> <mn>2</mn> </msup> </mfrac> <mo>]</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mo>;</mo> </mrow>⑤确定初始温度T0,降温速度C,迭代次数P;在温度T0下,各列阵元宽度的改变量Δdi=AT0random(N),A为常系数以控制该变量的大小,random(N)为一组数量为N,范围‑1~1的随机数,并满足∑random(N)=0,新生成的柱面透镜列阵元宽度di′=di+Δdi,且满足di′≥dlim,dlim为设计限定的最小宽度,对应的柱面透镜列阵焦距fi′=Fdi′/L;⑥新生成柱面透镜的位相函数: <mrow> <msup> <mi>&Phi;</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mi>&Sigma;</mi> <mi>i</mi> </munder> <mi>rect</mi> <mo>[</mo> <mrow> <mo>(</mo> <mn>2</mn> <mi>x</mi> <mo>+</mo> <mi>D</mi> <mo>-</mo> <msup> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>&prime;</mo> </msup> <mo>)</mo> </mrow> <mo>/</mo> <msup> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>&prime;</mo> </msup> <mo>]</mo> <mi>exp</mi> <mo>[</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mi>&pi;</mi> <mrow> <mi>&lambda;</mi> <msub> <mi>f</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>D</mi> <mo>/</mo> <mn>2</mn> <mo>-</mo> <msup> <msub> <mi>d</mi> <mi>i</mi> </msub> <mo>&prime;</mo> </msup> <mo>/</mo> <mn>2</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>]</mo> </mrow>靶面上焦线光强分布: <mrow> <msup> <mi>I</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mi>&lambda;</mi> <mn>2</mn> </msup> <msup> <mi>f</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <mo>&Integral;</mo> <mo>&Integral;</mo> <mi>A</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <msup> <mi>&Phi;</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mo>[</mo> <mo>-</mo> <mi>j</mi> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> <mi>&lambda;F</mi> </mfrac> <mi>xu</mi> <mo>]</mo> <mi>dxdy</mi> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mrow>计算初始温度下的目标函数rms’: <mrow> <msup> <mi>rms</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <msup> <mrow> <mo>[</mo> <mfrac> <mrow> <mo>&lang;</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>I</mi> <mo>&prime;</mo> </msup> <mo>-</mo> <mover> <msup> <mi>I</mi> <mo>&prime;</mo> </msup> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&rang;</mo> </mrow> <msup> <mover> <msup> <mi>I</mi> <mo>&prime;</mo> </msup> <mo>&OverBar;</mo> </mover> <mn>2</mn> </msup> </mfrac> <mo>]</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mo>;</mo> </mrow>⑦若rms′<rms,则接受该组新生成的列阵元宽度、焦距数据并接受rms=rms′;⑧若rms′≥rms,则以Metropolis准则以一定概率接受新生成的列阵元宽度、焦距数据并接受rms=rms′,接受概率p=exp[‑A′(rms′‑rms)/T0],A′为常系数调节接受概率;⑨系统退火冷却过程,Tm=CTm‑1,降温速度C控制冷却进度:在温度Tm下重复步骤6,7,8的操作,随着温度的降低,优化搜索范围Δdi缩小,同时Metropolis准则的接受概率p降低,最终在低温时收敛;⑩迭代P次,结束运算,选取最优解,获得最优的柱面透镜列阵元的宽度di,焦距fi。
地址 201800 上海市嘉定区800-211邮政信箱