主权项 |
一种基于二阶相干的热光测距方法,其特征在于包括以下步骤:步骤一、根据输入的时间序列标签进行循环的符合运算,得到曲线拟合所需要的一系列归一化的时间二阶相干函数的样本点;在某个延时τ下,两路时间序列的符合数和二阶相干函数g(2)(τ)有如下的关系: <mrow> <msup> <mi>g</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>τ</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>n</mi> <mrow> <mo>(</mo> <mi>τ</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>Tδ</mi> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>γ</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>γ</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>Tδ</mi> <msub> <mi>R</mi> <mn>1</mn> </msub> <msub> <mi>R</mi> <mn>2</mn> </msub> </mrow> <mrow> <mi>Tδ</mi> <msub> <mi>R</mi> <mn>1</mn> </msub> <msub> <mi>R</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>式中,T是总采样时间,δ为符合门宽,R1和R2为探测器1和探测器2的有用光子计数率,r1和r2为对应探测器的暗计数率和环境噪声引起的计数率之和,n(τ)为符合数;且当Ri>>ri(i=1,2)时,方程(1)简化为: <mrow> <msup> <mi>g</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>τ</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>n</mi> <mrow> <mo>(</mo> <mi>τ</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>Tδ</mi> <msub> <mi>R</mi> <mn>1</mn> </msub> <msub> <mi>R</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>在一路时间序列CH0标签中每个标签值i上加上某个延时τ;选取另外一路CH1的时间标签序列作为基准,将其中每一个时间标签作为一个时间区间的中点,而时间区间的长度即为符合门宽;在每一个时间区间内搜索有延时τ的那路的时间标签,若有时间标签落入区间内,则记符合数加一;搜索完信号路CH0所有的时间标签后,即得到某个延时τ下的符合数n(τ),并根据方程(2)计算这个延时τ下归一化的时间二阶相干函数g(2)(τ),完成一次符合计算;调整延时τ重新计算符合数,重新进行符合计算,直到满足曲线拟合所需的样本点数;步骤二、对于实验数据(xi,yi)(i=1,2,...,m),根据选择拟合函数w(x),使得误差ri=w(xi)‑yi平方和最小,使下式达到极小值: <mrow> <mi>u</mi> <mrow> <mo>(</mo> <msub> <mi>α</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mo>,</mo> <msub> <mi>α</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msubsup> <mi>r</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msup> <mrow> <mo>[</mo> <mi>w</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>]</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>式中,αi,i=1,…,k为w(x)中需要拟合的k个参数;由旋转毛玻璃所形成的赝热光场据具有如下的二阶相干的函数: <mrow> <msup> <mi>g</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>τ</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mi>f</mi> <mrow> <mn>1</mn> <mo>+</mo> <mi>f</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>τ</mi> <mo>-</mo> <mi>Δτ</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>/</mo> <msup> <mi>p</mi> <mn>2</mn> </msup> </mrow> </msup> <mo>+</mo> <mfrac> <mrow> <mn>2</mn> <mi>f</mi> </mrow> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>f</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>τ</mi> <mo>-</mo> <mi>Δτ</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>/</mo> <msup> <mi>q</mi> <mn>2</mn> </msup> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>式中,f为赝热光与相干光光强的比值,Δτ为二阶相干函数峰值所对应的值,p和q分别为相干光场和赝热光场的照度参数;根据符合方法输出一组二阶相干函数的样本值,通过曲线拟合方法,拟合出二阶相干函数峰值所对应的Δτ的值。 |