[发明专利]群时延失真的修正方法及其装置

权利要求书

1.一种群时延失真的修正方法，其特征在于，包括：

步骤一、接收训练序列h’；

步骤二、对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f；

步骤三、对频率特性grf_f进行计算，得到群时延系数grpts_est；

步骤四、对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV；

步骤五、根据修正后的群时延系数Grpts_est_INV对输入信号进行卷积，完成群时延失真的修正；

步骤二包括，根据频率相除法对接收到的训练序列h’进行计算，得到频率特性grf_f；频率相除法的计算公式为

其中，DFT为离散傅里叶变换，n为测试点，n的取值范围为1至N的正整数，h’为训练序列，real为matlab软件中的常用描述，表示取实部。

2.根据权利要求1所述的群时延失真的修正方法，其特征在于，还包括：在步骤一之前，对原始训练序列h进行信号转换生成训练序列h’。

3.根据权利要求2所述的群时延失真的修正方法，其特征在于，对原始训练序列h进行信号转换生成训练序列h’包括：

采集原始训练序列h；

将采集到的原始训练序列h进行数模转换，转换成模拟信号；

将转换成模拟信号的原始训练序列h依次进行射频转换和信号放大；

将经过射频转换和信号放大后的原始训练序列h进行模数转换，生成训练序列h’。

4.根据权利要求1～3任一项所述的群时延失真的修正方法，其特征在于，原始训练序列h和训练序列h均为数字信号，该数字信号包括：多音信号、时域脉冲信号、以及经过逆傅里叶变换之后的信号。

5.根据权利要求1所述的群时延失真的修正方法，其特征在于，步骤二包括，根据相位相关法对接收到的训练序列h’进行计算，得到频率特性grf_f；

相位相关法包括：

S21、根据群时延测试公式得到群时延有效测试值

S22、对N-1个测试点的群时延有效测试值进行频域平滑，

S23、根据信道内的带宽BW以及有效相位测试点的个数N，计算有效相位测试点之间的间距，群时延的孔径Δw＝BW/(N-1)；

S24、设参考相位点为根据参考相位点运用参考相位平均微分法，计算出参考相位点的群时延测试值参考相位点的群时延测试值为

S25、根据以下公式计算频率特性grf_f，

上述公式中，群时延的孔径Δω＝BW/(N-1)；BW为信道带宽，用Hz表示；i为信道带宽BW内有效相位测试点，i的取值范围为1到N的正整数；参考相位点的取值范围为2到N-1的正整数；为相位移；ω为群时延的角频率。

6.根据权利要求1所述的群时延失真的修正方法，其特征在于，步骤三包括，根据公式grpts_est＝IFFT(H(ω))＝IFFT(|H(ω)|·exp(jθ))对频率特性grf_f进行计算，得到群时延系数grpts_est；

其中，|H(ω)|为快速傅里叶变换后的幅频特性的值，θ为群时延滤波器系数的相位信息，ω为角频率，j为虚部表示符号，IFFT为逆傅里叶变换。

7.根据权利要求1所述的群时延失真的修正方法，其特征在于，步骤四包括，根据系数倒序法对对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV；

系数倒序法的计算公式为grpts_est_INV＝grpts_est(end:-1:1)；

其中，end:-1:1为matlab软件中的常用描述，表示按从尾到头的顺序往前数。