1.一种α非高斯噪声下chirp信号的参数估计方法,其特征在于包括下列步骤:步骤一:采集含有噪声信号含有加性噪声的chirp信号数学模型为:其中s(t)表示发射的chirp信号,其数学模型表示为:其中t表示时间,T表示时宽,f0表示初始频率,k表示调频率;
A表示chirp信号在传输过程中的衰减因子,n(t)表示α非高斯噪声;
步骤二:对信号进行限幅预处理限幅处理后的信号如式(3)所示:p是限幅值;
步骤三:量纲归一化
设信号X(t)的时域区间为 对应的频域区间为 将时域和频域都转换成量纲统一的域,引入一个量纲归一化因子S,其中T表示时宽,F表示带宽;
并定义量纲归一化坐标为
其中f表示频率;
新坐标系(t′,f′)实现了量纲归一化;
时域和频域2个区间都归一化为其中 △x为采样频率;
信号采样间隔变为
步骤四:FrFT
通过对含噪声chirp信号进行分数阶Fourier变换,可以估计出参数f0和k,具体算法如下:首先对含噪声chirp信号进行FrFT得到Xa(u),函数X(t)的a阶分数阶Fourier变换(FrFT)定义如下:a是分数阶阶次,u表示采样点,设 则核函数
其中, 为时频平面的旋转角度,δ(t)单位脉冲函数,n=1,2,...的正整数;
根据式(7)和式(8),FrFT的定义式改写为:采样型离散分数阶Fourier变换DFrFT可分解为以下步骤:步骤1:用chirp信号 与信号X(t)相乘,即步骤2:g(t)与chirp信号 做卷积,即步骤3:用chirp信号 与信号h(u)相乘,即然后搜索|Xa(u)|2最大值对应的 和旋转角通过式(13)可以估计出参数 和步骤五尺度变换,得到参数估计值 和和 就是α非高斯噪声下chirp信号的参数f0、k的估计值。
2.如权利要求1所述的一种α非高斯噪声下chirp信号的参数估计方法,其特征在于:所述步骤二中限幅值p是发射的chirp信号功率的1~1.5倍。