关于Matlab的IQ信号处理.pdfVIP

一、实习课题的具体内容 1、ASK、FSK 解调部分 (1)、从Tektronix 公司的RSA3303A 型号频谱仪中导出频谱仪I/Q 数据； (2)、在matlab 上编程实现基于I/Q 数据的多帧ASK 和FSK 解调； (3)、要求：Matlab 解调结果要和频谱仪所示吻合，并实现解调自动化。 二、实习课题完成的技术路线或实施方法 2.1、ASK、FSK 解调部分 基于Matlab 的I/Q 解调分析有两种方法，可以在计算机与监测设备联机的 情况下进行实时采集与分析，也可以利用设备提供的I/Q 解调功能采集数据，保 存在后缀名为“.mat”的文件中进行脱机分析处理。我采用的分析方法是后者， 利用 Tektronix 公司的 RSA 系列频谱分析仪采集数据，然后保存到后缀名为 “.mat”的文件中，再通过U 盘拷贝的方式导入PC 端，然后在Matlab 上分析处 理。 2.1.1、IQ 数据简介 极坐标参量与直角坐标参量间的转换示意如下图，极坐标图表示的矢量的参 数是幅值和瞬时相位，一般表示为A ，因为所有调制形式针对载波的改变都 是幅值、相位以及可以由相位推导的频率，极坐标图是分析调制技术矢量参数的 主要描述方式。I/Q 直角坐标图表示矢量的参数不是直接的幅值和瞬时相位，而 是把它们投影在I/Q 直角坐标轴上，采用I 轴和Q 轴的投影分量来确定矢量。 2.1.2、数据导入 我利用Tektronix频谱仪的I/Q解调功能采集射频数据，然后用MATLAB代码 “LOAD EXAMPLE.MAT”将数据导入。下图是利用RSA3303A采集数据后存储为 “.mat”文件的内容。 1 通过前期对泰克频谱仪用户手册和文献资料的查阅，我了解到频谱仪导出数 据的含义基本如下： 其中，Y变量是一个复变量，包含了采样数据点的同相I分量和正交Q分量。 InputCenter变量显示了采样时的中心频率，据此可计算出实际频率分量。 InputZoom变量值为1，表示采集数据已下变频到基带。XDelta为采样周期，由此 可以计算出采样率（fs=1/Xdelta）。 2.1.3、数据分析 根据“.mat”文件中的每个采样时刻的I/Q采样数据Y，可得到相应的I分量 和Q分量，用MATLAB代码表示为：I=real(Y),Q=imag(Y)。采集到的数据可用幅度 分析、频率分析、相位分析以及小波分析等方法进行处理。 输入信号的幅度 分析： 2 2 。根据I/Q信号的每一个采样值， Am Am I Q 可以计算出对应采样时刻的幅度值。在测量时间内的所有幅度数值构成I/Q幅度 数组，对应于“射频功率－时间曲线”。 输入信号的相对相位 分析： 。根据I/Q信号的每一个 m m arctan(Q/ I) 采样值，可以计算出对应采样时刻的相对相位值。在测量时间内的所有 数值 m 构成I/Q相对相位数组 。 m 输入信号的频率 计算： d (m) 。在测量时间内的所有 Fm数值构 Fm Fm Fm dt 成I/Q频率数组 ，对应于频率曲线。 Fm 频谱计算：Spectrum = FFT (Y)。对I/Q数据Y进行傅里叶变换，即可计算出 其频谱值。其中，Spectrum为输入信号中频带宽内的频谱。 2.1.4、解调实例分析 （1）FSK解调功能的实现 FSK 就是用数字信号去调制载波的频率，FSK 解调一般采用两种方式：非相 干解调和相干解调。相干解调需要同步信号，适用于窄带调频信号解调；而非相 干解调不需要同步信号，对于窄频带调频（NBFM）信号和宽频带调频（WBFM）信 号均适用