[发明专利]一种CPFSK调制识别方法

权利要求书

1.一种CPFSK调制识别方法，其特征在于，包括如下步骤：

S1、对接收信号进行采样，得到r(n)＝s(n)+w(n) (n＝1,2,...,N)，其中，s(n)为发送信号，所述发送信号的调制类型为2CPFSK、4CPFSK和8CPFSK中的任意一种，w(n)为加性高斯白噪声，N为采样点的个数；

S2、对S1所述r(n)进行希尔伯特变换，得到则接收信号对应的复信号为其中，*表示卷积运算，j为虚部；

S3、使用阶数为N_f，归一化通带截止频率为ω_p＝1.2Rb/f_s，归一化阻带截止频率为ω_s＝2ω_p的滤波器h(n)对S1所述接收信号进行滤波，滤除带外噪声，得到新的信号其中，Rb为S1所述接收信号符号速率，f_s为采样频率，所述滤波器h(n)采用等纹波逼近法设计，N_f的值取2的次方；

S4、提取S3所述信号的相位信息，即相位其中，为S3所述新的信号实部，为S3所述新的信号虚部；

S5、以滑窗的方式对S4所述相位φ(n)进行N-D+1次线性拟合，得到瞬时频率值f(n)，其中，D为所述滑窗的窗长，D≤接收信号符号的周期；

S6、令f(n)∈[f_min,f_max]，其中，f_min为S5所述瞬时频率值f(n)的最小值，f_max为S5所述瞬时频率值f(n)的最大值；

S7、将S6所述[f_min,f_max]分为10个等间隔区间，统计每个区间的f(n)个数，以直方图的形式呈现f(n)在所述10个等间隔区间的分布情况，使用H＝[h(1),h(2),...,h(k),...,h(10)]表示10个等间隔区间内分布的f(n)点个数，其中，h(k)表示第k个等间隔区间的f(n)点个数，u＝1,2,3,...,10；

S8、利用判决式对发送信号的调制类型，

若满足则判定发送信号的调制类型为2CPFSK，

若满足则判定发送信号的调制类型为4CPFSK，

若满足则判定发送信号的调制类型为8CPFSK，其中，0.25≤TH1≤0.4，0.15≤TH2≤0.25。

2.根据权利要求1所述的一种CPFSK调制识别方法，其特征在于：S4所述提取S3所述信号的相位信息的方法具体如下：使用解缠绕处理得到连续变化的瞬时相位值，即：

当φ(n)-φ(n-1)＞π时，φ(n)＝φ(n)-2π，

当φ(n)-φ(n-1)＜-π时，φ(n)＝φ(n)+2π。

3.根据权利要求1所述的一种CPFSK调制识别方法，其特征在于：S5所述得到瞬时频率值f(n)的具体步骤如下：

S51、采用最小二乘法的方式得到拟合直线其中，n_i是采样时间序列，即n_i＝i，i＝1,2,3,...,N，为每次处理滑窗内n_i的平均值，是每次处理滑窗内φ(i)的平均值，

S52、将S51所述所有以每D个点为一组，提取出每组的绝对值最大值作为瞬时频率值，即其中，表示向下取整运算，

4.根据权利要求1所述的一种CPFSK调制识别方法，其特征在于：S8所述TH1＝0.3，TH2＝0.2。