[发明专利]一种雷达回波动态估算噪声的方法和气象雷达系统

权利要求书

1.一种雷达回波动态估算噪声的方法，其特征在于，包括：

将雷达回波信号转换为IQ时序信号，对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱sp_M，M为正整数；

按升序对功率谱sp_M进行排序，记排序后的功率谱为sp_S；

设变量mc和b0，初始化时令mc＝N，b0＝0；

判断mc>2，且b0<1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：

b0＝b2/b1，spS(n)表示sp_S中的第n个功率谱系数，n为不大于N的正整数；

每一次循环计算后，判断是否有b1≤0，并记录当前的mc；

如果b1≤0，则采用当前的mc估算系统噪声，得到第M个距离库的系统噪声为

2.根据权利要求1所述的方法，其特征在于，所述将雷达回波信号转换为IQ时序信号包括：

将雷达回波信号相移90度得到IQ时序信号。

3.根据权利要求1所述的方法，其特征在于，所述对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱sp_M包括：

对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换，得到信号

其中，k＝0,1,…,N-1；

根据公式sp_M(k)＝|y_M(k)|²计算功率谱系数，得到功率谱sp_M。

4.根据权利要求1、2或3所述的方法，其特征在于，所述判断mc>2，且b0<1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及b0，其中：

b0＝b2/b1包括：

判断mc>2，且b0<1时，令mc＝mc-1，计算b1和b2及新的b0，其中：

b0＝b2/b1；

重复上述步骤，直到mc不大于2或者b0不小于1。

5.一种气象雷达系统，其特征在于，包括：

信号转换模块，用于将雷达回波信号转换为IQ时序信号；

功率谱模块，用于对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱sp_M ，M为正整数；按升序对功率谱sp_M 进行排序，记排序后的功率谱为sp_S；

循环计算模块，用于设变量mc和b0，初始化时令mc＝N，b0＝0；判断mc>2，且b0<1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：

b0＝b2/b1，sp_S(n)表示sp_S中的第n个功率谱系数，n为不大于N的正整数；

噪声计算模块，用于每一次循环计算后，判断是否有b1≤0，并记录当前的mc；如果b1≤0，则采用当前的mc估算系统噪声，得到第M个距离库的系统噪声为

6.根据权利要求5所述的系统，其特征在于，

所述信号转换模块，具体用于将雷达回波信号相移90度得到IQ时序信号。

7.根据权利要求5所述的系统，其特征在于，

所述功率谱模块，具体用于对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换，得到信号其中，k＝0,1,…,N-1；根据公式sp_M(k)＝|y_M(k)| ²计算功率谱系数，得到功率谱sp_M 。

8.根据权利要求5、6或7所述的系统，其特征在于，

所述循环计算模块，具体用于判断mc>2，且b0<1时，令mc＝mc-1，计算b1和b2及新的b0，其中：b0＝b2/b1；重复上述步骤，直到mc不大于2或者b0不小于1。