删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于数字非线性锁相环的相干态相位估计

本站小编 Free考研考试/2021-12-29

摘要:基于量子理论获取相位参数的导航机制, 理论上可以突破经典物理极限对导航精度的限制. 利用量子零拍探测对相干态光场相位进行测量时, 通常需要相位与之正交的本振光才能使测量精度达到量子标准极限. 由于导航信号相位的高非线性特点, 想要利用传统的线性锁相环获取完全满足条件的本振光具有一定的难度. 为此, 本文设计了一种基于容积准则的非线性锁相环, 实现了在非正交本振光的条件下对相干态相位进行精确测量的功能. 首先, 利用相干态的Wigner函数推导了其相位在量子零拍探测的输出结果, 设计了量子相位估计的非线性数字锁相环框架. 然后基于正交单纯形容积准则设计了非线性滤波算法实现锁相环功能, 该锁相环通过对本振相位进行多次状态更新, 最终实现非线性迭代估计. 实验结果表明, 本文方法突破了本振光相位需与相干态相位正交的局限性, 避免了传统量子锁相环方法引入的线性化误差, 实现了对相干态相位的准确、稳定估计.
关键词: 相干态相位估计/
量子零拍探测/
非线性锁相环/
正交单纯形容积卡尔曼滤波

English Abstract


--> --> -->
相位估计是导航控制、图像处理、大气波导、混沌通信等领域的关键问题之一[1?4]. 例如, 在测角定位系统中, 对相位信息估计的精确度直接决定了波达角等参数获取的精度. 然而, 经典相位的精度在理论上有不可逾越的界限, 成为制约各类参数获取精度的瓶颈.
近年来, 随着量子技术的发展, 尤其是对连续变量量子相位测量的研究, 利用相干态相位信息获取相关参数, 有望获得超越经典理论极限的精度[5?7]. 目前, 已有研究通常是将相干态通过量子零拍探测器, 并建立基于相干态光场强度与其相位相关统计量的函数关系来进行相位估计[8,9], 这能够有效提高相位估计精度[10,11]. Wiseman[12]提出量子反馈测量方案, 即对量子零拍探测器的本振相位进行反馈控制, 从而实现近似正则测量, 并得到了验证[13]. 随后, Berry等[14]提出了相干态光场和压缩态光场连续变化相位的自适应测量跟踪方法, 证明了在相位跟踪中反馈测量优于非反馈测量. Tsang等[15]设计了自适应零拍线性锁相环对实时相位和瞬时频率进行测量, 并且在此基础上将非线性方程进行一阶泰勒展开, 利用卡尔曼-布什滤波实现了相位的实时跟踪[16,17]. 上述方法均是在本振光相位与待测光场相位正交的前提下能够达到标准量子极限, 如果信号为相位压缩态光场, 其相位测试精度能够突破标准量子极限. 然而对于导航信号来说, 相干态相位是未知的, 即使利用自适应锁相环对本振相位进行调整, 也并不能保证本振相位与相干态相位完全正交, 因此利用上述方法进行线性化处理时, 虽然能简化相位估计过程, 但也会不可避免地引入线性化误差, 并且在调整本振相位时带来了时延, 难以保证波达角等参数的实时获取, 这对于导航系统来说显然不合理.
为此, 本文在传统自适应反馈线性化连续锁相环的基础上, 设计了一种基于正交单纯形容积卡尔曼滤波的数字非线性锁相环来实现相干态相位的精确估计. 首先分析相干态相位的量子零拍探测输出模型, 然后利用AD变换将测量信息进行数字化处理, 再引入数据处理模块, 得到相位估计的数字锁相环框架. 设计数字锁相环的原因在于: 一是导航领域中所需的各类参数信息不一定是连续的, 因而无需估计连续的相位信息; 二是进行数字化更利于锁相环的实现与算法迭代; 三是可以避免连续数据隐藏的缺陷对模型的影响, 使模型更加稳定. 数据处理中, 为避免线性化误差, 采用基于容积准则[18,19]的非线性滤波算法. 在标准容积准则的基础上, 利用正则单纯形[20]来构建容积公式, 提高后验密度的数值积分精度; 然后根据考虑非线性对滤波的影响, 进一步利用正交矩阵将单纯形容积点进行变换, 调整非线性高阶项的影响, 从而避免局部采样效应, 得到正交单纯形容积卡尔曼滤波(OSCKF). 一方面, 上述方法避免了因本振光相位与信号相位无法完全正交引起的误差, 更易于处理且稳定性更好; 另一方面, OSCKF能够有效避免各类误差, 提高量子相位估计精度.
2
2.1.量子零拍测相
-->本文利用Wigner函数推导相干态相位在量子零拍探测的输出模型. 为了得到相干态的Wigner函数, 先分析真空态的Wigner函数. 定义正交分量场$x \equiv \dfrac{1}{2}\left( {a + {a^ + }} \right)$, $y \equiv \dfrac{1}{{2{\rm{i}}}}\left( {a - {a^ + }} \right)$, 式中, $a$${a^ + }$分别是光子湮灭算符与光子产生算符. 由不确定度原理知, $x$$y$的标准偏差满足$\Delta x\Delta y \geqslant \displaystyle\frac{1}{4}$.
记真空态的正交分量分别为${x_0}$${y_0}$, 则其Wigner函数可表示为
$W\left( {{x_0},{y_0}} \right) = \frac{{\rm{2}}}{{\text{π}}}\exp \left( { - 2{x_0}^2 - 2{y_0}^2} \right).$
从(1)式可以看出, ${x_0}$${y_0}$的方差均为$\dfrac{1}{4}$, 标准偏差$\Delta {x_0}\Delta {y_0} = \dfrac{1}{4}$, 可见真空态是最小不确定度态. 具有平均相位$\phi $的相干态$\left| \alpha \right\rangle $在相空间中可通过沿相位调制${x_0}$方向平移真空态得到, 其Wigner函数具体可以写为:
$W\left( {x,y\left| \phi \right.} \right) = W\left[ {{x_0}\left( {x,y} \right),{y_0}\left( {x,y} \right)} \right],$
${x_0} = x\cos \phi + y\sin \phi - \left| \alpha \right|,$
${y_0} = - x\sin \phi + y\cos \phi .$
零拍输出的正交场分量为$\eta = \dfrac{1}{2}\left( {{\rm{e}}^{ - {\rm{i}}\varphi }}a +\right.$$ \left.{{\rm{e}}^{{\rm{i}}\varphi }}{a^ +} \right)$, 其中$a$是相干态$\left| \alpha \right\rangle $的湮灭算符. 用正交分量可表示为$a = x + {\rm{i}}y$, 根据(1)式和(2)式, 可用真空态的正交分量表示为
$a = \left( {{x_0} + {\rm{i}}{y_0} + \left| \alpha \right|} \right){{\rm{e}}^{{\rm{i}}\phi }},$
则零拍输出$\eta $可写为
$\eta = \left| \alpha \right|\cos \left( {\phi - \varphi } \right) + {z_0},$
式中${z_0} = {x_0}\cos \left( {\phi - \varphi } \right) - {y_0}\sin \left( {\phi - \varphi } \right)$. 由(1)式知, 真空态的正交场分量服从零均值的高斯分布, 其统计分布与相位$\phi - \varphi $不相关, 则${z_0}$的统计特性服从高斯分布.可以看出, 正是${z_0}$引入了不确定性, 真空波动被分离到了${z_0}$中. 将${z_0}$处理成高斯噪声, 则(6)式变成高斯噪声环境下的参数估计问题, 经典估计技术能够很好地解决该问题. 首先将本振光的相位做${\text{π}/2}$偏置, 令$\varphi = {\varphi _{\rm{f}}} + \dfrac{\text{π}}{2}$, 这里定义${\varphi _{\rm{f}}}$为本振相位, 则(6)式可重写为
$\eta = \left| \alpha \right|\sin \left( {\phi - {\varphi _{\rm{f}}}} \right) + {z_0}',$
式中${z_0}' = {x_0}\sin \left( {\phi - {\varphi _{\rm{f}}}} \right) + {y_0}\cos \left( {\phi - {\varphi _{\rm{f}}}} \right)$, 与${z_0}$具有相同的统计特性. 之后将不再区分${z_0}$${z_0}'$, 统一记作高斯噪声${z_0}$. 在文献[16,17]中, 对(7)式线性化前均做了一个假设, 即
$\left\langle {{{\left( {\phi - {\varphi _{\rm{f}}}} \right)}^2}} \right\rangle < < 1.$
由此, 将(7)式进行线性近似可以得到
$\eta = \left| \alpha \right|\left( {\phi - {\varphi _{\rm{f}}}} \right) + {z_0}'.$
如果本振相位${\varphi _{\rm{f}}}$与相干态相位$\phi $较为接近, 则该假设成立, 利用上述线性化模型进行相位估计能够得到较好的结果. 然而对于导航系统来说, 目标的位置并不确定, 也就是说难以找到相干态相位$\phi $的准确大小. 因此无法保证本振相位${\varphi _{\rm{f}}}$与相位$\phi $满足上述假设条件. 如果本振相位与相干态相位相差过大, 则会引入较大的线性化误差, 甚至导致相位估计失败.
2
2.2.非线性数字锁相环设计
-->在量子零拍探测基础上, 本节建立非线性锁相环框架, 与传统线性锁相环不同的是, 这里采用AD变换引入非线性数据处理方法. 该锁相环能够突破传统方法中本振相位需要接近相干态相位的限制要求, 并避免线性化带来的误差, 进而实现相位的精确估计.
本文将非线性锁相环进行离散化的原因在于: 在导航系统中, 对各类参数的需求不一定连续, 得到的量测信息通常也是离散的, 因此通过AD采样得到足够的量测信息, 即可通过估计的手段得到所需参数. 离散化还可以将数据中隐藏的缺陷得以解决, 使模型结果更加稳定. 例如, 数据中的极端值是影响模型效果的一个重要因素. 极端值导致模型参数过高或过低, 或导致模型被虚假现象“迷惑”, 把原来不存在的关系作为重要模式来处理[21]. 而离散化处理可以有效地减弱极端值和异常值的影响. 此外, 离散化还有利于数值计算与算法迭代, 增强了物理可实现性.
设计的锁相环模型如图1所示. 由真空态进行平移和相位调制得到相干态, 其中${ D}(\left| \alpha \right|)$表示位移算符, $\exp ({\rm{i}} \phi )$表示相位调制, 然后将相干态通过初始本振相位为$\varphi $的量子零拍探测器得到输出$\eta $, 接着进行AD采样, 将所得数据送入数据处理模块(SP)中, 在数据处理后将更新的本振相位反馈给量子零拍探测器, 再重复此过程, 最终能够得到更为精确的估计结果.
图 1 数字零拍锁相环
Figure1. Digital homodyne phase-lock loop.

在所得锁相环框架的基础上, 在SP模块设计一种针对非线性模型的算法进行数据处理, 从而实现锁相环对相干态相位估计的功能. 因此, 所设计算法的好坏也会影响锁相环的相位估计性能.
得到数字零拍锁相环框架后, 问题就转化为将零拍数据离散化, 进而利用SP模块实现相干态相位估计. 对于离散非线性估计问题, 通常将非线性模型线性化, 并采用扩展卡尔曼滤波(EKF), 然而从(9)式可知, 量测模型是三角函数, 其非线性程度较高, 因此采用EKF会违背局部线性假设, 从而造成较大的线性化误差, 甚至导致滤波发散. 为此, 这里引入正交单纯形容积卡尔曼滤波(OSCKF)方法, 利用正交单纯形三阶球面径向容积准则, 计算相位信息对应的后验密度, 从而在不进行非线性模型泰勒展开、不引入线性化误差的前提下实现对相位$\phi $的精确实时估计.
2
3.1.系统模型和容积卡尔曼滤波框架
-->在SP模块中, 分别建立状态方程和量测方程为
$\quad{{ \phi} _k} = {F}{{ \phi} _{k - 1}} + {{ v}_{k - 1}}, $
${\eta _k} = h\left( {{{ \phi} _k}} \right) + {w_k}, $
其中状态向量${{ \phi} _k}$$k$时刻的相位信息, 状态转移矩阵${F} = 1$, 这里假定不含过程噪声, 即${{ v}_k} = 0$. 量测方程中, ${\eta _k}$$k$时刻量子零拍探测的量测值$h\left( {{{ \phi} _k}} \right) = \left| \alpha \right|\sin \left( {{{ \phi} _k} - {\varphi _{\rm{f}}}} \right)$, ${w_k}$为量测噪声, 服从均值为零的高斯白噪声.
CKF利用三阶球面径向容积准则进行数值积分, 从而近似后验密度对应的高维积分[22]. 其容积点可表示为${\xi _j} = \sqrt {{n_\phi }} {\left[ {\begin{array}{*{20}{c}} {{{I}_{{n_\phi }}}}&{- {{I}_{{n_\phi }}}} \end{array}} \right]_j}$, 其中${{I}_{{n_\phi }}}$${n_\phi }$维单位矩阵, ${\left[ \cdot \right]_j}$代表矩阵$\left[ \cdot \right]$的第$j$列, $j = 1,2, \cdots ,2{n_\phi }$. 在本问题中, 滤波过程如下.
1) 时间更新
由于系统的状态模型是线性的, 因此其预测状态${\hat \phi _{k\left| {k - 1} \right.}}$和预测协方差${{P}_{k\left| {k - 1} \right.}}$可直接表示为:
${\hat \phi _{k\left| {k - 1} \right.}} = {F}{\hat \phi _{k - 1}}, $
${{P}_{k\left| {k - 1} \right.}} = {F}{P_{k - 1}}{{F}^{\rm{T}}}. $
2) 测量更新
${{P}_{k\left| {k - 1} \right.}}$进行柯西分解${{P}_{k\left| {k - 1} \right.}} = {{S}_{k\left| {k - 1} \right.}}{S}_{k\left| {k - 1} \right.}^{\rm{T}}$, 得到容积点${\phi _{j,\left. k \right|k - 1}}$
${\phi _{j,\left. k \right|k - 1}} = {{S}_{k\left| {k - 1} \right.}}{\xi _j} + {\hat \phi _{k\left| {k - 1} \right.}}. $
3) 根据观测方程传播容积点
${{\eta }_{j,\left. k \right|k - 1}} = h\left( {{{\phi }_{j,\left. k \right|k - 1}}} \right) $
分别计算预测测量值${\hat {{\eta }}_{k\left| {k - 1} \right.}}$, 信息协方差矩阵${P}_{\eta \eta ,\left. k \right|k - 1}^{}$, 互协方差矩阵${{P}_{\phi \eta ,\left. k \right|k - 1}}$和卡尔曼增益${K_k}$,
${\hat {{\eta }}_{k\left| {k - 1} \right.}} = \frac{1}{{2{n_\phi }}}\sum\limits_{j = 1}^{2{n_\phi }} {{{\eta }_{j,\left. k \right|k - 1}}} , $
$\begin{split}{P}_{\eta \eta ,\left. k \right|k - 1} = & \frac{1}{{2{n_\phi }}}\sum\limits_{j = 1}^{2{n_\phi }} {{{\eta }_{j,\left. k \right|k - 1}}{\eta }_{j,\left. k \right|k - 1}^{\rm{T}}} \\ & - {\hat {{\eta }}_{k\left| {k - 1} \right.}}\hat {{\eta }}_{k\left| {k - 1} \right.}^{\rm{T}} + {{R}_k}, \end{split}$
${P}_{\phi \eta ,\left. k \right|k - 1}^{} = \frac{1}{{2{n_\phi }}}\!\sum\limits_{j = 1}^{2{n_\phi }}\!{{{\phi }_{j,\left. k \right|k - 1}}{\eta }_{j,\left. k \right|k - 1}^{\rm{T}}} - {\hat {{\eta }}_{k\left| {k - 1} \right.}}\hat {{\eta }}_{k\left| {k - 1} \right.}^{\rm{T}}, $
${K_k} = {{P}_{\phi \eta ,\left. k \right|k - 1}}{P}_{\eta \eta ,\left. k \right|k - 1}^{ - 1}. $
状态向量更新为
${\hat {{\phi }}_k} = {\hat {{\phi }}_{k\left| {k - 1} \right.}} + {{K}_k}\left( {{{\eta }_k} - {{\hat {{\eta }}}_{k\left| {k - 1} \right.}}} \right). $
协方差矩阵更新为
${{P}_k} = {{P}_{k\left| {k - 1} \right.}} - {{K}_k}{P}_{\eta \eta }^{}{K}_k^{\rm{T}}. $

2
3.2.正交单纯形容积卡尔曼滤波
-->考虑${n_\phi }$维正则单纯形, 可建立三阶球面单纯形容积准则[22], 其容积点可表示为
$\begin{array}{l}{\alpha }{\rm{ = }} - \left[ {{{\alpha }_1},{{\alpha }_2}, \cdots ,{{\alpha }_{{n_\phi } + 1}}, - {{\alpha }_1}, - {{\alpha }_2}, \cdots , - {{\alpha }_{{n_\phi } + 1}}} \right], \\ \end{array}$
式中${{\alpha }_j} \!=\! {\left(\! {{\alpha _{j,1}},{\alpha _{j,2}}, \!\ldots \!,{\alpha _{j,{n_\phi }}}} \!\right)^{\rm{T}}}$($j \!=\! 1,2, \cdots ,$$ {n_\phi } \!+\! 1$), 其中
$\begin{array}{l}{\alpha _{j,m}} = \left\{\!\!\begin{aligned}& - \sqrt {\frac{{{n_\phi } + 1}}{{{n_\phi }\left( {{n_\phi } - m + 2} \right)\left( {{n_\phi } - m + 1} \right)}}},\, m{\rm{ < }}j, \\& \sqrt {\frac{{\left( {{n_\phi } + 1} \right)\left( {{n_\phi } - j + 1} \right)}}{{{n_\phi }\left( {{n_\phi } - j + 2} \right)}}},\quad \quad \quad \quad \quad\, m=j, \\ \\ & 0,\quad \quad\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad m > j. \\ \end{aligned} \right. \end{array}$
为进一步提高滤波精度, 需要分析非线性对滤波的影响. 无迹卡尔曼滤波 (UKF)中, 尺度因子可以用来调节非线性对滤波的影响, 而基于容积准则的滤波方法没有相应的可调参数, 且容积点对称而无中心点, 这样均值估计与容积点之间有一定距离, 容易导致非局部采样效应[23], 进而影响滤波精度.
事实上, 假设${A}$为正交矩阵, 则将容积采样点进行正交变换后, 其仍然符合数值积分公式. 也就是说, 可以利用正交变换来调节高阶项对滤波的影响. 于是问题就转化为构建合适的正交矩阵, 使得变换后的单纯形容积点能够有效降低高阶干扰项对滤波的影响.
构建${n_\phi } \times {n_\phi }$正交矩阵${B} = \left[ {{{B}_1},{{B}_2}, \cdots ,{{B}_{{n_\phi }}}} \right]$, ${{B}_j} = {\left( {{\beta _{j,1}},{\beta _{j,2}}, \cdots ,{\beta _{j,{n_\phi }}}} \right)^{\rm T}}$, 其中
$\begin{gathered} {\beta _{j,2r - 1}} = \sqrt {2/{n_\phi }} \cos \left[ {\left( {2r - 1} \right){\rm{j {\text{π}} }}/{n_\phi }} \right], \\ {\beta _{j,2r}} = \sqrt {2/{n_\phi }} \sin \left[ {\left( {2r - 1} \right){\rm{j {\text{π}} }}/{n_\phi }} \right], \\ \end{gathered} $
式中$r = 1,2, \cdots ,\max \left\{ {p \in \mathbb{Z}|p \leqslant {n_\phi }/2} \right\}$, 如果${n_\phi }$为奇数, 则${\beta _{j,{n_\phi }}} = {\left( { - 1} \right)^j}/\sqrt {{n_\phi }} $.
将单纯形容积点列利用正交矩阵${B}$进行变换, 得到正交单纯形容积点
${\gamma } = {B} \times {\alpha }, $
式中${\gamma }=\left( {{{\gamma }_1},{{\gamma }_2}, \ldots ,{{\gamma }_{2{n_\phi } + 2}}} \right)$, ${{\alpha }_j} =\left({\alpha _{j,1}},{\alpha _{j,2}}, \cdots,\right.$$\left.{\alpha _{j,{n_\phi }}} \right)^{\rm{T}}$, $j = 1,2, \ldots ,2{n_\phi } + 2$.
这样, 就能利用正交变换来减小非线性影响程度, 从而提高滤波精度. 将(25)式中新的容积点及其权重代入(12)—(21)式中, 即可得到OSCKF算法.
根据图1和(12)—(21)式, 得到完整的基于量子相位测量的非线性锁相环结构, 如图2所示.
图 2 基于OSCKF算法的量子非线性数字锁相环
Figure2. Digital quantum nonlinear phase-lock loop based on OSCKF algorithm.

图2所示, 首先输入一个本振相位, 通过OSCKF算法迭代解算出相干态相位${\phi _k}$, 然后将本振相位设为${\phi _k}$, 再解算相干态相位${\phi _{k + 1}}$, 不断重复该过程, 本振相位将会逐渐逼近相干态相位, 从而实现锁相环的功能. 需要强调的是, 该方法并没有要求本振相位初始值接近相干态相位, 也能实现相干态相位的精确估计.
由于相干态经过零拍探测后可看作是经典信号加上高斯噪声[16], 且主要考察锁相环的有效性, 为此, 在实验验证过程中, 将采用直接模拟经过零拍后的经典信号, 取代构建真实物理实验平台的方法. 假设系统中相干态平均相位$\phi = {{\rm{\pi }}/3}$, 量测噪声服从均值为0, 方差${\sigma ^2} = {10^{ - 6}}$的高斯分布. 则当初始本振相位$\varphi = 0$时, 连续对1000个采样数据进行观察, 由(9)式可得零拍数据, 如图3所示, 图中测量值的抖动是由于噪声即真空波动引起的.
图 3 零拍数据示意图
Figure3. Homodyne data schematic diagram.

为了验证所提方法即使在不满足初始本振相位接近相干态相位的条件下, 也能精确估计相干态相位, 分别设置7个不同的初始本振相位进行实验, 每次实验仿真次数为200次, 取每次仿真第100次本振相位更新值作为最终相干态相位估计值, 计算200次仿真结果的平均误差为
$\varDelta {\rm{ = }}\frac{{\rm{1}}}{L}\sum\limits_{i = 1}^L {\left( {\left| {{{\hat \phi }_{i,100}} - \phi } \right|} \right)} , $
其中L为实验仿真次数, ${\hat \phi _{i,100}}$为每次实验第100次迭代后的估计值. 实验结果如图4所示.
图 4 不同初始本振相位下所提方法的性能
Figure4. Performance comparison of phase-locked loop in different local oscillator phase.

图4所示, 在其他条件相同的情况下, 不同初始本振相位下的相位估计结果相差并不大, 其中最高平均误差和最低平均误差只相差0.0016, 这说明所提方法的相位估计精度并没有受到初始本振相位设置的影响, 突破了传统自适应反馈量子相位测量需要保证本振相位和真实相位相差不大的限制.
为验证基于OSCKF算法的数字非线性锁相环的精确性, 利用EKF以及CKF算法实现锁相环相位估计, 并将结果进行对比. 定义平均均方根误差为
${\rm{RMS}} = \sqrt {\frac{1}{L}\sum\limits_{j = 1}^L {{{({{\hat \phi }_{k,j}} - {\phi _{k,j}})}^2}} } , $
其中${\hat \phi _{k,j}}$为第$j$次仿真$k$次迭代的相位估计值, L为Monte-Carlo次数.
为保证比较的合理性, 假定各算法的初始条件相同, 将Monte-Carlo实验仿真次数L和每次实验迭代次数$k$设为100次. 令${w_k} \propto N(0,{10^6})$, 初始协方差$P = 2$, 本振相位$\varphi = {\rm{0}}$. 各算法的收敛曲线如图5所示.
图 5 各算法误差收敛曲线
Figure5. Convergence curve of various algorithms.

图5可以看出, 随着迭代次数的增加, 所有算法的相位估计误差都趋于收敛. 其中 EKF将非线性模型线性化, 虽然处理方式简单, 但会引入较大的线性化误差, 因此收敛速度和最终精度均低于基于容积准则的方法, 而OSCKF算法由于引入了单纯形容积准则, 且利用正交变换进一步调节了非线性对滤波的影响, 因此比传统CKF算法精度更高. 此外, 从仿真实验发现, 采用EKF进行相位估计结果容易发散, 这说明EKF方法稳定性较差. 本文所提方法稳定性相对更好, 且能更快收敛, 与EKF方法相比, 平均误差降低了60%. 此外, 本文的状态模型是线性的, 而所设计的锁相环对于高非线性状态模型, 尤其是导航中的高非线性运动状态, 有望获得更好的效果. 因此, 设计非线性数字锁相环, 并采用正交单纯形容积准则, 不仅能够突破本振相位需要接近相干态相位的限制要求, 更能有效避免线性化误差, 提高相干态相位估计精度.
在利用量子零拍探测对相干态相位进行估计时, 对本振相位和相干态相位有严格限制要求. 如果本振相位不满足接近相干态相位的假设条件, 则非线性测量模型在进行线性化时会造成较大的误差. 针对该问题, 设计了一种数字非线性锁相环, 并提出OSCKF算法实现了锁相功能. OSCKF算法利用球面径向容积准则做数值计算近似高维积分, 可以有效避免模型线性化带来的误差, 同时减小了由于量测模型的高非线性带来的局部采样效应对精度的影响. 从仿真实验结果可以得出该锁相环算法利用不满足假设条件的本振相位也能得到相干态相位, 且数值稳定, 误差更小, 能够实现相干态相位的精确估计, 实用性更好. 下一步, 将在此锁相环精确估计的基础上, 实现对连续变化量子相位的精确跟踪, 并对噪声为非高斯的情况做进一步的研究.
相关话题/测量 设计 数字 实验 信息

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于二维声子晶体的大尺寸夹心式换能器的优化设计
    摘要:夹心式换能器应用极为广泛,但当其横向尺寸过大时,存在耦合振动,影响其辐射面的位移分布.本文通过在大尺寸夹心式换能器的前盖板中加工周期排列的槽,来形成一种二维声子晶体结构.随后,采用有限元法对基于二维声子晶体的大尺寸夹心式换能器的振动传输特性、共振频率以及发射电压响应进行仿真模拟,讨论了开槽高度 ...
    本站小编 Free考研考试 2021-12-29
  • 甲胺分子的紫外光解离动力学实验研究
    摘要:在280287.5nm区域内,通过实验测定共振增强多光子电离-时间飞行质谱、碎片离子的分质量激发谱以及光强指数等对甲胺分子的光解离通道进行了研究.实验结果证实甲胺分子在单光子能量范围内存在一个电子排斥态,主要的光解离过程为甲胺分子共振吸收1个光子到达该电子排斥态后解离成中性碎片,然后是中性碎片 ...
    本站小编 Free考研考试 2021-12-29
  • 基于量子增强型光纤马赫-曾德尔干涉仪的低频信号测量
    摘要:利用低频光通信波段真空压缩态光场可实现基于光纤的量子精密测量.本文利用简并光学参量振荡器实验制备出1550nm低频真空压缩态光场.在分析频段10500kHz范围内压缩态光场的压缩度均达3dB.用实验制备的1550nm真空压缩态光场填补光纤马赫-曾德尔干涉仪的真空通道,实现了量子增强型光纤马赫- ...
    本站小编 Free考研考试 2021-12-29
  • 部分状粒子处理方法及其对云微物理参数测量的影响
    摘要:作为云微物理过程测量的重要利器,机载云降水粒子成像仪在云降水物理与人工影响天气研究中具有重要的作用.从采样结果来看,机载云降水粒子成像仪所测粒子图像中含有大量的粒子图像仅是粒子的一部分而已,即部分状粒子.因其数量较多,对该类粒子所选处理方法不同,会引起测量结果的很大差异.本文介绍并分析了现有部 ...
    本站小编 Free考研考试 2021-12-29
  • <sup>252</sup>Cf自发裂变中子发射率符合测量的回归分析
    摘要:252Cf同位素源具有标准的自发裂变中子能谱,但由于其半衰期较短,应用中常需要对源强进行标定修正.随着源年龄增加,来自源中250Cf和248Cm自发裂变的影响愈加凸显,不能简单按252Cf的衰变规律计算源中子发射率,而通过锰浴活化的间接测量方法周期较长,且在源强低于104n/s时误差较大.最近 ...
    本站小编 Free考研考试 2021-12-29
  • 探测器对量子增强马赫-曾德尔干涉仪相位测量灵敏度的影响
    摘要:研究了强度差测量方案下,探测器量子效率对光子数态、关联数态、压缩真空态三种量子光源注入的马赫-曾德尔干涉仪相位测量灵敏度的影响.获得了相位测量灵敏度与效率的定量关系,比较了探测效率对不同量子态注入的干涉仪相位灵敏度的影响.研究表明:光子数态注入时,相位测量灵敏度始终不能超越标准量子极限;关联数 ...
    本站小编 Free考研考试 2021-12-29
  • 实验条件不完美对薛定谔猫态制备的影响
    摘要:薛定谔猫态是一类重要的非经典光场,实验上可以通过真空压缩态减光子的方案获得.本文从理论上研究了实验条件对制备薛定谔猫态的影响,主要考虑了包括压缩态的压缩度和纯度、单光子探测器的效率及噪声以及零拍探测器的效率等诸多因素的影响.理想情况下通过减光子方案制备得到的薛定谔猫态为奇光子数态,其相空间原点 ...
    本站小编 Free考研考试 2021-12-29
  • 铁电存储器中高能质子引发的单粒子功能中断效应实验研究
    摘要:利用中国原子能科学研究院的中高能质子实验平台,针对两款商用铁电存储器开展了中高能质子单粒子效应实验研究,发现其中一款器件在质子辐照下发生了单粒子翻转和单粒子功能中断.本文主要针对单粒子功能中断效应展开了后续实验研究.首先通过改变质子能量对器件进行辐照,发现单粒子功能中断截面随质子能量的提高而增 ...
    本站小编 Free考研考试 2021-12-29
  • 利用气泡探测器测量激光快中子
    摘要:在利用超强激光驱动中子源的研究和应用研究中,中子源的产额及其角分布至关重要.我们在星光Ⅲ号激光装置上采用气泡探测器对强激光驱动的中子源的产额及其角分布进行了测量.利用超强皮秒激光与碳氘薄膜靶相互作用产生高能氘离子束撞击次级碳氘靶,通过氘-氘核反应产生准单能快中子.实验发现中子束的发射具有一定的 ...
    本站小编 Free考研考试 2021-12-29
  • 基于拉曼激光雷达的大气三相态水同步精细探测分光系统的设计与仿真分析
    摘要:水是惟一具有三相态的大气参数,三相态水的分布研究对认识云微物理、云降水物理以及人工影响天气过程具有重要的科学意义.在大气三相态水的拉曼激光雷达探测技术中,需首先解决三相态水的高光谱分光技术,以保证对回波信号的精细提取和高信噪比探测.考虑到水汽、液态水和固态水的拉曼光谱特性,本文首先通过理论仿真 ...
    本站小编 Free考研考试 2021-12-29