[发明专利]一种基于蒙特卡罗的窄带辐射源频差估计方法及装置

权利要求书

1.一种基于蒙特卡罗的窄带辐射源频差估计方法，其特征在于，包括如下步骤：

1)通过两路采样通道获取两路窄带辐射源信号；两路窄带辐射源信号的模型分别为：

其中，s(t)为信号的复包络，w₁(t)和w₂(t)为零均值的高斯白噪声且方差为和a₁和a₂为两路窄带辐射源信号的传播系数，f_d1和f_d2为两路窄带辐射源信号的多普勒频移，和为两路窄带辐射源信号的初始相位，j为虚数单位，t为时间，f_c为载频；

2)根据两路窄带辐射源信号获取两路窄带辐射源信号频差的似然函数，并根据所述似然函数构造频差的伪概率密度函数；所述频差的似然函数表示为：

其中，y₂＝[y₂(0),y₂(1),...,y₂(N-1)]^T

其中，f_d为两路窄带辐射源信号的多普勒频移的差值，N为采样点数，T_S为采样周期，j为虚数单位，y₁和y₂为两路窄带辐射源信号，n＝0,1,...,N-1；

3)利用蒙特卡罗方法生成服从所述伪概率密度函数的频差参数样本，根据生成的频差参数样本估算频差；

利用蒙特卡罗的接受-拒绝抽样方法生成服从所述伪概率密度函数的频差参数样本，频差参数样本的生成过程包括以下步骤：

(A)利用提议分布函数q(f_d)生成候选样本其中f_dmin,f_dmax分别表示频差参数取值的上下界；

(B)计算候选样本的接受概率即

(C)生成随机数u～U[0,1]；

(D)如果接受概率则接受候选样本作为的样本，否则拒绝该样本；

(E)重复步骤(A)～(D)，直到生成设定数量的频差参数样本。

2.根据权利要求1所述的基于蒙特卡罗的窄带辐射源频差估计方法，其特征在于，所述频差表示为：

其中，为频差，表示从伪概率密度函数中抽取的样本。