CN1682542A - 确定调制信号的包络线的方法 - Google Patents
确定调制信号的包络线的方法 Download PDFInfo
- Publication number
- CN1682542A CN1682542A CNA038213516A CN03821351A CN1682542A CN 1682542 A CN1682542 A CN 1682542A CN A038213516 A CNA038213516 A CN A038213516A CN 03821351 A CN03821351 A CN 03821351A CN 1682542 A CN1682542 A CN 1682542A
- Authority
- CN
- China
- Prior art keywords
- sample
- fourier transform
- frequency
- numerical value
- inverse
- 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.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/44—Receiver circuitry for the reception of television signals according to analogue transmission standards
- H04N5/455—Demodulation-circuits
-
- 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/20—Arrangements for detecting or preventing errors in the information received using signal quality detector
- H04L1/206—Arrangements for detecting or preventing errors in the information received using signal quality detector for modulated signals
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N17/00—Diagnosis, testing or measuring for television systems or their details
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/10—Monitoring; Testing of transmitters
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/38—Transmitter circuitry for the transmission of television signals according to analogue transmission standards
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Multimedia (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Quality & Reliability (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Measurement Of Current Or Voltage (AREA)
- Complex Calculations (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种用于确定调制输入信号(S)的包络线的方法,包括以下方法步骤:通过对输入信号(S)的数字采样(1)产生数字样本(An);通过对数字样本(An)的傅立叶变换产生傅立叶变换样本(Bn);通过从傅立叶变换样本(Bn)去除负频率范围或正频率范围,产生边频清除的傅立叶变换样本(B′n);通过对边频清除的傅立叶变换样本(B′n)的反傅立叶变换(4)产生反变换样本(Cn);以及形成反变换样本(Cn)的绝对值。
Description
技术领域
本发明涉及一种用于确定调制信号(modulated signal)的包络线(envelope curve)的方法,例如用于确定CCDF曲线的值。
背景技术
不但在确定互补累积分布函数(CCDF)中特别需要确定调制信号的包络线,而且在其它应用中也需要确定调制信号的包络线。CCDF曲线图表示所分析的信号的包络线的信号电平(signal level)超过特定电平值的概率。根据CCDF曲线轨迹,可以确定峰值因子(crest factor)等参数,峰值因子表示信号峰值处的功率与平均功率之比。峰值因子用于帮助高频调制发射机(modulated high frequency transmitter)的操作员确定发射机放大器的最佳调制。一方面,希望发射功率尽可能的高,以使接收机端的信噪比(signal-to-noise)尽可能的大。另一方面,发射功率不应过大,以避免由于发射放大器中的短功率峰值(short power peaks)而造成的损坏。如果与理想信号的轨迹一起来表示所检测的CCDF轨迹,则可以得到关于发射信号中非线性和极限效应(limitation effects)的结论。
根据德国专利DE 199 10 902 A1,可以了解用于CCDF曲线的测量值检测设备和显示设备。其中在确定调制信号的包络线或包络线的功率中包括信号处理步骤。从第10列第47行到第11列28行,提出了确定包络线功率的操作:用四倍符号频率对信号进行采样,将包含四个取样值的一组数字数值进行平方运算,求和后除以4。因此,对应于低通滤波(loss-pass filtering),生成了调制信号瞬时振幅的功率值的变化平均值(sliding average value)。但是,这种操作模式存在以下缺点:对采样得到的数字数值的必要的平方处理,导致更高频的频谱分量(higher-frequency spectral components)。随后的非理想的低通滤波导致CCDF测量的不精确。更准确地说,取样值的平方导致高频频谱分量,该高频频谱分量不能通过取平均的方式(相当于采用具有sin(x)/x的频率响应的滤波器进行滤波)正确去除。
发明内容
因此,本发明的主要目的在于提供一种用于确定调制信号的包络线的方法,能够获得相对高的计算精度。
上述目的是由权利要求1的技术特征来实现。
与现有技术相比,根据本发明,不是通过低通滤波,而是对数字样本在频域进行傅立叶变换来实现确定包络线。在频域中,就去除了频域中的正频率范围或负频率范围。然后在时域进行傅立叶反变换。接着才形成反变换样本的数值。本申请后面也指出反变换样本的绝对值表示经过调制的高频信号的包络线。
与数值形成和随后的低通滤波相比,根据本发明的方法具有如下优点:本发明方法的执行与低通滤波的质量、信号类型及其谱位置无关。另外,本发明的执行和将要测量的高频信号的同步状态也无关。而且,根据本发明的方法实质上比采用低通滤波的现有方法更加精确。
从属权利要求涉及对本发明的较佳发展。
除了去除负频率范围或正频率范围,本发明优点还在于在频域进行傅立叶变换后,还可以去除其中的直流(DC)频率0处的电平分量。因此,可以确信的是,非理想的模拟/数字转换器的直流电压偏移(direct voltageoffset)对根据本发明的方法没有影响。理想信号在中频平面内没有直接电压分量,使得直流电压分量的去除不会歪曲测量结果。
此外,对在时域中进行了反变换的样本只在有限范围进一步内进行处理是明智的,在该有限范围中,傅立叶变换和反傅立叶变换造成的信号的循环连续(cyclic continuation of the signal)得到了抑制。
权利要求6、7、8和9涉及基于本发明方法的相应数字存储介质、计算机程序或计算机程序产品。
附图说明
下面参考附图对本发明进行更详细的说明。其中:
图1为CCDF曲线的实例;
图2为根据本发明方法的方框图;
图3为解释根据本发明方法的操作模式的示意图;
图4为在频域进行傅立叶变换的采样样本;和
图5为在时域进行反变换的采样样本。
具体实施方式
下面通过应用实例说明根据本发明的方法,该应用实例用于为CCDF曲线确定包络线的瞬时功率。但是,如已描述的,根据本发明的方法不局限于该应用,而是适合于所有需要包络线的瞬时电平或如功率(即电平的平方)等信号数值的应用。
图2通过方框图的方式说明了根据本发明的方法。其中,经过调制信号调制的高频输入信号S,首先通过采样与保持电路1进行数字采样。这样就产生了输入信号S的数字样本An。然后,将样本An进行傅立叶变换,例如可用快速傅立叶变换(FFT)算法对样本An进行傅立叶变换,这样就产生了傅立叶变换样本(Fourier-transformed samples)Bn。傅立叶变换如图2的方框2所示。
由于经过采样的实信号(real signal)的傅立叶变换,如已知的那样产生了傅立叶变换样本,该样本覆盖整个负频率范围和整个正频率范围。根据本发明,负频率范围或者正频率范围都将从傅立叶变换样本Bn中去除。如果傅立叶变换样本Bn的指标n的范围是,比如从-2N/2到2N/2-1,N是自然数,那么负频率范围对应于样本Bn的n<0的部分,或者正频率范围对应于样本Bn的n>0的部分。
图2所示的B′n是去除操作后保留的样本,该样本或者只是正频率或者只是负频率。图2的方框3表示去除负频率范围的样本,方框3具有变换函数H(f),该函数只在正频率范围内与0不同。经过边频清除(sideband-cleaned)处理的傅立叶变换样本B′n,接着经过反傅立叶变换处理变换回到时域。同样地,可以使用如图2的方框4所示的快速数字傅立叶反变换(IFFT)。这样,时域中呈现反变换样本Cn,其数值最终在数值形成器5中形成。在时域中经过反变换的样本的数值如图2的Dm所示。
在应用于CCDF曲线的情况下,现在必须在方框6中建立数值样本的平方D2m与平均功率D2eff相比,在以分贝(dB)为单位的对数坐标上,超过的门限值x的相对频率,该数值样本的平方D2m对应于功率。优选地,不是在取对数处理之前而是在取对数处理之后进行平方运算,即不是用比例系数10,而是用比例系数20进行乘法运算,结果记作:
这样,CCDF曲线就可以显示在显示装置7,比如屏幕上。
如图5所示,先经过傅立叶变换然后经过时域反变换的信号包括数字样本Cn,由于最终的时域和频域采样,该信号是周期的,如图5的实例所示,该信号的周期长度为m2-m1-1。图5中,指标n的范围是从0到2N-1。因此,优选地,进一步只在有限范围13内对反变换样本Cn进行处理,以消除周期持续。也就是,采用公式Cm=Cn,m1≤m≤m2。数量数值只有由反变换样本的有限部分Cm进行计算,该反变换样本的有限部分Cm对应于图2的描述。这样,数值形成根据如下公式实现:
重复执行用于确定反变换样本的数值Dm的步骤,直到得到足够数量的数值Dm,这样可以根据已知的运算法则确定数值序列的有效值Deff。这样,该有效值的功率是用于表示CCDF曲线横轴(0dB)上的数值的参考数值。在CCDF曲线的纵轴上,属于各个功率电平的CCDF数值被绘制成曲线,也就是平均功率D2eff超过功率数值x的相对频率。这可以通过如下公式实现:
其中:
p为发生概率或者相对频率;
D为包络线的瞬时数值;
Deff为包络线的有效值。
这里,不用比较电平大小或电压大小,也自然可以知道对应的功率大小(瞬时功率D2和平均功率D2eff)是彼此直接相关的。但是,对数处理的预定系数(prefactor)已经从20变为10。
参考图3和图4更详细地说明根据本发明方法的功能。可以利用傅立叶序列对信号S进行因式分解,也就是,任何随机输入信号都可以由一系列具有不同信号电平和相位的余弦信号构成。下面仅考虑傅立叶分量中的一个,通常可以记作如下:
s1(t)=A(t)·cos(ω·t+) (4)
因此,这里所要确定的包络线为A(t)。发射信号涉及到的实信号可以复数表示如下:
这种关系可以通过向量曲线图形化表示,如图3所示。
信号s1(t)包括以角频率ω向左旋转的第一旋转向量8(first rotatingvector)和与此同步地以同样角频率ω向右旋转的第二旋转向量9。根据本发明,负频率范围的删除将导致旋转向量9的消除。反之,仅可能作为一种替换,正频率范围的删除将导致旋转向量8的消除。因此,频域中的滤波将导致公式(5)两项中的一项被删除。例如,如果公式(4)中的负频率分量,即图3中向左旋转的旋转向量9被删除,则进行数量形成之后将产生如下结果:
根据图3,该数值对应于保留向量的长度。当使用信号s2(t)确定CCDF曲线时,由于数值形成造成的s2(t)只能是正的事实就不重要了。在CCDF曲线中,只有正功率可以互相比较。除以系数2的运算也同样不会影响CCDF曲线的结果。
根据傅立叶分量得到的结论当然可以应用于表示高阶傅立叶分量的线性叠加的全部信号。为此,图4中表示了傅立叶变换样本Bn。这里,指标n从-2N/2到2N/2-1。可以看出,实输入信号S中,负频率范围10是正频率范围11的镜像。
如果在进一步的信号处理中,负频率范围10被删除,即:
B′n=0,对于n<0(for n<0)并且
Bn=B′n,对于n>0(for n>0),
或者正频率范围11被删除,即:
B′n=B′n,对于n<0(for n<0)并且
B′n=0,对于n>0(for n>0),
那么如前面图3所示的那样,在绝对值形成和时域反变换之后,就会自动生成包络线。
优选地,不仅将负频率范围10或者正频率范围11消除,而且也将0频率的电平分量12消除;这里使用指标中,即n=0时的B0。因此,可以消除可能出现的直流电压分量(DC偏移)。由于被评估的信号来自于中间频率平面,因此这些信号实际上并不包含直流电压分量。但是,如果出现直流电压分量,则这可能来自于模拟-数字转换器的直流电压偏移,这种直流电压分量的去除将增加测量的精度。
图1示出了CCDF曲线的一个例子,主要是利用本发明的方法获得其中的包络线。为此,如已描述的那样,在CCDF曲线中按照超过对数坐标上的特定电平D的方式,将相对频率p绘制成曲线。图3所示的实例中,根据8VSB的标准对输入信号进行数字调制,仍会出现有效功率超过3dB时,具有大约10%的相对频率,而出现有效功率超过6dB时,则具有比1%还小的相对频率。
如已经多次提到的那样,根据本发明的方法不局限于用于为CCDF曲线确定瞬时电平值或者瞬时功率值的应用,而一般也适用于确定调制信号的包络线。该方法可以采用数字硬件执行,比如使用现场可编程门阵列(FPGA),或者采用特殊处理器中的软件执行,理想地,采用数字信号处理器(DSP)。
Claims (9)
1、一种用于确定调制输入信号S的包络线的方法,包括以下方法步骤:
通过对输入信号S的数字采样产生数字样本An;
通过对数字样本An的傅立叶变换产生傅立叶变换样本Bn;
通过从傅立叶变换样本Bn去除负频率范围或正频率范围,产生边频清除的傅立叶变换样本B′n;
通过对边频清除的傅立叶变换样本B′n的反傅立叶变换产生反变换样本Cn;和
形成反变换样本Cn的绝对值Dm的数值。
2、根据权利要求1所述的方法,其特征在于,
为产生边频清除的傅立叶变换样本B′n,除了去除负频率或正频率范围之外,还去除0频率处的电平分量。
3、根据权利要求1或2所述的方法,其特征在于,
进一步仅在有限范围中处理反变换样本Cn,以使得由傅立叶变换和反傅立叶变换引起的周期持续被消除。
4、根据权利要求1至3任一项所述的方法,其特征在于,
绝对值Dm的数值相对于反变换样本的有效值Deff进行取对数处理。
5、根据权利要求4所述的方法,其特征在于,
取对数处理的数值的频率分布作为对数化电平CCDF曲线的函数显示。
6、一种具有电子可读控制信号的数字存储介质,其与可编程计算机或数字信号处理器配合,执行根据权利要求1至5任一项的方法。
7、一种具有程序编码工具的计算机程序产品,其存储在机器可读的载体中,当在计算机或数字信号处理器上运行该程序时,执行根据权利要求1至5任一项的所有步骤。
8、一种具有程序编码土具的计算机程序,当在计算机或数字信号处理器上运行该程序时,执行根据权利要求1至5任一项的所有步骤。
9、一种具有程序编码工具的计算机程序,该程序存储在机器可读的数据载体中,并且执行根据权利要求1至5任一项的所有步骤。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE10242333A DE10242333A1 (de) | 2002-09-12 | 2002-09-12 | Verfahren zum Bestimmen der Hüllkurve eines modulierten Signals |
DE10242333.4 | 2002-09-12 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1682542A true CN1682542A (zh) | 2005-10-12 |
CN100473175C CN100473175C (zh) | 2009-03-25 |
Family
ID=31724695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB038213516A Expired - Fee Related CN100473175C (zh) | 2002-09-12 | 2003-08-07 | 确定调制信号的包络线的方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US7424404B2 (zh) |
EP (1) | EP1537750B1 (zh) |
JP (1) | JP4242346B2 (zh) |
CN (1) | CN100473175C (zh) |
DE (2) | DE10242333A1 (zh) |
DK (1) | DK1537750T3 (zh) |
WO (1) | WO2004034712A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI411788B (zh) * | 2005-11-04 | 2013-10-11 | Tektronix Inc | 時間任意信號功率統計測量裝置及方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005002195A1 (de) * | 2005-01-17 | 2006-07-27 | Siemens Ag | Verfahren und Anordnung zur Regeneration eines optischen Datensignals |
US8130632B2 (en) | 2007-04-06 | 2012-03-06 | Panasonic Corporation | Transmitter and SSB signal generation method |
US9159325B2 (en) * | 2007-12-31 | 2015-10-13 | Adobe Systems Incorporated | Pitch shifting frequencies |
WO2014032698A1 (en) * | 2012-08-28 | 2014-03-06 | Rohde & Schwarz Gmbh & Co. Kg | Method and device for measuring weak inband interference |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2830837C2 (de) * | 1977-07-15 | 1983-06-09 | Tokyo Shibaura Denki K.K., Kawasaki, Kanagawa | Verfahren und Vorrichtung zur Messung der Kennwerte eines Lautsprechers |
DE3338193C2 (de) * | 1983-10-20 | 1987-05-14 | Institut für Rundfunktechnik GmbH, 8000 München | Schaltungsanordnung zum Messen des Quantisierungsrauschens oder des Leerkanalrauschens oder der Gesamtverzerrung eines digital arbeitenden Audiogerätes |
US4827516A (en) * | 1985-10-16 | 1989-05-02 | Toppan Printing Co., Ltd. | Method of analyzing input speech and speech analysis apparatus therefor |
US5469127A (en) * | 1992-08-04 | 1995-11-21 | Acrodyne Industries, Inc. | Amplification apparatus and method including modulator component |
US5956372A (en) * | 1994-03-17 | 1999-09-21 | Digital Compression Technology, L.P. | Coding system for digital transmission compression |
US20010055320A1 (en) * | 1994-12-15 | 2001-12-27 | Pierzga Wayne Francis | Multiplex communication |
US6272172B1 (en) * | 1998-03-31 | 2001-08-07 | Tektronix, Inc. | Measurement acquisition and display apparatus |
US6314146B1 (en) * | 1998-06-05 | 2001-11-06 | The Board Of Trustees Of The Leland Stanford Junior University | Peak to average power ratio reduction |
US6433558B1 (en) * | 1999-05-13 | 2002-08-13 | Microtest, Inc. | Method for diagnosing performance problems in cabling |
DE60027123T2 (de) * | 1999-12-28 | 2006-09-07 | Ntt Docomo Inc. | Schaltung zur Komprimierung des dynamischen Bereich eines Signals |
US6323412B1 (en) * | 2000-08-03 | 2001-11-27 | Mediadome, Inc. | Method and apparatus for real time tempo detection |
DE10130943A1 (de) * | 2001-05-18 | 2002-11-21 | Rohde & Schwarz | Signalgenerator mit Anzeigeeinrichtung |
-
2002
- 2002-09-12 DE DE10242333A patent/DE10242333A1/de not_active Withdrawn
-
2003
- 2003-08-07 WO PCT/EP2003/008792 patent/WO2004034712A1/de active IP Right Grant
- 2003-08-07 CN CNB038213516A patent/CN100473175C/zh not_active Expired - Fee Related
- 2003-08-07 US US10/525,638 patent/US7424404B2/en not_active Expired - Fee Related
- 2003-08-07 DE DE50301766T patent/DE50301766D1/de not_active Expired - Lifetime
- 2003-08-07 DK DK03807783T patent/DK1537750T3/da active
- 2003-08-07 EP EP03807783A patent/EP1537750B1/de not_active Expired - Lifetime
- 2003-08-07 JP JP2004542297A patent/JP4242346B2/ja not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI411788B (zh) * | 2005-11-04 | 2013-10-11 | Tektronix Inc | 時間任意信號功率統計測量裝置及方法 |
Also Published As
Publication number | Publication date |
---|---|
DE50301766D1 (de) | 2005-12-29 |
WO2004034712A1 (de) | 2004-04-22 |
DK1537750T3 (da) | 2006-04-03 |
JP4242346B2 (ja) | 2009-03-25 |
US7424404B2 (en) | 2008-09-09 |
CN100473175C (zh) | 2009-03-25 |
JP2006502631A (ja) | 2006-01-19 |
EP1537750B1 (de) | 2005-11-23 |
US20060155493A1 (en) | 2006-07-13 |
EP1537750A1 (de) | 2005-06-08 |
DE10242333A1 (de) | 2004-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Peetz | Dynamic testing of waveform recorders | |
US7161511B2 (en) | Linearization system and method | |
US8843335B2 (en) | Wavelet denoising for time-domain network analysis | |
US6760674B2 (en) | Audio spectrum analyzer implemented with a minimum number of multiply operations | |
CN1812295A (zh) | 一种3g终端的射频一致性测试方法 | |
CN111046791A (zh) | 一种基于含可变因子广义s变换的电流信号滤波去噪方法 | |
CN111693136B (zh) | 一种采用回波信号自相关相位谱的声表面波谐振器频率估计算法 | |
WO2024109031A1 (zh) | 一种基于hht变换的电压质量扰动检测方法 | |
CN1682542A (zh) | 确定调制信号的包络线的方法 | |
JPWO2002103377A1 (ja) | ジッタ測定装置、及びジッタ測定方法 | |
CN106375039B (zh) | 一种提高矢量网络分析仪接收机动态范围的方法 | |
CN117368774A (zh) | 一种基于阻抗谱重构技术的锂离子电池健康状态估计方法 | |
CN117110850A (zh) | 一种基于st2500测试机的mcu片上adc测试方法 | |
CN116046968A (zh) | 一种液相色谱工作站数据处理方法、系统及可存储介质 | |
CN110070886A (zh) | 噪声特征分析方法及噪声特征分析系统 | |
CN113114245B (zh) | Adc参数测试中超量程输入与非相干采样下的信号恢复方法 | |
CN2700889Y (zh) | 实时电能质量综合分析仪 | |
CN115290919A (zh) | 一种非常规的汽轮机转速测量方法 | |
CN100346171C (zh) | 一种核磁共振成像信号的接收方法 | |
Nunzi et al. | A procedure for highly reproducible measurements of ADC spectral parameters | |
CN107180636B (zh) | 一种基于安卓系统傅立叶变换的音频数据获取方法 | |
CN110442989A (zh) | 一种失真动态范围性能指标的计算方法 | |
US6873923B1 (en) | Systems and methods for performing analysis of a multi-tone signal | |
CN1844940A (zh) | 一种测试adsl线路差模噪声的装置及其方法 | |
Chornik et al. | Deconvolution in electron spectroscopy revisited: computational aspects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090325 Termination date: 20100807 |