CN103532906B - Digital modulation error parameter metering method and system based on waveform design - Google Patents

Digital modulation error parameter metering method and system based on waveform design Download PDF

Info

Publication number
CN103532906B
CN103532906B CN201310499303.8A CN201310499303A CN103532906B CN 103532906 B CN103532906 B CN 103532906B CN 201310499303 A CN201310499303 A CN 201310499303A CN 103532906 B CN103532906 B CN 103532906B
Authority
CN
China
Prior art keywords
sequence
symbol
baseband
waveform
error
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
CN201310499303.8A
Other languages
Chinese (zh)
Other versions
CN103532906A (en
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.)
China Academy of Information and Communications Technology CAICT
Original Assignee
Ministry Of Industry & Information Technology Telecommunication Metrology Center
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 Ministry Of Industry & Information Technology Telecommunication Metrology Center filed Critical Ministry Of Industry & Information Technology Telecommunication Metrology Center
Priority to CN201310499303.8A priority Critical patent/CN103532906B/en
Publication of CN103532906A publication Critical patent/CN103532906A/en
Application granted granted Critical
Publication of CN103532906B publication Critical patent/CN103532906B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses a digital modulation error parameter metering method and a system based on waveform design. The metering method comprises the following steps that a, beginning and end continuity of digital modulation signal baseband waveform with finite length is realized; b, setting of a digital modulation error is realized based on error symbol design; and c, an intact error symbol frame is formed. The method and the system achieve metering of digital modulation error parameters such as EVM (Error Vector Magnitude); an experiment indicates that the metering method and the system have the following characteristics that 1, continuous adjustability of a digital modulation error set value within a larger range can be realized; expected digital modulation error parameters can be set; 2, measured values of the digital modulation parameters such as the EVM are highly identical with expected calculated values; 3, readings of the measured values of the digital modulation parameters have very good repeatability; and 4, the digital modulation error parameter metering method and the system are suitable for almost all digital modulation modes, such as 256 QAM (Quadrature Amplitude Modulation) and higher-order modulation modes.

Description

Digital modulation error parameter metering method based on Waveform Design and system
Technical field
The present invention relates to digital communication technology field, integrated wireless communication tester or VSA are carried out Digital modulation error parameter metering method based on Waveform Design and system.
Background technology
Digital modulation signals are the important carriers of present information transmission, it has also become one of foundation stone of information-intensive society, such as no Line communication system just using multiple digital modulation, such as gmsk, qpsk, 16qam, 64qam, 256qam etc..Wherein digital modulation is missed Difference parameter includes Error Vector Magnitude root-mean-square value (hereinafter referred to as evmrms), Error Vector Magnitude peak value (hereinafter referred to as Evmpeak), range error root-mean-square value (hereinafter referred to as magerrrms), range error peak value be (hereinafter referred to as Magerrpeak), phase error root-mean-square value (hereinafter referred to as phaseerrrms), phase error peak value be (hereinafter referred to as Phaseerrpeak) it is the important errors parameter weighing digital modulation signals error.The same with all of physical quantity, these errors Parameter also has to pass through metering that is believable, can tracing to the source, and guarantee value is accurately and unified, and guarantee information system is just Often operate.However, still suffering from following problem in current digital modulation error parameter metering:
1st, closed loop is mutually measured it is difficult to trace to the source.This is primarily referred to as: according to digital modulation error parameter metering side both domestic and external at present Method (such as jjf1128-2004 VSA calibrating standard or jjf1174-2007 digital signal generator calibrating standard The digital modulation error parameter metering method of regulation), on the one hand, digital signal generator is with VSA (below Abbreviation vsa) come to calibrate, that is, each error parameter of the generated digital modulation signals of digital signal generator is not known , but measured and analyzed by vsa and draw;On the other hand, digital signal generator is then based on as standard to the calibration of vsa Device.Above-mentioned metering process is substantially the process that a closed loop is mutually measured, and is the process that a closed loop that can not trace to the source is mutually measured, and is In measuring system " castles in the air ", as shown in Figure 1.That is, the error that current digital modulation error parameter metering obtains Parameter can not be traced to the source, and this is it is difficult to ensure that the accuracy of metering process and credibility.
2nd, the digital modulation signals as calibration vsa standard lack the setting of error parameter.In systems in practice, due to many The interference of the factor of kind, often there is error in digital modulation signals.Therefore, numeral tune can accurately be measured in order to ensure vsa Each error parameter of signal processed, the error parameter as the digital modulation signals of calibration vsa standard within the specific limits should When can arrange.One of effect of vsa is exactly the evm value measuring in relative broad range, and such as td-lte standard requires 16qam The evmrms < 8% of evmrms < the 12.5%.64qam signal of signal, should arrange different evm values at least in the range of this. As the range of certain slide calliper rule is 10 centimetres, and only make the gauge block of 0.1 centimetre of length remove to calibrate slide calliper rule, be inadequate.China is studied Personnel have indicated Similar Problems in pertinent literature.However, the calibration steps of Contemporary Digital modulation error parameter is but not Error parameter is configured, this is an off real requirement.
Content of the invention
For solving technical problem present in above-mentioned prior art, realize the purpose of the present invention, the present invention proposes a kind of base Digital modulation error parameter metering method and system in Waveform Design.
The method comprises the following steps: a, realizes the continuous step of finite length digital modulation signals baseband waveform head and the tail; B, the basis based on error symbol design, the step realizing the setting of digital modulation error;C, form complete error symbolic frame Step.
Wherein, step a refers under the hardware structure of vector modulation signal generator, the iq of one finite duration of design Baseband waveform, is then introduced into the baseband portion of signal generator, and this limited will be held in the way of loop play by signal generator The baseband waveform of continuous time carries out head and the tail and connects, thus it is continuous to realize finite length digital modulation signals baseband waveform head and the tail.Enter One step ground, realizes step a, can be in three kinds of modes below including:
(1) step a comprises the steps:
Step 1: construction symbol sebolic addressing, it is l symbol period that the time domain impulse of baseband filter blocks length, and construction is not filtered The rear l symbol of the front l symbol of ripple symbol sebolic addressing and sequence constitutes symmetrical relationss;
Step 2: pretreatment, according to the baseband waveform sample rate setting and baseband signalling speed, by unfiltered iq symbol Sequence and baseband filter shock response waveform carry out suitable interpolation and process, and make two waveform corresponding time-domain sampling interval phases Deng;
Step 3: convolution, pretreated for step 2 two values sequence is done convolution, participates in the iq symbol of convolution The sampling number of number sequence is n, and corresponding baseband filter shock response waveform sampling points are g, the points of convolution results It is n+g-1;
Step 4: intercept, the convolution results of step 3 are intercepted, intercept the n numerical point in stage casing, intercepting result is Required baseband waveform sequence of values, this sequence has the continuous characteristic of head and the tail.
(2) step a comprises the steps:
Step 1: pretreatment, according to the baseband waveform sample rate setting and baseband signalling speed, by unfiltered iq symbol Sequence and baseband filter shock response waveform carry out suitable interpolation and process and null filling, make the corresponding time domain of two waveforms Sampling interval is equal, and counts and equal be n;
Step 2: fast Fourier transform, pretreated for step 1 two values sequence is in quick Fu of n point respectively Leaf transformation, forms the sequence of two new n ' points: xk1 and xk2, then xk1 with xk2 execution corresponding element order is multiplied one by one Computing, its result yk is also the sequence of values of a n point;
Step 3: an inverse fast fourier, an inverse fast fourier computing is executed to yk, forms the ripple of a time domain Figurate number value sequence s, this sequence is required baseband waveform sequence of values, and this sequence has the continuous characteristic of head and the tail.
(3) step a comprises the steps:
Step 1: sequence repeats, initial iq symbol sebolic addressing comprises m symbol, repeats this sequence, forms new sequence;
Step 2: pretreatment, according to the baseband waveform sample rate setting and baseband signalling speed, by new sequence and base band Filter impulse response wave shape carries out suitable interpolation and processes, or adjustment sample rate, makes the corresponding time-domain sampling of two waveforms Interval is equal;
Step 3: convolution, pretreated for step 1 two values sequence is done convolution computing, participates in convolution The sampling number of iq symbol sebolic addressing is 2n, and corresponding baseband filter shock response waveform sampling points are g, convolution results sequence Represented with ck1, its points are 2n+g-1;
Step 4: intercept, convolution results ck1 of step 2 are intercepted, from selected starting point, continuous n after intercepting The sequence that individual numerical point is constituted, this intercepting result ck2 is required baseband waveform sequence of values, and it is continuously special that this sequence has head and the tail Property, above-mentioned selected starting point in the selection range in ck1 is: from the beginning of the g+1 point, to the n-th -2 points, all optional.
Wherein, step c further includes following concrete steps:
Step 1: construction frame, list the sequence of complex numbers s of one group of standard symbol1、s2、s3…sk…sk, for each symbol sk, According to step b, design corresponding l error symbol sk1、sk2…skl, form complete error symbolic frame, this symbolic frame be framee;
Step 2: combined frames, according to the requirement of actual test degree symbol sebolic addressing length, by multiple complete error symbolic frames Framee combines, and forms a new sign matrix framem, sign matrix framem restructuring is become the plural number of single file Rearrangement is upset in the position of each symbol then by sign magnitude sequence at random, forms symbol sebolic addressing sframe
Step 3: construction waveform, according to step a, based on symbol sebolic addressing sframeGenerate the continuous baseband signal waveform of head and the tail wbframeIt is ensured that baseband signal waveform wbframeIn corresponding symbol sebolic addressing, each symbol that symbolic frame framee comprises is strict Equiprobability occurs;
Step 4: generate rf modulated signal, by baseband waveform wbframeCarry out proper treatment and be allowed to adaptive vector modulation letter The operation of number generator requires, and baseband waveform is imported the baseband portion of vector modulation signal generator, and then produces radio frequency and adjust Signal wr processedframe, VSA vsa tested for the input of this radiofrequency signal is demodulated;
Step 5: setting VSA, wb is setframeCorresponding symbol total number is p, by tested vector letter The symbolic analysis length of number analyser vsa is set to the integral multiple of p.
The present invention also proposes a kind of digital modulation error parameter metering system based on Waveform Design, comprising: for realizing The device of described step a;For realizing the device of described step b;For realizing the device of described step c.
Brief description
A kind of Fig. 1 flow process of the digital modulation error parameter metering method based on Waveform Design provided in an embodiment of the present invention Figure
The metering method of the digital modulation parameter of the closed loop that can not trace to the source in Fig. 2 prior art
Fig. 3 a unfiltered sign magnitude sequence
Fig. 3 b baseband filter shock response waveform
Fig. 4 do not take steps a process baseband filtering after the 16qam signal baseband waveform of finite duration head Tail junction pattern
The corresponding Error Vector Magnitude of each symbol (evm) after the discontinuous waveform loop play of Fig. 5 head and the tail
The construction of the iq symbol sebolic addressing of Fig. 6 the first implementation of step a
The initial part of Fig. 7 iq symbol sebolic addressing
The initial part of the new iq symbol sebolic addressing of Fig. 8 the third implementation of step a
Fig. 9 is using the head and the tail junction of the 64qam signal baseband waveform after the continuous scheme of waveform head and the tail
Figure 10 is using the corresponding error of each symbol after the 64qam baseband signal loop play after the continuous scheme of waveform head and the tail Amplitude of the vector (evm)
Digital modulation error coordinate points setting in the corresponding symbol judgement domain of Figure 11 a point
Axial equal proportion equilateral triangle design in Figure 12 iq plane
The complete error symbol frame of Figure 13 mono-: framee
Figure 14, according to 16qam modulation system, is provided with the vsa demodulated methed vectogram of the digital modulation signals of standard error
Figure 15 is according to expected setting evmrms value (expected evmrms) of 16qam modulation system, value of calculation (calculated) and measurement result (measured) relation;
Wherein, top figure: measure the evmrms obtaining;Middle figure: evmrms measured value and the difference calculating desired value, it is expressed as Measurement error (measurement error);Base map: the standard deviation (reading time interval 0.1 second) of survey measurementss, in figure evmrms refers to Error Vector Magnitude root-mean-square value;
Figure 16 is according to expected setting magerrrms value (expected magerrrms) of 16qam modulation system, value of calculation (calculated) and measurement result (measured) relation,
Wherein, top figure: measure the magerrrms obtaining;Middle figure: magerrrms measured value and the difference calculating desired value, It is expressed as measurement error (measurement error;Base map: the standard deviation (reading time interval 0.1 second) of survey measurementss, in figure Magerrrms fingerbreadth degree error mean square root;
Figure 17 is according to expected setting phaseerrrms value (expected phaseerrrms) of 16qam modulation system, meter Calculation value (calculated) and the relation of measurement result (measured),
Wherein, top figure: measure the phaseerrrms obtaining;Middle figure: phaseerrrms measured value and calculating desired value Difference, is expressed as measurement error (measurement error);Base map: the standard deviation (reading time interval 0.1 second) of survey measurementss, In figure phaseerrrms refers to phase error root-mean-square value;
Figure 18 is according to 64qam modulation system, " equilateral triangle radius " γ (linear dimension) and measurement in error setting model The relation of result,
Wherein, top figure: measure the evmrms obtaining;Middle figure: evmrms measured value (measured) and calculating desired value (calculated) difference, is expressed as measurement error (measurement error);Base map: the standard deviation of survey measurementss (is read Number interval 0.1 second), in figure evmrms refers to Error Vector Magnitude root-mean-square value;
Figure 19, according to 64qam modulation system, is provided with the vsa demodulated methed vector of the digital modulation signals of standard error Figure;
Figure 20 is according to 256qam modulation system, " equilateral triangle radius " γ (linear dimension) and measurement in error setting model The relation of result,
Wherein, top figure: measure the evmrms obtaining;Middle figure: evmrms measured value and the difference calculating desired value, it is expressed as Measurement error (measurement error);Base map: the standard deviation (reading time interval 0.1 second) of survey measurementss, in figure evmrms refers to Error Vector Magnitude root-mean-square value;
Figure 21 is according to 256qam modulation system, " equilateral triangle radius " γ (linear dimension) and measurement in error setting model The relation of result,
Wherein, top figure: measure the magerrrms obtaining;Middle figure: magerrrms measured value and the difference calculating desired value, It is expressed as measurement error (measurement error);Base map: the standard deviation (reading time interval 0.1 second) of survey measurementss, in figure Magerrrms fingerbreadth degree error mean square root;
Figure 22 is according to 256qam modulation system, " equilateral triangle radius " γ (linear dimension) and measurement in error setting model The relation of result,
Wherein, top figure: measure the phaseerrrms obtaining;Middle figure: phaseerrrms measured value and calculating desired value Difference, is expressed as measurement error (measurement error);Base map: the standard deviation (reading time interval 0.1 second) of survey measurementss, In figure phaseerrrms refers to phase error root-mean-square value.
Specific embodiment
The embodiment provides a kind of digital modulation error parameter metering method based on Waveform Design.As Fig. 1 Shown, the method comprises the following steps: a, realizes the continuous step of finite length digital modulation signals baseband waveform head and the tail;B, base In the basis of error symbol design, the step realizing the setting of digital modulation error;C, the step forming complete error symbolic frame.
First, finite length digital modulation signals baseband waveform head and the tail are realized continuous
The iq base band of a finite duration under the hardware structure of existing vector modulation signal generator, can be designed Waveform, is then introduced into the baseband portion of signal generator, and signal generator is by when in the way of loop play by this finite-duration Between baseband waveform carry out head and the tail connect, thus constantly modulate carrier wave, formed rf modulated signal.
Form the sequence of values of the iq baseband waveform of a finite duration using the mode of Digital Signal Processing, existing Three steps are typically had in technology:
1st, define the sequence of values of unfiltered symbol (i road and q road define respectively).Shown in (a) of following Fig. 3, in order to It is easy to represent, Fig. 3 (a) illustrate only the sequence on i road, but the symbol sebolic addressing participating in computing comprises i road and q road, typically with multiple The form of number real part and imaginary part represents, no longer individually illustrates below.
2nd, define the time domain impulse response sequence of values of baseband filter, such as conventional raised cosine filter, root liter are remaining String wave filter, Gaussian filter etc..
3rd, the time domain impulse response wave shape of unfiltered sign magnitude sequence and baseband filter does convolution in time domain, realizes filter Ripple, that is, define " the iq baseband waveform of finite duration ".
But in the case of not taking special handling, the waveform that above-mentioned convolution obtains has disadvantages in that waveform sequence Row first and last numerical value are discontinuous, so due to being to connect to come by head and the tail in the internal waveform of signal generator Realize loop play, then the waveform in head and the tail junction will be discontinuous, as shown in Figure 4.
As shown in figure 4, baseband waveform will lead to the extra error of signal in the discontinuous of head and the tail junction, such as measurement figure The Error Vector Magnitude (evm) on every symbol of signal shown in 4, result such as Fig. 5.
Analysis Fig. 5 understands, due to discontinuous in head and the tail junction, leads to head and the tail junction to characterize the evm of signal errors Value significantly increases, and obvious error and uprushes.Obviously, such signal waveform cannot function as calibrating the standard signal of vsa.
Only achieve baseband waveform head and the tail continuously, signal is possible to the metering for digital modulation error parameter, is Realize that baseband waveform head and the tail are continuous, the present invention proposes 3 kinds of technical schemes that can independently use.
Waveform head and the tail continuously realize technical scheme 1:
Step 1: construction symbol sebolic addressing.Frequency domain filtering is equivalent to convolution, and the time domain impulse of such as baseband filter blocks Length is l symbol period, then the rear l symbol of the front l symbol and sequence that construct unfiltered symbol sebolic addressing constitutes symmetrical closing System.As shown in Figure 6.
Step 2: pretreatment.According to the baseband waveform sample rate setting and baseband signalling speed, by unfiltered iq symbol Baseband filter shock response waveform shown in sequence and Fig. 3 (b) carries out suitable interpolation and processes, when making two waveforms corresponding The domain sampling interval is equal.
Step 3: convolution.Two values sequence pretreated described in step 2 is done convolution.If participating in volume The sampling number of long-pending iq symbol sebolic addressing is n, and corresponding baseband filter shock response waveform sampling points are g, then volume The points of long-pending result are n+g-1.
Step 4: intercept.Convolution results described in step 3 are intercepted, only intercepts the n numerical point in stage casing.Should Intercepting result is exactly the final baseband waveform sequence of values needing, and this sequence has the continuous characteristic of head and the tail.
Waveform head and the tail continuously realize technical scheme 2:
Step 1: pretreatment.According to the baseband waveform sample rate setting and baseband signalling speed, by unfiltered iq symbol Baseband filter shock response waveform shown in sequence and Fig. 3 (b) carries out suitable interpolation and processes and null filling, makes two ripples Shape corresponding time-domain sampling interval is equal, and counts and equal be n.
Step 2: fast Fourier transform.Pretreated for step 1 two values sequence is in quick Fu of n point respectively Leaf transformation (fft), this is a kind of digital signal processing algorithm of standard.Define the sequence of two new n points: xk1 and xk2, Then to xk1 and xk2 execution corresponding element order multiplication operation one by one, its result yk is also the sequence of values of a n point.
Step 3: an inverse fast fourier.To yk execution an inverse fast fourier computing (ifft), when forming one It is simply that the final baseband waveform sequence of values needing, this sequence has the continuous characteristic of head and the tail for wave numerics sequence s in domain.
Waveform head and the tail continuously realize technical scheme 3:
Step 1: sequence repeats.Iq symbol sebolic addressing as initial comprises m symbol, can be expressed as form shown in Fig. 6.So After repeat this sequence, form new sequence, as shown in Figure 8.
Step 2: pretreatment.According to the baseband waveform sample rate setting and baseband signalling speed, by the sequence of iq symbol shown in Fig. 8 Row and the baseband filter shock response waveform shown in Fig. 3 (b) carry out suitable interpolation and process, or adjustment sample rate, make two Individual waveform corresponding time-domain sampling interval is equal.
Step 3: convolution.Two values sequence pretreated described in step 2 is done convolution computing.If ginseng Sampling number with the iq symbol sebolic addressing of convolution is 2n, and corresponding baseband filter shock response waveform sampling points are g, that Convolution results sequence can be represented with ck1, and its points are 2n+g-1.
Step 4: intercept.Convolution results ck1 described in step 3 are intercepted, can connect after intercepting from selected starting point The sequence that n continuous numerical point is constituted, this intercepting result ck2 is exactly the final baseband waveform sequence of values needing, and this sequence has There is the continuous characteristic of head and the tail.Above-mentioned selected starting point in the selection range in ck1 is: from the beginning of the g+1 point, to the n-th -2 Point, all optional.
Above three scheme is independent, and is all effective, after the continuous technical scheme of head and the tail waveform, Ke Yida Shown in effect following Fig. 9, the Figure 10 arriving.
Contrast Fig. 9 and Fig. 4, Figure 10 and Fig. 5 is it is clear that have: the continuous scheme of waveform head and the tail is effective, thus determines this One basic technology of invention.
2nd, the digital modulation error setting on the basis based on error symbol design
The symbol sebolic addressing designing digital modulation error it is necessary to (undistorted error free) of establishing criteria goes design error to accord with Number.The sequence of complex numbers assuming the corresponding standard symbol of certain digital modulation mode is s1、s2、s3…sk…skIt is clear that each is multiple Coordinate points in the corresponding iq plane of numerical symbol.Generally, it is assumed that skCorresponding a point in iq plane is logical according to classical numeral Letter is theoretical, there is a region and be referred to as " symbol s in iq planekJudgement domain ", if the sample vector that obtains of judgement falls into This region, all adjudicates as symbol sk, as schematic diagram, " symbol skJudgement domain " following Figure 11 gray area shown in, then Point beyond several a points of design in judgement domain, these points are labeled as a1、a2…al, their corresponding complex values are followed successively by sk1、sk2…skl.So just can count in iq plane, using the definition of the digital modulation error parameter such as geometrical relationship and evm Calculate this setting corresponding digital modulation error parameter, thus realizing the setting of digital modulation error.L >=1, for different symbols Number, l can equal it is also possible to.
Iq coordinate points a1、a2…alGeometric pattern design can be varied, such as triangle, pentagon, hexagon etc. Deng.The present invention provides one kind and preferably designs example, and this method can be referred to as " axial equal proportion equilateral triangle design ", As shown in figure 12.
In Figure 12, a point is the corresponding iq coordinate points of standard symbol, if a point complex coordinates are: i+jq, vector's Amplitude is ma, phase angle is θa, then adjudicate a of setting in domain in a point1、a2、a33 points of phasor coordinate is respectively as formula (1)~(3) Shown.
oa 1 &rightarrow; = m a ( 1 - &gamma; ) exp ( j&theta; a ) - - - ( 1 )
oa 2 &rightarrow; = oa 1 &rightarrow; + 3 &gamma; m a exp [ j ( &theta; a - &pi; 6 ) ] - - - ( 2 )
oa 3 &rightarrow; = oa 1 &rightarrow; + 3 &gamma; m a exp [ j ( &theta; a + &pi; 6 ) ] - - - ( 3 )
Can the γ in (1)~(3) be called " equilateral triangle radius ", due to the geometrical property of equal proportion, can be by such as Lower simplified operation is calculating corresponding digital modulation error.If the corresponding complex coordinates of a point are 1, then have:
oa 1 &rightarrow; = 1 - &gamma; - - - ( 4 )
oa 2 &rightarrow; = 1 + &gamma; 2 + j 3 &gamma; 2 - - - ( 5 )
oa 3 &rightarrow; = 1 + &gamma; 2 - j 3 &gamma; 2 - - - ( 6 )
According to symmetry, reference constellation point a ' should be on the extended line of oa, and the computational methods according to root-mean-square value obtain:
oa &prime; &rightarrow; = m a &prime; = 1 + &gamma; 2 - - - ( 7 )
phaseerrpeak = arctan ( 3 &gamma; 2 1 + &gamma; 2 ) - - - ( 8 )
phaseerrrms = 2 phaseerrpeak 2 3 = 2 3 phaseerrpeak - - - ( 9 )
magerrpeak = evmpeak = 1 + &gamma; 2 - ( 1 - &gamma; ) 1 + &gamma; 2 - - - ( 10 )
Magerrrms and evmrms carries out vector calculus according to definition.
magerri=| ma'-|oai| | i=1,2,3 (11)
magerrrms = &sigma; i = 1 3 magerr i 2 3 m a ' &times; 100 % - - - ( 12 )
evm i = | oa i &rightarrow; - oa &prime; &rightarrow; | , i = 1,2,3 - - - ( 13 )
evmrms = &sigma; i = 1 3 evm i 2 3 m a &prime; &times; 100 % - - - ( 14 )
Formula (8) to (14) gives the functional relationship of modulation error parameter and " equilateral triangle radius " γ.In some works In Cheng Yingyong, give the digital modulation parameters such as the evmrms needing setting it is desirable to corresponding γ-value, that is, need to be similar to (15) formula Functional relationship:
γ=f (evmrms) (15)
For functional relationship shown in (15) formula, contrafunctional method can be analysed it is also possible to count by interpolation etc. by solving Value method is tried to achieve, and can reach very high precision.
3rd, form complete error symbolic frame
This step will ensure that strict equiprobability each error symbol occurs in the measurements, reduces symbol and randomness to metering The adverse effect of test is it is ensured that the strict unification of mathematical model and measurement mechanism.
Specifically comprise the following steps that
1st, construct frame.List the sequence of complex numbers s of certain one group of standard symbol of digital modulation mode1、s2、s3…sk…sk, right In each symbol sk, according to step b, design corresponding l error symbol sk1、sk2…skl.Then formed complete as shown in fig. 13 that Whole error symbol frame, this symbolic frame is named as framee.As shown in figure 13.
2nd, combined frames.According to the requirement of actual test degree symbol sebolic addressing length, by multiple complete error symbol frames Framee combines, and material is thus formed a new sign matrix framem, framem restructuring is become the plural number of single file Rearrangement is upset in the position of each symbol then by sign magnitude sequence at random, forms symbol sebolic addressing sframe.
3rd, construct waveform.Using step a, based on sframeGenerate continuous baseband signal waveform wb of head and the tailframeIt is ensured that wbframeIn corresponding symbol sebolic addressing, each symbol that framee comprises is that strict equiprobability occurs.
4th, generate rf modulated signal.By baseband waveform wbframeCarry out proper treatment and be allowed to adaptive vector modulated signal sending out The operation of raw device requires, and baseband waveform is imported the baseband portion of vector modulation signal generator, and then produces rf modulations letter Number wrframe, VSA (vsa) tested for the input of this radiofrequency signal is demodulated.
5th, setting VSA (vsa).If wbframeCorresponding symbol total number is p, then by tested vector The symbolic analysis length of signal analyzer (vsa) is set to the integral multiple of p, this reduces symbol and randomness to metering The adverse effect of test is it is ensured that the mathematical model of digital modulation error setting and measurement mechanism can strictly be unified.
Based on said method and coordinate respective vectors modulation signal generator hardware of the prior art, just construct " base Digital modulation error parameter metering device in head and the tail continuous wave ", if calculate thering is the preferable goodness of fit with actual measurement, and Digital modulation error measuring value has good repeatability and stability, then this device can serve as Transfer Standards or metering Standard uses.
Build experiment test system according to the method described above, produced and there is the digital modulation signals wr being provided with errorframe, Then input a signal into tested VSA (vsa) to measure, shown in following Figure 14~22 of measurement result.
Solve the measurement problem of the digital modulation error parameters such as evm by the method for the present invention and system, experiment shows this Metering method and system have characteristics that 1, can realize digital modulation error arranges value continuously may be used interior in a big way Adjust, expected digital modulation error parameter can be set out.2nd, the measured value of the digital modulation parameter such as evm and expected calculated value It is identical.3rd, the measured value reading of digital modulation parameter has very good repeatability, is embodied in the standard of reading Difference very little.4th, this digital modulation error parameter metering method and system are applied to almost all of digital modulation mode, such as 256qam modulation and its modulation system of higher order.
Described above is only presently preferred embodiments of the present invention it is impossible to limit the interest field of the present invention with this, according to The equal change that scope of the present invention patent is done, still belongs to the scope that the present invention is covered.

Claims (6)

1. a kind of digital modulation error parameter metering method based on Waveform Design is it is characterised in that the method includes following step Rapid:
A, realize the continuous step of finite length digital modulation signals baseband waveform head and the tail;
B, the basis based on error symbol design, the step realizing the setting of digital modulation error;
C, the step forming complete error symbolic frame;
Wherein, step a refers under the hardware structure of vector modulation signal generator, the baseband wave of one finite duration of design Shape, is then introduced into the baseband portion of signal generator, and described vector modulation signal generator was by should in the way of loop play The baseband waveform of finite duration carries out head and the tail and connects, thus realize finite length digital modulation signals baseband waveform head and the tail connecting Continuous.
2. metering method according to claim 1, wherein, step a comprises the steps:
Step 1: construction symbol sebolic addressing, it is l symbol period that the time domain impulse of baseband filter blocks length, constructs unfiltered symbol The rear l symbol of the front l symbol of number sequence and sequence constitutes symmetrical relationss;
Step 2: pretreatment, according to the baseband waveform sample rate setting and baseband signalling speed, by unfiltered symbol sebolic addressing and base Band filter shock response waveform carries out suitable interpolation and processes, and makes the corresponding time-domain sampling interval of two waveforms equal;
Step 3: convolution, pretreated for step 2 two values sequence is done convolution, participates in the iq symbol sequence of convolution The sampling number of row is n, and corresponding baseband filter shock response waveform sampling points are g, and the points of convolution results are n+ g-1;
Step 4: intercept, the convolution results of step 3 are intercepted, intercept the n numerical point in stage casing, it is required for intercepting result Baseband waveform sequence of values, this sequence has the continuous characteristic of head and the tail.
3. metering method according to claim 1, wherein, step a comprises the steps:
Step 1: pretreatment, according to the baseband waveform sample rate setting and baseband signalling speed, by unfiltered iq symbol sebolic addressing Carry out suitable interpolation with baseband filter shock response waveform to process and null filling, make the corresponding time-domain sampling of two waveforms Interval is equal, and counts and equal be n;
Step 2: fast Fourier transform, the fast Fourier that pretreated for step 1 two values sequence is done n point respectively becomes Change, form the sequence that two new n ' put: xk1 and xk2, then xk1 with xk2 is executed with the fortune that corresponding element order is multiplied one by one Calculate, its result yk is also the sequence of values of a n point;
Step 3: an inverse fast fourier, an inverse fast fourier computing is executed to yk, forms the waveform number of a time domain Value sequence s, this sequence is required baseband waveform sequence of values, and this sequence has the continuous characteristic of head and the tail.
4. metering method according to claim 1, wherein, step a comprises the steps:
Step 1: sequence repeats, initial iq symbol sebolic addressing comprises m symbol, repeats this sequence, forms new sequence;
Step 2: pretreatment, according to the baseband waveform sample rate setting and baseband signalling speed, by new sequence and baseband filtering Device shock response waveform carries out suitable interpolation and processes, or adjustment sample rate, makes the corresponding time-domain sampling interval of two waveforms Equal;
Step 3: convolution, pretreated for step 2 two values sequence is done convolution computing, participates in the iq symbol of convolution The sampling number of number sequence is 2n, and corresponding baseband filter shock response waveform sampling points are g, and convolution results sequence is used Ck1 represents, its points are 2n+g-1;
Step 4: intercept, convolution results ck1 of step 3 are intercepted, from selected starting point, continuous n number after intercepting The sequence that value point is constituted, this intercepting result ck2 is required baseband waveform sequence of values, and this sequence has the continuous characteristic of head and the tail, Above-mentioned selected starting point in the selection range in ck1 is: from the beginning of the g+1 point, to the n-th -2 points, all optional.
5. metering method according to claim 1, wherein, step c comprises the steps:
Step 1: construction frame, list the sequence of complex numbers s of one group of standard symbol1、s2、s3…sk, for each symbol sk, according to step B, designs corresponding l error symbol sk1、sk2…skl, form complete error symbolic frame, this complete error symbolic frame named For framee;
Step 2: combined frames, according to the requirement of actual test degree symbol sebolic addressing length, by multiple complete error symbolic frame framee Combine, form a new sign matrix framem, sign matrix framem restructuring is become the complex symbol number of single file Rearrangement is upset in the position of each symbol then by value sequence at random, forms symbol sebolic addressing sframe
Step 3: construction waveform, according to step a, based on symbol sebolic addressing sframeGenerate the continuous baseband waveform wb of head and the tailframeIt is ensured that Baseband waveform wbframeIn corresponding symbol sebolic addressing, each symbol that complete error symbolic frame framee comprises is strict equiprobability Occur;
Step 4: generate rf modulated signal, by baseband waveform wbframeCarry out proper treatment and be allowed to adaptive vector modulated signal sending out The operation of raw device requires, and baseband waveform is imported the baseband portion of vector modulation signal generator, and then produces rf modulations letter Number, VSA vsa tested for the input of this rf modulated signal is demodulated;
Step 5: setting VSA, baseband waveform wb is setframeCorresponding symbol total number is p, by tested arrow The symbolic analysis length of amount signal analyzer vsa is set to the integral multiple of p.
6. a kind of metering system of the metering method realized described in one of claim 1-5, comprising:
For realizing the device of described step a;
For realizing the device of described step b;
For realizing the device of described step c.
CN201310499303.8A 2013-06-20 2013-10-22 Digital modulation error parameter metering method and system based on waveform design Active CN103532906B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310499303.8A CN103532906B (en) 2013-06-20 2013-10-22 Digital modulation error parameter metering method and system based on waveform design

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN201310247509.1 2013-06-20
CN201310247509 2013-06-20
CN2013102475091 2013-06-20
CN201310499303.8A CN103532906B (en) 2013-06-20 2013-10-22 Digital modulation error parameter metering method and system based on waveform design

Publications (2)

Publication Number Publication Date
CN103532906A CN103532906A (en) 2014-01-22
CN103532906B true CN103532906B (en) 2017-02-01

Family

ID=49934591

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310499303.8A Active CN103532906B (en) 2013-06-20 2013-10-22 Digital modulation error parameter metering method and system based on waveform design

Country Status (1)

Country Link
CN (1) CN103532906B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105302935B (en) * 2015-08-10 2018-03-27 工业和信息化部电信研究院 Digital demodulation and measurement analysis method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101499861A (en) * 2008-02-02 2009-08-05 大唐移动通信设备有限公司 Measuring method and apparatus for error vector amplitude
CN101641931A (en) * 2006-12-18 2010-02-03 艾利森电话股份有限公司 Pulse width modulator
CN102916753A (en) * 2012-08-08 2013-02-06 工业和信息化部通信计量中心 Measurement and calculation methods of digital modulation error

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7519383B2 (en) * 2005-07-05 2009-04-14 Litepoint Corp. Method for efficient calibration of EVM using compression characteristics
CN101694503B (en) * 2009-10-10 2011-10-26 工业和信息化部通信计量中心 Method and system for metering digital demodulating error parameter based on amplitude modulation method or phase modulation method
CN101702697B (en) * 2009-11-04 2012-05-23 工业和信息化部通信计量中心 Method and system for metering digital demodulation error parameters based on continuous wave combination method
CN103051585A (en) * 2012-07-27 2013-04-17 工业和信息化部通信计量中心 Metering method and device for digital modulation parameter based on IQ (Intelligence Quotient) plane elliptical orbit

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101641931A (en) * 2006-12-18 2010-02-03 艾利森电话股份有限公司 Pulse width modulator
CN101499861A (en) * 2008-02-02 2009-08-05 大唐移动通信设备有限公司 Measuring method and apparatus for error vector amplitude
CN102916753A (en) * 2012-08-08 2013-02-06 工业和信息化部通信计量中心 Measurement and calculation methods of digital modulation error

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
数字调制参数校准仿真研究;王希东,孔路平;《电子测量与仪器学报》;20090831;第23卷(第8期);全文 *

Also Published As

Publication number Publication date
CN103532906A (en) 2014-01-22

Similar Documents

Publication Publication Date Title
CN102882818B (en) Amending method and amending system directing at unbalanced zero intermediate frequency feedback
CN105978643B (en) The measuring method of quadrature phase shift keyed modulators delay testing
CN103647554A (en) Method for measuring bias and gain errors of time alternative ADC acquisition system
CN106018955B (en) The low rate end frequency estimating methods of fast convolution tunable filter group
CN103454494A (en) High-precision harmonic analysis method
CN105302935A (en) Digital demodulating and measurement analysis method
CN102387098A (en) Frequency deviation and phase deviation combination measurement method and apparatus thereof
CN105510706A (en) High-precision undersampling frequency measurement method
CN109274623A (en) A kind of implementation method based on ultra long FFT amendment carrier wave frequency deviation
CN103532906B (en) Digital modulation error parameter metering method and system based on waveform design
CN104185271A (en) Identification and positioning method for multiple passive intermodulation generation points
CN107454032A (en) OFDM frequency offset estimation method based on amplitude product between subcarriers
CN105491585B (en) It is a kind of based on non-signaling TD_LTE_Advanced base station signal analytical equipment and method
CN106031046A (en) Compressive sensing-based signal processing method and device
CN103873416A (en) EVM (Error Vector Magnitude) phase estimating and compensating method
CN111999587B (en) Iterative vector error correction method for modulation domain network parameter test
CN105471525A (en) Signal processing method for four-channel compressed sensing digital receiver of vector network analyzer
CN107947899B (en) Multi-user signal-to-noise ratio blind estimation method in single-carrier interleaved frequency division multiple access system
CN109740225A (en) A method of it calculates and assesses for broad-band TWT
CN103616566B (en) A kind of linear average cymoscope implementation method based on software radio
CN103051585A (en) Metering method and device for digital modulation parameter based on IQ (Intelligence Quotient) plane elliptical orbit
CN102944885B (en) Method and device for simulating multi-channel consistency by satellite navigation signal
CN109030942B (en) Harmonic phase angle analysis method
CN103344930B (en) A kind of metering method of FSK digital demodulating error parameter and system
CN101672877B (en) Amplitude-modulated signal accompanying phase modulation measurement method based on vector analysis

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211229

Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee after: CHINA ACADEMY OF INFORMATION AND COMMUNICATIONS

Address before: 100083 No. 52 Garden North Road, Beijing, Haidian District

Patentee before: TELECOMMUNICATION METROLOGY CENTER OF MINISTRY OF INDUSTRY AND INFORMATION TECHNOLOGY