发明名称 一种用线阵CCD测量光束中心位置的方法
摘要 一种用线阵CCD测量光束中心位置的方法,涉及一种测量入射到线阵CCD表面的光束中心位置的光电测量方法及运用该方法测位移和角位移的具体方法。本发明使CCD表面光斑的辐射功率密度呈高斯分布或准高斯分布;由CCD、相关电路特性等定出数字信号输出的误差限并选定截尾阈值;对截尾准高斯分布的有效数字信号用加权回归法计算出光束中心位置估值;用蒙特卡罗法数值模拟误差分布定出光束宽度范围。本发明方法可使线阵CCD测光束中心位置的有效分辨率和精密度从1W左右降到0.03W以下(W为像元间距),以实现亚微米精密度位移及1秒有效分辨率角度的测量,并使其能被运用于冲击电流计转镜转角、线状光谱波长等其它相关量的测量领域。
申请公布号 CN1333231C 申请公布日期 2007.08.22
申请号 CN200510012066.3 申请日期 2005.07.01
申请人 清华大学 发明人 朱鹤年
分类号 G01B11/02(2006.01);G01B11/26(2006.01) 主分类号 G01B11/02(2006.01)
代理机构 代理人
主权项 1.一种用线阵CCD测量光束中心位置的方法,其特征在于:1)使CCD探测器表面的光斑沿线阵方向的幅射功率线密度呈高斯分布或准高斯分布;2)求出截尾阈值:设CCD探测器第j个中心坐标为x<sub>j</sub>的像元的数字信号输出为V<sub>Dj</sub>,所述V<sub>Dj</sub>为LSB的整数倍;根据CCD探测器及相关电路特性、输出比特数N及噪声幅度,用公式<maths num="001"><![CDATA[ <math> <mrow> <mo>&PlusMinus;</mo> <msub> <mi>U</mi> <msub> <mi>V</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <mo>&PlusMinus;</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>c</mi> <mo>%</mo> <msub> <mi>V</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow> </math> ]]></maths>表示出V<sub>Dj</sub>的误差限,式中a为正常数,c%为下比例系数;用(a+c%×2<sup>N</sup>)到2<sup>N</sup>/20之间的某一整数作截尾阈值V<sub>th</sub>,使截尾后不小于阈值V<sub>th</sub>的像元的有效数字输出信号V<sub>Dj</sub>呈截尾准高斯分布;3)对中心坐标为x<sub>j</sub>的像元的截尾准高斯分布的有效数字信号输出V<sub>Dj</sub>,采用回归方程模型<maths num="002"><![CDATA[ <math> <mrow> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>=</mo> <msub> <mi>b</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mn>2</mn> </msub> <msubsup> <mi>x</mi> <mi>j</mi> <mn>2</mn> </msubsup> </mrow> </math> ]]></maths>,分别以x<sub>j</sub>和x<sub>j</sub><sup>2</sup>为自变量、以y<sub>j</sub>=ln(V<sub>Dj</sub>)为因变量,作权因子为<maths num="003"><![CDATA[ <math> <mrow> <msub> <mi>w</mi> <mi>yj</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>V</mi> <mi>Dj</mi> </msub> <mo>/</mo> <msub> <mi>U</mi> <msub> <mi>v</mi> <mi>j</mi> </msub> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </math> ]]></maths>的加权回归求b<sub>1</sub>和b<sub>2</sub>,再算出光束中心位置估值<maths num="004"><![CDATA[ <math> <mrow> <msub> <mover> <mi>x</mi> <mo>^</mo> </mover> <mi>c</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>/</mo> <mrow> <mo>(</mo> <msub> <mrow> <mn>2</mn> <mi>b</mi> </mrow> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </math> ]]></maths>;4)用蒙特卡罗法数值模拟误差分布定出光束宽度范围,其具体步骤是:A).像元间隔W已知时,对于光束中心位置x<sub>c</sub>,按幅射功率线密度的高斯分布E<sub>l</sub>(x<sub>c</sub>,x)=E<sub>lm</sub>exp(-(x-x<sub>c</sub>)<sup>2</sup>/(2σ<sup>2</sup>)),在某一分布标准差σ的取值下、对固定的采样时间间隔和像元内的位置变量x积分算出中心坐标为x<sub>j</sub>的第j个像元的模拟信号V<sub>j</sub>=V<sub>j</sub>(x<sub>j</sub>,x<sub>c</sub>),作为模数转换前的光电测量信号初值;V<sub>j</sub>是相对值,积分时选取比例常量使V<sub>j</sub>的峰值为0.8×2<sup>N</sup>;B).按照误差限为<maths num="005"><![CDATA[ <math> <mrow> <mo>&PlusMinus;</mo> <msub> <mi>U</mi> <msub> <mi>V</mi> <mi>j</mi> </msub> </msub> <mo>&PlusMinus;</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>c</mi> <mo>%</mo> <msub> <mi>V</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow> </math> ]]></maths>的规律数值模拟一组随机误差分布<maths num="006"><![CDATA[ <math> <mrow> <msub> <mi>&epsiv;</mi> <msub> <mi>V</mi> <mi>j</mi> </msub> </msub> <mo>=</mo> <msub> <mi>r</mi> <mi>aj</mi> </msub> <mo>+</mo> <msub> <mi>r</mi> <mi>cj</mi> </msub> <mo>%</mo> <msub> <mi>V</mi> <mi>j</mi> </msub> </mrow> </math> ]]></maths>,这早r<sub>aj</sub>和r<sub>cj</sub>分别是误差限为±a和±c的均匀分布随机数;取(V<sub>j</sub>+ε<sub>Vj</sub>)的整数部分作为模数转换后的数字信号输出V<sub>Dj</sub>;C).截取V<sub>Dj</sub>≥V<sub>th</sub>的有效数字信号;D).对方程<maths num="007"><![CDATA[ <math> <mrow> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>=</mo> <mn>1</mn> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>V</mi> <mi>Dj</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>b</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mn>2</mn> </msub> <msubsup> <mi>x</mi> <mi>j</mi> <mn>2</mn> </msubsup> </mrow> </math> ]]></maths>作权因子为<maths num="008"><![CDATA[ <math> <mrow> <msub> <mi>w</mi> <mi>yj</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>V</mi> <mi>Dj</mi> </msub> <mo>/</mo> <msub> <mi>U</mi> <msub> <mi>v</mi> <mi>j</mi> </msub> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </math> ]]></maths>的加权回归,求出<maths num="009"><![CDATA[ <math> <mrow> <msub> <mover> <mi>x</mi> <mo>^</mo> </mover> <mi>c</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>b</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </math> ]]></maths>并算出误差<maths num="010"><![CDATA[ <math> <mrow> <msub> <mi>&epsiv;</mi> <mi>xc</mi> </msub> <mo>=</mo> <msub> <mover> <mi>x</mi> <mo>^</mo> </mover> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>c</mi> </msub> </mrow> </math> ]]></maths>来;E).对多组相同规律的模拟误差分布分别计算出各组的误差ε<sub>xc</sub>,进而算出误差的标准差s<sub>xc</sub>,此即与光束分布标差σ对应的光束中心位置x<sub>c</sub>的测量标准差;F).对不同的σ值重复上述A到E的步骤计算对应的标准差s<sub>xc</sub>,得出s<sub>xc</sub>/W和σ/W之间的关系曲线,进而找出使s<sub>xc</sub>/W不大于0.08的光束分布参量σ/W的取值范围。
地址 100084北京市北京100084-82信箱