CN103543426B - A kind of Network Analyzer subrane calibration interpolation compensation method - Google Patents

A kind of Network Analyzer subrane calibration interpolation compensation method Download PDF

Info

Publication number
CN103543426B
CN103543426B CN201310512957.XA CN201310512957A CN103543426B CN 103543426 B CN103543426 B CN 103543426B CN 201310512957 A CN201310512957 A CN 201310512957A CN 103543426 B CN103543426 B CN 103543426B
Authority
CN
China
Prior art keywords
error coefficient
calibration
interpolation
error
data
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
CN201310512957.XA
Other languages
Chinese (zh)
Other versions
CN103543426A (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.)
CLP Kesiyi Technology Co Ltd
Original Assignee
CETC 41 Institute
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 CETC 41 Institute filed Critical CETC 41 Institute
Priority to CN201310512957.XA priority Critical patent/CN103543426B/en
Publication of CN103543426A publication Critical patent/CN103543426A/en
Application granted granted Critical
Publication of CN103543426B publication Critical patent/CN103543426B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Resistance Or Impedance (AREA)
  • Complex Calculations (AREA)

Abstract

The present invention provides a kind of Network Analyzer subrane calibration interpolation compensation method, comprises the following steps: step one: preserve the current setting state of instrumental effects calibration effect;Step 2: the measurement carrying out standard component performs calibration;Step 3: be calculated error coefficient;Step 4: judge whether measuring state changes, if yes then enter calibration interpolative data handling process, generates new error coefficient matrix;If being otherwise directly entered step 5;Step 5: measure measured piece, connects measured piece, performs measurement;Step 6: error correction, extracts error coefficient from error coefficient matrix, and according to the suitable error model of calibrating mode different choice, by measurement data and error coefficient, solving equation obtains the actual value of measured piece;Step 7: judge whether to be measured, terminates program in this way;As otherwise returned step 4.Use such scheme, utilize suitable interpolation algorithm to calculate the error coefficient of non-calibrated point, thus complete error correction.

Description

A kind of Network Analyzer subrane calibration interpolation compensation method
Technical field
The invention belongs to technical field of measurement and test, in particular a kind of Network Analyzer partial wave Section calibration interpolation compensation method.
Background technology
The most complicated relative to Other Instruments, the most complete place of vector network analyzer is that Having error correction function, the method for its application software makes up the deficiency of hardware design, logical Cross and set up suitable error model, use calibrating device known to characteristic to complete to calibrate, thus pole The earth improves the certainty of measurement of system.
Calibration is a process that requirement is the strictest, even connects the dynamics difference of calibrating device all Final test precision can be affected, be in general not allow to change instrument state after calibration, Must recalibrate after change state.But for the particular/special requirement of certain user, such as, think Measure the measured piece of different frequency range simultaneously, measurement result is required it is not the highest simultaneously, but Being the highest to measuring efficiency requirements, be not suitable for repeatedly calibrating, current solution is to utilize The error coefficient of phase near point completes error correction or calibration data is carried out simple interpolation. Close point is utilized to carry out error correction probably because phase near point causes interpolation error across wave band; And only consider that calibration data does not considers to cause the simple interpolation of every factor of interpolation error, Also there will be interpolation result the most incorrect undesirable situation.
Therefore, prior art existing defects, need to improve.
Summary of the invention
The technical problem to be solved is for the deficiencies in the prior art, it is provided that a kind of The calibration interpolation compensation method of Network Analyzer subrane.
Technical scheme is as follows:
A kind of Network Analyzer subrane calibration interpolation compensation method, wherein, comprises the following steps:
Step one: preserve the current setting state of instrumental effects calibration effect, described currently set Configuration state is frequency and power;
Step 2: the measurement carrying out standard component performs calibration;
Step 3: be calculated error coefficient, is saved in error coefficient matrix;
Step 4: judge whether measuring state changes, if yes then enter calibration interpolative data Handling process;If otherwise entering step 5;
Step 5: measure measured piece, connects measured piece, performs measurement;
Step 6: error correction, extracts error coefficient, according to school from error coefficient matrix The quasi-suitable error model of mode different choice, by measurement data and error coefficient, solves Equation obtains the actual value of measured piece;
Step 7: judge whether to be measured, terminates program in this way;As otherwise returned step Four.
Described Network Analyzer subrane calibration interpolation compensation method, wherein, described step 4 In, concretely comprising the following steps of calibration interpolative data handling process:
Step 401: extract initial error coefficient, by the true error under actual alignment state Coefficient extracts;
Step 402: plural number error coefficient is converted to amplitude data and phase data;
Step 403: phase data is converted to inverse folding linear phase data;
Step 404: subrane performs amplitude data, the interpolation of phase data, interpolation algorithm Choose Lagrangian Arithmetic, perform interpolation according to currently practical measuring state, obtain all surveys The error coefficient of amount point;
Step 405: amplitude, phase data are re-converted into complex data, generates new Error coefficient, returns error coefficient matrix.
Described Network Analyzer subrane calibration interpolation compensation method, wherein, described step In 401, initial error coefficient is plural number.
Described Network Analyzer subrane calibration interpolation compensation method, wherein, described step In 404, when subrane performs the interpolation of amplitude data, phase data, do not do in wave band Insert.
The implementation using above-mentioned calibration interpolation is relatively suitable for the measurement of vector network analyzer Feature, so interpolation is more satisfactory, utilizes suitable interpolation algorithm to calculate non-calibrated point Error coefficient, thus complete error correction.Present invention operation is very simple, can be whole Individual algorithm is embedded in vector net analysis instrument host software, it is also possible to utilize outside instrument Program control mode realizes, and practicality, autgmentability are the best.
Accompanying drawing explanation
Fig. 1 is the flow chart of interpolation compensation method of the present invention.
Fig. 2 is that the present invention calibrates interpolative data process chart.
Fig. 3 is different modes interpolation results schematic diagram.
Detailed description of the invention
Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.
Embodiment 1
Fig. 1 is that vector network analyzer of the present invention measures process flow diagram flow chart.
The whole measurement process of vector network analyzer substantially can be divided into following step:
It is first begin to, starts program and run;
Step one: preserving instrument current setting state, main preservation affects calibration effect here Instrument current setting state such as such as frequency, power etc.;
Step 2: perform calibration, the calibration of vector network analyzer mainly include mechanically calibrated, Two kinds of calibration modes of electric calibration, can also be subdivided into the multiple calibration side such as SOLT, TRL simultaneously Formula, main process is just by the measurement of standard component;
Step 3: utilize the program curing of vector network analyzer to complete solving of equation group, Finally give error coefficient, be saved in error coefficient matrix;
Step 4: judge whether align mode changes, system to be monitored before formal measurement State changes, if the instrument state of the current state of vector network analyzer and preservation differs Cause it is necessary to enter calibration interpolation flow process.Primary operational carries exactly from error coefficient matrix Take initial error coefficient, use suitable interpolation algorithm, generate new error coefficient, passback In error coefficient matrix;
Step 5: measure measured piece, connects measured piece, performs measurement;
Step 6: error correction, extracts error coefficient, according to school from error coefficient matrix The quasi-suitable error model of mode different choice, by measurement data and error coefficient, solves Equation obtains the true response of measured piece;
Step 7: judge whether to be measured, is measured end program;Measurement is not fully complete Return to step 4 continue to measure.
Present invention is generally directed to step 4 alignment interpolation and perform method, first explain lower calibration The concept of interpolation, vector network analyzer belongs to sweep measurement instrument, it is assumed that we are initial The calibration performed under the state that frequency 1GHz, termination frequency 2GHz, number of scan points are at 101, If number of scan points to be changed into 201 points, then be to newly increase by the frequency having 100 points Frequency, do not have correspondence error coefficient, at this moment will utilize suitable interpolation algorithm, from Existing 101 point data expand the error coefficient matrix of 201.
Fig. 2 is that the present invention performs vector network analyzer calibration interpolative data process chart.
The Railway Project that need to process when specifically introducing the implementation method of interpolation in the present invention.
The method of curve matching is a lot, through lot of experiments, analysis, we determined that glug is bright Day interpolation algorithm is best suitable for the calibration interpolation of vector network analyzer, and the calibration result obtained is Ideal, therefore we use Lagrange's interpolation algorithm as the processing method of calibration interpolation. Being embodied in a lot of data of Lagrange's interpolation algorithm is all described, the most superfluous at this State.
After have selected vector network analyzer calibration interpolation algorithm, also have and what time should be noted that:
(1) measurement data of vector network analyzer is plural number, and corresponding error coefficient is also Plural number, comprises real part, imaginary part, and current various interpolation algorithms is difficult to insert plural number Value, so needing to be separated into real number to carry out interpolation arithmetic.Here amplitude, phase place are used Interpolation method, needs that complex data is converted into amplitude, phase place carries out interpolation.
(2) interpolation for phase place has any to pay special attention to, usual vector network analyzer Phase Processing general provision-180 degree to 180 degree intervals, if processed in such a way, It is easy to the trip point that phase place occurs, causes big error to interpolation, so phase place is being inserted Before value, phase place to be changed into non-collapsible linear phase.For example, if two known frequencies The phase place of rate point f1 and f2 is-90 degree and+90 degree respectively, then central point (f1+f2)/2 Phase place will cause unusual, be 0 degree or 180 degree on earth, if employing non-folding Folded linear phase, the change direction of phase place is fixing, would not produce the puzzlement of interpolation. So-called linear phase, it is simply that phase place was not limited in the cycle, does not specify at-180 degree To+180 degree, the most do not specify at 0 degree to 360 degree, but according to the variation tendency line of phase place Property cumulative, namely arbitrary neighborhood point-to-point transmission phase difference is not more than 180 degree.
(3) due to the special technical requirement of vector network analyzer own, it is difficulty with at full range Continuous scanning in section, is typically necessary tape swapping and processes, and particularly frequency span is bigger Vector network analyzer, there may be tens tape swapping points especially.At tape swapping point, due to Wave band switches, and causes whole hardware corridor to change so that the characteristic of system creates very Big saltus step, this also causes the error coefficient of two wave bands at tape swapping point notable difference occur. Therefore, calibration interpolation is not suitable for operating continuously in whole frequency range, otherwise can be attached at tape swapping point Near there is bigger error, carries out so our interpolation calculates in being limited in wave band, so may be used To obtain reasonable interpolation.
Integrate, as in figure 2 it is shown, the present invention calibration interpolation perform process can be divided into Under several steps:
Step 401: extract initial error coefficient, by the true error under actual alignment state Coefficient extracts;
Step 402: plural number error coefficient is converted to amplitude, phase data;
Step 403: phase data carries out inverse folding linearization process, forms linear phase Data;
Step 404 subrane performs amplitude, phase interpolation, and interpolation algorithm chooses Lagrange Algorithm.Perform interpolation according to currently practical measuring state, obtain all error systems measuring point Number;
Step 405: amplitude, phase data are re-converted into complex data, generates new Error coefficient, returns error coefficient matrix, completes follow-up error correction.
Fig. 3 represents the concrete effect for plural number interpolation in different ways.
For plural number interpolation, be for real part, the direct interpolation of imaginary part or be converted into amplitude, It is more suitable that phase place carries out interpolation?It is concrete that Fig. 3 clearly gives that two kinds of interpolation methods obtain Result, it is assumed that A, B 2 are the true error coefficients that actual alignment obtains, and A point frequency Rate be f1, B dot frequency be f2, we respectively symbolically be Af1 and Bf2, now Require the error coefficient at f=(f1+f2)/2, if using amplitude, phase-interpolation, I The result that obtains be Cf;If using real part, imaginary part interpolation, the result that we obtain is Df.Therefrom we can be clearly seen that, it is more preferable that use amplitude, phase place carry out interpolation, The data obtained are more accurate, and therefore, we use amplitude, phase-interpolation method.
Embodiment 2
On the basis of above-described embodiment, further, a kind of Network Analyzer partial wave of the present invention Section calibration interpolation compensation method, wherein, comprises the following steps:
Step one: preserve the current setting state of instrumental effects calibration effect, described currently set Configuration state is frequency and power;
Step 2: the measurement carrying out standard component performs calibration;
Step 3: be calculated error coefficient, is saved in error coefficient matrix;
Step 4: judge whether measuring state changes, if yes then enter calibration interpolative data Handling process;If otherwise entering step 5;
Step 5: measure measured piece, connects measured piece, performs measurement;
Step 6: error correction, extracts error coefficient, according to school from error coefficient matrix The quasi-suitable error model of mode different choice, by measurement data and error coefficient, solves Equation obtains the actual value of measured piece;
Step 7: judge whether to be measured, terminates program in this way;As otherwise returned step Four.
Furthermore, in described step 4, the concrete steps of calibration interpolative data handling process For:
Step 401: extract initial error coefficient, by the true error under actual alignment state Coefficient extracts;
Step 402: plural number error coefficient is converted to amplitude data and phase data;
Step 403: phase data is converted to inverse folding linear phase data;
Step 404: subrane performs amplitude data, the interpolation of phase data, interpolation algorithm Choose Lagrangian Arithmetic, perform interpolation according to currently practical measuring state, obtain all surveys The error coefficient of amount point;
Step 405: amplitude, phase data are re-converted into complex data, generates new Error coefficient, returns error coefficient matrix.
Furthermore, in described step 401, initial error coefficient is plural number.
Furthermore, in described step 404, subrane performs amplitude data, phase data Interpolation time, do not do across wave band interpolation.
The implementation using above-mentioned calibration interpolation is relatively suitable for the measurement of vector network analyzer Feature, so interpolation is more satisfactory, utilizes suitable interpolation algorithm to calculate non-calibrated point Error coefficient, thus complete error correction.Present invention operation is very simple, can be whole Individual algorithm is embedded in vector net analysis instrument host software, it is also possible to utilize outside instrument Program control mode realizes, and practicality, autgmentability are the best.
It should be appreciated that for those of ordinary skills, can according on state Bright improved or converted, and all these modifications and variations all should be belonged to the appended power of the present invention The protection domain that profit requires.

Claims (1)

1. a Network Analyzer subrane calibration interpolation compensation method, it is characterised in that include Following steps:
Step one: preserve the current setting state of instrumental effects calibration effect, described current setting State is frequency and power;
Step 2: the measurement carrying out standard component performs calibration;
Step 3: be calculated error coefficient, is saved in error coefficient matrix;
Step 4: judge whether measuring state changes, at calibration interpolative data Reason flow process;If otherwise entering step 5;
Step 5: measure measured piece, connects measured piece, performs measurement;
Step 6: error correction, extracts error coefficient, according to calibration from error coefficient matrix The suitable error model of mode different choice, by measurement data and error coefficient, solving equation Obtain the actual value of measured piece;
Step 7: judge whether to be measured, terminates program in this way;As otherwise returned step 4;
In described step 4, concretely comprising the following steps of calibration interpolative data handling process:
Step 401: extract initial error coefficient, by the true error system under actual alignment state Number extracts;
Step 402: plural number error coefficient is converted to amplitude data and phase data;
Step 403: phase data is converted to inverse folding linear phase data;
Step 404: subrane performs amplitude data, the interpolation of phase data, interpolation algorithm selects Take Lagrangian Arithmetic, perform interpolation according to currently practical measuring state, obtain all measurements point Error coefficient;
Step 405: amplitude, phase data are re-converted into complex data, generates new mistake Difference coefficient, returns error coefficient matrix;
In described step 401, initial error coefficient is plural number, in described step 404, partial wave Duan Zhihang amplitude data, phase data interpolation time, do not do across wave band interpolation.
CN201310512957.XA 2013-10-28 2013-10-28 A kind of Network Analyzer subrane calibration interpolation compensation method Active CN103543426B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310512957.XA CN103543426B (en) 2013-10-28 2013-10-28 A kind of Network Analyzer subrane calibration interpolation compensation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310512957.XA CN103543426B (en) 2013-10-28 2013-10-28 A kind of Network Analyzer subrane calibration interpolation compensation method

Publications (2)

Publication Number Publication Date
CN103543426A CN103543426A (en) 2014-01-29
CN103543426B true CN103543426B (en) 2016-06-29

Family

ID=49967050

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310512957.XA Active CN103543426B (en) 2013-10-28 2013-10-28 A kind of Network Analyzer subrane calibration interpolation compensation method

Country Status (1)

Country Link
CN (1) CN103543426B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772174A (en) * 2016-11-15 2017-05-31 中国电子科技集团公司第四十研究所 A kind of Electronic Calibration part automatic marking method
CN107860973B (en) * 2017-11-10 2020-05-12 中国电子科技集团公司第四十一研究所 Frequency response compensation method and system applied to spectrum analyzer
CN107976646B (en) * 2017-11-20 2020-02-07 中国电子科技集团公司第四十一研究所 Signal power characteristic compensation method and device based on vector network analyzer
CN108871406B (en) * 2018-04-28 2021-05-25 常州同惠电子股份有限公司 Algorithm for searching excellent calibration points
CN109782208A (en) * 2019-03-04 2019-05-21 上海精密计量测试研究所 Vector network analyzer auto-calibration device and calibration method
CN109870622A (en) * 2019-03-14 2019-06-11 中国电子科技集团公司第二十九研究所 The measurement method of efficient process array channel radio frequency delay inequality
CN112003651A (en) * 2020-08-18 2020-11-27 深圳宇宙桥无线通信技术有限公司 Network division parameter compensation method, system, storage medium and network division compensation management platform
CN111988063B (en) * 2020-08-31 2021-10-08 中电科思仪科技股份有限公司 Cable loss compensation method and system for comprehensive tester
CN112152732A (en) * 2020-09-27 2020-12-29 深圳宇宙桥无线通信技术有限公司 Curve adjusting method, curve adjusting system, curve adjusting medium, debugging preparation method, debugging preparation system and debugging preparation medium
CN112231985B (en) * 2020-11-04 2023-03-24 中国电子科技集团公司第二十九研究所 Radio frequency filter modeling method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1442704A (en) * 2002-03-05 2003-09-17 特克特朗尼克公司 Improved calibration of vector network analyzer
CN1719791A (en) * 2004-07-07 2006-01-11 华为技术有限公司 Measuring method and apparatus based on network analyzer
EP2101182A2 (en) * 2008-03-13 2009-09-16 Spinner GmbH Vector network analyser calibration system
CN103364751A (en) * 2013-07-11 2013-10-23 中国电子科技集团公司第四十一研究所 Electronic calibration part of vector network analyzer and calibration method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7124049B2 (en) * 2005-01-03 2006-10-17 Agilent Technologies, Inc. Method for implementing TRL calibration in VNA

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1442704A (en) * 2002-03-05 2003-09-17 特克特朗尼克公司 Improved calibration of vector network analyzer
CN1719791A (en) * 2004-07-07 2006-01-11 华为技术有限公司 Measuring method and apparatus based on network analyzer
EP2101182A2 (en) * 2008-03-13 2009-09-16 Spinner GmbH Vector network analyser calibration system
CN103364751A (en) * 2013-07-11 2013-10-23 中国电子科技集团公司第四十一研究所 Electronic calibration part of vector network analyzer and calibration method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
矢量网络分析仪中非***器件的校准方法;郭永瑞等;《仪器仪表学报》;20100831;第31卷;文章第1节 *

Also Published As

Publication number Publication date
CN103543426A (en) 2014-01-29

Similar Documents

Publication Publication Date Title
CN103543426B (en) A kind of Network Analyzer subrane calibration interpolation compensation method
Crotti et al. Frequency response of MV voltage transformer under actual waveforms
CN103364751B (en) A kind of vector network analyzer Electronic Calibration part and calibration steps
CN107505499A (en) For the current transformer and the method and apparatus of the phase error of power-measuring device for learning to include current transformer error correction
CN106918795B (en) Precision resister calibration system based on FPGA and the resistance calibrating method using system realization
CN106053899B (en) A kind of reference signal generating system and method
CA3074115C (en) Measurement system and calibration method with wideband modulation
Kaczmarek et al. Measuring system for testing the transformation accuracy of harmonics of distorted voltage by medium voltage instrument transformers
CN108809447A (en) A kind of calibration of amplitude and phase method of phased array antenna
CN103595404A (en) Frequency mixer group delay measuring circuit and method
Crotti et al. A new industry-oriented technique for the wideband characterization of voltage transformers
CN105717389A (en) A Mobile Device Tester For Precise Inductive Power Measurement And A Calibration Unit Therefor
CN107860973B (en) Frequency response compensation method and system applied to spectrum analyzer
CN109283391A (en) A kind of synchronized phasor method for measurement based on nonlinear fitting
CN103454489B (en) The loss power scaling method of matching network and system
Jebroni et al. Gain correction and phase compensation of a smart electrical energy meter
CN109088675A (en) A kind of the access calibration method and device of radio-frequency signal source
Gu et al. NIST-traceable microwave power measurement in a waveguide calorimeter with correlated uncertainties
US10422846B2 (en) Method for calibrating a radio frequency test instrument and radio frequency test instrument
US8169211B2 (en) Method for verifying the bandwidth and phase of a digital power control system
CN109100672A (en) A kind of adjustment system and adjusting process of the voltage and current sampling module suitable for non-intrusion type load monitoring
Crotti et al. Frequency calibration of voltage transformers by digital capacitance bridge
CN101227240B (en) Characterization of a frequency response for a frequency translation device
US10830859B2 (en) Validating fundamental-only energy measurement
Gallo et al. Evaluation of metrological performance of electromagnetic current measurement transformers in non-sinusofdal conditions

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for 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: 20190225

Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Patentee after: China Electronics Technology Instrument and Meter Co., Ltd.

Address before: 266555 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong

Patentee before: The 41st Institute of CETC

CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Patentee after: CLP kesiyi Technology Co.,Ltd.

Address before: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Patentee before: CHINA ELECTRONIC TECHNOLOGY INSTRUMENTS Co.,Ltd.