CN106850103B - 群时延失真的修正方法及其装置 - Google Patents

群时延失真的修正方法及其装置 Download PDF

Info

Publication number
CN106850103B
CN106850103B CN201710110430.2A CN201710110430A CN106850103B CN 106850103 B CN106850103 B CN 106850103B CN 201710110430 A CN201710110430 A CN 201710110430A CN 106850103 B CN106850103 B CN 106850103B
Authority
CN
China
Prior art keywords
group delay
training sequence
grpts
est
coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710110430.2A
Other languages
English (en)
Other versions
CN106850103A (zh
Inventor
熊军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Yu Fei Electronic Technology Co., Ltd.
Original Assignee
Xi'an Yufei Electronic Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xi'an Yufei Electronic Technology Co ltd filed Critical Xi'an Yufei Electronic Technology Co ltd
Priority to CN201710110430.2A priority Critical patent/CN106850103B/zh
Publication of CN106850103A publication Critical patent/CN106850103A/zh
Application granted granted Critical
Publication of CN106850103B publication Critical patent/CN106850103B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • H04B17/364—Delay profiles
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/24—Testing correct operation
    • H04L1/248—Distortion measuring systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Nonlinear Science (AREA)
  • Quality & Reliability (AREA)
  • Electromagnetism (AREA)
  • Amplifiers (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明公开一种群时延失真的修正方法及其装置，涉及卫星通信技术领域，用于有效地缩短群时延波动范围。群时延失真的修正方法包括：接收训练序列h’；对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f；对频率特性grf_f进行计算，得到群时延系数grpts_est；对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV；根据修正后的群时延系数Grpts_est_INV对输入信号进行卷积，完成群时延失真的修正。本发明用于模拟卫星信道中数据的传输。

Description

群时延失真的修正方法及其装置
技术领域
本发明涉及卫星通信技术领域，尤其涉及一种群时延滤波器失真的修正方法及其装置。
背景技术
宽带卫星通信系统的高速数据传输达到了几百Mbit/s，在高速数据传输的过程中，群时延(系统在某频率处的相位对于频率的变化率)在宽带卫星通信的通道(以下简称信道)中的传播会影响宽带卫星通信的传输，例如，若群时延波动过大，会使传输信号产生畸变，进一步地，以数字信号为例，由于数字信号相邻码元之间不同频率分量之间的时延不同，导致接收信号会在相邻码元之间产生串扰，增加误码率。群时延的理想状态为一条水平直线，其值理想状态下为一个常数。
目前，对于群时延造成的传输的高速数据的失真(以下简称群时延失真)可通过群时延滤波器进行调整。高速数据通过群时延滤波器后允许有一定的时域延时，但不希望产生在频域上产生波形的失真。为了保持波形，可在群时延滤波器中添加时延均衡器。时延均衡器虽然能够使群时延波动相对变小，但是时延均衡器的参数是固定的,只能被动地相对减小群时延波动，对于波动较大的波形则效果不显著。
发明内容
本发明的目的在于提供一种群时延失真的修正方法及其装置，用于有效地缩短群时延波动范围。
为达到上述目的，本发明采用如下技术方案：
本发明的第一方面提供一种群时延失真的修正方法，该群时延失真的修正方法包括：
步骤一、接收训练序列h’。
步骤二、对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f。
步骤三、对频率特性grf_f进行计算，得到群时延系数grpts_est。
步骤四、对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV。
步骤五、根据修正后的群时延系数Grpts_est_INV对输入信号进行卷积，完成群时延失真的修正。
本发明的第一方面的有益效果如下：根据群时延失真的修正方法，通过训练序列h’得到所有需要频率点的群时延信息，与现有技术中通过单个频率点扫描相比，一次训练就可得到训练序列h’，操作简单快捷，并且节省资源。接着，将收集到训练序列h’，经过三步计算，即可得到修正后的群时延系数。另外，该群时延失真的修正方法可周期性的修正群时延系数，通过不断修正最佳值，将群时延限制在一定范围内，进一步地，可将群时延波动范围控制在0.5ns。
群时延失真的修正方法还包括：在步骤一之前，对原始训练序列h进行信号转换生成训练序列h’。
对原始训练序列h进行信号转换生成训练序列h’包括：采集原始训练序列h；将采集到的原始训练序列h进行数模转换，转换成模拟信号；将转换成模拟信号的原始训练序列h依次进行射频转换和信号放大；将经过射频转换和信号放大后的原始训练序列h进行模数转换，生成训练序列h’。
原始训练序列h和训练序列h均为数字信号，该数字信号包括：多音信号、时域脉冲信号、以及经过逆傅里叶变换之后的信号。
步骤二包括，根据频率相除法对接收到的训练序列h’进行计算，得到频率特性grf_f；频率相除法的计算公式为
Figure BDA0001234264120000031
其中，DFT为离散傅里叶变换，n为测试点，n的取值范围为1至N的正整数，h’为训练序列。
步骤二包括，根据相位相关法对接收到的训练序列h’进行计算，得到频率特性grf_f；相位相关法包括：
S21、根据群时延测试公式
Figure BDA0001234264120000032
得到群时延有效测试值
Figure BDA0001234264120000033
S22、对N-1个测试点的群时延有效测试值进行频域平滑，
Figure BDA0001234264120000034
S23、根据信道内的带宽BW以及有效相位测试点的个数N，计算有效相位测试点之间的间距，有效相位测试点之间的间距bw＝BW/(N-1)。
S24、设参考相位点为
Figure BDA0001234264120000035
根据参考相位点
Figure BDA0001234264120000036
运用参考相位平均微分法，计算出参考相位点的群时延测试值
Figure BDA0001234264120000037
参考相位点的群时延测试值为
Figure BDA0001234264120000038
S25、根据以下公式计算频率特性grf_f，
Figure BDA0001234264120000039
上述公式中，群时延的孔径Δω＝BW/(N-1)；BW为信道带宽，用Hz表示；i为信道带宽BW内有效相位测试点，i的取值范围为1到N的正整数；参考相位点
Figure BDA00012342641200000310
的取值范围为2到N-1的正整数；
Figure BDA00012342641200000311
为相位移；ω为群时延的角频率。
步骤三包括，根据公式grpts_est＝IFFT(H(ω))＝IFFT(H(ω)|·exp(jθ))对频率特性grf_f进行计算，得到群时延系数grpts_est；其中，|H(ω)|为快速傅里叶变换后的幅频特性的值，θ为群时延滤波器系数的相位信息，ω为角频率，j为虚部表示符号，IFFT为逆傅里叶变换。
步骤四包括，根据系数倒序法对对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV；系数倒序法的计算公式为grpts_est_INV＝grpts_est(end:-1:1)。
本发明的第二方面提供一种群时延失真的修正装置，该群时延失真的修正装置用于对群时延失真的修正，群时延失真的修正装置适用通过上述本发明的第一方面所述的群时延失真的修正方法计算得到的修正后的群时延系数。群时延失真的修正装置包括：
接收模块，用于接收训练序列h’；与接收模块相连的频率特性计算模块，用于对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f；与频率特性计算模块相连的群时延系数计算模块，用于对频率特性grf_f进行计算，得到群时延系数grpts_est；与群时延系数计算模块相连的逆函数计算模块，用于对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV，其中，逆函数计算模块通过高速群时延滤波模块与数模转换模块相连。
该群时延失真的修正装置还包括：
采集模块，用于采集原始训练序列h；与采集模块相连的数模转换模块，用于将采集到的原始训练序列h进行数模转换，转换成模拟信号；与数模转换模块相连的射频转换模块，用于将转换成模拟信号的原始训练序列h进行射频转换，转换成射频信号；与射频转换模块相连的功率放大模块，用于将转换成射频信号的原始训练序列h进行信号放大；与功率放大模块相连的射频转换接收模块，用于接收信号放大后的原始训练序列h；与射频转换接收模块相连的模数转换模块，用于将经过射频转换和信号放大后的原始训练序列h转换成训练序列h’，其中，模数转换模块还与接收计算模块相连。
根据本发明的第二方面所述的群时延失真的修正装置，可得到修正后的群时延系数，操作简单，实时性强，该群时延失真的修正装置可周期性的修正群时延系数，通过不断修正最佳值，将群时延限制在一定范围内，进一步地，可将群时延波动范围控制在0.5ns。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为现有技术中的群时延-频率曲线图；
图2为本发明实施例中群时延失真的修正方法流程图；
图3为本发明实施例中群时延失真的修正装置示意图；
图4为本发明实施例中修正后的群时延系数波动图；
图5为本发明实施例中群时延-载波信号频率曲线图一；
图6为本发明实施例中群时延-载波信号频率曲线图二；
图7为本发明实施例中群时延-载波信号频率曲线图三；
图8为本发明实施例中群时延-载波信号频率曲线图四；
图9为本发明实施例中群时延-载波信号频率曲线图五；
图10为本发明实施例中群时延-载波信号频率曲线图六。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为了本领域技术人员更好的理解“群时延”，以举例的方式给出群时延的定义：
宽带卫星通信系统中电离层的充电粒子的减缓了无线电信号在信道中的传播，超过自由空间传播时间的时间延迟可称为群时延，通常以t表示。当频率不同时，时延不同，此时一群频率的延时就是群时延。在MSS系统中，它是必须考虑的重要因素。现有技术中可通过以下公式计算群时延的数量：t＝1.345NT/f2×10-7，其中，t表示与真空中传播相比的时延，单位为s；f表示频率，单位为Hz；NT表示总电子含量(以下简称TEC)，单位为e/m2，TEC可由倾斜的传播路径决定。
示例性地，当TEC在1016～1019e/m2区间变化时，在1600MHz频带附近，信号的群时延大概在0.5ns到500ns区间变化。
示例性地，如图1所示，当TEC在1017e/m2，频带在1.6GHz～1.8GHz时，群时延随着频率的增大而减小。
实施例一
本实施例提供一种群时延失真的修正方法，如图2所示，该群时延失真的修正方法包括：
步骤一、接收训练序列h’。
步骤二、对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f。
步骤三、对频率特性grf_f进行计算，得到群时延系数grpts_est。
步骤四、对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV。
步骤五、根据修正后的群时延系数Grpts_est_INV对输入信号进行卷积，完成群时延失真的修正。
示例性地，输入信号可为用户的业务数据信号。
根据群时延失真的修正方法，通过训练序列h’得到所有需要频率点的群时延信息，与现有技术中通过单个频率点扫描相比，一次训练就可得到训练序列h’，操作简单快捷，并且节省资源。接着，将收集到训练序列h’，经过三步计算，即可得到修正后的群时延系数。另外，该群时延失真的修正方法可周期性的修正群时延系数，通过不断修正最佳值，将群时延限制在一定范围内，进一步地，可将群时延波动范围控制在0.5ns。
群时延失真的修正方法还包括：在步骤一之前，对原始训练序列h进行信号转换生成训练序列h’。
进一步地，对原始训练序列h进行信号转换生成训练序列h’包括：
采集原始训练序列h；
将采集到的原始训练序列h进行数模转换，转换成模拟信号；
将转换成模拟信号的原始训练序列h依次进行射频转换和信号放大；
将经过射频转换和信号放大后的原始训练序列h进行模数转换，生成训练序列h’。
另外，本实施例中涉及到的原始训练序列h和训练序列h均为数字信号，示例性地，该数字信号可包括：多音信号、时域脉冲信号、以及经过逆傅里叶变换之后的信号。
值得一提的是，本实施例中的原始训练序列h是群时延时域上的一个脉冲，可以发送时域脉冲信号，其代表的频率为全通带。
为了得到频率特性grf_f，示例性地，可通过如下两种方式：
第一种方式，步骤二包括，根据频率相除法对接收到的训练序列h’进行计算，得到频率特性grf_f；频率相除法的计算公式为
Figure BDA0001234264120000071
其中，DFT为离散傅里叶变换，n为测试点，n的取值范围为1至N的正整数，h’为训练序列。
第二种方式步骤二包括，根据相位相关法对接收到的训练序列h’进行计算，得到频率特性grf_f。
利用该相位相关法计算得到频率特性grf_f的步骤包括：
S21、根据群时延测试公式
Figure BDA0001234264120000081
得到群时延有效测试值
Figure BDA0001234264120000082
S22、对N-1个测试点的群时延有效测试值进行频域平滑，
Figure BDA0001234264120000083
S23、根据信道内的带宽BW以及有效相位测试点的个数N，计算有效相位测试点之间的间距，有效相位测试点之间的间距bw＝BW/(N-1)。
S24、设参考相位点为
Figure BDA0001234264120000084
根据参考相位点
Figure BDA0001234264120000085
运用参考相位平均微分法，计算出参考相位点的群时延测试值
Figure BDA0001234264120000086
参考相位点的群时延测试值为
Figure BDA0001234264120000087
S25、根据以下公式计算频率特性grf_f，频率特性
Figure BDA0001234264120000088
上述公式中，群时延的孔径Δω＝BW/(N-1)；BW为信道带宽，用Hz表示；i为信道带宽BW内有效相位测试点，i的取值范围为1到N的正整数；参考相位点
Figure BDA0001234264120000089
的取值范围为2到N-1的正整数；
Figure BDA00012342641200000810
为相位移；ω为群时延的角频率。
步骤三包括，根据公式grpts_est＝IFFT(H(ω))＝IFFT(H(ω)|·exp(jθ))对频率特性grf_f进行计算，得到群时延系数grpts_est；其中，|H(ω)|为快速傅里叶变换后的幅频特性的值，θ为群时延滤波器系数的相位信息，ω为角频率，j为虚部表示符号，IFFT为逆傅里叶变换。
进一步地，为了更好的理解群时延系数的计算过程，下面将对群时延系数的计算过程做详细的描述：
首先，对幅频特性|H(ω)|以及相频特性θ(ω)分别进行快速傅里叶变换变换，具体如下：
对幅频特性进行快速傅里叶变换包括：使用公式
Figure BDA0001234264120000091
对幅频特性进行快速傅里叶变换，其中，|H(ω)|为快速傅里叶变换后的幅频特性的值，n的取值范围为1至N的正整数，N为快速傅里叶变换的长度，τS为群时延在各个频率点的时延值。
对相频特性进行快速傅里叶变换包括：使用公式
Figure BDA0001234264120000092
对相频特性进行快速傅里叶变换，其中，θ(ω)为快速傅里叶变换后的相频特性的值，n的取值范围为1至N的正整数，N为快速傅里叶变换的长度，τS为群时延在各个频率点的时延值。
更进一步地，根据计算得到的幅频特性的值以及相频特性的值，得出群时延系数。
步骤四包括，根据系数倒序法对对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV；系数倒序法的计算公式为grpts_est_INV＝grpts_est(end:-1:1)。
实施例二
本实施例提供一种群时延失真的修正装置，如图3所示，该群时延失真的修正装置用于对群时延失真的修正，群时延失真的修正装置适用上述实施例一中通过群时延失真的修正方法计算得到的修正后的群时延系数。群时延失真的修正装置包括：
接收模块7，用于接收训练序列h’。
与接收模块7相连的频率特性计算模块8，用于对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f。
与频率特性计算模块8相连的群时延系数计算模块9，用于对频率特性grf_f进行计算，得到群时延系数grpts_est。
与群时延系数计算模块9相连的逆函数计算模块10，用于对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV，其中，逆函数计算模块10通过高速群时延滤波模块11与数模转换模块2相连。
继续参见图3，该群时延失真的修正装置还包括：
采集模块1，用于采集原始训练序列h。
与采集模块1相连的数模转换模块2，用于将采集到的原始训练序列h进行数模转换，转换成模拟信号。
与数模转换模块2相连的射频转换模块3，用于将转换成模拟信号的原始训练序列h进行射频转换，转换成射频信号。
与射频转换模块3相连的功率放大模块4，用于将转换成射频信号的原始训练序列h进行信号放大。
与功率放大模块4相连的射频转换接收模块5，用于接收信号放大后的原始训练序列h。
与射频转换接收模块3相连的模数转换模块6，与射频转换接收模块相连的模数转换模块，用于将经过射频转换和信号放大后的原始训练序列h转换成训练序列h’，其中，模数转换模块还与接收计算模块相连。
需要说明的是，上述各模块具体执行其功能的方式参见实施例一即可，此处不再进行赘述。
根据本实施例中的群时延失真的修正装置，可得到修正后的群时延系数，操作简单，实时性强，该群时延失真的修正装置可周期性的修正群时延系数，通过不断修正最佳值，将群时延限制在一定范围内，进一步地，可将群时延波动范围控制在0.5ns。
实施例三
为了更好的阐述上述实施例一的群时延失真的修正方法以及实施例二中群时延失真的修正装置，本实施例提供一种验证方法，下面对验证方法进行详细的描述：
本实施例中采用MATLAB进行验证，首先对信道中的群时延进行模拟，在信道中设置参数如下：
开始参数为：frf＝1.6*109；1.6GHz。
扫描带宽参数为：BW＝200*106；200MHz。
扫频间隔参数为：detf＝1*106；1MHz。
带宽参数为：f＝[frf+(detf:detf:bw)]；160MHz，精度为1MHz，扫频为200MHz。
其中，Nt＝[1017]，表示总电子含量(TEC，m-2)。垂直柱状体的TEC在1016到1018el/m2的范围内变化；群时延的单位是秒(s)，本实施例用用纳秒(ns)表示；tao_freq＝1.345*Nt/f2*10-7；tao_ns＝tao.*10-9。
接着，设置训练序列h，训练序列h是群时延时域上的一个脉冲，可以发送时域脉冲信号，代表的频率为全通带。如图4所示，全通道(信道)信号中，不同频率点之间波动限制在0.002DBC之间。进一步地，从图4可知，根据群时延失真的修正方法对群时延失真进行修正后，群时延从原来的+-4ns，限制到+-0.15ns，提高了约27倍。
参见图5，可从中直观的得出，通过MATLAB仿真计算得到的群时延(仿真设置频率下的群时延)与通过群时延失真的修正方法得到的群时延(采集反馈信号h分析得到的群时延)，彼此重叠在一起，说明利用上述群时延失真的修正方法得到的修正后的群时延系数是真实有效的。
进一步地，对上述图5做进一步地放大得到图6，如图6所示，利用群时延失真的修正方法得到的修正后的群时延系数的波动较小。
上述图5以及图6中采集反馈信号h即为本实施例中所述的训练序列h’。
参见图7，利用群时延失真的修正方法得到的修正后的群时延系数，与不经过修正的群时延系数的对比图，从图7中可直观的得知，修正后的群时延系数接近一条笔直的横线，上下波动范围在0.2ns，而不经过修正的群时延系数(修正前的群时延)是一条斜率较大的线，这样的斜线，容易造成较大的误码率。
本发明中提供的群时延失真的修正方法，不仅适用于波形较为平整群时延信号，也适用波形较复杂的群时延信号，如图8所示，选取任意波形的群时延信号，从而利用MATLAB模拟计算得到的群时延系数(仿真设置频率下的群时延)与利用群时延失真的修正方法计算得到的修正后的群时延系数，彼此重叠，说明群时延失真的修正方法同样适用于任意波形。
需要说明的是，图8中的采集反馈信号h即为本实施例中所述的训练序列h’。
参见图9，对任意波形的群时延信号，利用群时延失真的修正方法得到的修正后的群时延系数(修正后的群时延)为一条比较平直的水平线。而没有经过修正的群时延系数(修正前的群时延)斜率较大，且波动较大，这样的曲线会使信道中传输信号失真，增大误码率。
进一步地，对图9进行放大得到图10，从图10中可直接看出，利用群时延失真的修正方法修正的群时延系数(修正后的群时延)波动在0.02ns之间，约为一条直线。而没有修正的群时延系数(修正前的群时延)波动范围很大，呈现出两条竖直线，说明误码率大。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (9)

1.一种群时延失真的修正方法，其特征在于，包括：
步骤一、接收训练序列h’；
步骤二、对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f；
步骤三、对频率特性grf_f进行计算，得到群时延系数grpts_est；
步骤四、对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV；
步骤五、根据修正后的群时延系数Grpts_est_INV对输入信号进行卷积，完成群时延失真的修正；
步骤二包括，根据频率相除法对接收到的训练序列h’进行计算，得到频率特性grf_f；频率相除法的计算公式为
Figure FDA0002482763550000011
其中，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、根据群时延测试公式
Figure FDA0002482763550000021
得到群时延有效测试值
Figure FDA0002482763550000022
S22、对N-1个测试点的群时延有效测试值进行频域平滑，
Figure FDA0002482763550000023
S23、根据信道内的带宽BW以及有效相位测试点的个数N，计算有效相位测试点之间的间距，群时延的孔径Δw＝BW/(N-1)；
S24、设参考相位点为
Figure FDA0002482763550000024
根据参考相位点
Figure FDA0002482763550000025
运用参考相位平均微分法，计算出参考相位点的群时延测试值
Figure FDA0002482763550000026
参考相位点的群时延测试值为
Figure FDA0002482763550000027
S25、根据以下公式计算频率特性grf_f，
Figure FDA0002482763550000028
上述公式中，群时延的孔径Δω＝BW/(N-1)；BW为信道带宽，用Hz表示；i为信道带宽BW内有效相位测试点，i的取值范围为1到N的正整数；参考相位点
Figure FDA0002482763550000029
的取值范围为2到N-1的正整数；
Figure FDA00024827635500000210
为相位移；ω为群时延的角频率。
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软件中的常用描述，表示按从尾到头的顺序往前数。
8.一种群时延失真的修正装置，该群时延失真的修正装置用于对群时延失真的修正，其特征在于，群时延失真的修正装置适用通过上述权利要求1～7任一项所述的群时延失真的修正方法计算得到的修正后的群时延系数；
群时延失真的修正装置包括：
接收模块，用于接收训练序列h’；
与接收模块相连的频率特性计算模块，用于对接收到的训练序列h’进行计算，得到群时延的频率特性grf_f；
与频率特性计算模块相连的群时延系数计算模块，用于对频率特性grf_f进行计算，得到群时延系数grpts_est；
与群时延系数计算模块相连的逆函数计算模块，用于对群时延系数grpts_est进行逆函数计算，得到修正后的群时延系数grpts_est_INV，其中，逆函数计算模块通过高速群时延滤波模块与数模转换模块相连。
9.根据权利要求8所述的群时延失真的修正装置，其特征在于，该群时延失真的修正装置还包括：
采集模块，用于采集原始训练序列h；
与采集模块相连的数模转换模块，用于将采集到的原始训练序列h进行数模转换，转换成模拟信号；
与数模转换模块相连的射频转换模块，用于将转换成模拟信号的原始训练序列h进行射频转换，转换成射频信号；
与射频转换模块相连的功率放大模块，用于将转换成射频信号的原始训练序列h进行信号放大；
与功率放大模块相连的射频转换接收模块，用于接收信号放大后的原始训练序列h；
与射频转换接收模块相连的模数转换模块，用于将经过射频转换和信号放大后的原始训练序列h转换成训练序列h’，其中，模数转换模块还与接收计算模块相连。
CN201710110430.2A 2017-02-28 2017-02-28 群时延失真的修正方法及其装置 Active CN106850103B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710110430.2A CN106850103B (zh) 2017-02-28 2017-02-28 群时延失真的修正方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710110430.2A CN106850103B (zh) 2017-02-28 2017-02-28 群时延失真的修正方法及其装置

Publications (2)

Publication Number Publication Date
CN106850103A CN106850103A (zh) 2017-06-13
CN106850103B true CN106850103B (zh) 2020-09-29

Family

ID=59133651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710110430.2A Active CN106850103B (zh) 2017-02-28 2017-02-28 群时延失真的修正方法及其装置

Country Status (1)

Country Link
CN (1) CN106850103B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111490785B (zh) * 2020-04-26 2021-11-09 华中科技大学 一种数据采集系统的时延测量装置及方法
CN111580137B (zh) * 2020-05-18 2022-04-22 中国人民解放军国防科技大学 一种高精度导航接收机射频通道群时延特性的拟合方法
CN113708853B (zh) * 2021-08-09 2024-03-26 中国计量科学研究院 一种针对天线群时延测量的阻抗失配修正方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102325058A (zh) * 2011-09-06 2012-01-18 北京空间飞行器总体设计部 一种变频系统群时延测试方法
CN103281268A (zh) * 2013-04-25 2013-09-04 西安空间无线电技术研究所 一种卫星导航信号生成系统带内群时延波动的预补偿方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105959064A (zh) * 2016-04-25 2016-09-21 武汉元昊科技有限公司 实现高速信号传输的强度调制与直接检测系统及方法
CN106341101B (zh) * 2016-08-31 2018-12-14 北京睿信丰科技有限公司 群时延滤波器系数的计算方法以及滤波器的模拟系统
CN106357323A (zh) * 2016-10-08 2017-01-25 航天恒星科技有限公司 卫星信道群延时补偿方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102325058A (zh) * 2011-09-06 2012-01-18 北京空间飞行器总体设计部 一种变频系统群时延测试方法
CN103281268A (zh) * 2013-04-25 2013-09-04 西安空间无线电技术研究所 一种卫星导航信号生成系统带内群时延波动的预补偿方法

Also Published As

Publication number Publication date
CN106850103A (zh) 2017-06-13

Similar Documents

Publication Publication Date Title
US9362981B2 (en) Mechanisms for low-amplitude echo estimation
CN106850103B (zh) 群时延失真的修正方法及其装置
US6344749B1 (en) Test system for measuring frequency response and dynamic range on cable plant
CN103281268B (zh) 一种卫星导航信号生成系统带内群时延波动的预补偿方法
CN103888209A (zh) 一种宽带接收阵列天线通道幅相误差时域校正方法
WO2004010603A2 (en) Frequency domain equalization of communication signals
CN101729461B (zh) 消除单频干扰及多频干扰的系统及方法
US9641286B2 (en) Transmitter noise injection
CN110221268B (zh) 一种相位-幅度校准滤波器以及提升短脉冲非相参雷达信号处理性能的方法
CN107797099B (zh) 一种多通道数字接收机实时内定标处理方法及装置
CN104486267B (zh) 一种短波信道下基于小波去噪的sc‑fde信道估计方法
JP2001332999A (ja) 適応等化器のトレーニング回路及びモデム装置並びに通信装置
US20210266096A1 (en) Systems and methods for multi-carrier signal echo management using pseudo-extensions
US20240250858A1 (en) Transmitter and receiver for and method of transmitting and receiving symbols over an orthogonal time frequency space communication channel subject to doppler spread and multipath
CN110927750A (zh) 一种基于格型滤波Burg谱估计算法的低轨卫星多普勒频偏捕获方法
CN103763225A (zh) 一种基于时间反演的信号处理装置及处理方法
CN114189293B (zh) 一种宽带接收阵列天线通道幅相校准方法及系统
CN112383492B (zh) 应用于短波ofdm双选天波信道估计的递归压缩感知方法及系统
Ivanov et al. Mitigation of Dispersion Distortions of Transionospheric Communication Channels when Total Electron Content Measurements are Corrupted with Stochastic Error
CN110808750A (zh) 一种基于逆滤波的邻道干扰抑制方法和装置
JP3910453B2 (ja) 受信装置および伝送路推定方法
CN106911409B (zh) 群时延信道系统的优化方法及其装置
CN101682367B (zh) 线缆网络中的信号质量确定
CN110798176A (zh) 一种任意波宽带信号预失真滤波器构建方法及数字滤波器
CN101436876A (zh) 一种对消测试信号的方法和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20180111

Address after: 710075 Shaanxi city of Xi'an province high tech Zone Jinye Road Green Central Square - Blue Ocean fourth unit 1 building 7 layer 10713.

Applicant after: Xi'an Yu Fei Electronic Technology Co., Ltd.

Address before: 100193 room 1, 17 floor, Zhongguancun Software Park, No. 8, No. 8, Wanxi Road, northeast Haidian District, Beijing

Applicant before: BEIJING RINFON TECHNOLOGY CO., LTD.

GR01 Patent grant
GR01 Patent grant