一种弹道目标RCS时间序列的进动周期提取方法

一种弹道目标rcs时间序列的进动周期提取方法
技术领域
1.本发明涉及计算机技术领域，尤其是一种弹道目标rcs时间序列的进动周期提取方法。


背景技术：

2.空间目标识别是指包括飞机，导弹及临近空间飞行器等目标在内的跟踪与识别技术，该技术以雷达为主要传感器，利用雷达获取空间目标的回波信号，从中提取目标的位置、速度、结构等特征信息，进而实现对空间目标的类型或属性的识别。随着现代雷达技术的进步，雷达对目标的精细刻画能力越来越高，可用于识别的主要特征包括但不限于目标等效散射面积(rcs)，微动，极化，一维像二维像等。微动作为目标的独特特征，相比于平动包含了更多的细节信息，对于空间目标探测与识别具有重要意义。
3.研究表明，在较长的时间段内，rcs序列具有典型的非平稳性、相关性和拟周期性的特点，因此常规的谱分析周期估计方法并不能得到满意的效果，且精度通常与观测数据长度成正比，这对于实时性要求很高的目标识别防御系统来说往往相互矛盾，同时，当进动角较大时，rcs序列在姿态变化范围内起伏程度增大，这对姿态角近似余弦函数的形状破坏十分严重，因此使得虚假周期分量增大。基于rcs周期提取的另一个思路是基于序列的估计，如目前所提出的二重循环自相关法，正余弦函数线性组合拟合法，方差分析法等。然而上述方法均只利用了单通道极化信息进行处理，存在特征提取精度不够高，极化信息利用不充分等问题，同时由于各极化通道处理效果不同，通常难以确定某一具体目标极化通道。


技术实现要素：

4.有鉴于此，本发明实施例提供一种精度高且稳定性高的，弹道目标rcs时间序列的进动周期提取方法。
5.本发明的一方面提供了一种弹道目标rcs时间序列的进动周期提取方法，包括：
6.计算目标的运动引起的电波入射角变化信息后，通过对所述目标的暗室静态测量数据进行插值，得到所述目标的动态数据序列；
7.对所述动态数据序列进行平滑处理后，采用小波变换阈值去噪法进行去噪处理，得到初始数据；
8.对所述初始数据进行全极化通道融合处理，得到进动周期信息。
9.可选地，所述计算目标的运动引起的电波入射角变化信息后，通过对所述目标的暗室静态测量数据进行插值，得到所述目标的动态数据序列，包括：
10.建立质心运动模型，根据预先设定的导弹参数模拟计算出质心的位置矢量和速度矢量信息；
11.建立姿态运动模型，模拟计算出目标的姿态角的变化信息；
12.根据极化基变换关系，将所述姿态角和目标静态散射测量数据相结合，插值得到特定姿态下目标的散射特性，获取目标的动态数据序列。
13.可选地，所述对所述动态数据序列进行平滑处理后，采用小波变换阈值去噪法进行去噪处理，得到初始数据，具体为：
14.通过savitzky-golay滤波器对所述动态数据序列进行平滑处理；
15.采用小波变换阈值去噪法进行去噪处理，得到初始数据。
16.可选地，所述通过savitzky-golay滤波器对所述动态数据序列进行平滑处理，包括：
17.通过多阶多项式来拟合一个窗口内的一组数据，确定拟合数据点与原数据点之间的残差平方和；其中，当所述残差平方和对所述多阶多项式的系数求偏导结果为零时，残差平方和最小，拟合效果最好；
18.根据需要拟合的单边点数、多项式阶次以及待拟合数据，计算得到所述多阶多项式；
19.根据所述多阶多项式求取窗口内的中心点估计值，完成对所述动态数据序列的平滑处理。
20.可选地，所述采用小波变换阈值去噪法进行去噪处理，得到初始数据，包括：
21.对去噪处理后的信号进行小波分解处理，得到各尺度下的近似系数和细节系数；
22.根据预设的处理规则对所述近似系数和所述细节系数进行阈值处理，得到小波系数；
23.对所述小波系数进行信号重构，得到去噪后的初始数据。
24.可选地，所述阈值处理中的阈值采用固定式阈值，所述阈值处理根据每一层小波分解的噪声水平估计进行调整。
25.可选地，所述对所述初始数据进行全极化通道融合处理，得到进动周期信息，包括：
26.假设以平面电磁波照射一个有限尺寸目标，获取雷达发射天线辐射波在目标处的电场矢量；
27.用平面波近似散射波，并根据目标散射的远场区中的给定观测点以及目标与接受天线之间的距离，确定散射电场；
28.定义一个复的极化散射矩阵来描述入射波和远区散射波各极化分量间的变换关系，构建极化散射矩阵；
29.定义目标的功率矩阵；
30.根据所述功率矩阵定义功率散射矩阵以及功率散射矩阵的迹；
31.根据所述功率散射矩阵以及所述功率散射矩阵的迹，定义去极化系数；
32.确定所述去极化系数的时间序列图以及功率谱密度结果；
33.根据所述时间序列图以及所述功率谱密度结果，得到进动周期信息。
34.本发明实施例的另一方面还提供了一种弹道目标rcs时间序列的进动周期提取装置，包括：
35.第一模块，用于计算目标的运动引起的电波入射角变化信息后，通过对所述目标的暗室静态测量数据进行插值，得到所述目标的动态数据序列；
36.第二模块，用于对所述动态数据序列进行平滑处理后，采用小波变换阈值去噪法进行去噪处理，得到初始数据；
37.第三模块，用于对所述初始数据进行全极化通道融合处理，得到进动周期信息。
38.本发明实施例的另一方面还提供了一种电子设备，包括处理器以及存储器；
39.所述存储器用于存储程序；
40.所述处理器执行所述程序实现如前面所述的方法。
41.本发明实施例的另一方面还提供了一种计算机可读存储介质，所述存储介质存储有程序，所述程序被处理器执行实现如前面所述的方法。
42.本发明实施例还公开了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行前面的方法。
43.本发明的实施例计算目标的运动引起的电波入射角变化信息后，通过对所述目标的暗室静态测量数据进行插值，得到所述目标的动态数据序列；对所述动态数据序列进行平滑处理后，采用小波变换阈值去噪法进行去噪处理，得到初始数据；对所述初始数据进行全极化通道融合处理，得到进动周期信息。本发明提高了精度和稳定性。
附图说明
44.为了更清楚地说明本技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
45.图1为本发明的整体步骤流程图；
46.图2为本发明实施例提供的仿真数据获取过程的步骤流程图；
47.图3为本发明实施例提供的模拟弹头类目标的形状示意图；
48.图4为仿真目标在中段的各通道rcs时间序列示意图；
49.图5为趋势项提取结果及剩余项的示意图；
50.图6为小波分解和重构过程的示意图；
51.图7为小波阈值去噪结果的示意图；
52.图8为预处理前后的示意图；
53.图9为去极化系数d时频域图；
54.图10为圆锥体的提取误差对比示意图；
55.图11为锥球体的提取误差对比示意图。
具体实施方式
56.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
57.术语说明：
58.rcs：雷达截面积，(radar cross section,rcs)是目标在雷达接收方向上反射雷达信号能力的度量，一个目标的rcs等于单位立体角目标在雷达接收天线方向上反射的功
率(每单独立体角)与人射到目标处的功率密度(每平方米)之比。
59.进动：在弹道导弹攻防系统中，导弹离开大气层后，会分离动力装置，并完成姿态调整，从而进入中段飞行。为了保持弹头在大气层外飞行的稳定性和提高命中精度，弹头在中段一般会进行姿态控制。姿态控制的方式有多种，其中自旋稳定是弹头在中段最常用的姿态控制方式。在弹头与母舱的分离过程中，弹头会受到冲击力矩的作用，力矩消失后对称轴将在平衡位置做椭圆锥运动，即进动。
60.极化：对各种矢量波来说，极化是指用一个场矢量来描述空间某一固定点所观测到的矢量波(电场，应变，自旋)随时间变化的特性，对电磁波而言，极化描述了电场矢量端点作为时间的函数所形成的空间轨迹的形状和旋向。根据电场强度矢量端点运动轨迹的不同，电磁波可以被分为完全极化波，部分极化波和完全非极化波。
61.针对现有技术存在的问题，本发明实施例提供了一种弹道目标rcs时间序列的进动周期提取方法，包括：
62.计算目标的运动引起的电波入射角变化信息后，通过对所述目标的暗室静态测量数据进行插值，得到所述目标的动态数据序列；
63.对所述动态数据序列进行平滑处理后，采用小波变换阈值去噪法进行去噪处理，得到初始数据；
64.对所述初始数据进行全极化通道融合处理，得到进动周期信息。
65.可选地，所述计算目标的运动引起的电波入射角变化信息后，通过对所述目标的暗室静态测量数据进行插值，得到所述目标的动态数据序列，包括：
66.建立质心运动模型，根据预先设定的导弹参数模拟计算出质心的位置矢量和速度矢量信息；
67.建立姿态运动模型，模拟计算出目标的姿态角的变化信息；
68.根据极化基变换关系，将所述姿态角和目标静态散射测量数据相结合，插值得到特定姿态下目标的散射特性，获取目标的动态数据序列。
69.可选地，所述对所述动态数据序列进行平滑处理后，采用小波变换阈值去噪法进行去噪处理，得到初始数据，具体为：
70.通过savitzky-golay滤波器对所述动态数据序列进行平滑处理；
71.采用小波变换阈值去噪法进行去噪处理，得到初始数据。
72.可选地，所述通过savitzky-golay滤波器对所述动态数据序列进行平滑处理，包括：
73.通过多阶多项式来拟合一个窗口内的一组数据，确定拟合数据点与原数据点之间的残差平方和；其中，当所述残差平方和对所述多阶多项式的系数求偏导结果为零时，残差平方和最小，拟合效果最好；
74.根据需要拟合的单边点数、多项式阶次以及待拟合数据，计算得到所述多阶多项式；
75.根据所述多阶多项式求取窗口内的中心点估计值，完成对所述动态数据序列的平滑处理。
76.可选地，所述采用小波变换阈值去噪法进行去噪处理，得到初始数据，包括：
77.对去噪处理后的信号进行小波分解处理，得到各尺度下的近似系数和细节系数；
78.根据预设的处理规则对所述近似系数和所述细节系数进行阈值处理，得到小波系数；
79.对所述小波系数进行信号重构，得到去噪后的初始数据。
80.可选地，所述阈值处理中的阈值采用固定式阈值，所述阈值处理根据每一层小波分解的噪声水平估计进行调整。
81.可选地，所述对所述初始数据进行全极化通道融合处理，得到进动周期信息，包括：
82.假设以平面电磁波照射一个有限尺寸目标，获取雷达发射天线辐射波在目标处的电场矢量；
83.用平面波近似散射波，并根据目标散射的远场区中的给定观测点以及目标与接受天线之间的距离，确定散射电场；
84.定义一个复的极化散射矩阵来描述入射波和远区散射波各极化分量间的变换关系，构建极化散射矩阵；
85.定义目标的功率矩阵；
86.根据所述功率矩阵定义功率散射矩阵以及功率散射矩阵的迹；
87.根据所述功率散射矩阵以及所述功率散射矩阵的迹，定义去极化系数；
88.确定所述去极化系数的时间序列图以及功率谱密度结果；
89.根据所述时间序列图以及所述功率谱密度结果，得到进动周期信息。
90.本发明实施例的另一方面还提供了一种弹道目标rcs时间序列的进动周期提取装置，包括：
91.第一模块，用于计算目标的运动引起的电波入射角变化信息后，通过对所述目标的暗室静态测量数据进行插值，得到所述目标的动态数据序列；
92.第二模块，用于对所述动态数据序列进行平滑处理后，采用小波变换阈值去噪法进行去噪处理，得到初始数据；
93.第三模块，用于对所述初始数据进行全极化通道融合处理，得到进动周期信息。
94.本发明实施例的另一方面还提供了一种电子设备，包括处理器以及存储器；
95.所述存储器用于存储程序；
96.所述处理器执行所述程序实现如前面所述的方法。
97.本发明实施例的另一方面还提供了一种计算机可读存储介质，所述存储介质存储有程序，所述程序被处理器执行实现如前面所述的方法。
98.本发明实施例还公开了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行前面的方法。
99.下面结合说明书附图，对本发明的具体实现过程进行详细描述：
100.本发明立足于空间目标识别背景，针对弹道导弹微动特征提取精度不够高的问题，对rcs时间序列的预处理环节进行了详细的算法设计和仿真研究，并利用雷达回波中蕴含的丰富极化信息结合多源数据统计分析方法来进行全极化通道数据融合，从而对目标微动特征进行更有效的进动参数估计，使得特征提取精度更高，效果更稳健，具有重要应用价
值。如图1所示，本发明的整体步骤包括以下三大步骤：1、仿真数据获取；2、rcs时间序列预处理；3、全极化通道融合。
101.下面详细描述各个步骤的具体实现原理：
102.1、仿真数据获取：
103.在缺乏目标运动动态实测数据的情况下，一种可行的研究方法是通过计算目标的运动引起的电波入射角变化，用目标的暗室静态测量数据插值得到目标的动态数据序列。暗室仿真测量的逼真性是后期数据的准确获取及分析评估的基础，具体的仿真步骤如图2所示：首先建立质心运动模型，根据预先设定的导弹参数模拟计算出质心的位置矢量和速度矢量等信息，接着建立姿态运动模型，模拟推算出目标姿态角的变化，借助极化基变换关系，将计算得到的姿态角和目标静态散射测量数据相结合，插值得到特定姿态下目标的散射特性，进而推算并获取目标的动态回波。
104.本发明所使用的模拟弹头类目标有两类，模型具体形状及参数如下图3所示：分别为圆锥体和无缝隙锥球体。测量时目标的俯仰角和横滚角均为0
°
，方位角为0
°
～180
°
，雷达采用步进频测量方式，频扫范围8.75～10.75ghz，步长20mhz。中段总飞行时间370s，跟踪数据率为10hz，中段总数据长度为3701个数据点，极化状态为水平极化和垂直极化。
105.图4为仿真目标在中段的各通道rcs时间序列(以圆锥体为例)，由图可知该序列的时间起伏特性剧烈，趋势项明显，基本不具备平稳性，且四通道的数据具有很大差别。
106.2、rcs时间序列预处理过程：
107.savitzky-golay滤波器是一种用于平滑噪声数据的低通滤波器，通过在时域上不断移动窗口位置利用最小二乘法进行多项式拟合。该滤波器在平滑曲线的同时可以确保信号的形状、宽度不变，计算简便迅速，可操作性强，被广泛地运用于数据流平滑除噪。假定用一个n阶多项式f(i)来拟合一个窗口内的一组数据x[i]，(i＝-m,
…
m)，则拟合数据点与原数据点的残差平方和为：
[0108][0109]
当e对该多项式的系数bnk求偏导结果为零时，残差平方和最小，拟合效果最好：
[0110][0111]
解得：
[0112][0113]
当需要拟合的单边点数m，多项式阶次n及待拟合数据x[i]已知后，即可求出该多项式。拟合出来的多项式用于求取该窗口内的中心点估计值，然后不断移动窗口重复操作即可。savitzky-golay平滑公式为：
[0114]
[0115]
其中hi/h为平滑系数，由最小二乘法拟合多项式求得。以圆锥体数据为例，图5分别为hh通道趋势项提取结果及剩余项。可以看出该方法的趋势拟合效果很好，所提取出的趋势项可以反映rcs序列的宏观演化行为。
[0116]
虽然目前的rcs序列已接近平滑，但仍含有较强的噪声，因此要进行降噪处理。小波变换具有良好的时频域多分辨率处理能力，目前已被广泛应用于信号，图像处理等研究。本发明采用小波变换阈值去噪方法进行处理。小波的分解重构算法如下图6所示，具体包括以下步骤：
[0117]
1)对原始信号进行小波分解，得到各尺度近似系数(低频)和细节系数(高频)；
[0118]
2)按照一定规则对各尺度系数进行阈值处理；
[0119]
3)利用处理后的小波系数进行信号重构，得到去噪后信号；
[0120]
为了尽量避免实际信号中由不规则旋转体自旋所带来的峰值影响并提升去噪效果，采用软阈值处理，阈值选择方法采用无偏似然估计阈值，阈值标准采用固定式阈值，阈值处理根据每一层小波分解的噪声水平估计进行调整，去噪前后信号的结果如图7所示。
[0121]
可以看出，进行小波去噪处理后的rcs时间序列在波形上的周期性更加明显，波峰和波谷处的噪声波动被有效去除(95s-110s)，且由自旋所引起的幅度峰值也得到了抑制，因此具有很好的效果。图8是rcs数据经过整个预处理前后的效果。
[0122]
3、全极化通道融合处理过程：
[0123]
极化是与目标物理特性密切相关的特征量。理论上，当选定了散射空间坐标系以及相应的极化基后，目标的变极化效应就可以用极化散射矩阵的形式来表示，它代表了特定姿态和观测频率下目标的全部极化信息。目标的极化散射矩阵不仅取决于目标本身的物理结构和材料属性，也取决于观测条件。同时，收发天线坐标系及极化基的选取也直接决定着散射矩阵的取值。然而，对雷达观测者而言，在特定的观测状态下，坐标系或者极化基的改变并不会增添关于目标的新信息。“极化不变量”作为目标的特征信号，能够在给定的条件下直观的反映出目标的某些物理特性，而与雷达极化基选取无关。假设以平面电磁波照射一个有限尺寸目标，雷达发射天线辐射波在目标处的电场矢量具有如下形式：
[0124][0125]
式中e
ix
和e
iy
分别为入射场在x,y方向的复振幅，ωt-kz为入射波的瞬时相位，k代表波数。在目标散射的远场区某一给定观测点，散射波可用平面波近似，令r为目标与接受天线的距离，此时的散射电场为：
[0126][0127]
定义一个复的极化散射矩阵来描述入射波和远区散射波各极化分量间的变换关系：
[0128]esc
＝g(r)se
inc
#(3.1.3)
[0129]
式中g(r)称为球面波因子，s称为目标的极化散射矩阵，将e
sc
和e
inc
用矢量形式表示为：
[0130]
[0131][0132]
则极化散射矩阵可以写为：
[0133][0134]
其中s
hv
物理意义对应着以垂直极化照射目标时后向散射波的水平极化分量，其他类似。目标的graves功率矩阵是一个2阶hermite矩阵，定义为：
[0135][0136]
功率散射矩阵的迹p1定义为：
[0137]
p1＝|g
11
|2 |g
12
|2 |g
2１
|
２
|ｇ
22
|2#(3.1.8)
[0138]
p1表示天线所接收到的总功率，在不考虑目标表面条件的情况下大致可以表示目标的大小。令tr(s)代表极化散射矩阵的迹，则去极化系数d的定义为：
[0139][0140]
去极化系数d的取值范围为0～1，在一定程度上可以表示目标散射中心的数量。
[0141]
以锥球体数据为例，取仿真频率freq＝0.5hz，图9分别为去极化系数d的时间序列图及功率谱密度结果，可以看出，由锥球体数据全极化时间序列推导得到的去极化系数d在时域上具有很强的周期性，因而其功率谱密度图在仿真频率freq＝0.5hz时有很高的辨识度。
[0142]
最后，本发明实施例设置仿真进动频率freq＝0.4～1.0hz，分别取圆锥体和锥球体各通道极化数据做提取误差效果对比如图10和图11所示。其中，各通道数据均通过预处理，hh单通道处理方法为采用二重循环自相关法得到的时间序列频率估计，去极化系数d则由极化散射矩阵推导得出，各方法得到的序列由周期图法求归一化功率谱密度图，最终频率提取值为功率谱密度峰值点对应频率，频率提取误差则取自提取结果与仿真频率的距离。由误差对比图可以看出，在无观测噪声条件下，本发明所提方法对目标进动周期的提取效果显然要优于采用二重循环自相关法的单通道处理结果，且提取效果更稳定，平均误差由单通道的0.01减少至0.002。
[0143]
相比较于目前所提出的基于单通道处理的rcs时间序列进动周期提取方法，本发明所提出的方法从极化散射矩阵的角度出发，导出了极化不变量的概念，并进一步利用去极化系数d作为融合参数来进行进动特征提取。该方法充分融合了全极化通道数据，有效克服了因目标起伏带来的特征不明显和寄生分量问题，对目标微动特征进行了更有效更稳健的估计，进一步提高了周期提取精度与目标识别效率，填补了雷达目标识别领域中对目标极化信息利用较少的空缺，丰富了基于雷达回波rcs时间序列提取进动特征的方法，具有重要应用价值。
[0144]
综上所述，本发明立足于空间目标识别背景，针对弹道导弹微动特征提取精度不够高，极化信息利用不充分等问题，结合现代典型非平稳信号处理手段与多元统计数据融合方法，采用理论分析及暗室静态反演数据进行仿真研究，具有如下优点：
[0145]
1、将雷达rcs回波序列看做具有一定周期性的非平稳时间信号，对其进行趋势项
剔除与降噪等预处理，使其成为更易于观测及分析的时间序列，为后期微动特征提取的精确性提供帮助；
[0146]
2、提出了利用雷达回波极化信息进行全极化通道数据融合的方法，并利用最终融合结果精确提取了目标微动特征。经过仿真分析和误差对比表明，该融合方法在不同进动频率及进动角设置下周期提取结果均优于单通道处理结果，且提取效果更稳健。
[0147]
在一些可选择的实施例中，在方框图中提到的功能/操作可以不按照操作示图提到的顺序发生。例如，取决于所涉及的功能/操作，连续示出的两个方框实际上可以被大体上同时地执行或所述方框有时能以相反顺序被执行。此外，在本发明的流程图中所呈现和描述的实施例以示例的方式被提供，目的在于提供对技术更全面的理解。所公开的方法不限于本文所呈现的操作和逻辑流程。可选择的实施例是可预期的，其中各种操作的顺序被改变以及其中被描述为较大操作的一部分的子操作被独立地执行。
[0148]
此外，虽然在功能性模块的背景下描述了本发明，但应当理解的是，除非另有相反说明，所述的功能和/或特征中的一个或多个可以被集成在单个物理装置和/或软件模块中，或者一个或多个功能和/或特征可以在单独的物理装置或软件模块中被实现。还可以理解的是，有关每个模块的实际实现的详细讨论对于理解本发明是不必要的。更确切地说，考虑到在本文中公开的装置中各种功能模块的属性、功能和内部关系的情况下，在工程师的常规技术内将会了解该模块的实际实现。因此，本领域技术人员运用普通技术就能够在无需过度试验的情况下实现在权利要求书中所阐明的本发明。还可以理解的是，所公开的特定概念仅仅是说明性的，并不意在限制本发明的范围，本发明的范围由所附权利要求书及其等同方案的全部范围来决定。
[0149]
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0150]
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。
[0151]
计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(ram)，只读存储器(rom)，可擦除可编辑只读存储器(eprom或闪速存储器)，光纤装置，以及便携式光盘只读存储器(cdrom)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其
他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
[0152]
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
[0153]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0154]
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
[0155]
以上是对本发明的较佳实施进行了具体说明，但本发明并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换，这些等同的变形或替换均包含在本技术权利要求所限定的范围内。