主权项 |
1.用于检测非球面的高精度位相型计算全息图的编码方法,其特征是:它由以下步骤实现:步骤一、根据被测旋转对称非球面的二次曲面偏心率e,确定旋转对称非球面的类型;然后根据公式:<maths num="0001"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msup><mi>cS</mi><mn>2</mn></msup><mrow><mn>1</mn><mo>+</mo><msqrt><mn>1</mn><mo>-</mo><mrow><mo>(</mo><mi>L</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><msup><mi>c</mi><mn>2</mn></msup><msup><mi>S</mi><mn>2</mn></msup></msqrt></mrow></mfrac><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><msup><mi>S</mi><mn>4</mn></msup><mo>+</mo><msub><mi>A</mi><mn>2</mn></msub><msup><mi>S</mi><mn>6</mn></msup><mo>+</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>+</mo><msub><mi>A</mi><mi>n</mi></msub><msup><mi>S</mi><mrow><mn>2</mn><mi>n</mi></mrow></msup></mrow></math>]]></maths>生成被测旋转对称非球面模型W(x,y),单位为微米(um);式中:S<sup>2</sup>=x<sup>2</sup>+y<sup>2</sup>;(x,y)为被测旋转对称非球面上任意一点在底面投影的坐标;<img file="FDA0000468970870000012.GIF" wi="196" he="160" />R<sub>0</sub>为被测旋转对称非球面的顶点曲率半径;A<sub>1</sub>、A<sub>2</sub>、LA<sub>n</sub>为高次非球面各高次项系数;L=-e<sup>2</sup>;步骤二、对被测旋转对称非球面模型W(x,y)采用位相型计算全息编码;具体为:首先对被测旋转对称非球面模型W(x,y)进行抽样离散化处理,获得W(i,j);<maths num="0002"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><mi>W</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>·</mo><mi>comb</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>N</mi></mfrac><mo>)</mo></mrow><mo>·</mo><mi>comb</mi><mfrac><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mi>N</mi></mfrac><mo>·</mo><mi>rect</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>D</mi></mfrac><mo>,</mo><mfrac><mi>y</mi><mi>D</mi></mfrac><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>式中:i,j=1,2,3LN,N为采样数,D为被测旋转对称非球面的口径;然后以2π为模对W(i,j)进行周期量化,获得位相图,通过放大取整后,将位相图转换为灰阶图,灰阶分布范围为0-255,生成位相型计算全息图W<sub>CGH</sub>(i,j);<maths num="0003"><![CDATA[<math><mrow><msub><mi>W</mi><mi>CGH</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>round</mi><mo>[</mo><mn>255</mn><mo>·</mo><mfrac><mn>1</mn><mi>λh</mi></mfrac><mo>·</mo><mi>mod</mi><msub><mrow><mo>[</mo><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>·</mo><mi>h</mi><mo>]</mo></mrow><mrow><mn>2</mn><mi>π</mi></mrow></msub><mo>]</mo><mrow><mo>(</mo><mi>h</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>λ</mi></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>式中:λ为波长,单位为微米(um);步骤三、提取步骤二获得的位相型计算全息图W<sub>CGH</sub>(i,j)上包含对称中心点的一维数组B<sub>i</sub>,即:B<sub>1</sub>、B<sub>2</sub>……B<sub>N</sub>,从一维数组B<sub>i</sub>中提取每个量化周期的最小边界值b<sub>k</sub>,建立线性方程组:z<sub>k</sub>=b<sub>k</sub>-b<sub>0</sub>并求解最小边界值偏差z<sub>k</sub>;式中:b<sub>0</sub>为期望值,k为量化周期数,k=1,2,3…C;C是最大量化周期数;步骤四、建立线性方程组:<maths num="0004"><![CDATA[<math><mrow><msub><mi>β</mi><mi>k</mi></msub><mo>=</mo><mfrac><msub><mi>a</mi><mn>0</mn></msub><mrow><msub><mi>a</mi><mi>k</mi></msub><mo>-</mo><msub><mi>z</mi><mi>k</mi></msub></mrow></mfrac><mo>-</mo><msub><mi>b</mi><mn>0</mn></msub></mrow></math>]]></maths>求取:每个量化周期的调整系数β<sub>k</sub>;式中:a<sub>k</sub>为各周期最大边界值,a<sub>0</sub>为期望值;步骤五、通过非球面量化周期的边界定义式:<img file="FDA0000468970870000022.GIF" wi="1423" he="230" />提取位相型计算全息图W<sub>CGH</sub>(i,j)每个量化周期数据W<sub>k</sub>(i,j);式中:λ为波长,单位为微米(um);步骤六、用一维数组调整系数β<sub>k</sub>通过公式:R<sub>k</sub>(i,j)=(1+β<sub>k</sub>)[W<sub>k</sub>(i,j)-z<sub>k</sub>(i,j)]调整对应量化周期数据W<sub>k</sub>(i,j),得到调整后量化周期数据R<sub>k</sub>(i,j);其中:最小边界值偏差矩阵z<sub>k</sub>(i,j)由一维数组z<sub>k</sub>构造:<img file="FDA0000468970870000023.GIF" wi="1134" he="198" />再累加调整后每个量化周期数据:<maths num="0005"><![CDATA[<math><mrow><msub><mi>W</mi><mrow><mi>CGH</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>round</mi><mrow><mo>(</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>C</mi></munderover><msub><mi>R</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>生成校正后位相型计算全息图W<sub>CGH1</sub>(i,j);步骤七、判断步骤六生成的校正后位相型计算全息图W<sub>CGH1</sub>(i,j)的最小边界值偏差是否满足|z<sub>k</sub>|<z<sub>min</sub>,z<sub>min</sub>为最小边界值偏差期望值,如果判断结果为否,则返回执行步骤三;如果判断结果为是,则完成高精度位相型计算全息图的编码。 |