CN108414001B - 非均匀采样正弦波形失真度的确定方法 - Google Patents
非均匀采样正弦波形失真度的确定方法 Download PDFInfo
- Publication number
- CN108414001B CN108414001B CN201711101376.1A CN201711101376A CN108414001B CN 108414001 B CN108414001 B CN 108414001B CN 201711101376 A CN201711101376 A CN 201711101376A CN 108414001 B CN108414001 B CN 108414001B
- Authority
- CN
- China
- Prior art keywords
- sampling
- waveform
- sine
- fitting
- sequence
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D11/00—Component parts of measuring arrangements not specially adapted for a specific variable
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R13/00—Arrangements for displaying electric variables or waveforms
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Current Or Voltage (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明属于无线电计量测试技术领域,涉及一种非均匀采样正弦波形失真度的确定方法。
背景技术
有许多物理现象可以最终归结为正弦现象,例如振动、摆动、波动,涉及到声、光、电、无线电、力等众多学科领域物理量的信号波形,因而正弦波形是计量测试中常用的基本信号波形,在无线电测量、振动测量、声学计量、电学计量、力学计量中均有应用。
理想的正弦波形是符合数学关系表述的二维信号波形,通常一个维度是均匀流逝的时间量,而另一维度为瞬时变化的幅度信息。实际测量过程中所涉及到的正弦波形,均是含有噪声和失真的非理想波形,因而,有关正弦波失真度的测量评价成为正弦波形参量评价中的一个基本问题,失真度测量的意义和价值,不仅仅体现在失真度本身,更在于它会影响幅度、频率、相位、直流分量等另外几个正弦波形参量的测量准确度。
通常的正弦波形失真度测量有频域法和时域法两种。频域法是建立在理想的正弦波形曲线在频域是一条单一的频率谱的基础上,对所测量的正弦波形进行频域分解,将基波以外的所有频谱分量皆归结为失真,从而获得正弦波形的失真度。时域法通常是在时域对正弦波形测量序列进行最小二乘拟合,将拟合残差归结为失真,从而获得正弦波形的失真度。
以上方法仅仅局限于均匀采样的技术条件下,失真也是连续波形的情况。此时,只要量化误差足够小、采样间隔足够细密,从原理上就可以获得接近于定义的正弦波形失真。
参见图1至图3,高速采样技术的发展,使得用多个低速A/D“拼接”成单个高速A/D技术逐渐成熟并广泛应用,使得高速采样***可以获得几乎任意高的实时采样速率,而目前的绝大多数高速采样***均使用了合成A/D技术。图2采用并行合成A/D技术,图3采用串行合成A/D技术,图1是上述两种合成A/D技术的典型采样序列的示意图。其缺点是:由于拼接不完善导致的非均匀采样效应出现,使得采样序列的失真增大,用其进行正弦波形失真度分析遇到了困难。由此可见,当正弦波采样序列为非均匀采样序列时,使用以均匀采样为前提条件的频域分析方法便无法再用了,若强行使用,会带来具有不确定性的方法误差,甚至给出错误结果。而以均匀采样为前提的时域最优估计方式也将无法直接使用,由于失真在时域表现并非是“均匀”分布的,直接使用最优估计方法带来的失真测量结果也将受采样均匀程度的影响而与失真定义存在差异,从而不能进行测量应用。
发明内容
本发明的目的是:提出一种非均匀采样正弦波形失真度的确定方法,以便消除由于拼接不完善导致的非均匀采样效应造成的采样序列的失真增大,保证正弦波形失真度分析的准确度。
本发明的技术方案是:非均匀采样正弦波形失真度的确定方法,通过非均匀采样采集实际正弦波采样序列yi,实际正弦波采样序列yi中每个采样值的设定采样时刻为ti;其特征在于,波形失真度补偿的步骤如下:
1、整数周期截取:将实际正弦波采样序列yi进行数据序列截取,使其恰好含有整数个正弦周期,i是采样值的序列号,i=0,1,…,n-1;
2、标定采样时刻:对截取后的实际正弦波采样序列yi中的每个采样值的设定采样时刻ti进行精确标定,即:用实测采样时刻t′i取代每个采样值的设定采样时刻ti;将不含有实测采样时刻的实际采样序列yi转化成同时含有幅度信息和实测采样时刻的双坐标采样序列(yi,t′i);
3、确定采样正弦波形的理论函数:
采样正弦波形的理论函数为:
y(t)=Ecos(2πft+Φ)+Q [1]
式中,E是正弦波形幅度;Φ是初始相位角;f为正弦波形频率;Q是叠加到正弦波形信号上的直流分量
4、非均匀采样序列波形拟合:使用同时含有幅度信息和实测采样时刻的双坐标采样序列(yi,t′i)进行正弦曲线拟合,使公式[2]所述残差平方和ε最小,获得采样正弦波形的包括幅度、频率、相位、直流分量及拟合回归残差;
式中,A为正弦波形幅度E的拟合值,f为正弦波形频率拟合值,θ为正弦波形初始相位角Φ的拟合值,C为叠加到正弦波形信号上的直流分量Q的拟合值;
拟合残差有效值ρ为:
拟合残差εi为:
至此,获得拟合参数A、f、θ、C,得到正弦波形周期为T=1/f;
5、计算εi的加权值wi:依照实际正弦波采样序列yi中每个采样点所在实测采样时刻t′i的前后采样时刻差,按照下式计算出该采样点的加权值wi;
7、计算非均匀采样正弦波形失真度η:
至此,完成非均匀采样正弦波形失真度的确定。
本发明的优点是:提出了一种非均匀采样正弦波形失真度的确定方法,能够消除由于拼接不完善导致的非均匀采样效应造成的采样序列的失真增大,保证了正弦波形失真度分析的准确性。
附图说明
图1是合成A/D技术的典型采样序列的示意图。
图2是并行合成A/D技术的原理示意图。由于各个延迟环节的非线性规律将产生周期性非均匀采样结果。
图3是串行合成A/D技术的原理示意图。由于各个延迟环节的不均匀规律将产生随机性非均匀采样结果。
具体实施方式
下面对本发明做进一步详细说明。参见图1至图3,非均匀采样正弦波形失真度的确定方法,通过非均匀采样采集实际正弦波采样序列yi,实际正弦波采样序列yi中每个采样值的设定采样时刻为ti;其特征在于,波形失真度补偿的步骤如下:
1、整数周期截取:将实际正弦波采样序列yi进行数据序列截取,使其恰好含有整数个正弦周期,i是采样值的序列号,i=0,1,…,n-1;
2、标定采样时刻:对截取后的实际正弦波采样序列yi中的每个采样值的设定采样时刻ti进行精确标定,即:用实测采样时刻t′i取代每个采样值的设定采样时刻ti;将不含有实测采样时刻的实际采样序列yi转化成同时含有幅度信息和实测采样时刻的双坐标采样序列(yi,t′i);
3、确定采样正弦波形的理论函数:
采样正弦波形的理论函数为:
y(t)=Ecos(2πft+Φ)+Q [1]
式中,E是正弦波形幅度;Φ是初始相位角;f为正弦波形频率;Q是叠加到正弦波形信号上的直流分量
4、非均匀采样序列波形拟合:使用同时含有幅度信息和实测采样时刻的双坐标采样序列(yi,t′i)进行正弦曲线拟合,使公式[2]所述残差平方和ε最小,获得采样正弦波形的包括幅度、频率、相位、直流分量及拟合回归残差;
式中,A为正弦波形幅度E的拟合值,f为正弦波形频率拟合值,θ为正弦波形初始相位角Φ的拟合值,C为叠加到正弦波形信号上的直流分量Q的拟合值;
拟合残差有效值ρ为:
拟合残差εi为:
至此,获得拟合参数A、f、θ、C,得到正弦波形周期为T=1/f;
5、计算εi的加权值wi:依照实际正弦波采样序列yi中每个采样点所在实测采样时刻t′i的前后采样时刻差,按照下式计算出该采样点的加权值wi;
7、计算非均匀采样正弦波形失真度η:
至此,完成非均匀采样正弦波形失真度的确定。
本发明的工作原理是:针对非均匀采样序列波形,以与相邻点采样间隔在总采样长度中所占比重作为该点拟合残差在总波形失真中所占比重的加权值,由此计算非均匀采样波形序列的总失真度。技术过程如下:
针对非均匀采样正弦波形序列,首先,对其实际采样相对时刻进行标定,形成如图1所示的同时包含幅度和时刻信息的非均匀采样序列;其次,进行整数波形周期的截取,形成图1所示的包含整数个波形周期的非均匀采样序列;然后,对该波形序列进行非均匀采样最小二乘拟合,获得各个采样点的拟合残差;其后,根据各个采样点与相邻点的间隔,计算各个采样点拟合残差的加权值,并计算加权修正的拟合残差有效值最终,计算非均匀采样正弦波形失真度η。
实施例
针对如图2所示的非均匀采样序列波形,当仅有2个ADC进行采样,最终合成获得波形序列情况,其中,正弦信号峰值幅度4V,频率11.000Hz,失真度为0,每个ADC1和ADC2的采样速率均为200Sa/s,ADC1的延迟τ1=0s,ADC2的延迟τ2=1μs,由此生成非均匀采样序列。
原正弦信号为24位A/D的理想信号,失真度小于-145dB。经此非均匀采样后,若未进行非均匀采样补偿,仍然按照均匀采样序列处理,将获得其总失真度为4.8×10-5,为-68.12dB>-145dB。数据处理方法带来了额外失真。
按照本发明的非均匀采样序列处理方式,将获得其总失真度为9.3×10-8,为-154.96dB<-145dB。数据处理方法没有带来额外失真。
本发明针对非均匀采样方式的正弦波形序列,提出非均匀采样方式的正弦波形失真度的一种测量评价方法。它是一种在非均匀采样序列正弦曲线拟合基础上,寻找出一种权值计算方式,然后通过时域加权计算方法获得非均匀采样正弦波形失真度,可以在非均匀采样条件下进行正弦波形失真度的测量与评价。并且,当该非均匀采样序列变成均匀采样序列后,本发明所述方法即变成均匀采样条件下正弦波形失真度测量方法,与现有方法相一致。现有方法成为本发明所述方法的特例。
通过本发明方法,可以实现非均匀采样正弦波形失真度的测量评价,结束非均匀采样正弦波形失真度一直无法测量的问题。
本发明所述方法为软件构造式加权计算测量方法,对于硬件***没有任何额外要求,因而可以适用于任何一种波形采样***,包括均匀采样和非均匀采样***,为非均匀采样测量序列的失真测量提供了具有普适意义的方法,并且,从理论上可以推广应用到任意确定波形的非均匀采样序列失真度测量中,拥有非常广泛的应用前景。
特别说明,由人工观测,并记录观测时刻和观测结果的正弦测量数据序列,属于最常用的随机性非均匀采样序列,可以使用本发明上述方法计算其正弦波形失真。
以上所述为本发明的较佳实施例而已,本发明不应该局限于该实施例和附图所公开的内容。凡是不脱离本发明所公开的精神下完成的等效或修改,都落入本发明保护的范围。
Claims (1)
1.非均匀采样正弦波形失真度的确定方法,通过非均匀采样采集实际正弦波采样序列yi,实际正弦波采样序列yi中每个采样值的设定采样时刻为ti;其特征在于,波形失真度补偿的步骤如下:
1.1、整数周期截取:将实际正弦波采样序列yi进行数据序列截取,使其恰好含有整数个正弦周期,i是采样值的序列号,i=0,1,...,n-1;
1.2、标定采样时刻:对截取后的实际正弦波采样序列yi中的每个采样值的设定采样时刻ti进行精确标定,即:用实测采样时刻t′i取代每个采样值的设定采样时刻ti;将不含有实测采样时刻的实际采样序列yi转化成同时含有幅度信息和实测采样时刻的双坐标采样序列(yi,t′i);
1.3、确定采样正弦波形的理论函数:
采样正弦波形的理论函数为:
y(t)=Ecos(2πft+Φ)+Q [1]
式中,E是正弦波形幅度;Φ是初始相位角;f为正弦波形频率;Q是叠加到正弦波形信号上的直流分量
1.4、非均匀采样序列波形拟合:使用同时含有幅度信息和实测采样时刻的双坐标采样序列(yi,t′i)进行正弦曲线拟合,使公式[2]残差平方和ε最小,获得采样正弦波形的包括幅度、频率、相位、直流分量及拟合回归残差;
式中,A为正弦波形幅度E的拟合值,f为正弦波形频率拟合值,θ为正弦波形初始相位角Φ的拟合值,C为叠加到正弦波形信号上的直流分量Q的拟合值;
拟合残差有效值ρ为:
拟合残差εi为:
至此,获得拟合参数A、f、θ、C,得到正弦波形周期为T=1/f;
1.5、计算εi的加权值wi:依照实际正弦波采样序列yi中每个采样点所在实测采样时刻t′i的前后采样时刻差,按照下式计算出该采样点的加权值wi;
1.7、计算非均匀采样正弦波形失真度η:
至此,完成非均匀采样正弦波形失真度的确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711101376.1A CN108414001B (zh) | 2017-11-09 | 2017-11-09 | 非均匀采样正弦波形失真度的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711101376.1A CN108414001B (zh) | 2017-11-09 | 2017-11-09 | 非均匀采样正弦波形失真度的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108414001A CN108414001A (zh) | 2018-08-17 |
CN108414001B true CN108414001B (zh) | 2020-12-29 |
Family
ID=63125276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711101376.1A Active CN108414001B (zh) | 2017-11-09 | 2017-11-09 | 非均匀采样正弦波形失真度的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108414001B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115598416B (zh) * | 2022-09-16 | 2024-01-30 | 珠海多创科技有限公司 | 台区采样信号的处理方法、***、存储介质及计算机设备 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101082666B (zh) * | 2006-06-02 | 2011-07-27 | 中国科学院电子学研究所 | 基于自动测试***实现对雷达脉冲信号高精度调制的方法 |
CN102882517B (zh) * | 2012-03-30 | 2014-12-31 | 杭州电子科技大学 | 一种生成低失真度的低频正弦信号的装置及方法 |
CN102645585A (zh) * | 2012-05-14 | 2012-08-22 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种超低频正弦信号相位差的同步快速测量方法及装置 |
CN103823120A (zh) * | 2014-03-10 | 2014-05-28 | 北京航空航天大学 | 一种任意波形失真度的确定方法 |
CN104166793A (zh) * | 2014-08-06 | 2014-11-26 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种基于二次曲线拟合的冲击波形峰值测量方法 |
CN104796116B (zh) * | 2014-08-06 | 2017-08-01 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种基于参数拟合的冲击波形峰值测量方法 |
-
2017
- 2017-11-09 CN CN201711101376.1A patent/CN108414001B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108414001A (zh) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8519719B2 (en) | System for complex impedance measurement | |
Vucijak et al. | A simple algorithm for the estimation of phase difference between two sinusoidal voltages | |
US7933942B2 (en) | Low cost, high purity sign wave generator | |
Giaquinto et al. | Fast and accurate ADC testing via an enhanced sine wave fitting algorithm | |
CN106124033B (zh) | 一种激光测振校准用大触发延迟的累积校准方法 | |
US7378854B2 (en) | Dual sine-wave time stamp method and apparatus | |
CN108414001B (zh) | 非均匀采样正弦波形失真度的确定方法 | |
Pálfi | An improved sine wave histogram test method for ADC characterization | |
CN108196217B (zh) | 一种用于非车载充电机现校仪的直流计量方法及*** | |
CN108572277B (zh) | 多频信号测量方法及*** | |
JP5553334B2 (ja) | 正弦波パラメータ推定方法 | |
Bertocco et al. | Sine-fit versus discrete Fourier transform-based algorithms in SNR testing of waveform digitizers | |
Radil et al. | DSP based portable impedance measurement instrument using sine-fitting algorithms | |
Ireland et al. | A method for using Josephson voltage standards for direct characterization of high performance digitizers to establish AC voltage and current traceability to SI | |
Belega et al. | Estimation of the effective number of bits of ADCs using the interpolated DFT method | |
CN108982966A (zh) | 基于线性修正算法的谐相角分析方法 | |
Belega et al. | Choice of the window used in the interpolated discrete Fourier transform method | |
Lapuh et al. | Digital oscilloscope calibration using asynchronously sampled signal estimation | |
Raze et al. | Non coherent spectral analysis of ADC using FFT windows: An alternative approach | |
Chávez et al. | Comparison between sine wave fitting and zero-crossing methods applied to QCM impedance measurements | |
Baccigalupi et al. | Least square procedures to improve the results of the three-parameter sine fitting algorithm | |
Lapuh et al. | Measurement of repetitive arbitrary waveform RMS value | |
CN117928707A (zh) | 一种激光测振校准用大触发延迟的双频激励测量方法 | |
Radil et al. | Frequency domain parameter estimation of two common frequency single-tone signals | |
Dallet et al. | ADCs dynamic testing by multiharmonic sine fitting algorithms |
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 |