无线电信号分析技术

无线电信号分析技术现在是1页\一共有51页\编辑于星期一2内容总览第一部分：数字无线电监测系统第二部分：无线电信号频谱分析第三部分：无线电信号时频分析第四部分：无线电信号参数分析第五部分：信号分离技术现在是2页\一共有51页\编辑于星期一3第一部分：数字无线电监测系统

——下一代无线电监测系统

一、概述

二、数字无线电监测系统组成

三、数字信号采样四、数字接收机关键指标五、I、Q输出现在是3页\一共有51页\编辑于星期一4一、概述

模拟无线电监测系统的演进；数字信号处理技术发展的结晶；采集和存储信号的数字采样信息；通过数字信号处理技术实现信号再生、复原以及参数测量。现在是4页\一共有51页\编辑于星期一5二、数字无线电监测系统组成现在是5页\一共有51页\编辑于星期一6三、数字信号采样

采样方式：1、基本采样采样原则：采样数值能够准确地确定原信号。适用范围：基带信号采样速率：fs≥2fmax

，fmax为信号最高频率

2、带通采样适用范围：中频信号采样速率：

n为满足上述条件的最大正整数。现在是6页\一共有51页\编辑于星期一7四、数字接收机关键指标1、A/D位数确定A/D转换灵敏度、信噪比

n为A/D位数，B为模拟信号带宽，fs为采样信号频率

2、采样速率fs决定数字接收机可以处理的模拟信号带宽

现在是7页\一共有51页\编辑于星期一8四、数字接收机关键指标(续)3、接收机动态范围其中，IP3为放大器三阶截点的输出功率；

G为放大器增益；

NF为噪声系数（含放大器、下变频器和A/D噪声系数）；

BF为为滤波器带宽。从中看出：IP3是关键指标

现在是8页\一共有51页\编辑于星期一9五、I、Q输出许多数字接收机提供I和Q两路输出。数学含义：其中：不含射频信息，但含有基带信息。称为信号的同相分量和正交分量。现在是9页\一共有51页\编辑于星期一10第二部分：无线电信号频谱分析

一、概述

二、离散傅立叶变换（DFT）

三、常用算法：FFT现在是10页\一共有51页\编辑于星期一11一、概述

频谱管理是无线电管理工作的重点频谱分析是无线电管理的工具频谱分析仪是常用的频谱分析设备频谱分析工具——离散傅立叶变换现在是11页\一共有51页\编辑于星期一12二、离散傅立叶变换（DFT）采集得到数字信号x(n),n=0,1,2,3,…….DFT:,k=0,1,2,3…N-1;

其中X(k)为复数，含有幅度相位信息。IDFT:n=0,1,2,3…,N-1;现在是12页\一共有51页\编辑于星期一13三、加窗问题计算离散时间信号需要用到信号在全体时间上的值，在实际应用中我们仅能在有限时宽内观测信号，相当于对信号进行加窗。加窗的结果导致信号能量的泄漏。由于泄漏，加窗不仅扭曲了谱估计，也降低了谱的分辨率。现在是13页\一共有51页\编辑于星期一14三、加窗问题（续1）令{X(n)}为被分析信号，如果序列时宽为区间0≤n≤l-1内的l个样本，这相当于用长度为l的矩形窗乘，即：

令

为窗函数序列的傅立叶变换则W(ω)的傅立叶变换频率现在是14页\一共有51页\编辑于星期一15三、加窗问题（续2）还可以给信号加其它的窗加窗遵循规律带宽准确性高：主瓣宽度越大，旁瓣数目越少。分辨率准确性高：主瓣宽度越小，旁瓣数目越多。现在是15页\一共有51页\编辑于星期一16三、常用算法：FFTFFT

是DFT的快速计算方法

Matlab

语句：fft(x(n),N)

x(n)为数字采样信号；

N为离散傅里叶变换数据x(n)的长度。特点：使用简单、经济方便现在是16页\一共有51页\编辑于星期一17FFT的关键指标1、频率分辨率：,fs为信号采样速率2、实际频率：,n为FFT输出序列3、对称性：输出数据关于n/2对称，相位相反。因此n/2以上数据无意义。N(512)点fft输出系列n/2以上数据无意义现在是17页\一共有51页\编辑于星期一18第三部分：无线电信号时频分析

一、目的

二、短时傅立叶变换(STFT)

三、应用现在是18页\一共有51页\编辑于星期一19一、时频分析目的频谱分析缺点：不能反映信号频谱随时间的变化；时频分析目的：分析信号频率随时间的变换规律，为信号分析提供更多的信息。现在是19页\一共有51页\编辑于星期一20二、短时傅立叶变换STFT数学表达式：

n为数据长度，其取值大于N。类似于离散傅立叶变换，是信号序列的DFT。Matlab语句：tfrstft(x,t,N,h)

其中：t时刻，N为长度，h为归一化的频率平滑窗。现在是20页\一共有51页\编辑于星期一21三、时频分析应用

对于ASK和FSK信号调制方式识别非常容易。举例一：ASK和FSK信号。举例二：大功率无绳电话开机密码解码。现在是21页\一共有51页\编辑于星期一22举例一：2FSK（实际信号）现在是22页\一共有51页\编辑于星期一23举例二：

大功率无绳电话开机密码解码现在是23页\一共有51页\编辑于星期一24第四部分：无线电信号参数分析

一、无线电信号二、信号参数分析目的三、希尔伯特变换四、工程实现方法五、应用现在是24页\一共有51页\编辑于星期一25一、无线电信号无线电信号表达式：三个参数瞬时幅度：a(t)瞬时相位：瞬时频率：现在是25页\一共有51页\编辑于星期一26二、信号参数分析目的手段：获取瞬时幅度、瞬时频率、瞬时相位；从而全面掌握信号特征，识别信号类型；完成信号参数测量和识别解调等任务。现在是26页\一共有51页\编辑于星期一27三、希尔伯特变换1、解析信号实信号：解析信号：相位：虚部：现在是27页\一共有51页\编辑于星期一28问题：为什么使用解析信号？回答：解析信号易于获取信号的三个参数，为信号分析提供手段。现在是28页\一共有51页\编辑于星期一292、希尔伯特变换实信号：希尔伯特变换：物理解释：对信号进行一次正交相移，从而产生正交信号。正交含义：现在是29页\一共有51页\编辑于星期一30频谱解释实信号：——共轭对称性物理含义：正负频率分量对称，相位相反。解析信号

实信号的频谱正、负频率对称解析信号的频谱f<0,频谱为0f>0,频谱为2U(f)现在是30页\一共有51页\编辑于星期一313、参数计算三个特征参数：瞬时幅度瞬时相位瞬时频率

现在是31页\一共有51页\编辑于星期一32四、工程实现方法

信号采样——模拟信号数字化

u(n),n=0,1,2,…

数字信号处理计算步骤：（1）（2）现在是32页\一共有51页\编辑于星期一33四、工程实现方法（续）数字信号处理计算步骤（续）：（3）（4）现在是33页\一共有51页\编辑于星期一34五、应用1、调制方式识别

2、参数测量 数字接收机I/Q输出即为正交信号，可以从I/Q输出完成信号参数的测量。调制方式瞬时幅度瞬时相位特征瞬时频率特征调幅随时间变化以2/f0为周期变化固定值调频固定值随时间变化随时间变化调相固定值随时间变化固定值现在是34页\一共有51页\编辑于星期一35第五部分：信号分离技术

一、循环谱技术二、独立分量（ICA）技术

现在是35页\一共有51页\编辑于星期一36一、循环谱技术

1、循环谱相关函数

如果信号的自相关函数具有周期特性，则可以通过Fourier变换求得信号的循环谱相关函数，即：

（1）

这里：[-T,T]为感兴趣的时间。

现在是36页\一共有51页\编辑于星期一37一、循环谱技术（续）2、循环谱密度函数对循环谱相关函数求Fourier变换即可得到循环谱密度函数，即（2）

当=0时，（1）和（2）退化为通常的自相关函数和谱密度函数，可见循环谱相关函数和循环谱密度函数是传统自相关函数和功率谱密度函数的推广。由此我们可以判断从其周期频率方面可以提供更多的信息量。现在是37页\一共有51页\编辑于星期一38一、循环谱技术（续）3、谱相关特点：1）具有相同功率谱密度（如BPSK、QPSK）而调制类型不同的信号具有不同的谱相关函数；2）平稳噪声和干扰没有谱相关特性，或者说谱相关函数值恒等于0；3）谱相关函数包含与调制信号时间参数相关联的相位和频率信息；4）信号中谱相关分量相互关联，可以利用信号中某些谱元素分量估计其它分量。现在是38页\一共有51页\编辑于星期一39一、循环谱技术（续）4、应用：信号分离微弱信号提取TDOA定位调制方式识别和信号参数提取现在是39页\一共有51页\编辑于星期一40一、循环谱技术（续）5、举例:BPSK＋AM信号分离BPSK信号循环谱现在是40页\一共有51页\编辑于星期一41AM信号循环谱现在是41页\一共有51页\编辑于星期一42BPSK＋AM信号的循环谱现在是42页\一共有51页\编辑于星期一43普通功率谱现在是43页\一共有51页\编辑于星期一44BPSK与BPSK＋AM（1/Tc）循环谱对比现在是44页\一共有51页\编辑于星期一45二、独立分量（ICA）技术问题：目的：ICA数学模型现在是45页\一共有51页\编辑于星期一46目标：寻找解混矩阵W使输出信号之间互相独立

手段：ICA算法=目标函数+优化算法目标函数评价：

非高斯测度信息论优化算法：梯度下降法等ICA算法现在是46页\一共有51页\编辑于星期一47非高斯测度：中心极限定理，一般情况下，多个相互独立的随机变量的线性混合比其中任何一个随机变量都更接近高斯分布。所以，可以把非高斯性作为代价函数，通过非高斯性的最大化来达到提取独立分量的目的。信息论：利用输入或输出之间的信息关系，比如互信息最小化是利用平均互信息来衡量输出分量之间的独立性，当输出分量之间相互独立时，他们的互信息为零。现在是47页\一共有51页\编辑于星期一48ICA假设条件源信号为平稳随机过程，且相互独立。混合矩阵是列