本发明属于信号处理参技术领域,具体涉及一种超宽带信号参数估计方法,可用于雷达、无源侦察等高速数据采集和处理。
背景技术:
参数估计是现代信号处理的重要研究内容之一,一直受到国内外学者的广泛关注。近年来,随着电子技术不断发展及高速大带宽电子设备需求的不断增长,超宽带信号参数估计已经成为参数估计研究的热点。
目前,超宽带信号参数估计的研究主要有两个方面:基于压缩感知的信号参数估计和基于光电复合的信号参数估计。其中:
d.donoho、j.romberg、e.candes和t.taot等科学家于2004年提出基于信号稀疏特性的压缩感知理论,该理论的提出使得宽带频谱感知技术可以突破奈奎斯特采样定理的束缚,在远小于奈奎斯特采样率的条件下,以较低的adc采样速率对超宽带信号进行采样,然后通过非线性恢复算法重建信号。但是这种估计有其限制条件,如信号必须具有稀疏性,同时存在恢复算法复杂度高和信噪比要求较高的问题。
基于光电复合的超宽带信号参数估计方法:其利用高精度、低时间抖动的光脉冲对超宽带微波信号进行调制,并将调制后携带了微波信号的高速光脉冲序列经过解复用器并行输出为n路的低速信号,每一路低速信号经过光电转换后采用低速率、高精度的adc进行保持和量化,进而转换为数字信号,最终利用电域的信号处理算法进行信号恢复,来获取信号参数。这种方法虽兼具光域和电域各自的优势,是研究超宽带信号参数估计的重要途径,但其信号恢复的复杂度很高,不易在实际工程中应用。
技术实现要素:
本发明的目的在于针对上述现有技术的不足,提出一种基于光电复合的超宽带信号参数估计方法,以降低信号参数估计的复杂度,便于实际工程应用。
实现本发明目的的技术方案是采用无恢复算法的压缩感知技术进行信号参数的直接估计,步骤包括如下:
1)用脉冲激光光源产生超短光脉冲,将超短光脉冲通过第一单模光纤拉伸,变为一次拉伸的光脉冲e。
2)将射频信号x(t)通过马赫曾德尔调制器调制到一次拉伸的光脉冲e上,得到调制信号eout。
3)将调制信号eout通过第二单模光纤进行拉伸,得到二次拉伸信号e'out。
4)用光电探测器将二次拉伸信号e′out转变为电信号,即光域拉伸信号x'(t)。
5)建立光域拉伸信号x'(t)下的两个基函数:u1[n]=cos(2πfktn),u2[n]=sin(2πfktn)
其中,
6)用gold码设计观测矩阵φ:φ=dhr,其中d为低速ad采样,h为低通滤波器,r为n×n的对角矩阵;将观测矩阵写成列向量的形式,并根据该列向量得到感知矩阵
7)根据感知矩阵
8)根据基扩展系数
9)将在fk下计算出的wk按照从小到大顺序排列,取wk为最大值下的fk即为估计的信号载频
本发明与现有技术相比具有以下优点:
本发明由于利用时域拉伸技术对超宽带信号进行光域压缩,利用压缩采样技术对拉伸后的信号进行第二次压缩,大大减小了adc的采样速率;同时由于在不恢复信号波形的情况下直接进行信号参数的估计,大大减小压缩感知算法的参数估计运算量。
同时仿真结果表明,相比传统的超宽带信号参数估计方法,本发明具有更大的信号带宽、更高的估计精度,同时参数估计运算量更小。
附图说明
图1是本发明的实现流程图。
图2是用本发明多频信号进行仿真实验的频谱分布估计图。
图3是用本发明对多频信号估计的幅度误差图。
图4是用本发明对混合信号进行仿真实验的频谱分布估计图。
图5是用本发明对混合信号的幅度估计误差图。
具体实施方式
参照图1,本发明的实施步骤如下:
步骤1,对超短光脉冲进行拉伸。
光纤对光脉冲进行拉伸的原理如下:
在不考虑非线性效应和高阶色散的条件下,光脉冲在单模光纤中传输的频域通解表示为
其中,ω为光波角频率相对于脉冲中心角频率的偏离;z是光脉冲在光线中的传输距离;
一般高斯光脉冲光场归一化幅度的时域通解为
对于入射光场为无初始啁啾的高斯光脉冲,其时域通解为
式中τe为脉冲振幅1/e处的脉冲半宽度。由式<2>和式<3>,得到沿光纤方向任意一点z处的光脉冲光场归一化幅度为
由式<4>可知,光脉冲在传输过程中形状不变,但宽度增加。
将式<4>写成
定义色散长度
根据上述原理,本步骤的具体实现如下:
1.1)用脉冲激光光源随机产生超短光脉冲;
1.2)设置单模光纤:即选用群速度色散参数β2=20ps2/km,长度l1=1km的单模光纤作为第一单模光纤,选用群速度色散参数β2=20ps2/km、长度l2=4km的单模光纤作为第二单模光纤;
1.3)将超短脉冲通过第一单模光纤拉伸,变为第一路一次拉伸的光脉冲e。
步骤2,用马赫曾德尔调制器将射频信号调制到光脉冲上。
马赫曾德尔调制器将射频信号调制光脉冲上的原理如下:
输入马赫曾德尔调制器的光脉冲的表达式为
其中
对于理想的调制器,其消光比δ为无穷大,这时γ=1,所以,经过理想的调制器后,调制信号eout为:
光脉冲的相位变化为:
其中,v(t)为调制器的输入电压信号,vπ为半波电压,为调制器产生附加相位为π时的电压,对应为整个光波相位周期的一半;
由式<7>和式<8>可知,调制信号eout可写成如下形式:
其中,v1(t)为调制器的上交流电压v1rf和上直流电压v1dc之和,v2(t)为调制器的下交流电压v2rf和下直流电压v2dc之和。
根据上述原理,本步骤的具体实现如下:
2.1)生成射频信号x(t);
2.2)将射频信号x(t)通过马赫曾德尔调制器调制到一次拉伸的光脉冲e上,得到调制信号eout。
步骤3,对调制信号进行二次拉伸。
二次拉伸的原理如下:
将射频信号在时域上进行展宽的过程主要是由群速度色散参数β2决定,然而光线中还存在着三阶、四阶等更高阶等更高阶的非线性色散现象,光纤的模传输特性方程β(ω)在信号频谱中心频率ω0附近展开为泰勒级数,其表示如下:
其中
群速度色散参数β2在光纤中可由色散参数d(λ)表示,其色散参数d(λ)可表示为
式中s0为λ=λ0时的色散斜率,λ为光脉冲的波长,λ0为零色散波长,又由
可得
其中λr代表参考波长的值,τg(λ)为单位长度的群时延。
因此,在经过第一单模光纤后脉冲展宽时延t1为
t1=l1×τg1(λ),<14>
其中l1为第一单模光纤的长度;
同理,在经过第二单模光纤进行进一步展宽后,脉冲展宽时延t为
t=t1+t2=l1×τg1(λ)+l2×τg2(λ)<15>
其中l2为第二单模光纤的长度。
通过式<14>和式<15>可知,只要两段光纤有相等的色散参数,即τg1(λ)=τg2(λ),则从光脉冲进入第一单模光纤到出第二单模光纤的过程中,时域拉伸倍数为
从式<16>中可知,采用两段色散参数相同的光纤不会引起展宽系数的变化,拉伸倍数只与光纤的长度有关系。
所以,当采用两段色散参数相同的光纤时,拉伸倍数rst为:
由以上分析可知,调制信号经过第二单模光纤后被拉伸了rst倍,得到了二次拉伸的信号e′out。
根据上述原理,本步骤将调制信号eout通过第二单模光纤拉伸,得到了二次拉伸的信号e′out。
步骤4,用光电探测器将二次拉伸信号e′out转变为电信号,即光域拉伸信号x'(t):
x'(t)=cg2(t,τ2)x(t/m)<18>
其中c是一个与信号调制深度、链路损耗、光电探测器灵敏度有关的常数。
步骤5,建立光域拉伸信号x'(t)下的两个基函数。
选择n个采样点数,奈奎斯特采样频率fnyq,计算采样点的时间间隔
将原信号x'[t]经过采样后表示为x'[n]=a0cos(2πf0tn+θ0),其中a0,f0,θ0分别为信号的幅度,频率和相位,它们是未知的,需要对其估计。
引入正弦和余弦函数函数,定义估计信号为:
x[n]=a1cos(2πfktn)+a2sin(2πfktn)<19>
令u1=u1[n]=cos(2πfktn),u2=u2[n]=sin(2πfktn),
定义u1=u1[n]=cos(2πfktn),u2=u2[n]=sin(2πfktn)为两组正交基。
步骤6,用gold码设计观测矩阵φ,并根据该矩阵得到感知矩阵
6.1)用二进制线性反馈寄存器产生m序列,并对该m序列进行模2加操作,得到伪随机二进制序列,得到gold码序列;
6.2)用gold码序列的每一个元素ri构建对角矩阵:
r=diag(ri)<20>
6.3)将gold码序列采样信号x'[n]相乘实现随机混频,得到随机混频信号z(t);
6.4)将随机混频信号z(t)与低通滤波器冲激响应h[n]进行卷积,实现低通滤波器h,其中h为n×n的矩阵;
6.5)根据电域压缩倍数rc,计算采样矩阵d中的元素:
d=δ(i-j/rc),i=1,2...m,j=1,2,...n,<21>
其中m=n/rc;
6.6)根据6.2),6.4)和6.5)的结果,计算得到观测矩阵φ:
φ=dhr。<22>
6.7)将观测矩阵φ写成列向量的形式
并根据该列向量得到感知矩阵vfk,表示如下:
其中u1=u1[n],u2=u2[n],<φi,u1>表示为φi与u1的内积,<φi,u2>表示为φi与u2的内积。
步骤7,用最小二乘法求得基扩展系数
7.1)不考虑噪声的影响,利用观测矩阵φ对原采样信号x'[n]进行调制,得到观测向量y=φx'[n];
7.2)利用观测矩阵φ对估计信号x[n]进行调制,得到观测结果为y=φx[n];
7.3)定义
根据式<23>知
其中
步骤8,根据基扩展系数
8.1)根据奈奎斯特采样定理,估计信号的载频fk在(0,fnyq/2)内,故是在0和fnyq/2之间定义估计信号载频fk;
8.2)通过求解上式<24>可以得到fk对应的基扩展系数
8.3)根据8.1)和8.2)结果建立代价函数
步骤9,获得光域拉伸信号x'(t)信号参数频率
9.1)将在fk下计算出的wk按照从小到大顺序排列;
9.2)取wk为最大值下的fk即为估计的信号载频
9.3)根据式<25>得到此信号载频fk下的基扩展系数
本发明的效果通过以下仿真实验进一步说明:
1.仿真条件:
仿真参数如下:光脉冲功率为20dbm,频率为1550nm,马赫-曾德尔(mzm)调制在双边带调制模式,并且传输功率最大,设置第一段光纤长度为l1=1km,第二段光纤长度为l2=4km,即光域拉伸倍数为5倍,cs压缩倍数为rc=6,加性噪声20db。
2.仿真内容:
仿真1.用本发明方法对采样点数n=2048,频率为:f1=4ghz、f2=6ghz、f3=8ghz、f4=12ghz的多频,幅度为0.1的正弦信号进行频率估计仿真,得到的频谱分布估计结果如图2。其中图2(a)表示原始信号频率图,图2(b)表示时域拉伸信号频率图,图2(c)表示估计出的信号频率估计图,从图2中可以看出,本发明能够准确地估计出信号的频率。
仿真2.用本发明方法对采样点数n=2048,频率为:f1=4ghz、f2=6ghz、f3=8ghz、f4=12ghz的多频,幅度为0.1的正弦信号进行幅度估计仿真,得到的幅度分布估计结果如图3,其中四种曲线分别表示4ghz,6ghz,,8ghz,12ghz频率下的幅度误差,从图3可以看出,本发明估计信号的幅度效果一般。
仿真3.用本发明方法对采样点数为n=2048,频率分别为f1=15ghz、f2=10ghz、f3=6ghz,幅度均为0.1的2psk、2ask与正弦信号相叠加的混合信号进行频率估计仿真,得到的频谱分布估计,结果如图4。其中图4(a)表示包括原始信号频率图,图4(b)表示时域拉伸信号频率图,图4(c)表示估计出的信号频率估计图,从本发明可以看出,此方法能够准确地估计出信号的频率。
仿真4.用本发明方法对采样点数为n=2048,频率分别为f1=15ghz、f2=10ghz、f3=6ghz,幅度均为0.1的2psk、2ask与正弦信号相叠加的混合信号进行幅度估计仿真,得到的频谱分布估计结果如图5。其中三种曲线分别表示2psk,2ask,和正弦信号下的幅度误差,从图5可以看出,本发明能够较好地地估计出信号的幅度。