雷达细微特征分析(希尔伯特).docxVIP

雷达信号细微特征希尔伯特分析法姓名： 超哥 学号： 09001 完成时间：2012年11月7日一．实验目的1、掌握用希尔波特分析分析雷达信号细微特征法的分析方法；2、掌握雷达信号包络及瞬时频率的提取；3、了解用希尔波特分析分析雷达信号细微特征法的分析方法特点。二．实验原理雷达信号的上升沿时间、下降沿时间、脉宽是雷达信号具有个体特征的重要参数，精确测量出该参数对雷达的个体识别有着重要的意义。我们采用希尔伯特（Hilbert）变换法进行分析，提取分析雷达信号的脉冲包络、上升沿时间、下降沿时间、脉宽以及脉幅等时域特征和瞬时频率特征，并能据此识别信号脉内调制类型。对于因果实信号，它的傅里叶变换的实部和虚部以及幅频响应和相频响应之间存在希尔伯特变换关系。利用此关系，我们就能够构成原始的解析信号，并通过求解析信号的模值提取原始信号的包络特征，然后提取上升、下降沿的时间、脉宽和脉幅。同时，可以利用解析信号的相位特性来提取雷达信号的瞬时频率特征，并据此识别信号脉内调制类型。1.希尔伯特变化的定义以及原理希尔伯特变换定义为： 其中，是的希尔伯特变换，由上式可以看出，希尔伯特变换可以看作是信号通过冲击响应为的线性滤波器的结果。该滤波器的频率响应为： 可见即希尔伯特变换可以看作是幅频特性为1的全通滤波器。当信号通过变换器是，其负频率成分做+90度的相移，而正频率成分做—90度的相移，因此变换后的信号和是正交的。2.雷达信号的解析表示设某个雷达信号克表示为：其中，微信号包络，为信号的瞬时角度调制，为信号载频。的解析信号为： 它的频域如下：可见，解析信号的频谱不再包含负频分量，其正频分量幅度变化为的2倍。3.利用解析信号提取雷达信号的瞬时频率利用信号解析化理论，我们得到构造解析信号的方法：分别以作为实部，以其希尔伯特变换作为虚部来构成解析信号。解析信号与原始信号有相同的幅度和频率范围，并包含了原始信号中的相位信息。因此，可以通过计算的模值来获得原始信号的包络，通过解析信号的相位可以获取原始信号的瞬时频率信息。解析信号的时域表达式为： 其模值为：即为原始信号的幅度调制率，由此我们可以提取其时域特征（如上升沿、下降沿时间、幅度、脉宽等）。的相位为：的瞬时频率为: 由上式我们可以提取雷达信号的瞬时频率特征。三．实验条件实验参数设定值线性调频信号点数t_f（画图）1000采样频率fs20000线性调频信号载频fc500+0.5t_f01随机序列长度bit_len10码元传播速率RB10调相信号点数t_P（画图）1000采样频率Fs1000调相信号载频Fc30四．实验内容（一）线性调频信号1.模拟产生线性调频信号 首先，要产生线性调频信号，载频要随时间变化而变化，在实验绘图中，采样点共为1000个，令线性调频信号载频fc=500+0.5t_f。 然后，画一个1000点长度的具有衰减特点的人为包络m，在前100点让其以斜率为1上升，后100点以斜率为-1下降。为了使模拟信号更真实，上升下降沿用包络衰减来表示，对生成包络在加入高斯白噪声，在实验中取信噪比为40。代码如下： 最后，把线性变换的载频带入到余弦函数中，得到线性调频信号，在乘以加入噪声的包络得到模拟真实的线性调频信号。产生流程见下图2.对信号做希尔伯特变换并构造解析信号 因为希尔伯特变换可以看作是信号通过冲击响应为的线性滤波器的结果。该滤波器的频率响应为： 由于为让后面点乘的点数一样，转为频域时都做1000点dft。是1000点，相当于从-500点到500点，由于没有负频，所以要将前面的搬移到后面，所以0到500点取-i，500到1000点取+i。写时如下：3.对解析信号求模值得原信号的时域特征对解析信号直接abs得出它的包络，与原模拟真实信号包络做对比，比较结果还是很理想的，甚至从解析信号得出的包络比原信号的包络更平滑。对比图如下：4.由解析信号相位提取雷达信号的瞬时频率特征用angle函数对解析信号运算，得出解析信号的相位，由于angle函数得出的相位都在-pi到pi之间，所以相位的图像出现下图这种锯齿状，但是由图可看出来，随着点数的向后，相位增长的斜率也越来越大，这也说明了载频在增加。对angle后的解析信号相位求微分，由于相位有向下的突然跳变，所以，求导时会有向下的冲激，但是从瞬时频率图看来，除了头尾有波动外，还是很标准的线性增长，这也验证的书上的利用希尔伯特变换求瞬时频率的原理。若相位画图时没有跳变，瞬时频率的图将会是更标准的增长的直线。代码如下：（二）调相信号1.模拟产生调相信号，步骤同线性调频信号，过程如下：2.对信号做希尔伯特变换并构造解析信号 原理步骤同线性调频信号。3.对解析信号求模值得原信号的时域特征，可以看出来调相信号的效果没有线性调频信号好4.由解析信号相位提取雷达信号的瞬