CN102183753A - 一种混沌寻优雷达海杂波预报系统及方法 - Google Patents

Abstract

一种混沌寻优雷达海杂波预报系统，包括雷达、数据库以及上位机；雷达、数据库和上位机依次相连，所述雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库，所述的上位机包括数据预处理模块、预报模型建模模块、混沌寻优模块、海杂波预报模块、判别模型更新模块以及结果显示模块。以及提出了一种混沌寻优雷达海杂波预报方法。本发明提供一种能快速寻优且避免人为因素影响、高精度的雷达海杂波预报系统及方法。

Description

一种混沌寻优雷达海杂波预报系统及方法

技术领域

本发明涉及雷达数据处理领域，特别地，涉及一种混沌寻优雷达海杂波预报系统及方法。

背景技术

海杂波，即来自于被雷达发射信号照射的一片海面的后向散射回波。由于海杂波对来自海面或接近海面的“点”目标，如航海浮标及漂浮在海上的冰块等目标的雷达回波的可检测性形成严重制约，因此海杂波的研究对海洋背景中轮船等目标的探测性能具有十分重要的影响从而具有重要理论意义和实用价值。

习惯上海杂波被视为单一随机过程，如对数正态分布、K分布等。然而这些模型在实际应用中都有其特定的局限性，其中一个重要原因是海杂波看似随机的波形，实际上并不具有随机分布特性。

发明内容

为了克服传统雷达海杂波预报方法精度不高、易受人为因素影响的不足，本发明提供一种避免人为因素影响、高精度的混沌寻优雷达海杂波预报系统及方法。

本发明解决其技术问题所采用的技术方案是：

一种混沌寻优雷达海杂波预报系统，包括雷达、数据库以及上位机，雷达、数据库和上位机依次相连，所述雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库，所述的上位机包括：

数据预处理模块，用以进行雷达海杂波数据预处理，采用如下过程完成：

1)从数据库中采集N个雷达海杂波回波信号幅值x_i作为训练样本，i＝1，...，N；

2)对训练样本进行归一化处理，得到归一化幅值

${\stackrel{\‾}{x}}_i=\frac{x_i-\min x}{\max x-\min x}$

其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；

3)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：

$X=\left[\begin{array}{cccc}{\stackrel{\‾}{x}}_1& {\stackrel{\‾}{x}}_2& L& {\stackrel{\‾}{x}}_D\\ {\stackrel{\‾}{x}}_2& {\stackrel{\‾}{x}}_3& L& {\stackrel{\‾}{x}}_{D+1}\\ M& M& O& M\\ {\stackrel{\‾}{x}}_{N-D}& {\stackrel{\‾}{x}}_{N-D+1}& L& {\stackrel{\‾}{x}}_{N-1}\end{array}\right]$

$Y=\left[\begin{array}{c}{\stackrel{\‾}{x}}_{D+1}\\ {\stackrel{\‾}{x}}_{D+2}\\ M\\ {\stackrel{\‾}{x}}_N\end{array}\right]$

其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；

预报模型建模模块，用以建立预报模型，采用如下过程完成：

将得到的X、Y代入如下线性方程：

$\left[\begin{array}{cc}0& 1_v^T\\ 1_v& K+V_{\mathrm{\γ}}\end{array}\right]\left[\begin{array}{c}{b}^{*}\\ {\mathrm{\α}}^{*}\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]$

其中 $V_{\mathrm{\γ}}=\mathrm{diag}\{\frac{1}{\mathrm{\γ}{v}_1},K,\frac{1}{\mathrm{\γ}{v}_M}\}$

权重因子v_i由下式计算：

其中

是误差变量ξ_i标准差的估计，c₁，c₂为常量；

求解得待估计函数f(x)：

$f\left(x\right)=\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{{\mathrm{\α}}_i}^{*}\exp (-||x-x_i||/{\mathrm{\θ}}^2)+b^{*}$

其中，M是支持向量的数目，1_v＝[1，...，1]^T，

K＝exp(-||x_i-x_j||/θ²)，上标T表示矩阵的转置，是拉格朗日乘子，其中，i＝1，...，M，j＝1，…，M，b^*是偏置量，

和exp(-||x-x_i||/θ²)均为支持向量机的核函数，x_j为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；

混沌寻优模块，用以采用混沌遗传算法对预报模型的核参数θ和惩罚系数γ进行优化，采用如下过程完成：

5.1)采用实数编码方法对θ和γ进行编码；

5.2)随机产生初始种群；

5.3)计算每个个体的适应度，并判断是否符合算法终止条件，若符合，输出最佳个体及其代表的最优解，并结束计算，否则继续迭代；

5.4)采用正态分布概率选择个体；

5.5)通过单点线性交叉产生新个体；

5.6)对适应度最小的10％个体用Logistic映射进行混沌扰动；

5.7)产生的新一代种群，返回5.3)进行迭代；

其中，初始种群大小为50-200，最大代数50-300，最佳个体选择概率为0.05-0.1，交叉概率为0.5-0.9，个体适应度选择模型的泛化均方根误差，终止条件为达到最大代数或者连续五代适应度不变；

海杂波预报模块，用以进行海杂波预测，采用如下过程完成：

1)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[x_t-D+1，K，x_t]，x_t-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，x_t表示第t采样时刻的海杂波回波信号幅值；

2)进行归一化处理；

$\stackrel{\‾}{\mathrm{TX}}=\frac{\mathrm{TX}-\min x}{\max x-\min x}$

3)代入预报模型建模模块得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。

作为优选的一种方案：所述上位机还包括：判别模型更新模块，用以按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型；

作为优选的另一种方案：所述上位机还包括：结果显示模块，用以将海杂波预报模块计算得到的预报值在上位机显示。

一种混沌寻优雷达海杂波预报系统所使用的雷达海杂波预报方法，所述的方法包括以下步骤：

(1)从数据库中采集N个雷达海杂波回波信号幅值x_i作为训练样本，i＝1，...，N；

(2)对训练样本进行归一化处理，得到归一化幅值

${\stackrel{\‾}{x}}_i=\frac{x_i-\min x}{\max x-\min x}$

其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；

(3)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：

$X=\left[\begin{array}{cccc}{\stackrel{\‾}{x}}_1& {\stackrel{\‾}{x}}_2& L& {\stackrel{\‾}{x}}_D\\ {\stackrel{\‾}{x}}_2& {\stackrel{\‾}{x}}_3& L& {\stackrel{\‾}{x}}_{D+1}\\ M& M& O& M\\ {\stackrel{\‾}{x}}_{N-D}& {\stackrel{\‾}{x}}_{N-D+1}& L& {\stackrel{\‾}{x}}_{N-1}\end{array}\right]$

$Y=\left[\begin{array}{c}{\stackrel{\‾}{x}}_{D+1}\\ {\stackrel{\‾}{x}}_{D+2}\\ M\\ {\stackrel{\‾}{x}}_N\end{array}\right]$

其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；

(4)将得到的X、Y代入如下线性方程：

$\left[\begin{array}{cc}0& 1_v^T\\ 1_v& K+V_{\mathrm{\γ}}\end{array}\right]\left[\begin{array}{c}{b}^{*}\\ {\mathrm{\α}}^{*}\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]$

其中 $V_{\mathrm{\γ}}=\mathrm{diag}\{\frac{1}{\mathrm{\γ}{v}_1},K,\frac{1}{\mathrm{\γ}{v}_M}\}$

权重因子v_i由下式计算：

其中是误差变量ξ_i标准差的估计，c₁，c₂为常量；

求解得待估计函数f(x)：

$f\left(x\right)=\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{{\mathrm{\α}}_i}^{*}\exp (-||x-x_i||/{\mathrm{\θ}}^2)+b^{*}$

其中，M是支持向量的数目，1_v＝[1，...，1]^T，K＝exp(-||x_i-x_j||/θ²)，上标T表示矩阵的转置，

是拉格朗日乘子，其中，i＝1，...，M，j＝1，…，M，b^*是偏置量，和exp(-||x-x_i||/θ²)均为支持向量机的核函数，x_j为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；

(5)用混沌遗传算法对步骤(4)的核参数θ和惩罚系数γ进行优化，

5.1)采用实数编码方法对θ和γ进行编码；

5.2)随机产生初始种群；

5.3)计算每个个体的适应度，并判断是否符合算法终止条件，若符合，输出最佳个体及其代表的最优解，并结束计算，否则继续迭代；

5.4)采用正态分布概率选择个体；

5.5)通过单点线性交叉产生新个体；

5.6)对适应度最小的10％个体用Logistic映射进行混沌扰动；

5.7)产生的新一代种群，返回5.3)进行迭代；

其中，初始种群大小为50-200，最大代数50-300，最佳个体选择概率为0.05-0.1，交叉概率为0.5-0.9，个体适应度选择模型的泛化均方根误差，终止条件为达到最大代数或者连续五代适应度不变；

(6)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[x_t-D+1，K，x_t]，x_t-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，x_t表示第t采样时刻的海杂波回波信号幅值；

(7)进行归一化处理；

$\stackrel{\‾}{\mathrm{TX}}=\frac{\mathrm{TX}-\min x}{\max x-\min x}$

(8)代入步骤(4)得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。

作为优选的一种方案：所述的方法还包括：

(9)、按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。

作为优选的另一种方案：在所述的步骤(8)中，将计算得到的海杂波预报值在上位机显示。

本发明的技术构思为：本发明针对雷达海杂波的混沌特性，对雷达海杂波数据进行重构，并对重构后的数据进行非线性拟合，引入混沌寻优方法，从而建立雷达海杂波的混沌最优预报模型。

本发明的有益效果主要表现在：1、建立了雷达海杂波预报模型，可以在线预测雷达海杂波；2、所用的建模方法只需较少样本即可、响应速度快；3、精度高、避免了人为因素的影响。

附图说明

图1是本发明所提出的系统的硬件结构图；

图2是本发明所提出的上位机的功能模块图。

具体实施方式

下面结合附图对本发明作进一步描述。本发明实施例用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范围。

实施例1

参照图1、图2，一种混沌寻优雷达海杂波预报系统，包括雷达1连接的数据库2、及上位机3，雷达1、数据库2和上位机3依次相连，所述雷达1对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库2，所述的上位机3包括：

数据预处理模块4，用以进行雷达海杂波数据预处理，采用如下过程完成：

1)从数据库中采集N个雷达海杂波回波信号幅值x_i作为训练样本，i＝1，...，N；

2)对训练样本进行归一化处理，得到归一化幅值

${\stackrel{\‾}{x}}_i=\frac{x_i-\min x}{\max x-\min x}$

其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；

3)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：

$X=\left[\begin{array}{cccc}{\stackrel{\‾}{x}}_1& {\stackrel{\‾}{x}}_2& L& {\stackrel{\‾}{x}}_D\\ {\stackrel{\‾}{x}}_2& {\stackrel{\‾}{x}}_3& L& {\stackrel{\‾}{x}}_{D+1}\\ M& M& O& M\\ {\stackrel{\‾}{x}}_{N-D}& {\stackrel{\‾}{x}}_{N-D+1}& L& {\stackrel{\‾}{x}}_{N-1}\end{array}\right]$

$Y=\left[\begin{array}{c}{\stackrel{\‾}{x}}_{D+1}\\ {\stackrel{\‾}{x}}_{D+2}\\ M\\ {\stackrel{\‾}{x}}_N\end{array}\right]$

其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；

预报模型建模模块5，用以建立预报模型，采用如下过程完成：

将得到的X、Y代入如下线性方程：

$\left[\begin{array}{cc}0& 1_v^T\\ 1_v& K+V_{\mathrm{\γ}}\end{array}\right]\left[\begin{array}{c}{b}^{*}\\ {\mathrm{\α}}^{*}\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]$

其中 $V_{\mathrm{\γ}}=\mathrm{diag}\{\frac{1}{\mathrm{\γ}{v}_1},K,\frac{1}{\mathrm{\γ}{v}_M}\}$

权重因子v_i由下式计算：

其中

是误差变量ξ_i标准差的估计，c₁，c₂为常量；

求解得待估计函数f(x)：

$f\left(x\right)=\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{{\mathrm{\α}}_i}^{*}\exp (-||x-x_i||/{\mathrm{\θ}}^2)+b^{*}$

其中，M是支持向量的数目，1_v＝[1，...，1]^T，

K＝exp(-||x_i-x_j||/θ²)，上标T表示矩阵的转置，

是拉格朗日乘子，其中，i＝1，...，M，j＝1，…，M，b^*是偏置量，

和exp(-||x-x_i||/θ²)均为支持向量机的核函数，x_j为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；

混沌寻优模块6，用以采用混沌遗传算法对预报模型的核参数θ和惩罚系数γ进行优化，采用如下过程完成：

5.1)采用实数编码方法对θ和γ进行编码；

5.2)随机产生初始种群；

5.3)计算每个个体的适应度，并判断是否符合算法终止条件，若符合，输出最佳个体及其代表的最优解，并结束计算，否则继续迭代；

5.4)采用正态分布概率选择个体；

5.5)通过单点线性交叉产生新个体；

5.6)对适应度最小的10％个体用Logistic映射进行混沌扰动；

5.7)产生的新一代种群，返回5.3)进行迭代；

其中，初始种群大小为50-200，最大代数50-300，最佳个体选择概率为0.05-0.1，交叉概率为0.5-0.9，个体适应度选择模型的泛化均方根误差，终止条件为达到最大代数或者连续五代适应度不变；

海杂波预报模块7，用以进行海杂波预测，采用如下过程完成：

1)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[x_t-D+1，K，x_t]，x_t-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，x_t表示第t采样时刻的海杂波回波信号幅值；

2)进行归一化处理；

$\stackrel{\‾}{\mathrm{TX}}=\frac{\mathrm{TX}-\min x}{\max x-\min x}$

3)代入预报模型建模模块得到的待估计函数f(x)得到采样时刻(t+1)的海杂波预报值

所述的上位机3还包括：判别模型更新模块8，按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。

所述上位机3还包括：结果显示模块9，用于将海杂波预报模块计算得到的预报值在上位机显示。

所述上位机3的硬件部分包括：I/O元件，用于数据的采集和信息的传递；数据存储器，存储运行所需的数据样本和运行参数等；程序存储器，存储实现功能模块的软件程序；运算器，执行程序，实现指定的功能；显示模块，显示设置的参数和运行结果。

实施例2

参照图1、图2，一种混沌寻优雷达海杂波预报方法，所述的方法包括以下步骤：

(1)从数据库中采集N个雷达海杂波回波信号幅值x_i作为训练样本，i＝1，...，N；

(2)对训练样本进行归一化处理，得到归一化幅值

${\stackrel{\‾}{x}}_i=\frac{x_i-\min x}{\max x-\min x}$

其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；

(3)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：

$X=\left[\begin{array}{cccc}{\stackrel{\‾}{x}}_1& {\stackrel{\‾}{x}}_2& L& {\stackrel{\‾}{x}}_D\\ {\stackrel{\‾}{x}}_2& {\stackrel{\‾}{x}}_3& L& {\stackrel{\‾}{x}}_{D+1}\\ M& M& O& M\\ {\stackrel{\‾}{x}}_{N-D}& {\stackrel{\‾}{x}}_{N-D+1}& L& {\stackrel{\‾}{x}}_{N-1}\end{array}\right]$

$Y=\left[\begin{array}{c}{\stackrel{\‾}{x}}_{D+1}\\ {\stackrel{\‾}{x}}_{D+2}\\ M\\ {\stackrel{\‾}{x}}_N\end{array}\right]$

其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；

(4)将得到的X、Y代入如下线性方程：

$\left[\begin{array}{cc}0& 1_v^T\\ 1_v& K+V_{\mathrm{\γ}}\end{array}\right]\left[\begin{array}{c}{b}^{*}\\ {\mathrm{\α}}^{*}\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]$

其中 $V_{\mathrm{\γ}}=\mathrm{diag}\{\frac{1}{\mathrm{\γ}{v}_1},K,\frac{1}{\mathrm{\γ}{v}_M}\}$

权重因子v_i由下式计算：

其中

是误差变量ξ_i标准差的估计，c₁，c₂为常量；

求解得待估计函数f(x)：

$f\left(x\right)=\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{{\mathrm{\α}}_i}^{*}\exp (-||x-x_i||/{\mathrm{\θ}}^2)+b^{*}$

其中，M是支持向量的数目，1_v＝[1，...，1]^T，

K＝exp(-||x_i-x_j||/θ²)，上标T表示矩阵的转置，

是拉格朗日乘子，其中，i＝1，...，M，j＝1，…，M，b^*是偏置量，和exp(-||x-x_i||/θ²)均为支持向量机的核函数，x_j为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；

(5)用混沌遗传算法对步骤(4)的核参数θ和惩罚系数γ进行优化，

5.1)采用实数编码方法对θ和γ进行编码；

5.2)随机产生初始种群；

5.3)计算每个个体的适应度，并判断是否符合算法终止条件，若符合，输出最佳个体及其代表的最优解，并结束计算，否则继续迭代；

5.4)采用正态分布概率选择个体；

5.5)通过单点线性交叉产生新个体；

5.6)对适应度最小的10％个体用Logistic映射进行混沌扰动；

5.7)产生的新一代种群，返回5.3)进行迭代；

其中，初始种群大小为50-200，最大代数50-300，最佳个体选择概率为0.05-0.1，交叉概率为0.5-0.9，个体适应度选择模型的泛化均方根误差，终止条件为达到最大代数或者连续五代适应度不变；

(6)在采样时刻t采集D个海杂波回波信号幅值得到，x_t-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，x_t表示第t采样时刻的海杂波回波信号幅值；

(7)进行归一化处理；

$\stackrel{\‾}{\mathrm{TX}}=\frac{\mathrm{TX}-\min x}{\max x-\min x}$

(8)代入步骤(4)得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。

所述的方法还包括：(9)、按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。

所述的方法还包括：在所述的步骤(8)中将计算得到的海杂波预报值在上位机显示。

Claims (6)

1.一种混沌寻优雷达海杂波预报系统，包括雷达、数据库以及上位机，雷达、数据库和上位机依次相连，其特征在于：所述雷达对所检测海域进行照射，并将雷达海杂波数据储存到所述的数据库，所述的上位机包括：

数据预处理模块，用以进行雷达海杂波数据预处理，采用如下过程完成：

1)从数据库中采集N个雷达海杂波回波信号幅值x_i作为训练样本，i＝1，...，N；

2)对训练样本进行归一化处理，得到归一化幅值

${\stackrel{\‾}{x}}_i=\frac{x_i-\min x}{\max x-\min x}$

其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；

3)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：

$X=\left[\begin{array}{cccc}{\stackrel{\‾}{x}}_1& {\stackrel{\‾}{x}}_2& L& {\stackrel{\‾}{x}}_D\\ {\stackrel{\‾}{x}}_2& {\stackrel{\‾}{x}}_3& L& {\stackrel{\‾}{x}}_{D+1}\\ M& M& O& M\\ {\stackrel{\‾}{x}}_{N-D}& {\stackrel{\‾}{x}}_{N-D+1}& L& {\stackrel{\‾}{x}}_{N-1}\end{array}\right]$

$Y=\left[\begin{array}{c}{\stackrel{\‾}{x}}_{D+1}\\ {\stackrel{\‾}{x}}_{D+2}\\ M\\ {\stackrel{\‾}{x}}_N\end{array}\right]$

其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；

预报模型建模模块，用以建立预报模型，采用如下过程完成：

将得到的X、Y代入如下线性方程：

$\left[\begin{array}{cc}0& 1_v^T\\ 1_v& K+V_{\mathrm{\γ}}\end{array}\right]\left[\begin{array}{c}{b}^{*}\\ {\mathrm{\α}}^{*}\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]$

其中 $V_{\mathrm{\γ}}=\mathrm{diag}\{\frac{1}{\mathrm{\γ}{v}_1},K,\frac{1}{\mathrm{\γ}{v}_M}\}$

权重因子v_i由下式计算：

其中

是误差变量ξ_i标准差的估计，c₁，c₂为常量；

求解得待估计函数f(x)：

$f\left(x\right)=\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{{\mathrm{\α}}_i}^{*}\exp (-||x-x_i||/{\mathrm{\θ}}^2)+b^{*}$

其中，M是支持向量的数目，1_v＝[1，...，1]^T，

K＝exp(-||x_i-x_j||/θ²)，上标T表示矩阵的转置，

是拉格朗日乘子，其中，i＝1，...，M，j＝1，…，M，b^*是偏置量，

和exp(-||x-x_i||/θ²)均为支持向量机的核函数，x_j为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；

混沌寻优模块，用以采用混沌遗传算法对预报模型的核参数θ和惩罚系数γ进行优化，采用如下过程完成：

5.1)采用实数编码方法对θ和γ进行编码；

5.2)随机产生初始种群；

5.3)计算每个个体的适应度，并判断是否符合算法终止条件，若符合，输出最佳个体及其代表的最优解，并结束计算，否则继续迭代；

5.4)采用正态分布概率选择个体；

5.5)通过单点线性交叉产生新个体；

5.6)对适应度最小的10％个体用Logistic映射进行混沌扰动；

5.7)产生的新一代种群，返回5.3)进行迭代；

其中，初始种群大小为50-200，最大代数50-300，最佳个体选择概率为0.05-0.1，交叉概率为0.5-0.9，个体适应度选择模型的泛化均方根误差，终止条件为达到最大代数或者连续五代适应度不变；

海杂波预报模块，用以进行海杂波预测，采用如下过程完成：

1)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[x_t-D+1，K，x_t]，x_t-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，x_t表示第t采样时刻的海杂波回波信号幅值；

2)进行归一化处理；

$\stackrel{\‾}{\mathrm{TX}}=\frac{\mathrm{TX}-\min x}{\max x-\min x}$

3)代入预报模型建模模块得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。

2.如权利要求1所述的混沌寻优雷达海杂波预报系统，其特征在于：所述上位机还包括：判别模型更新模块，用以按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型；

3.如权利要求1或2所述的混沌寻优雷达海杂波预报系统，其特征在于：所述上位机还包括：结果显示模块，用以将海杂波预报模块计算得到的预报值在上位机显示。

4.一种如权利要求1所述的混沌寻优雷达海杂波预报系统所使用的雷达海杂波预报方法，其特征在于：所述的方法包括以下步骤：

(1)从数据库中采集N个雷达海杂波回波信号幅值x_i作为训练样本，i＝1，...，N；

(2)对训练样本进行归一化处理，得到归一化幅值

${\stackrel{\‾}{x}}_i=\frac{x_i-\min x}{\max x-\min x}$

其中，minx表示训练样本中的最小值，maxx表示训练样本中的最大值；

(3)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y：

$X=\left[\begin{array}{cccc}{\stackrel{\‾}{x}}_1& {\stackrel{\‾}{x}}_2& L& {\stackrel{\‾}{x}}_D\\ {\stackrel{\‾}{x}}_2& {\stackrel{\‾}{x}}_3& L& {\stackrel{\‾}{x}}_{D+1}\\ M& M& O& M\\ {\stackrel{\‾}{x}}_{N-D}& {\stackrel{\‾}{x}}_{N-D+1}& L& {\stackrel{\‾}{x}}_{N-1}\end{array}\right]$

$Y=\left[\begin{array}{c}{\stackrel{\‾}{x}}_{D+1}\\ {\stackrel{\‾}{x}}_{D+2}\\ M\\ {\stackrel{\‾}{x}}_N\end{array}\right]$

其中，D表示重构维数，D为自然数，且D＜N，D的取值范围为50-70；

(4)将得到的X、Y代入如下线性方程：

$\left[\begin{array}{cc}0& 1_v^T\\ 1_v& K+V_{\mathrm{\γ}}\end{array}\right]\left[\begin{array}{c}{b}^{*}\\ {\mathrm{\α}}^{*}\end{array}\right]=\left[\begin{array}{c}0\\ Y\end{array}\right]$

其中 $V_{\mathrm{\γ}}=\mathrm{diag}\{\frac{1}{\mathrm{\γ}{v}_1},K,\frac{1}{\mathrm{\γ}{v}_M}\}$

权重因子v_i由下式计算：

其中

是误差变量ξ_i标准差的估计，c₁，c₂为常量；

求解得待估计函数f(x)：

$f\left(x\right)=\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{{\mathrm{\α}}_i}^{*}\exp (-||x-x_i||/{\mathrm{\θ}}^2)+b^{*}$

其中，M是支持向量的数目，1_v＝[1，...，1]^T，

K＝exp(-||x_i-x_j||/θ²)，上标T表示矩阵的转置，

是拉格朗日乘子，其中，i＝1，...，M，j＝1，…，M，b^*是偏置量，

和exp(-||x-x_i||/θ²)均为支持向量机的核函数，x_j为第j个雷达海杂波回波信号幅值，θ是核参数，x表示输入变量，γ是惩罚系数；

(5)用混沌遗传算法对步骤(4)的核参数θ和惩罚系数γ进行优化，

5.1)采用实数编码方法对θ和γ进行编码；

5.2)随机产生初始种群；

5.3)计算每个个体的适应度，并判断是否符合算法终止条件，若符合，输出最佳个体及其代表的最优解，并结束计算，否则继续迭代；

5.4)采用正态分布概率选择个体；

5.5)通过单点线性交叉产生新个体；

5.6)对适应度最小的10％个体用Logistic映射进行混沌扰动；

5.7)产生的新一代种群，返回5.3)进行迭代；

其中，初始种群大小为50-200，最大代数50-300，最佳个体选择概率为0.05-0.1，交叉概率为0.5-0.9，个体适应度选择模型的泛化均方根误差，终止条件为达到最大代数或者连续五代适应度不变；

(6)在采样时刻t采集D个海杂波回波信号幅值得到TX＝[x_t-D+1，K，x_t]，x_t-D+1表示第t-D+1采样时刻的海杂波回波信号幅值，x_t表示第t采样时刻的海杂波回波信号幅值；

(7)进行归一化处理；

$\stackrel{\‾}{\mathrm{TX}}=\frac{\mathrm{TX}-\min x}{\max x-\min x}$

(8)代入步骤(4)得到的待估计函数f(x)计算得到采样时刻(t+1)的海杂波预报值。

5.如权利要求4所述的雷达海杂波预报方法，其特征在于：所述的方法还包括：

(9)、按设定的采样时间间隔采集数据，将得到的实测数据与模型预报值比较，如果相对误差大于10％，则将新数据加入训练样本数据，更新预报模型。

6.如权利要求4或5所述的雷达海杂波预报方法，其特征在于：在所述的步骤(8)中，将计算得到的海杂波预报值在上位机显示。

Images