发明名称 一种基于粒子群优化的跳频信号参数估计方法
摘要 本发明公开了一种基于粒子群优化的跳频信号参数估计方法,包括如下步骤:基于跳频信号的多分量模型,采用多峰函数粒子群算法从原始跳频信号中并行搜索跳频分量,寻找各个跳频周期内的时间和频率中心位置,实现粒子在跳频分量时频中心位置的聚集;通过粒子群算法输出的粒子个体最优值获取跳频信号参数。本发明具有简单方便、调整参数少,搜素能力强,尤其是适应宽信噪比范围内估计等优势。
申请公布号 CN105429722A 申请公布日期 2016.03.23
申请号 CN201510712920.0 申请日期 2015.10.20
申请人 信阳师范学院 发明人 郭建涛;涂友超;王林;易伟;刘梦果
分类号 H04B17/391(2015.01)I;H04B17/336(2015.01)I;H04B1/713(2011.01)I 主分类号 H04B17/391(2015.01)I
代理机构 代理人
主权项 一种基于粒子群优化的跳频信号参数估计方法,其特征在于,包括如下步骤:S1、基于跳频信号的多分量模型,采用多峰函数粒子群算法从原始跳频信号中并行搜索跳频分量,寻找各个跳频周期内的时间和频率中心位置,实现粒子在跳频分量时频中心位置的聚集:S11、采用粒子群局域模型:假设在D维空间有M个粒子,i=1,2,…M,其每个粒子属性由位置x<sub>i</sub>=(x<sub>i1</sub>,x<sub>i2</sub>,…,x<sub>iD</sub>)和速度v<sub>i</sub>=(v<sub>i1</sub>,v<sub>i2</sub>,…v<sub>iD</sub>)决定;同时记录其在搜索过程中所经历过的个体最优位置P<sub>i</sub>=(P<sub>i1</sub>,P<sub>i2</sub>,…,P<sub>iD</sub>),邻域最优位置P<sub>n</sub>=(P<sub>n1</sub>,P<sub>n2</sub>,…,P<sub>nD</sub>);对第t+1代的第i个粒子第d维方向上,通过以下公式更新粒子速度和位置:v<sub>id</sub>(t+1)=wv<sub>id</sub>(t)+c<sub>1</sub>r<sub>1d</sub>(t)(P<sub>id</sub>(t)‑x<sub>id</sub>(t))+c<sub>2</sub>r<sub>2d</sub>(t)(P<sub>nd</sub>(t)‑x<sub>id</sub>(t));x<sub>id</sub>(t+1)=x<sub>id</sub>(t)+v<sub>id</sub>(t+1);i=1,2,…,M and d=1,2,…,D;式中,w、c<sub>1</sub>和c<sub>2</sub>分别表示惯性权重和学习因子,r<sub>1</sub>和r<sub>2</sub>是0~1内的随机数,经过多次迭代,实现粒子在适应度函数的不同峰值位置聚集;S12、粒子邻域确定:定义粒子i和j的测度如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>m</mi><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>P</mi><mrow><mi>j</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mrow></mfrac><mo>|</mo></mrow>]]></math><img file="FSA0000122613440000011.GIF" wi="491" he="160" /></maths>式中,P<sub>j</sub>表示第j个粒子当前个体最优位置,f(P<sub>j</sub>)和f(x<sub>i</sub>)分别表示第j个粒子的个体最优适应度和当前粒子适应度;测度中的分母维度选择2,表示时间;S13、原子库构建:选用实型高斯函数作为过完备库原子的原型函数,定义为:<img file="FSA0000122613440000021.GIF" wi="720" he="78" />其中,尺度因子s<sub>k</sub>控制原子时域方向所占宽度,参数t<sub>k</sub>、f<sub>k</sub>和<img file="FSA0000122613440000022.GIF" wi="47" he="44" />分别表示高斯原子的时间中心和频率中心及其相位,<img file="FSA0000122613440000023.GIF" wi="67" he="75" />是原子的规范化系数,k是原子索引,表示粒子个数;因此,D=4,第二维表示时间,原子由一个参数集<img file="FSA0000122613440000024.GIF" wi="275" he="53" />完全表示,分别潜在地表示跳频周期、每一跳频周期内的时间和频率中心及其相位;S2、通过粒子群算法输出的粒子个体最优值获取跳频信号参数:S21、对原子的频率参数按照上升排序,得到在0到f<sub>s</sub>/2之间的频率值f<sub>i</sub>,i=1,…,M及其索引ind1;S22、以最小频率分辨率或最小跳频频率间隔(若已知的情况下)为阈值th,通过下式求取f<sub>i</sub>的跳变时刻对应的索引ind2′,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>i</mi><mi>n</mi><mi>d</mi><msup><mn>2</mn><mo>&prime;</mo></msup><mo>=</mo><munder><mi>arg</mi><mi>i</mi></munder><mrow><mo>(</mo><mo>(</mo><mrow><msub><mi>f</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>f</mi><mi>i</mi></msub></mrow><mo>)</mo><mo>&gt;</mo><mi>t</mi><mi>h</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FSA0000122613440000025.GIF" wi="581" he="105" /></maths>并加上粒子首尾两个值构成排序后的频率近似区间划分,即ind2=[1;ind2′;M];S23、结合索引ind1和ind2找到这些近似频率对应的粒子,对所有区间内的频率参数值求平均获取跳频频率f<sub>k</sub>,k=1,…,K,其中K表示跳变时刻个数;对所有区间内的时间参数值求平均获取跳频分量的时间中心位置t<sub>k</sub>,k=1,…,K;f<sub>k</sub>和t<sub>k</sub>为每个跳频分量的时频中心参数;S24、对估计的t<sub>k</sub>,k=1,…,K估计参数,按照上升排序,差分后求取平均值,得到跳频周期的估计值<img file="FSA0000122613440000026.GIF" wi="72" he="66" />进一步,由<img file="FSA0000122613440000027.GIF" wi="370" he="123" />获取各个跳变时刻估计值。
地址 464000 河南省信阳市南湖路237号