无线电信号分析技术课件

无线电信号分析技术无线电信号分析技术内容总览第一部分：数字无线电监测系统第二部分：无线电信号频谱分析第三部分：无线电信号时频分析第四部分：无线电信号参数分析第五部分：信号分离技术内容总览第一部分：数字无线电监测系统第一部分：数字无线电监测系统

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

一、概述

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

三、数字信号采样四、数字接收机关键指标五、I、Q输出第一部分：数字无线电监测系统

——下一代无线电监测系统一、概述

模拟无线电监测系统的演进；数字信号处理技术发展的结晶；采集和存储信号的数字采样信息；通过数字信号处理技术实现信号再生、复原以及参数测量。一、概述模拟无线电监测系统的演进；二、数字无线电监测系统组成二、数字无线电监测系统组成三、数字信号采样

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

，fmax为信号最高频率

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

n为满足上述条件的最大正整数。三、数字信号采样采样方式：四、数字接收机关键指标1、A/D位数确定A/D转换灵敏度、信噪比

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

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

四、数字接收机关键指标1、A/D位数四、数字接收机关键指标(续)3、接收机动态范围其中，IP3为放大器三阶截点的输出功率；

G为放大器增益；

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

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

四、数字接收机关键指标(续)3、接收机动态范围五、I、Q输出许多数字接收机提供I和Q两路输出。数学含义：其中：不含射频信息，但含有基带信息。称为信号的同相分量和正交分量。五、I、Q输出许多数字接收机提供I和Q两路输出。第二部分：无线电信号频谱分析

一、概述

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

三、常用算法：FFT第二部分：无线电信号频谱分析一、概述

频谱管理是无线电管理工作的重点频谱分析是无线电管理的工具频谱分析仪是常用的频谱分析设备频谱分析工具——离散傅立叶变换一、概述频谱管理是无线电管理工作的重点二、离散傅立叶变换（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;二、离散傅立叶变换（DFT）采集得到数字信号x(n),三、加窗问题计算离散时间信号需要用到信号在全体时间上的值，在实际应用中我们仅能在有限时宽内观测信号，相当于对信号进行加窗。加窗的结果导致信号能量的泄漏。由于泄漏，加窗不仅扭曲了谱估计，也降低了谱的分辨率。三、加窗问题计算离散时间信号需要用到信号在全体时间上的值，在三、加窗问题（续1）令{X(n)}为被分析信号，如果序列时宽为区间0≤n≤l-1内的l个样本，这相当于用长度为l的矩形窗乘，即：

令

为窗函数序列的傅立叶变换则W(ω)的傅立叶变换频率三、加窗问题（续1）令{X(n)}为被分析信号，如果序列时宽三、加窗问题（续2）还可以给信号加其它的窗加窗遵循规律带宽准确性高：主瓣宽度越大，旁瓣数目越少。分辨率准确性高：主瓣宽度越小，旁瓣数目越多。三、加窗问题（续2）还可以给信号加其它的窗三、常用算法：FFTFFT

是DFT的快速计算方法

Matlab

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

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

N为离散傅里叶变换数据x(n)的长度。特点：使用简单、经济方便三、常用算法：FFTFFT是DFT的快速计算方法FFT的关键指标1、频率分辨率：,fs为信号采样速率2、实际频率：,n为FFT输出序列3、对称性：输出数据关于n/2对称，相位相反。因此n/2以上数据无意义。N(512)点fft输出系列n/2以上数据无意义FFT的关键指标1、频率分辨率：,fs为信号采第三部分：无线电信号时频分析

一、目的

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

三、应用第三部分：无线电信号时频分析一、目的一、时频分析目的频谱分析缺点：不能反映信号频谱随时间的变化；时频分析目的：分析信号频率随时间的变换规律，为信号分析提供更多的信息。一、时频分析目的频谱分析缺点：不能反映信号频谱二、短时傅立叶变换STFT数学表达式：

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

其中：t时刻，N为长度，h为归一化的频率平滑窗。二、短时傅立叶变换STFT数学表达式：三、时频分析应用

对于ASK和FSK信号调制方式识别非常容易。举例一：ASK和FSK信号。举例二：大功率无绳电话开机密码解码。三、时频分析应用对于ASK和FSK信号调制方式识别非常容易举例一：2FSK（实际信号）举例一：2FSK（实际信号）举例二：

大功率无绳电话开机密码解码举例二：

大功率无绳电话开机密码解码第四部分：无线电信号参数分析

一、无线电信号二、信号参数分析目的三、希尔伯特变换四、工程实现方法五、应用第四部分：无线电信号参数分析一、无线电信号一、无线电信号无线电信号表达式：三个参数瞬时幅度：a(t)瞬时相位：瞬时频率：一、无线电信号无线电信号表达式：二、信号参数分析目的手段：获取瞬时幅度、瞬时频率、瞬时相位；从而全面掌握信号特征，识别信号类型；完成信号参数测量和识别解调等任务。二、信号参数分析目的手段：获取瞬时幅度、瞬时频率、瞬时相位第四部分：无线电信号参数分析令{X(n)}为被分析信号，如果序列时宽为区间0≤n≤l-1内的l个样本，这相当于用长度为l的矩形窗乘，即：采集和存储信号的数字采样信息；2）平稳噪声和干扰没有谱相关特性，或者说谱相关函数值恒等于0；5、举例:BPSK＋AM信号分离2）平稳噪声和干扰没有谱相关特性，或者说谱相关函数值恒等于0；加窗的结果导致信号能量的泄漏。第四部分：无线电信号参数分析计算离散时间信号需要用到信号在全体时间上的值，在实际应用中我们仅能在有限时宽内观测信号，相当于对信号进行加窗。二、短时傅立叶变换STFT计算步骤（续）：物理含义：正负频率分量对三、希尔伯特变换1、解析信号实信号：解析信号：相位：虚部：第四部分：无线电信号参数分析三、希尔伯特变换1、解析信号问题：为什么使用解析信号？回答：解析信号易于获取信号的三个参数，为信号分析提供手段。问题：为什么使用解析信号？回答：2、希尔伯特变换实信号：希尔伯特变换：物理解释：对信号进行一次正交相移，从而产生正交信号。正交含义：2、希尔伯特变换实信号：频谱解释实信号：——共轭对称性物理含义：正负频率分量对称，相位相反。解析信号

实信号的频谱正、负频率对称解析信号的频谱f<0,频谱为0f>0,频谱为2U(f)频谱解释实信号：实信号的频谱正、负频率对称解析信号的频谱f3、参数计算三个特征参数：瞬时幅度瞬时相位瞬时频率

3、参数计算三个特征参数：四、工程实现方法

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

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

数字信号处理计算步骤：（1）（2）四、工程实现方法信号采样——模拟信号数字化四、工程实现方法（续）数字信号处理计算步骤（续）：（3）（4）四、工程实现方法（续）数字信号处理五、应用1、调制方式识别

2、参数测量 数字接收机I/Q输出即为正交信号，可以从I/Q输出完成信号参数的测量。调制方式瞬时幅度瞬时相位特征瞬时频率特征调幅随时间变化以2/f0为周期变化固定值调频固定值随时间变化随时间变化调相固定值随时间变化固定值五、应用1、调制方式识别调制方式瞬时幅度瞬时相位特征第五部分：信号分离技术

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

第五部分：信号分离技术一、循环谱技术一、循环谱技术

1、循环谱相关函数

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

（1）

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

一、循环谱技术1、循环谱相关函数一、循环谱技术（续）2、循环谱密度函数对循环谱相关函数求Fourier变换即可得到循环谱密度函数，即（2）

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

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

非高斯测度信息论优化算法：梯度下降法等ICA算法目标：寻找解混矩阵W使输出信号之间互相独立IC非高斯测度：中心极限定理，一般情况下，多个相互独立的随机变量的线性混合比其中任何一个随机变量都更接近高斯分布。所以，可以把非高斯性作为代价函数，通过非高斯性的最大化来达到提取独立分量的目的。信息论：利用输入或输出之间的信息关系，比如互信息最小化是利用平均互信息来衡量输出分量之间的独立性，当输出分量之间相互独立时，他们的互信息为零。非高斯测度：中心极限定理，一般情况下，多个相互独立的随机ICA假设条件源信号为平稳随机过程，且相互独立。混合矩阵是列满秩的源信号的混合是无噪的最多只允许一个源信号服从高斯分布源信号的数目不大于传感器的数目ICA假设条件源信号为平稳随机过程，且相互独立。FastICA方法基于非高斯测度的经典算法目标函数：峭度函数优化算法：梯度下降法流程：白化=〉初始化=〉权值迭代FastICA方法基于非高斯测度的经典算法实验仿真及结果实验仿真及结果与源信号对比r1=-0.9999r2=-0.9999r3=-1.0000与源信号对比r1=-0.9999r2=-0.9999无线电信号分析技术无线电信号分析技术内容总览第一部分：数字无线电监测系统第二部分：无线电信号频谱分析第三部分：无线电信号时频分析第四部分：无线电信号参数分析第五部分：信号分离技术内容总览第一部分：数字无线电监测系统第一部分：数字无线电监测系统

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

一、概述

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

三、数字信号采样四、数字接收机关键指标五、I、Q输出第一部分：数字无线电监测系统

——下一代无线电监测系统一、概述

模拟无线电监测系统的演进；数字信号处理技术发展的结晶；采集和存储信号的数字采样信息；通过数字信号处理技术实现信号再生、复原以及参数测量。一、概述模拟无线电监测系统的演进；二、数字无线电监测系统组成二、数字无线电监测系统组成三、数字信号采样

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

，fmax为信号最高频率

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

n为满足上述条件的最大正整数。三、数字信号采样采样方式：四、数字接收机关键指标1、A/D位数确定A/D转换灵敏度、信噪比

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

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

四、数字接收机关键指标1、A/D位数四、数字接收机关键指标(续)3、接收机动态范围其中，IP3为放大器三阶截点的输出功率；

G为放大器增益；

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

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

四、数字接收机关键指标(续)3、接收机动态范围五、I、Q输出许多数字接收机提供I和Q两路输出。数学含义：其中：不含射频信息，但含有基带信息。称为信号的同相分量和正交分量。五、I、Q输出许多数字接收机提供I和Q两路输出。第二部分：无线电信号频谱分析

一、概述

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

三、常用算法：FFT第二部分：无线电信号频谱分析一、概述

频谱管理是无线电管理工作的重点频谱分析是无线电管理的工具频谱分析仪是常用的频谱分析设备频谱分析工具——离散傅立叶变换一、概述频谱管理是无线电管理工作的重点二、离散傅立叶变换（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;二、离散傅立叶变换（DFT）采集得到数字信号x(n),三、加窗问题计算离散时间信号需要用到信号在全体时间上的值，在实际应用中我们仅能在有限时宽内观测信号，相当于对信号进行加窗。加窗的结果导致信号能量的泄漏。由于泄漏，加窗不仅扭曲了谱估计，也降低了谱的分辨率。三、加窗问题计算离散时间信号需要用到信号在全体时间上的值，在三、加窗问题（续1）令{X(n)}为被分析信号，如果序列时宽为区间0≤n≤l-1内的l个样本，这相当于用长度为l的矩形窗乘，即：

令

为窗函数序列的傅立叶变换则W(ω)的傅立叶变换频率三、加窗问题（续1）令{X(n)}为被分析信号，如果序列时宽三、加窗问题（续2）还可以给信号加其它的窗加窗遵循规律带宽准确性高：主瓣宽度越大，旁瓣数目越少。分辨率准确性高：主瓣宽度越小，旁瓣数目越多。三、加窗问题（续2）还可以给信号加其它的窗三、常用算法：FFTFFT

是DFT的快速计算方法

Matlab

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

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

N为离散傅里叶变换数据x(n)的长度。特点：使用简单、经济方便三、常用算法：FFTFFT是DFT的快速计算方法FFT的关键指标1、频率分辨率：,fs为信号采样速率2、实际频率：,n为FFT输出序列3、对称性：输出数据关于n/2对称，相位相反。因此n/2以上数据无意义。N(512)点fft输出系列n/2以上数据无意义FFT的关键指标1、频率分辨率：,fs为信号采第三部分：无线电信号时频分析

一、目的

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

三、应用第三部分：无线电信号时频分析一、目的一、时频分析目的频谱分析缺点：不能反映信号频谱随时间的变化；时频分析目的：分析信号频率随时间的变换规律，为信号分析提供更多的信息。一、时频分析目的频谱分析缺点：不能反映信号频谱二、短时傅立叶变换STFT数学表达式：

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

其中：t时刻，N为长度，h为归一化的频率平滑窗。二、短时傅立叶变换STFT数学表达式：三、时频分析应用

对于ASK和FSK信号调制方式识别非常容易。举例一：ASK和FSK信号。举例二：大功率无绳电话开机密码解码。三、时频分析应用对于ASK和FSK信号调制方式识别非常容易举例一：2FSK（实际信号）举例一：2FSK（实际信号）举例二：

大功率无绳电话开机密码解码举例二：

大功率无绳电话开机密码解码第四部分：无线电信号参数分析

一、无线电信号二、信号参数分析目的三、希尔伯特变换四、工程实现方法五、应用第四部分：无线电信号参数分析一、无线电信号一、无线电信号无线电信号表达式：三个参数瞬时幅度：a(t)瞬时相位：瞬时频率：一、无线电信号无线电信号表达式：二、信号参数分析目的手段：获取瞬时幅度、瞬时频率、瞬时相位；从而全面掌握信号特征，识别信号类型；完成信号参数测量和识别解调等任务。二、信号参数分析目的手段：获取瞬时幅度、瞬时频率、瞬时相位第四部分：无线电信号参数分析令{X(n)}为被分析信号，如果序列时宽为区间0≤n≤l-1内的l个样本，这相当于用长度为l的矩形窗乘，即：采集和存储信号的数字采样信息；2）平稳噪声和干扰没有谱相关特性，或者说谱相关函数值恒等于0；5、举例:BPSK＋AM信号分离2）平稳噪声和干扰没有谱相关特性，或者说谱相关函数值恒等于0；加窗的结果导致信号能量的泄漏。第四部分：无线电信号参数分析计算离散时间信号需要用到信号在全体时间上的值，在实际应用中我们仅能在有限时宽内观测信号，相当于对信号进行加窗。二、短时傅立叶变换STFT计算步骤（续）：物理含义：正负频率分量对三、希尔伯特变换1、解析信号实信号：解析信号：相位：虚部：第四部分：无线电信号参数分析三、希尔伯特变换1、解析信号问题：为什么使用解析信号？回答：解析信号易于获取信号的三个参数，为信号分析提供手段。问题：为什么使用解析信号？回答：2、希尔伯特变换实信号：希尔伯特变换：物理解释：对信号进行一次正交相移，从而产生正交信号。正交含义：2、希尔伯特变换实信号：频谱解释实信号：——共轭对称性物理含义：正负频率分量对称，相位相反。解析信号

实信号的频谱正、负频率对称解析信号的频谱f<0,频谱为0f>0,频谱为2U(f)频谱解释实信号：实信号的频谱正、负频率对称解析信号的频谱f3、参数计算三个特征参数：瞬时幅度瞬时相位瞬时频率

3、参数计算三个特征参数：四、工程实现方法

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

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

数字信号处理计算步骤：（1）（2）四、工程实现方法信号采样——模拟信号数字化四、工程实现方法（续）数字信号处理计算步骤（续）：（3）（4）四、工程实现方法（续）数字信号处理五、应用1、调制方式识别

2、参数测量 数字接收机I/Q输出即为正交信号，可以从I/Q输出完成信号参数的测量。调制方式瞬时幅度瞬时相位特征瞬时频率特征调幅随时间变化以2/f0为周期变化固定值调频固定值随时间变化随时间变化调相固定值随时间变化固定值五、应用1、调制方式识别调制方式瞬时幅度瞬时相位特征第五部分：信号分离技术

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

第五部分：信号分离技术一、循环谱技术一、循环谱技术

1、循环谱相关函数

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

（1）

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

一、循环谱技术1、循环谱相关函数一、循环谱技术（续）2、循环谱密度函数对循环谱相关函数求Fourier变换即可得到循环谱密度函数，即（2）

当=0时，（1）和（2）退化为通常的自相关函数和谱密度函数，可见循环谱相关函数和循环谱密度函数是传统自相关函数和功率谱密度函数的推广。由此我们可以判断从其周期频率方面可以提供更多的信息量。一、循环谱技术（续）2、循环谱密度函数一、循环谱技术（续）3、