专利名称：基于模糊关系熵的混沌扩频序列复杂度测度方法
技术领域：
本发明属于无线通信技术领域，特别是一种混沌伪随机扩频序列复杂度的测度方 法，可用于对混沌伪随机扩频序列的复杂度的测度。
背景技术：
混沌伪随机序列是通过对混沌映射产生的序列进行量化而得到的。混沌伪随机序 列具有类噪声、对初始值敏感以及非周期等特性，因此，可以取代传统的伪随机序列应用于 保密扩频通信之中。由于混沌伪随机扩频序列的复杂度是指它与随机序列的相似程度，是对利用序列 的部分恢复出整体的难易程度的度量，因此它是衡量保密通信系统中扩频序列抗干扰和截 获能力的重要指标。在序列复杂度测度研究中，众所周知，由于线性复杂度分析算法 (Berlekamp-Massey)是针对常规伪随机扩频序列提出的测度方法，不能够胜任混沌伪随机 扩频序列复杂度的准确测度。而另外两种已经提出的混沌伪随机扩频序列复杂度测度方法 为近似熵(ApEn)的方法和符号熵方法。其中，ApEn方法能够较为准确地测度不同混沌伪 随机扩频序列的复杂度，但是存在测度结果受参数选取影响较大的问题，从而限制了其测 度实用性。符号熵方法是采用符号动力学方法测度混沌伪随机扩频序列复杂度的方法。该 方法能够有效地测度不同复杂度的混沌伪随机扩频序列，并且进一步降低了参数选取对测 度结果的影响。但是，由于符号动力学方法中固有的粗粒化标示特点，该方法在序列符号空 间元素较多或者窗口向量长度较大的情况下，计算复杂度巨大，无法有效完成混沌伪随机 扩频序列的复杂度测度；此外，该方法在测度前必须已知待测扩频序列的符号空间大小，而 在实际应用中要得到这个先验知识是很困难的。
发明内容
本发明的目的在于克服上述现有技术的不足，将模糊理论中关于模糊关系的概念 引入到混沌伪随机扩频序列复杂度的统计学测度方法之中，提出一种基于模糊关系熵的混 沌扩频序列复杂度测度方法，以在无需已知待测扩频序列的符号空间大小的情况下，扩展 混沌伪随机扩频序列复杂度测度的实用性，提高混沌伪随机扩频序列复杂度测度的准确 性。实现本发明目的的技术方案是根据混沌伪随机扩频序列的特性，构造出相适应 的模糊关系，将模糊关系度量模块嵌入统计度量方法之中，具体步骤如下(1)将观察长度为N的混沌伪随机扩频序列，分割成为长度为m的窗口向量 sm(i)sm(i) = [u(i)，u(i+l)，…，u(i+m-l)]，i = 1，2，…，N—m+1其中，m= 1,2,…，N_2，u(i)是给定的长度为N的非负整数混沌伪随机扩频序列 [u(l)，U(2)，-,u(N)]中的元素，SPsm(i)为待测混沌伪随机扩频序列中按照原有序列元素顺序，从第i个元素开始取出的m个元素构成的向量；(2)计算m维的第i个窗口向量SnXi)中的元素u(i+k)与第j个窗口向量sm(j) 中的元素u(j+k)之间的模糊关系度量F^=R{u{i + k),u{j + k)) = e ()( - )式中，i，j = l，2，…，N-m+l，k = 0，1，…，m_l ;r表示分辨率参数，为正实数，用
于确定不同窗口向量之间距离的模糊关系度量精度，R( ·)是为计算元素之间的模糊关系 度量<而定义的模糊隶属函数；(3)根据步骤(2)中确定的元素之间的模糊关系度量if计算得到窗口向量sm(i) 和sm(j)之间的模糊关系度量^T
m k=0 m k=0 rn k=0
(4)通过计算m+1维的第i个窗口向量sm+1(i)和第j个窗口向量sm+1 (j)之间的
最大距离，得到sm+1⑴和sm+1 (j)之间的模糊关系度量
2
_ 6] ψ = R Xsm+1 (0,U)) =式中，i，j = l，2，...,N-m,R' (·)是用来计算窗口向量之间模糊关系度量的模 糊隶属函数，^+1表示m+1维窗口向量sm+1⑴禾口 sm+1(j)的最大距离= max (u(i + k)-u(j + k))；(5)通过求和平均的方法计算模糊关系度量《在j上的统计平均C/"，以及在 j上的统计平均Qk+1 ；(6)通过求和平均的方法计算步骤(5)中的Cf在i上的统计平均Φ"1，以及Q"+l在 i上的统计平均Φ"1+1 ；(7)根据步骤(6)中计算得到的Φ"1和Φ"1+1值，定义并计算混沌伪随机扩频序列 的模糊关系熵F-REn F-REn = lnOm-lnOm+10
图1是本发明的实现步骤框图；图2是本发明仿真使用的现有混沌序列的相空间结构图；图3是本发明与现有的ApEn测度方法及符号熵测度方法的测度值随测量维度参 数的变化情况仿真结果；图4是本发明与ApEn测度方法的测度值随分辨率参数变化的曲线。
具体实施例方式参照图1，本发明包括以下步骤
步骤1，将观察长度为N混沌伪随机扩频序列分割成为长度为m的窗口向量Sm(i)。令正整数N是混沌伪随机扩频序列的观测长度，对于正整数m，m彡N_2，给定非负 整数序列[U(I)，u (2)，…，u(N)]，其中，u(n) e {0,1,…，2^1}，1彡η彡N，q是混沌伪随 机扩频序列中每个序列值的二进制比特数，形成m维非负整数序列空间sm(l)，sm(2),…， sm(N-m+l) e Nm，其中 sm(i)表示为sm(i) = [u(i)，u(i+l)，...，u(i+m-l)]，l 彡 i <N-m+l。1)步骤2，计算第i个窗口向量sm(i)中的元素u(i+k)与第j个窗口向量sm(j)中 的元素u(j+k)之间的模糊关系度量if。(2. 1)设长度为 m 的窗口 向量为{sm(i) = [u(i)，u(i+l)，…，u(i+m_l) ]}和 {sm(j) = [u(j)，u(j+l)，...，u(j+m-l)]}，其中 i，j = 1，2，…N—m+1(2.2)定义一个模糊隶属函数R(*)Jfsm(i) ^P sm (j)中的元素u(i+k)和u(j+k)
的值代入R(·)中，计算得到if:Fk ^ Κ{μ{ + k)Mj + k)) = e-ln(2)KM^O±M)22)式中，i，j = l，2，…，N_m+l，k = 0，1，…，m_l ;r表示分辨率参数，为正实数，用
于确定不同窗口向量之间距离的模糊关系度量精度。步骤3，通过对元素之间的模糊关系度量if进行求和平均，计算m维的第i个窗口
向量sm(i)和第j个窗口向量sm(j)之间的模糊关系度量^
1 m-\ι m-1ι m-1 ^n,Mi+k)-u{j+k)L2
4 二去Σ^Ν去Σ所—·+从‘+幻)=去Σ^ _r^) ° 3)
m k=0 m k=0 m k=0步骤4，计算m+1维窗口向量sm+1(i)禾Psm+1(j)的最大距离^+1,得到sm+1 (i)和 sm+1(j)之间的模糊关系度量％+1。(4. 1)令d『+i表示m+1维窗口向量sm+1⑴和sm+1 (j)的最大距离 ψχ = kmd^j\u{i + k)-u{j + A:)|)； 4)(4.2)定义一个计算窗口向量之间模糊关系度量的模糊隶属函数为R' (·)，将 代入R' (·)中，从而得到m+1维的第i个窗口向量sm+1⑴和第j个窗口向量sm+1 (j)
之间的模糊关系度量Bm+X ^ R ^ sm+l ^ =，^ = 1,2’...斤-饥。5)步骤5，利用求和平均的方法，计算模糊关系度量在j上的统计平均，以及 β『+1在·]·上的统计平均C/"+1 步骤6，利用求和平均的方法，计算步骤6中(^在i上的统计平均Φ"1，以及Qn+1在 i上的统计平均Φ"1+1 步骤7，根据步骤6中计算得到的Φ-和Φ"1+1值，计算混沌伪随机扩频序列的模糊 关系熵(F-REn)F-REn = Fre (m，r, N) = InOm-InO111+110) 其中，Fre ( ·)表示计算混沌伪随机扩频序列的模糊关系熵的函数，参数为m、r与
N0本发明的效果可以通过以下理论推导与分析进一步说明(A)经过推导和证明得到以下两个重要性质性质1 纽 Fre(m，r,N) = 0性质2 当被测序列为常数序列时，Fre (m, r，N) = 0。(B)本发明与现有的ApEn测度方法及符号熵测度方法的性能分析为了对比分析方便，以logistic映射和耦合映射格子(Coupled map lattice) 系统迭代产生的伪随机序列为例，说明本发明的F-REn测度方法作为混沌伪随机扩频序列 复杂度测度具有更加优越的对序列符号空间的适用性，对测量维度的敏感性和对分辨率的 鲁棒性。图2表示了 logistic序列和耦合映射格子序列的相空间结构，其中图2(a)表示 logistic序列的相空间结构，图2(b)表示耦合映射格子序列的相空间结构，由图2可知，耦 合映射格子序列比logistic映射序列具有更高的复杂度。Bi)本发明与现有的ApEn测度方法及符号熵测度方法对序列符号空间的适用性 仿真分析序列符号空间是指待测混沌伪随机扩频序列的可选取元素集的大小，用参数q表 示，复杂度测度方法对序列符号空间的适用性定义为具有不同符号空间q的混沌伪随机 扩频序列对序列复杂度测度结果的影响程度。这种影响越大，序列复杂度测度方法对序列 符号空间的适用性越差。序列复杂度测度方法对序列符号空间的适用性主要体现在两个方面，一是待测序 列的符号空间是不是复杂度测度方法必须已知的信息；二是不同符号空间的序列是否影响 不同测度方法的测度结果。对于第一个问题，在现有的符号熵方法中，待测度序列的符号空间q作为先验知 识，必须是测度之前的已知信息。而现有的ApEn测度方法和本发明的F-REn测度方法并不 需要已知待测序列的符号空间大小，因此可以对任意未知符号空间大小的混沌伪随机序列 进行复杂度测度。对于第二个问题，不同的复杂度测度方法对不同符号空间的序列的测度适用性不同，用本发明的F-REn测度方法与现有的ApEn测度方法及符号熵测度方法对不同符号空间 (q)的序列进行复杂度测度仿真，得出表1所示结果。 表1不同测度方法对不同符号空间(q)的序列的复杂度测度仿真结果 从表1中可以看出，在现有的ApEn测度方法中，仅仅当q = 4时，两种序列的复杂 度能够被正确测度，随着q > 4，出现耦合映射格子序列测度值小于logistic映射序列测度 值的情况，这与实际序列复杂度相反，因此无法正确进行混沌伪随机扩频序列复杂度的正 确测度。在现有的符号熵测度方法中，当q = 4时，序列复杂度能够被正确测度，但是随着 q的不断增大，当q > 9时，由于序列空间较大，造成计算量急剧增大，仿真结果出现数据溢 出现象，因此在该仿真参数条件下，现有的符号熵测度方法对q > 9的混沌伪随机扩频序列 的复杂度测度的实用性不强。在本发明提出的F-REn测度方法中，随着参数q的增大，两种序列的测度值结果始 终与其实际序列复杂度相匹配，因此，序列符号空间对本发明的F-REn测度方法影响很小， 说明本发明的F-REn测度方法具有很好的序列适用性。B2)本发明与现有的ApEn测度方法及符号熵测度方法对测量维度的敏感性仿真 分析测量维度参数m是指待测序列的窗口向量长度。该参数的变化会对最终测度结果 产生影响。在实际测度中，我们总是希望测量维度参数的变化对最终测度结果的影响越小 越好。为了进一步说明本发明提出的F-REn测度方法的优越性，图3表示了随着测量维度参 数m的变化，利用本发明的F-REn测度方法与现有的ApEn测度方法及符号熵测度方法得到 的序列复杂度测度仿真结果。其中，图3(a)、图3(b)和图3(c)为q = 6情况下，本发明的 F-REn测度方法、现有的ApEn测度方法和现有的符号熵测度方法随着m变化所得到的测度 结果，图3(d)、图3(e)和图3(f)为q= 12情况下，本发明的F-REn测度方法、现有的ApEn 测度方法和现有的符号熵测度方法随着m变化所得到的测度结果。从图3中可以看出，对q = 6的序列进行复杂度测度时，在现有的ApEn测度方法 中，当m > 2时就出现复杂度测度错误的现象，而对q = 12的序列进行复杂度测度时，由于 q增大，能够准确测度的m取值进一步变小，因此，只有m= 1的情况下能够正确测度两种序 列的复杂度。
对于现有的符号熵测度方法，在q = 6的情况下，它在me [2,4]范围内能够很好 的完成序列复杂度的测度，但是随着m的逐渐增大，当m > 4时，由于计算量太大，符号熵测 度仿真出现数据溢出而无法得到测度值。并且在q = 12时，由于序列状态空间进一步增大， 造成m = 2时符号熵测度值出现错误，m > 2时符号熵测度仿真出现数据溢出无测度值。而对于本发明的F-REn测度方法来说，不仅对q = 6的序列，在m e [1，10]范围 内均能够完成复杂度的准确测度，而且对于q = 12的大符号空间序列，在m e [1,8]范围 内也能够完成复杂度的准确测度。因此，本发明的F-REn测度方法对测量维度参数的敏感 性非常小，测量维度参数的变化对最终测度结果的影响很小。B3)本发明与现有的ApEn测度方法对分辨率的鲁棒性仿真分析由于现有的符号熵测度方法中没有分辨率参数r，所以在分辨率参数鲁棒性的仿 真中仅对本发明的F-REn测度方法和现有的ApEn测度方法进行对比分析。在基于ApEn和 F-REn的FH/SS序列复杂度测度方法中，分辨率性能是衡量待测序列不同截取窗口之间差 异的门限水平，它确定了复杂度测度过程的精度。图4为本发明的F-REn测度方法和现有的ApEn测度方法对序列的复杂度测度值 随分辨率参数r的变化情况，其中，图4(a)、图4(b)表示在待测序列采样点数取N = 2000 的情况下，本发明的F-REn测度方法和现有的ApEn测度方法的测度值随分辨率参数的变化 情况，图4 (c)、图4 (d)表示在待测序列采样点数取N = 500的情况下，本发明的F-REn测度 方法和现有的ApEn测度方法的测度值随分辨率参数的变化情况。从图4中可知，在待测序列采样点数取N = 2000的情况下，现有的ApEn测度方法 在r = 0. 5和r = 5时出现序列复杂度测度错误的情况；而本发明的F-REn测度分辨率参 数在re
时均能够正确的完成复杂度测度。随着待测度序列采样点数的降低，由 于统计样本点减少，造成测度结果的出现不准确情况。从图4中可以看出，在待测序列采样 点数降低到N = 500的情况下，现有的ApEn测度方法能够准确测度的分辨率参数选取范围 进一步缩小，只有r e [2,4. 5]的情况下能够正确测度；而本发明的F-REn测度方法中分辨 率参数在re
时仍然能够正确的完成复杂度测度。因此，可以看出，与现有的ApEn 测度方法相比，本发明的F-REn测度方法对分辨率参数的变化具有很好的鲁棒性。综上，本发明所提出的基于模糊关系熵(F-REn)计算混沌伪随机扩频序列复杂度 的方法，在相同数据要求条件下，与其它混沌伪随机扩频序列复杂度测度方法相比，F-REn 测度方法具有更加良好的对序列符号空间的适用性、对测量维度的敏感性和对分辨率参数 的鲁棒性。
权利要求
一种基于模糊关系熵的混沌扩频序列复杂度测度方法，包括如下步骤(1)将观察长度为N的混沌伪随机扩频序列，分割成为长度为m的窗口向量sm(i)sm(i)＝[u(i)，u(i+1)，…，u(i+m 1)]，i＝1，2，…，N m+1其中，m＝1，2，…，N 2，u(i)是给定的长度为N的非负整数混沌伪随机扩频序列[u(1)，u(2)，…，u(N)]中的元素，即sm(i)为待测混沌伪随机扩频序列中按照原有序列元素顺序，从第i个元素开始取出的m个元素构成的向量；(2)计算m维的第i个窗口向量sm(i)中的元素u(i+k)与第j个窗口向量sm(j)中的元素u(j+k)之间的模糊关系度量 <mrow><msubsup> <mi>F</mi> <mi>ij</mi> <mi>k</mi></msubsup><mo>=</mo><mi>R</mi><mrow> <mo>(</mo> <mi>u</mi> <mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>k</mi><mo>)</mo> </mrow> <mo>,</mo> <mi>u</mi> <mrow><mo>(</mo><mi>j</mi><mo>+</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><msup> <mi>e</mi> <mrow><mo>-</mo><mi>ln</mi><mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo></mrow><mo>·</mo><msup> <mrow><mo>(</mo><mfrac> <mrow><mo>|</mo><mi>u</mi><mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><mi>u</mi><mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mi>k</mi> <mo>)</mo></mrow><mo>|</mo> </mrow> <mi>r</mi></mfrac><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow></msup> </mrow>式中，i，j＝1，2，…，N m+1，k＝0，1，…，m 1；r表示分辨率参数，为正实数，用于确定不同窗口向量之间距离的模糊关系度量精度，R(·)是为计算元素之间的模糊关系度量而定义的模糊隶属函数；(3)根据步骤(2)中确定的元素之间的模糊关系度量计算得到窗口向量sm(i)和sm(j)之间的模糊关系度量 <mrow><msubsup> <mi>A</mi> <mi>ij</mi> <mi>m</mi></msubsup><mo>=</mo><mfrac> <mn>1</mn> <mi>m</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>k</mi><mo>=</mo><mn>0</mn> </mrow> <mrow><mi>m</mi><mo>-</mo><mn>1</mn> </mrow></munderover><msubsup> <mi>F</mi> <mi>ij</mi> <mi>k</mi></msubsup><mo>=</mo><mfrac> <mn>1</mn> <mi>m</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>k</mi><mo>=</mo><mn>0</mn> </mrow> <mrow><mi>m</mi><mo>-</mo><mn>1</mn> </mrow></munderover><mi>R</mi><mrow> <mo>(</mo> <mi>u</mi> <mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>k</mi><mo>)</mo> </mrow> <mo>,</mo> <mi>u</mi> <mrow><mo>(</mo><mi>j</mi><mo>+</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>1</mn> <mi>m</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>k</mi><mo>=</mo><mn>0</mn> </mrow> <mrow><mi>m</mi><mo>-</mo><mn>1</mn> </mrow></munderover><msup> <mi>e</mi> <mrow><mo>-</mo><mi>ln</mi><mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo></mrow><mo>·</mo><msup> <mrow><mo>(</mo><mfrac> <mrow><mo>|</mo><mi>u</mi><mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><mi>u</mi><mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mi>k</mi> <mo>)</mo></mrow><mo>|</mo> </mrow> <mi>r</mi></mfrac><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow></msup><mo>;</mo> </mrow>(4)通过计算m+1维的第i个窗口向量sm+1(i)和第j个窗口向量sm+1(j)之间的最大距离得到sm+1(i)和sm+1(j)之间的模糊关系度量 <mrow><msubsup> <mi>B</mi> <mi>ij</mi> <mrow><mi>m</mi><mo>+</mo><mn>1</mn> </mrow></msubsup><mo>=</mo><msup> <mi>R</mi> <mo>′</mo></msup><mrow> <mo>(</mo> <msup><mi>s</mi><mrow> <mi>m</mi> <mo>+</mo> <mn>1</mn></mrow> </msup> <mrow><mo>(</mo><mi>i</mi><mo>)</mo> </mrow> <mo>,</mo> <msup><mi>s</mi><mrow> <mi>m</mi> <mo>+</mo> <mn>1</mn></mrow> </msup> <mrow><mo>(</mo><mi>j</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><msup> <mi>e</mi> <mrow><mo>-</mo><mi>ln</mi><mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo></mrow><mo>·</mo><msup> <mrow><mo>(</mo><mfrac> <msubsup><mi>d</mi><mi>ij</mi><mrow> <mi>m</mi> <mo>+</mo> <mn>1</mn></mrow> </msubsup> <mi>r</mi></mfrac><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow></msup> </mrow>式中，i，j＝1，2，…，N m，R′(·)是用来计算窗口向量之间模糊关系度量的模糊隶属函数，表示m+1维窗口向量sm+1(i)和sm+1(j)的最大距离 <mrow><msubsup> <mi>d</mi> <mi>ij</mi> <mrow><mi>m</mi><mo>+</mo><mn>1</mn> </mrow></msubsup><mo>=</mo><munder> <mi>max</mi> <mrow><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi> </mrow></munder><mrow> <mo>(</mo> <mo>|</mo> <mi>u</mi> <mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>k</mi><mo>)</mo> </mrow> <mo>-</mo> <mi>u</mi> <mrow><mo>(</mo><mi>j</mi><mo>+</mo><mi>k</mi><mo>)</mo> </mrow> <mo>|</mo> <mo>)</mo></mrow><mo>;</mo> </mrow>(5)通过求和平均的方法计算模糊关系度量在j上的统计平均以及在j上的统计平均(6)通过求和平均的方法计算步骤(5)中的在i上的统计平均Φm，以及在i上的统计平均Φm+1；(7)根据步骤(6)中计算得到的Φm和Φm+1值，定义并计算混沌伪随机扩频序列的模糊关系熵F REnF REn＝lnΦm 1nΦm+1。FSA00000186142100011.tif,FSA00000186142100013.tif,FSA00000186142100014.tif,FSA00000186142100015.tif,FSA00000186142100017.tif,FSA00000186142100018.tif,FSA00000186142100021.tif,FSA00000186142100023.tif,FSA00000186142100024.tif,FSA00000186142100025.tif,FSA00000186142100026.tif,FSA00000186142100027.tif,FSA00000186142100028.tif
全文摘要
本发明公开了一种基于模糊关系熵的混沌扩频序列复杂度测度方法，主要解决现有的序列复杂度测度方法实用性低和准确性差的问题。其实现步骤是(1)将待测序列分割成一组长度为m的窗口向量；(2)计算m维窗口向量之间的模糊关系度量；(3)计算m+1维窗口向量之间的模糊关系度量；(4)对步骤(2)、(3)中得出的m维和m+1维窗口向量的模糊关系度量分别进行统计平均；(5)将步骤(4)中统计平均的结果进行对数差运算，得到混沌伪随机扩频序列的模糊关系熵(F-REn)。本发明具有测度实用性强和准确性高的优点，可用于对任意未知符号空间大小的混沌伪随机扩频序列的复杂度测度。
文档编号H04B1/69GK101902323SQ20101022522
公开日2010年12月1日 申请日期2010年7月12日 优先权日2010年7月12日
发明者关磊, 司江勃, 姚磊, 李赞, 杜军朝, 蔡觉平, 谢梦, 郝本建, 陈小军, 高锐 申请人:西安电子科技大学