本发明涉及数字接收机领域，具体涉及一种接收机带宽测量方法。
背景技术：
在数字接收机领域，尤其是非协作通信领域，接收机对于当前所接收信号的具体参数未知。以数字超外差接收机为例，信号经过射频模块被搬移至中频，中频信号通过高速adc采样转化数字信号，便于后续的数字信号处理。数字下变频模块ddc将一路实信号进行正交混频，通过混频将信号搬移至零中频，混频之后的信号成为iq两路，二者同相位且相互正交。后续的频谱测量模块、解调模块等通过输入的iq信号进行相应的处理。
在数字信号处理过程中，需要多次进行数字滤波，因此滤波器带宽的选择直接影响后续解调的效果；如果滤波器选择过宽，则引入更多的噪声；如果滤波器选择过窄，则会丢失有用信号；因此在接收机工作时需要根据当前的频谱情况不断地手动进行修正参数设置，常见的接收机只提供电平测量和频偏测量，而实际上还可以通过运算从信号中进一步提取带宽信息，通过提取到的信息，可以在界面上提示使用者当前最合适的工作带宽等。
技术实现要素：
针对现有技术的上述不足，本发明提供了一种针对常用模拟调制及数字调制信号的数字接收机中信号带宽测量方法。
为达到上述发明目的，本发明所采用的技术方案为：
提供一种接收机带宽测量方法，其包括以下步骤：
s1：设置有效信号大于噪声基底的门限值delta，设置窗长度为winlen；
s2：对于输入的iq信号加窗，并将每个窗中的数据看做1包；
s3：对当前1包数据进行fft运算得到fft信号，再对fft信号进行fftshift运算得到fftshift信号；
s4：计算fftshift信号的功率谱密度值psd；
s5：对功率谱密度值psd计算平均，得到噪声基底估计值noisetmp，并采用多次迭代去强信号留弱信号的方式，计算信号的噪声基底值noise，以该值作为信号的实际噪声基底值；
s6：将噪声基底值noise加上门限值delta，得到有效信号的阈值threshold，并以功率谱密度值psd中大于阈值threshold的部分作为有效信号的范围，计算有效信号范围内信号的总功率power；
s7：计算总功率power的99％部分为所测信号的带宽范围。
进一步地，fft运算的计算公式为：
其中，x为x的频域幅度值，j为虚数单位；e为自然常数；i为每包时域信号x的下标；n为傅里叶变换的点数，并且n等于winlen。
进一步地，噪声基底值noisetmp的计算公式为：
其中，p(k)表示功率谱密度值psd；
采用多次迭代去强信号留弱信号的方式，计算信号真实的噪声基底值noise的计算公式为：
进一步地，步骤s6中计算有效信号的范围内信号的总功率power包括：
s61：在0～winlen-1的范围内遍历功率谱密度值psd，得到psd最大值的位置posmax；
s62：从位置posmax往winlen进行递加，且每次递加1，直到求得点per，per满足在范围[per，per+m]内，至少有90％×m个点的值低于阈值threshold，其中m＝winlen/10；
s63：从位置posmax往0进行递减，且每次递减1，直到求得某个点pel，满足在范围[pel-m，pel]内，至少有90％×m个点的值低于阈值threshold，其中m＝winlen/10；
s64：计算有效信号的总功率power：
其中，p(i)为功率谱密度值psd。
进一步地，步骤s7进一步包括：
s71：在[pel，per]的范围内，计算位置pbl：
s72：在[pel，per]的范围内，计算位置pbr：
s73：信号的带宽部分bw为：
其中，fs为信号采样率。
进一步地，门限值delta表示在信号的功率谱中，信号功率大于噪声基底值noise与门限值delta之和的信号视为有效信号。
本发明的有益效果为：本发明提出的接收机带宽测量方法，可通过设置不同的有效信号门限值确定有效信号的界定条件，通过设置不同的窗长度来控制信号测量的精度。在整个算法过程中的时间复杂度为o(nlogn)，其运行效率有利于在实际设备中进行使用。
通过测量信号的带宽信息，可以提示使用者最佳的参数设定，不仅有利于提升工作效率，而且便于进一步实现接收机自适应工作。
附图说明
图1为接收机带宽测量方法的流程图。
图2为采样率为1.6×106hz的带宽测量范围的仿真图。
图3为采样率为3.2×106hz的带宽测量范围的仿真图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，接收机带宽测量方法包括以下步骤：
s1：设置有效信号大于噪声基底的门限值delta，单位为db；该值表示在信号的功率谱中，信号功率大于噪声基底与门限值delta之和的信号视为有效信号；并设置窗长度为winlen；
s2：选用适当的窗，并使用对应的窗函数，将每个窗中的数据看做1包。
本方案优选使用矩形窗：
x(i)＝r(i)×w(i)i＝0,1......winlen-1
x为加窗后的iq值，r为iq时域数据，w为窗函数值；
s3：进行fft运算：
其中，x为x的频域幅度值，n为fft的步长，即为傅里叶变换的点数，n＝winlen。
因为信号fft运算后的范围为[0，fs]，而在后续的运算中，需要对应fft运算在
s4：将每个样点值x′(k)进行对数运算，并乘以10，得到信号功率谱密度值psd：
p(k)＝10×log10(x'(k))
s5：采用多次迭代去强信号留弱信号的方式，估算信号的噪声基底。本方案进行二次迭代。步骤如下：
s51：对功率谱密度值psd数据求和，再计算平均，得到一个粗略的估计值noisetmp：
s52：对于p(k)中大于noisetmp的值用当前的noisetmp值替代进行累加，然后求平均，得到均值noise，将其作为噪声基底值：
s6：令threshold＝noise+delta，将其作为有效信号的阈值，大于该值的部分视为有效信号部分。计算有效信号的总功率，具体流程如下：
s61：从0～winlen-1遍历psd值，得到psd取得最大值时的位置posmax；
s62：从位置posmax往winlen进行递加，每次递加1，直到求得某个点per，使得在范围[per，per+m]内，至少有90％×m个点的值低于阈值threshold，其中
s63：从位置posmax往0进行递减，每次递减1，直到求得某个点pel，满足在范围[pel-m，pel]内，至少有90％×m个点的值低于阈值threshold，其中
s64：pel和per的中间部分就是有效信号范围。累加从pel到per的范围内阈值threshold以上的功率，即累加时要减去有效信号的阈值threshold，得到有效信号的总功率power。
s7：求在[pel,per]范围内，99％×power的部分，将该部分所占带宽作为所需测量的信号带宽；步骤如下：
s71：从pel开始累加psd值，累加时要减去有效信号的阈值threshold，pel递增，直到累加值刚好超过0.5％×power，此时的位置为pbl，即：
s72：从per开始累加psd值，累加时要减去有效信号的阈值threshold，per递减，直到累加值刚好超过0.5％×power，此时的位置为pbr，即：
s73：从pbl到pbr的部分即为信号的带宽部分；带宽bw的计算为：
其中，fs为信号采样率。
本发明的仿真效果如图2和图3所示，仿真条件为：已调信号为fm信号，其载波频率为12.8mhz，经过模拟ddc，将信号搬移到零中频，并进行降采样，并且图中的采样率为降采样后的采样率。图中的两条竖线表示测量过程中的pbl和pbr，即测量出的带宽范围。
本发明提出的接收机带宽测量方法，可通过设置不同的有效信号门限值确定有效信号的界定条件，通过设置不同的窗长度来控制信号测量的精度。在整个算法过程中的时间复杂度为o(nlogn)，其运行效率有利于在实际设备中进行使用。
通过测量信号的带宽信息，可以提示使用者最佳的参数设定，不仅有利于提升工作效率，而且便于进一步实现接收机自适应工作。