DSP怎么实现调幅广播信号的监测系统？

　　随着通信与广播电视业务的发展，无线电频谱迅速、大量的被占用，频道拥挤和相互间干扰日趋严重，为了能有效地利用无线电频谱，减少相互间的干扰，信号监测业务随之成为必要。调幅广播信号监测系统是用于实时监测短波调幅广播信号的调幅度、载波频率的专用系统。
　　图1为调幅广播信号质量监测系统的系统框图。本系统由数据采集模块、总线控制模块、数据处理模块、上位机通信模块组成。其工作方式为：输入信号通过线性数控增益放大器后由A/D转换器采样，采样后的数据由FPGA送入DSP进行数据处理，所得到的监测结果由FPGA通过PCI接口送入上位机。同时由DSP对采样所得信号大小进行监测，通过FPGA对线性数控增益放大器的增益进行调整，使其输出信号满足系统测量要求。
　　
　　图1：监测系统框图
　　算法和数字处理软件
　　调幅广播信号的特点是载波频率相对稳定而调幅度实时变化，所以系统每监测一百次调幅度再监测一次载波频率同样可以达到监测信号的目的。系统所要监测的信号的带宽小于10kHz，频率范围为1.5MHz~30MHz。若对信号进行过采样，为保证测量精度则采样频率要达到240MHz，后续数据处理难度较大，所以系统选择对信号进行欠采样。采样点数为N=4096。
　　调幅信号调幅度的计算公式如下：
　　
　　计算公式
　　其中m为调幅度，A为包络上、下两正峰间的幅度，B为包络上、下两负峰间的幅度。
　　本系统中先用频域移相法求得A/D采样后信号的Hilbert变换，然后取绝对值即可得到AM信号的上包络。因为单音调幅信号的上、下包络关于零轴对称，所以只要得到上、下包络其中之一就能根据式（1）求得调幅度。
　　设输入信号为f（t），则由采样定理可知采样前信号的频谱F（w）和采样后信号的频谱FS（w）有如下关系：
　　
　　计算公式
　　输入信号最高频率为30MHz，系统选取最小公倍数为60MHz的两个采样频率分别为fs1、fs2（fs1=60/264MHz=0.228MHz、fs2=60/261MHz=0.230MHz）对信号进行采样，对采样后的信号作FFT变换，与fs1和fs2相对应的频率分别为f1=fs1×n1/N和f2=fs2×n2/N。由此可粗略的得到信号的载波频率为：
　　
　　计算公式
　　进而为了达到1Hz的测量精度，必须进行频谱细化。若要达到0.5Hz分辨率，则分析200Hz谱宽需要计算400点。在粗估计频率的前后100Hz内进行频率细化分析，分别对f1和f2作n1和n2周围400点的DFT变换。变换后得到频谱F1和F2，分别求出其幅度最大值对应的位置nF1和nF2，从而得到精估计的信号频率值和：
　　
　　计算公式
　　根据上述算法可得到DSP中数据处理软件的流程图（图2）。
　　
　　DSP中数据处理软件的流程图

　　FPGA逻辑设计
　　本系统中FPGA主要用来协调各个模块间的数据传输，分别为A/D采样数据到DSP的传输、DSP计算结果到PCI接口的传输以及数控增益放大器的增益控制。同时FPGA还为系统工作提供了必要的时钟、复位信号、控制信号（图3）。
　　
　　FPGA逻辑图
　　器件选择
　　A/D转换器是整个监测系统的关键部件，它的性能往往直接影响整个监测系统的技术指标。当A/D有效位数大于12位时量化损失为0.0055dB，其对测量精度的影响可忽略不计。系统选用的A/D转换器为ADI公司的AD9433。输入AD9433的信号幅度要控制在一定的范围内，否则会造成失真，甚至烧毁芯片，所以要在AD9433之前用运放对信号幅度进行调控。同时根据调幅广播信号幅度实时变化的特点，要求所选择的运放增益可变。基于上述要求系统选用ADI公司的线性数控增益放大器AD8320。
　　系统对信号采样点数为N=4096，算法采用Hilbert变换解调求调幅度和欠采样求载波频率，所以每计算100次调幅度和1次载波频率所需要的运算量大概为：
　　
　　计算公式
　　本系统选用ADI公司SHARC系列的ADSP-21262作为数据处理芯片。
　　根据ADSP-21262性能可估算出系统完成一次调幅度测量所需要的时间大概为800μs，完成一次载波频率测量所需要的时间大概为10ms，可以满足系统实时性要求。
　　在总线控制模块中，系统选用Altera公司Cyclone II系列中的EP2C8Q208C8 FPGA芯片。
　　PCI接口模块选用PLX公司的PCI总线控制芯片PC19054。