CN110830081B - 一种信号同轴电缆传输畸变补偿方法 - Google Patents

一种信号同轴电缆传输畸变补偿方法 Download PDF

Info

Publication number
CN110830081B
CN110830081B CN201911105097.1A CN201911105097A CN110830081B CN 110830081 B CN110830081 B CN 110830081B CN 201911105097 A CN201911105097 A CN 201911105097A CN 110830081 B CN110830081 B CN 110830081B
Authority
CN
China
Prior art keywords
signal
coaxial cable
compensated
distortion compensation
compensation
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
CN201911105097.1A
Other languages
English (en)
Other versions
CN110830081A (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.)
Institute of Applied Electronics of CAEP
Original Assignee
Institute of Applied Electronics of CAEP
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 Institute of Applied Electronics of CAEP filed Critical Institute of Applied Electronics of CAEP
Priority to CN201911105097.1A priority Critical patent/CN110830081B/zh
Publication of CN110830081A publication Critical patent/CN110830081A/zh
Application granted granted Critical
Publication of CN110830081B publication Critical patent/CN110830081B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/02Details
    • H04B3/04Control of transmission; Equalising
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/28Measuring attenuation, gain, phase shift or derived characteristics of electric four pole networks, i.e. two-port networks; Measuring transient response
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/02Details
    • H04B3/46Monitoring; Testing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mathematical Physics (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

本发明公开了一种信号同轴电缆传输畸变补偿方法,属于电磁场测试的技术领域,该方法包括以下步骤:(1)测量待补偿信号并计算关键参数;(2)根据步骤(1)所获取的结果,计算同轴电缆冲击响应函数h;(3)对待补偿信号进行预处理以获取处理后的信号z;(4)对所述信号z进行畸变补偿并获取畸变补偿后的信号x,以达到可同时实现信号的幅度畸变补偿和相位畸变补偿,弥补传统补偿方法的不足,实现更好的测量信号补偿效果,确保测试结果的准确性的目的。

Description

一种信号同轴电缆传输畸变补偿方法
技术领域
本发明属于电磁场测试的技术领域,具体而言,涉及一种信号同轴电缆传输畸变补偿方法。
背景技术
同轴电缆,作为一种高效的信号传输媒介,广泛用于电信、广播电视、雷达通信等电子信息领域。IEC等国际标准对同轴电缆的相关参数进行了规范,使信号可沿同轴电缆有效地传输。然而,由于同轴电缆存在衰减、色散特性,导致信号在电缆传输过程中出现畸变。这种畸变现象在高频、宽带信号远距离传输中尤其严重,如大区域电磁环境态势感知中,感知节点与信息处理中心之间的传输距离通常较远,信号在电缆传输过程中,存在严重的损耗衰减;特别的,对于高频、宽带信号而言,在电缆传输中,还会产生较大的非线性相移,进而产生严重的色散问题,电缆输出端信号无法正确表述被测电磁环境。因此,在信号电缆传输过程中,特别是在高频、宽带、远距离传输等测试环境下,必须考虑同轴电缆的衰减、色散等问题,需对测量信号进行畸变补偿,以提高测试结果的准确性。
传统上,同轴电缆传输信号补偿方式主要有硬件网络补偿方法和基于同轴电缆衰减因子补偿的软件补偿法。硬件网络补偿法不仅需要针对每根电缆进行设计,并且难以对信号的相位进行补偿,同时还会引入了新的畸变因素;衰减因子补偿法只适用于窄带高频信号,对于宽带信号或者低频信号而言,由于衰减曲线拟合精度不高,导致补偿效果差,普适性差。此外,国内外学者还提出了直接卷积法、一阶差分法等补偿方法,但由于存在矩阵病态解、求解精度差等问题,难以推广使用。
综上所述,目前的信号补偿方法难以实现被测信号的无畸变补偿。
发明内容
鉴于此,为了解决现有技术存在的上述问题,本发明的目的在于提供一种信号同轴电缆传输畸变补偿方法以达到可同时实现信号的幅度畸变补偿和相位畸变补偿,弥补传统补偿方法的不足,实现更好的测量信号补偿效果,确保测试结果的准确性的目的。
本发明所采用的技术方案为:一种信号同轴电缆传输畸变补偿方法,该方法包括以下步骤:
(1)测量待补偿信号并计算关键参数;
(2)根据步骤(1)所获取的结果,计算同轴电缆冲击响应函数h;
(3)对待补偿信号进行预处理以获取处理后的信号z;
(4)对所述信号z进行畸变补偿并获取畸变补偿后的信号x。
进一步地,所述步骤(4)具体包括:
A.根据同轴电缆冲击响应函数h,获取计算矩阵H;
B.设定畸变补偿精度ε,中间参数α(k)、β(k)、τ(k),畸变补偿后的信号X(k);初始化待补偿信号向量Y,k=0时的中间参数α(0)、β(0)、τ(0),以及k=0时,畸变补偿后的信号X(0)
C.计算得到第k+1次畸变补偿后的信号X(k+1)
Figure BDA0002271050980000021
其中,HT是矩阵H的转置,I是单位矩阵;
D.依次计算第k+1次的α(k+1)、第k+1次的β(k+1)、第k+1次的τ(k+1),并带入至式(1)中;
E.根据下式,判断Xk+1是否满足要求
Figure BDA0002271050980000031
其中,ε为计算精度;若不满足,重复执行上述C、D步骤;若满足,则停止迭代计算,进入F步骤:
F.计算畸变补偿后的信号x
x=X(k+1)-D (3)
其中,D为设定的偏置信号。
进一步地,所述式(1)采用Gauss-Seidel迭代法求解。
进一步地,根据以下式(4)计算第k+1次的α(k+1),式(4)如下:
Figure BDA0002271050980000032
其中,Γ为Gamma函数,η1、η2、η3、η4为调节参数,N为采样点数;
根据下式(5),计算得到第k+1次的βk+1
Figure BDA0002271050980000033
其中,η1、η4为调节参数,N为采样点数;
根据下式(6),计算得到第k+1次的τk+1
Figure BDA0002271050980000041
其中,η5、η6为调节参数,N为采样点数。
进一步地,所述式(4)采用二分法求解。
进一步地,所述步骤(1)包括:
1)通过频谱测量设备测量同轴电缆输出端信号的频谱,确定同轴电缆输出端信号的最高频率Fmax和3dB带宽BW;
2)通过信号测量设备测量同轴电缆输出端信号的时域波形y,且该信号测量设备的采样率Fs由式(7)确定
Fs≥4Fmax (7)
3)根据时域波形y得到输出端波形最小值ymin
4)根据时域波形y的持续时间T和采样率Fs,由式(8)确定采样点数N
N=T×Fs (8)
进一步地,所述频谱测量设备采用频谱分析仪;所述信号测量设备采用示波器。
进一步地,计算同轴电缆冲击响应函数的具体方法如下:
a.根据BW和Fmax确定同轴电缆冲击响应函数计算的上限频率:
fmax=Fmax+3BW (9)
b.根据BW和Fmax确定同轴电缆冲击响应函数计算的下限频率:
fmin=Fmax-3BW (10)
c.通过矢量网络分析仪测量同轴电缆S21参数曲线,包括衰减曲线Sa和相频曲线Sp,测量点数为N+1,测量下限频率为fmin,上限频率为fmax
d.根据所测的同轴电缆S21参数曲线,计算同轴电缆冲击响应函数h
Figure BDA0002271050980000051
进一步地,待补偿信号的预处理方法具体如下:
(a)通过矢量网络分析仪测量同轴电缆S21参数曲线,包括衰减曲线Sa和相频曲线Sp,测量点数为N+1,测量下限频率为fmin,上限频率为fmax,计算同轴电缆的最大衰减量:
At=max(|Sa|) (12)
(b)设定偏置信号D,其中,D满足以下条件:
min(D)≥|ymin|×At (13)
(c)对待补偿信号进行偏置处理:
z=y+D*h (14)
其中,z为处理后的信号。
本发明的有益效果为:
1.采用本发明所提供的信号同轴电缆传输畸变补偿方法,通过计算同轴电缆冲击响应函数h,充分考虑了同轴电缆传输特性对信号的影响,兼顾了信号的幅度特性与相位特性,可同时对信号的幅度畸变和相位畸变进行有效补偿,无须对衰减曲线进行近似处理,相比于现有技术中常用的衰减因子补偿法,补偿效果更好,且补偿效果大幅提升。
2.采用本发明所提供的信号同轴电缆传输畸变补偿方法,不受信号带宽和频率的限制,适用范围更广,且相比于传统方法,本发明的信号同轴电缆传输畸变补偿方法的计算结果更稳定且计算精度更高。
3.采用本发明所提供的信号同轴电缆传输畸变补偿方法,可有效解决信号沿同轴电缆传输带来的畸变问题,提高信号测试准确度,特别地,对大区域电磁环境高准确度测试,具有十分重要的意义。
附图说明
图1为本发明提供的信号同轴电缆传输畸变补偿方法的整体流程图;
图2为本发明提供的信号同轴电缆传输畸变补偿方法中信号传输畸变补偿的流程图;
图3为实施例1中方波调制信号经10米同轴电缆传输后的波形图;
图4为实施例1中10米同轴电缆输出信号采用本发明方法补偿后的波形图;
图5为实施例1中10米同轴电缆输入端波形与本发明方法补偿波形对比图;
图6为实施例1中10米同轴电缆输入端波形与衰减因子补偿法补偿波形对比图;
图7为实时例1中本发明方法补偿误差与衰减因子补偿法补偿误差对比图;
图8为实施例2中双指数信号经10米同轴电缆传输后的波形图;
图9为实施例2中10米同轴电缆输出信号采用本发明方法补偿后的波形图;
图10为实施例2中10米同轴电缆输入端波形与本发明方法补偿波形对比图;
图11为实施例2中10米同轴电缆输入端波形与衰减因子补偿法补偿波形对比图;
图12为实时例2中本发明方法补偿误差与衰减因子补偿法补偿误差对比图;
图13为实施例3中阻尼振荡信号经10米同轴电缆传输后的波形图;
图14为实施例3中10米同轴电缆输出信号采用本发明方法补偿后的波形图;
图15为实施例3中10米同轴电缆输入端波形与本发明方法补偿波形对比图;
图16为实施例3中10米同轴电缆输入端波形与衰减因子补偿法补偿波形对比图;
图17为实施例3中本发明方法补偿误差与衰减因子补偿法补偿误差对比图。
具体实施方式
下面详细描述本申请的实施例,所述实施例的示例在附图中示出。下面通过参考附图描述的实施例是示例性的,仅用于解释本申请,而不能理解为对本申请的限制。相反,本申请的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
本发明提供了一种信号同轴电缆传输畸变补偿方法,如图1所示,该方法包括以下步骤:
(1)测量待补偿信号并计算关键参数,具体如下:
1)通过频谱测量设备测量同轴电缆输出端信号的频谱,确定同轴电缆输出端信号的最高频率Fmax和3dB带宽BW;
2)通过信号测量设备测量同轴电缆输出端信号的时域波形y,所述信号测量设备采用示波器且该信号测量设备的采样率Fs由式(7)确定:
Fs≥4Fmax (7)
3)根据时域波形y得到输出端波形最小值ymin
4)根据时域波形y的持续时间T和采样率Fs,由式(8)确定采样点数N
N=T×Fs (8)
(2)根据步骤(1)所获取的结果,计算同轴电缆冲击响应函数h;具体如下:
a.根据BW和Fmax确定同轴电缆冲击响应函数计算的上限频率:
fmax=Fmax+3BW (9)
b.根据BW和Fmax确定同轴电缆冲击响应函数计算的下限频率:
fmin=Fmax-3BW (10)
c.通过矢量网络分析仪测量同轴电缆S21参数曲线,包括衰减曲线(幅频曲线)Sa和相频曲线Sp,测量点数为N+1,测量下限频率为fmin,上限频率为fmax
d.根据所测的同轴电缆S21参数曲线,计算同轴电缆冲击响应函数h
Figure BDA0002271050980000081
其中,ft=f×t,即代表频率与时间之间的积;
(3)对待补偿信号进行预处理以获取处理后的信号z,具体如下:
(a)通过矢量网络分析仪测量同轴电缆S21参数曲线,包括衰减曲线Sa和相频曲线Sp,测量点数为N+1,测量下限频率为fmin,上限频率为fmax,计算同轴电缆的最大衰减量:
At=max(|Sa|) (12)
(b)设定偏置信号D,其中,D满足以下条件:
min(D)≥|ymin|×At (13)
(c)对待补偿信号进行偏置处理:
z=y+D*h (14)
其中,z即为处理后的信号。
(4)对所述信号z进行畸变补偿并获取畸变补偿后的信号x;如图2所示,该步骤(4)具体包括:
A.根据同轴电缆冲击响应函数h,获取计算矩阵H;
Figure BDA0002271050980000091
其中,h1,h2,…,hN为冲击响应函数h各离散时刻点对应的值;
B.设定畸变补偿精度ε,中间参数α(k)、β(k)、τ(k),畸变补偿后的信号X(k);初始化待补偿信号向量Y,k=0时的中间参数α(0)、β(0)、τ(0),以及k=0时,畸变补偿后的信号X(0)
Y=[z1,z2,…,zN] (16)
X(0)=[X1 (0),X2 (0),...,XN (0)] (17)
C.计算得到第k+1次畸变补偿后的信号X(k+1)
Figure BDA0002271050980000092
其中,HT是矩阵H的转置,I是单位矩阵;式(1)采用Gauss-Seidel迭代法求解,具体如下:
C1.给定初值L(0)=X(k),设置计算精度ε1
C2.计算
Figure BDA0002271050980000093
Figure BDA0002271050980000094
其中:
Figure BDA0002271050980000095
C3.判断是否满足条件:
Figure BDA0002271050980000096
若不满足,则重复执行步骤C2;若满足,则执行步骤C4;
C4.计算X(k+1)
X(k+1)=L(m+1) (21)
迭代计算终止。
D.依次计算第k+1次的α(k+1)、第k+1次的β(k+1)、第k+1次的τ(k+1),并带入至式(1)中;具体如下:
根据以下式(4)计算第k+1次的α(k+1),式(4)如下:
Figure BDA0002271050980000101
其中,Γ为Gamma函数,η1、η2、η3、η4为调节参数,N为采样点数;式(4)采用二分法求解,具体如下:
D1.令:
Figure BDA0002271050980000102
D2.给定计算精度为ε2,中间参数a,b=α(k)
D3.计算中间参数c:
Figure BDA0002271050980000103
D4.判断是否满足条件:
f(c)<ε2 (24)
若满足,则执行步骤D5;若不满足,则执行步骤D6;
D5.计算α(k+1)
α(k+1)=c (25)
迭代计算终止。
D6.如果f(c)×f(a)<0,则令:
b=c (26)
如果f(c)×f(b)<0,则令:
a=c (27)
D7.重复执行D4、D5、D6,直至迭代计算终止。
根据下式(5),计算得到第k+1次的βk+1
Figure BDA0002271050980000111
其中,η1、η4为调节参数,N为采样点数;
根据下式(6),计算得到第k+1次的τk+1
Figure BDA0002271050980000112
其中,η5、η6为调节参数,N为采样点数。
E.根据下式,判断Xk+1是否满足要求:
Figure BDA0002271050980000113
其中,ε为计算精度;若不满足,重复执行上述C、D步骤;若满足,则停止迭代计算,进入F步骤:
F.计算畸变补偿后的信号x
x=X(k+1)-D (3)
其中,D为设定的偏置信号;
通过上述的补偿方法,最终得到补偿后的信号x。
实施例1
基于本发明所提供的信号同轴电缆传输畸变补偿方法,以方波调制信号在10米同轴电缆中的传输补偿为例,计算过程如下:
(1)通过频谱分析仪测量同轴电缆输出端信号的频谱,确定同轴电缆输出端信号的最高频率2.62GHz和3dB带宽320MHz;
(2)通过信号测量设备测量同轴电缆输出端信号时域波形y,如图3所示,其中测量设备采样率为12GS/s;
(3)根据时域波形y得到输出端波形最小值为-1.6V;
(4)时域波形y的持续时间为37.5ns,计算得采样点数为N=450;
(5)确定10米同轴电缆冲击响应函数计算的上限频率为3.58GHz;
(6)确定10米同轴电缆冲击响应函数计算的下限频率为1.66GHz;
(7)通过矢量网络分析仪测量10米同轴电缆S21参数曲线,包括衰减曲线(幅频曲线)Sa和相频曲线Sp,其中测量上限频率为3.58GHz、下限频率为1.66GHz、测试点数为451;
(8)根据所测的同轴电缆S21参数曲线,计算同轴电缆冲击响应函数h;
(9)设置偏置信号D,D=2;
(10)对待补偿信号进行偏置处理,得到处理后的信号z;
(11)根据同轴电缆冲击响应函数h,获取计算矩阵H;
(12)初始化待补偿信号向量Y,畸变补偿精度ε,k=0时的中间参数α(k)、β(k)、τ(k),以及k=0时,畸变补偿后的信号X(k),取值如下:
Y=[z1,z2,…,zN]
ε=1e-5
α(0)=10
β(0)=1
τ(0)=1
X(0)=[X1 (0),X2 (0),...,X450 (0)]=[1,1,…,1]
(13)根据图2算法流程,对偏置处理后的信号z进行补偿处理,其中,畸变补偿精度ε=1e-5,k=0时的中间参数α(k)、β(k)、τ(k),以及k=0时,畸变补偿后的信号X(k)初始化取值如下:
α(0)=10
β(0)=1
τ(0)=1
X(0)=[X1 (0),X2 (0),…,X500 (0)]=[1,1,…,1]
采用Gauss-Seidel迭代法求解方程(1)计算X(k+1),计算精度为ε1=1e-5;采用二分法求解方程(4),计算α(k+1),计算精度为ε2=1e-5,η1、η2、η3、η4取值如下:
η1=1
η2=1.2
η3=100
η4=0.01
计算τ(k+1)时,η5、η6的取值如下:η5=1、η6=1e-10;
最终,得到补偿后的信号x,如图4所示。
同轴电缆输入端波形与本发明方法补偿波形对比图、同轴电缆输入端波形与衰减因子补偿法补偿波形对比图分别如图5、图6所示,可以看出,本发明方法的补偿效果更好。进一步地,本发明方法补偿后的信号与线缆输入信号之间的误差、衰减因子补偿法补偿后的信号与线缆输入信号之间的误差对比如图7所示,从图中可以更加直观地看出本发明方法的补偿效果更优。
实施例2
基于本发明所提供的信号同轴电缆传输畸变补偿方法,以双指数信号在10米同轴电缆中的传输补偿为例,计算过程如下:
(1)通过频谱分析仪测量同轴电缆输出端信号的频谱,确定同轴电缆输出端信号的最高频率94MHz和3dB带宽31.32MHz;
(2)通过信号测量设备测量同轴电缆输出端信号时域波形y,如图8所示,其中测量设备采样率为5GS/s;
(3)根据时域波形y得到输出端波形最小值为-0.2V;
(4)时域波形y的持续时间为100ns,计算得采样点数为N=500;
(5)确定10米同轴电缆冲击响应函数计算的上限频率为40kHz;
(6)确定10米同轴电缆冲击响应函数计算的下限频率为187.96MHz;
(7)通过矢量网络分析仪测量10米同轴电缆S21参数曲线,包括衰减曲线(幅频曲线)Sa和相频曲线Sp,其中测量上限频率为40kHz、下限频率为187.96MHz、测试点数为501;
(8)根据所测的同轴电缆S21参数曲线,计算同轴电缆冲击响应函数h;
(9)设置偏置信号D,D=1;
(10)对待补偿信号进行偏置处理,得到处理后的信号z;
(11)根据同轴电缆冲击响应函数h,获取计算矩阵H;
(12)初始化待补偿信号向量Y,畸变补偿精度ε,k=0时的中间参数α(k)、β(k)、τ(k),以及k=0时,畸变补偿后的信号X(k),取值如下:
Y=[z1,z2,…,zN]
ε=1e-5
α(0)=10
β(0)=1
τ(0)=1
X(0)=[X1 (0),X2 (0),…,X450 (0)]=[1,1,…,1]
(13)根据图2算法流程,对偏置处理后的信号z进行补偿处理,其中,畸变补偿精度ε=1e-5,k=0时的中间参数α(k)、β(k)、τ(k),以及k=0时,畸变补偿后的信号X(k)初始化取值如下:
α(0)=10
β(0)=1
τ(0)=1
X(0)=[X1 (0),X2 (0),…,X500 (0)]=[1,1,…,1]
采用Gauss-Seidel迭代法求解方程(1)计算X(k+1),计算精度为ε1=1e-5;采用二分法求解方程(4),计算α(k+1),计算精度为ε2=1e-5,η1、η2、η3、η4取值如下:
η1=1
η2=1.2
η3=100
η4=0.01
计算τ(k+1)时,η5、η6的取值如下:η5=1、η6=1e-10;
最终,得到补偿后的信号x,如图9所示。
同轴电缆输入端波形与本发明方法补偿波形对比图、同轴电缆输入端波形与衰减因子补偿法补偿波形对比图分别如图10、图11所示,可以看出,本发明方法的补偿效果更好。进一步地,本发明方法补偿后的信号与线缆输入信号之间的误差、衰减因子补偿法补偿后的信号与线缆输入信号之间的误差对比如图12所示,从图中可以更加直观地看出本发明方法的补偿效果更优。
实施例3
基于本发明所提供的信号同轴电缆传输畸变补偿方法,以阻尼振荡信号在10米同轴电缆中的传输补偿为例,计算过程如下:
(1)通过频谱分析仪测量同轴电缆输出端信号的频谱,确定同轴电缆输出端信号的最高频率156.36MHz和3dB带宽12.55MHz;
(2)通过信号测量设备测量同轴电缆输出端信号时域波形y,如图13所示,其中测量设备采样率为4GS/s;
(3)根据时域波形y得到输出端波形最小值为-1.42V;
(4)时域波形y的持续时间为160ns,计算得采样点数为N=640;
(5)确定10米同轴电缆冲击响应函数计算的上限频率为194.01kHz;
(6)确定10米同轴电缆冲击响应函数计算的下限频率为118.71MHz;
(7)通过矢量网络分析仪测量10米同轴电缆S21参数曲线,包括衰减曲线(幅频曲线)Sa和相频曲线Sp,其中测量上限频率为194.01kHz、下限频率为118.71MHz、测试点数为641;
(8)根据所测的同轴电缆S21参数曲线,计算同轴电缆冲击响应函数h;
(9)设置偏置信号D,D=3;
(10)对待补偿信号进行偏置处理,得到处理后的信号z;
(11)根据同轴电缆冲击响应函数h,获取计算矩阵H;
(12)初始化待补偿信号向量Y,畸变补偿精度ε,k=0时的中间参数α(k)、β(k)、τ(k),以及k=0时,畸变补偿后的信号X(k),取值如下:
Y=[z1,z2,…,zN]
ε=1e-5
α(0)=10
β(0)=1
τ(0)=1
X(0)=[X1 (0),X2 (0),…,X450 (0)]=[1,1,…,1]
(13)根据图2算法流程,对偏置处理后的信号z进行补偿处理,其中,畸变补偿精度ε=1e-5,k=0时的中间参数α(k)、β(k)、τ(k),以及k=0时,畸变补偿后的信号X(k)初始化取值如下:
α(0)=10
β(0)=1
τ(0)=1
X(0)=[X1 (0),X2 (0),…,X500 (0)]=[1,1,…,1]
采用Gauss-Seidel迭代法求解方程(1)计算X(k+1),计算精度为ε1=1e-5;采用二分法求解方程(4),计算α(k+1),计算精度为ε2=1e-5,η1、η2、η3、η4取值如下:
η1=1
η2=1.2
η3=100
η4=0.01
计算τ(k+1)时,η5、η6的取值如下:η5=1、η6=1e-10;
最终,得到补偿后的信号x,如图14所示。
同轴电缆输入端波形与本发明方法补偿波形对比图、同轴电缆输入端波形与衰减因子补偿法补偿波形对比图分别如图15、图16所示,可以看出,本发明方法的补偿效果更好。进一步地,本发明方法补偿后的信号与线缆输入信号之间的误差、衰减因子补偿法补偿后的信号与线缆输入信号之间的误差对比如图17所示,从图中可以更加直观地看出本发明方法的补偿效果更优。
综合上述实施例1、实施例2以及实施例3的计算结果,建立如下表1,表1中出示了补偿前的信号相对误差、衰减因子补偿法补偿的相对误差以及采用信号同轴电缆传输畸变补偿方法补偿的相对误差,如下:
表1
实施例 RE1(%) RE2(%) RE3(%)
1 18.38 8.50 3.42
2 11.49 8.98 0.03
3 18.60 16.68 0.051
其中:1.RE1表示同轴电缆输出与输入信号之间的相对误差;RE2表示衰减因子补偿法补偿后的信号与输入信号之间的相对误差,衰减因子补偿法为《基于衰减特性的电缆传输畸变软件补偿》;RE3表示信号同轴电缆传输畸变补偿方法补偿后的信号与输入信号之间的相对误差。
2.相对误差=RE;表示如下:
Figure BDA0002271050980000191
结合上述表1对比可知,采用信号同轴电缆传输畸变补偿方法补偿后的信号与输入信号之间的相对误差最低,该方法充分考虑了同轴电缆传输特性对信号的影响,兼顾了信号的幅度特性与相位特性,可同时对信号的幅度畸变和相位畸变进行有效补偿,相比于传统的畸变补偿方法,补偿效果大幅提升;同时,有效解决了信号沿同轴电缆传输带来的畸变问题,提高信号测试准确度,特别地,对大区域电磁环境高准确度测试,具有十分重要的意义,相比于传统方法,计算结果更稳定、计算精度更高。
尽管上面已经示出和描述了本申请的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本申请的限制,本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (3)

1.一种信号同轴电缆传输畸变补偿方法,其特征在于,该方法包括以下步骤:
(1)测量待补偿信号并计算关键参数;其包括:
1)通过频谱测量设备测量同轴电缆输出端信号的频谱,确定同轴电缆输出端信号的最高频率Fmax和3dB带宽BW;
2)通过信号测量设备测量同轴电缆输出端信号的时域波形y,且该信号测量设备的采样率Fs由式(1)确定:
Fs≥4Fmax (1)
3)根据时域波形y得到输出端波形最小值ymin
4)根据时域波形y的持续时间T和采样率Fs,由式(2)确定采样点数N
N=T×Fs (2)
(2)根据步骤(1)所获取的结果,计算同轴电缆冲击响应函数h;
(3)对待补偿信号进行预处理以获取处理后的信号z;
(4)对所述信号z进行畸变补偿并获取畸变补偿后的信号x;其包括:
A.根据同轴电缆冲击响应函数h,获取计算矩阵H;
B.设定畸变补偿精度ε,中间参数α(k)、β(k)、τ(k),畸变补偿后的信号X(k);初始化待补偿信号向量Y,k=0时的中间参数α(0)、β(0)、τ(0),以及k=0时,畸变补偿后的信号X(0)
C.计算得到第k+1次畸变补偿后的信号X(k+1)
Figure FDA0003186030500000011
其中,HT是矩阵H的转置,I是单位矩阵;
D.依次计算第k+1次的α(k+1)、第k+1次的β(k+1)、第k+1次的τ(k+1),并带入至式(3)中;
E.根据下式,判断Xk+1是否满足要求
Figure FDA0003186030500000021
其中,ε为计算精度;若不满足,重复执行上述C、D步骤;若满足,则停止迭代计算,进入F步骤:
F.计算畸变补偿后的信号x
x=X(k+1)-D (5)
其中,D为设定的偏置信号。
2.根据权利要求1所述的信号同轴电缆传输畸变补偿方法,其特征在于,所述式(3)采用Gauss-Seidel迭代法求解。
3.根据权利要求1所述的信号同轴电缆传输畸变补偿方法,其特征在于,所述频谱测量设备采用频谱分析仪;所述信号测量设备采用示波器。
CN201911105097.1A 2019-11-13 2019-11-13 一种信号同轴电缆传输畸变补偿方法 Active CN110830081B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911105097.1A CN110830081B (zh) 2019-11-13 2019-11-13 一种信号同轴电缆传输畸变补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911105097.1A CN110830081B (zh) 2019-11-13 2019-11-13 一种信号同轴电缆传输畸变补偿方法

Publications (2)

Publication Number Publication Date
CN110830081A CN110830081A (zh) 2020-02-21
CN110830081B true CN110830081B (zh) 2021-11-09

Family

ID=69554470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911105097.1A Active CN110830081B (zh) 2019-11-13 2019-11-13 一种信号同轴电缆传输畸变补偿方法

Country Status (1)

Country Link
CN (1) CN110830081B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105785225A (zh) * 2016-03-01 2016-07-20 南方电网科学研究院有限责任公司 接收端信号时延估算误差的补偿方法
CN107817375A (zh) * 2017-10-30 2018-03-20 中国舰船研究设计中心 同轴电缆线内电磁脉冲耦合电压测试装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3000805A1 (fr) * 2013-01-04 2014-07-11 Commissariat Energie Atomique Procede d'analyse d'un cable par compensation de la dispersion subie par un signal lors de sa propagation au sein dudit cable

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105785225A (zh) * 2016-03-01 2016-07-20 南方电网科学研究院有限责任公司 接收端信号时延估算误差的补偿方法
CN107817375A (zh) * 2017-10-30 2018-03-20 中国舰船研究设计中心 同轴电缆线内电磁脉冲耦合电压测试装置

Also Published As

Publication number Publication date
CN110830081A (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN107592134B (zh) 一种面向同时同频全双工的多抽头模拟自干扰抵消方法
CN109617560A (zh) 一种iq信号校准补偿方法
CN106569074B (zh) 基于耦合装置的传导敏感度测试方法
CN106453170B (zh) 信号非线性时域测量和模拟方法及应用
US7161511B2 (en) Linearization system and method
CN111200896A (zh) 基于宽带反射系数曲线曲率分析的等离子体参数诊断方法
CN113949350A (zh) 基于基带-射频联合优化的数字预失真方法及系统
CN112213560A (zh) 一种基于z-adaline的高精度电网宽频信号测量方法
CN112415298B (zh) 一种基于时域测量的传导电磁干扰分离的实现方法
Costa et al. Access impedance in Brazilian in-home, broadband and low-voltage electric power grids
CN109583575B (zh) 基于深度学习提高仪器矢量信号分析性能的处理方法
Zietz et al. A general calibration procedure for measuring RF voltages and currents applied to the EMC analysis of automotive high-voltage power networks
CN110988490A (zh) 一种电源滤波器差损时域测量系统及方法
CN111614371B (zh) 一种无源互调电路干扰抑制设计方法
CN110830081B (zh) 一种信号同轴电缆传输畸变补偿方法
CN102095942A (zh) V型线阻抗稳定网络拓扑结构的特性参数分析及校准方法
CN111934694A (zh) 一种宽带零中频收发系统失真补偿装置
US10379162B1 (en) System for performing modulation analysis without using a modulated signal
US20020123869A1 (en) Simulation method and test arrangement for determining nonlinear signal distortion
CN111274534B (zh) 一种改进全泄漏抑制的短时dft插值方法
CN114926035A (zh) 一种电磁辐射多频钝感阻塞效应评估方法及系统
CN107831478A (zh) 一种风廓线雷达回波信号的自动定标方法和系统
CN111999587A (zh) 一种调制域网络参数测试的迭代矢量误差修正方法
CN106405341A (zh) 一种还原局部放电超高频信号的方法
CN110806555A (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
GR01 Patent grant
GR01 Patent grant