CN103616566B - A kind of linear average cymoscope implementation method based on software radio - Google Patents
A kind of linear average cymoscope implementation method based on software radio Download PDFInfo
- Publication number
- CN103616566B CN103616566B CN201310587171.4A CN201310587171A CN103616566B CN 103616566 B CN103616566 B CN 103616566B CN 201310587171 A CN201310587171 A CN 201310587171A CN 103616566 B CN103616566 B CN 103616566B
- Authority
- CN
- China
- Prior art keywords
- signal
- cymoscope
- linear average
- digital
- component signal
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 238000005070 sampling Methods 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000012546 transfer Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 6
- 238000012360 testing method Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 7
- 101100388212 Arabidopsis thaliana DSP3 gene Proteins 0.000 description 4
- 230000005611 electricity Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Landscapes
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The present invention provides a kind of linear average cymoscope implementation method based on software radio, comprises the following steps: step one: analog intermediate frequency is carried out signal sampling；Step 2: the signal after sampling is divided into in-phase component signal after Digital IF Processing and orthogonal component signal sends；Step 3: recover intermediate frequency envelope signal from in-phase component signal and orthogonal component signal；Step 4: intermediate frequency envelope signal is simulated artificial network computing；Step 5: the output of analog simulation network is obtained maximum detection value, as the output voltage of linear average cymoscope.Compare the analog circuit implementation method of routine, use such scheme, stability and concordance that linear average cymoscope is measured can be improved, reduce measuring instrument volume and weight, obtain more excellent measuring speed, particularly can realize the reconstruct of different time constant linear average cymoscope the most easily, reduce design and debugging cost.
Description
Technical field
The invention belongs to linear average reading detector technical field in electromagnetic compatibility test, in particular a kind of based on
The linear average cymoscope implementation method of software radio.
Background technology
The quickening come into the market along with various electronic products, valued degree is also in product test for electromagnetic compatibility test
More and more higher, not only at authentication test laboratory, research and develop in each big Electronic Products Corporation, test and the department such as quality inspection, electromagnetism
Compatibility test is also becoming an important link.
When carrying out electromagnetic interference emission (EMI) and measuring, it is required for disturbance Source Type and circulation way, uses not
Same test layouts, method of testing and subtest equipment.But it is identical for having a bit, this class testing is required for using highly sensitive
The electromagnetic interference test receiver (also known as EMI receiver) of degree.EMI receiver is that voltage amplitude is changed by one with frequency
The instrument that shows of rule, its display receives the frequency spectrum of signal, thus will appreciate that amplitude and the occurring source of electromagnetic interference.
In decades, along with the development of every field technology, EMI receiver towards wide-band, wide-range, high accuracy,
At high speed, digitized, miniaturization, intellectuality, modularity and systematized direction are developed.The most in recent years, microelectronics
The new techniques such as technology, Digital Signal Processing and software and radio technique are gradually applied to setting of novel EMI receiver
In meter exploitation, receiver performance is made to produce qualitative leap.Wherein, the central idea of software and radio technique is: construct a tool
There are opening, standardization, modular general hardware platform, various radio analysis functions, such as modulation /demodulation, detection etc. are used
Software completes, and make broadband A/D and D/A conversion as far as possible near antenna, with develop high degree of flexibility, open new one
The wireless communication system in generation.
In electromagnetic compatibility test field, CISPR (CISPR) from 1933 set up since, one
One of straight important maker as international EMC testing standard.Company follows the standard of CISPR all the time in the industry, as
The global specification of civilian EMC test.In recent years, some corrections that CISPR standard is made about cymoscope type,
In amendment A1:2002 of CISPR16-1:1999 (second edition) for test uses average reading detector be made that and repair
Change, use new linear average (CISPR Average) cymoscope to instead of original average reading detector；And define
The time constant of cymoscope.The main purpose of this amendment is able to measure the instability that intermittence occurs more accurately
Harassing and wrecking signal.According to new CISPR standard, EMI receiver is needed to realize the test of linear average cymoscope in testing.
New CISPR standard describes the conventional implementation method of linear average cymoscope as it is shown in figure 1, in envelope detection
Employ a pointer instrument artificial network 102 after device 101 and between A/D converter 103, be used for simulating one traditional
The measuring state of pointer instrument.By the method it can be seen that it mainly uses analog circuit (such as one piece of special line of design
Property average reading detector analog circuit board) realize envelope detector 101 and instrument simulation network 102, then by using
A/D converter 103 and microprocessor 104 carry out the output of continuous monitor table, it is thus achieved that peakreading is as linear average detection
The output reading of device.
Find, owing to analog circuit exists temperature not during the conventional implementation method of linear average reading detector uses
The factors such as stability, often there is bigger deviation with actual value in measurement result.And, the new standard of CISPR at least needs 2 kinds
The linear average cymoscope of different time constant (affecting the important parameter of geophone property), this just increases the design of circuit
Difficulty, too increases associated production and debugging cost simultaneously.If the board of this analog circuit being installed in EMI receiver, also
The volume and weight of equipment can be increased, be unfavorable for that the carrying of instrument and user use.It addition, the linear averaging that analog circuit realizes
The process time of value cymoscope is relatively long, makes troubles to user.
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 based on software radio
Linear average cymoscope implementation method.
Technical scheme is as follows:
A kind of linear average cymoscope implementation method based on software radio, wherein, comprises the following steps:
Step one: analog intermediate frequency is carried out signal sampling；
Step 2: the signal after sampling is divided into in-phase component signal and orthogonal component signal after Digital IF Processing
Send；
Step 3: recover intermediate frequency envelope signal from in-phase component signal and orthogonal component signal；
Step 4: intermediate frequency envelope signal is simulated artificial network computing；
Step 5: the output of analog simulation network is obtained maximum detection value, as the output electricity of linear average cymoscope
Pressure.
Described linear average cymoscope implementation method, wherein, described step 2: Digital IF Processing, including numeral
Down coversion, low pass filtered involve extraction, and first, the signal that IF input signals produces with digital local oscillator together enters digital mixer
Complete Digital Down Convert；Then, after Digital Down Convert, obtain zero intermediate frequency digital signal, including in-phase component signal and orthogonal point
Amount signal two paths of signals altogether, carries out digital low-pass filtering respectively and extraction processes, to reduce the speed of signal and to ensure signal letter
Breath is complete.
Described linear average cymoscope implementation method, wherein, in described step 3, described in recover intermediate frequency envelope letter
Number concretely comprise the following steps: by the most squared to in-phase component signal and orthogonal component signal and open radical sign computing again and obtain.
Described linear average cymoscope implementation method, wherein, in described step 4, described analog simulation network operations
Being realized by second order iir digital filter, the transfer function of this wave filter is:
Wherein, S is transmission signal, TMFor time constant, H (S) is transfer function.
Described linear average cymoscope implementation method, wherein, the transfer function etc. of described second order iir digital filter
Valency is the cascade system that is multiplied of two identical first order IIR filtering devices:
Described first order IIR filtering device transfer function is converted into difference equation obtain:
y(n)=b*(x(n)+(x(n-1))-a*y(n-1)
Wherein, x (n) is the input value of wave filter, and y (n) is the output valve of wave filter, b and a is the coefficient of wave filter, its
Value is by time constant TMSampling period T with wave filter input datasDetermining, specific formula for calculation is as follows:
Wherein, specify according to CISPR standard, determine time constant T according to the wave band of required measurementM, simultaneously according to same
The sampling period of branch signal and orthogonal component signal determines TsFinally obtain the value of coefficient a and b of wave filter.
Described linear average cymoscope implementation method, wherein, in described step 5, the described maximum detection value that obtains has
Body step is: be the maximum within the scheduled measurement time in recording step four.
Employing such scheme:
1, temperature stability is high, measures concordance good；With DSP program sequential operation, no matter at a temperature of which kind of, if front
The circuit stability of end, numerical calculation result will not change, and stability and concordance are the best.And use analog circuit, electricity therein
Sensibility reciprocal is relatively big, frequently with line around mode design, relatively big with temperature fluctuation, affect test result.
2, volume is little；Digital detector can calculate with digital intermediate frequency, FFT computing etc. and Digital Signal Processing related operation
Share a digital processing platform, it is not necessary to individually increase circuit.And analog linearity average reading detector needs individually designed electricity
Road.
3, processing speed is fast；Using analog detector, single is effectively measured at least to be needed in 2s, during if less than this
Between, just affect the verity of measurement result.Being measured by reality, the process time of the cymoscope of the present invention can be less than 1s,
Some process time measuring frequency range cymoscope is about 500ms.So compare, use the digital detection of the present invention to measure and significantly save
About measure the time, time-consuming cost.
4, achieve the reconstruct of cymoscope, decrease design and debugging cost；Realize linear flat by the way of DSP program
Average detector, main frame passes through dsp interface, resets different time constants, it is possible to realize different cymoscopes, pass through program
Mode achieve the reconstruct of cymoscope, decrease design and debugging cost.
Accompanying drawing explanation
Fig. 1 be in prior art linear average reading detector realize schematic diagram.
Fig. 2 is present invention linear average based on software radio cymoscope 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
As in figure 2 it is shown, first, A/D converter 1 carries out signal sampling to analog intermediate frequency, and the signal after sampling is fed to
FPGA field programmable logic array 2, with NCO (i.e. digital local oscillator) 201 produce signal together enter digital mixer 202 and
212 complete Digital Down Convert.By Digital Down Convert, analog if signal has been converted to zero intermediate frequency digital signal, includes I/Q altogether
(inphase/orthogonal) two paths of signals, respectively through low-pass filtering 203 and 213, extraction 204 and 214 etc. (low-pass filtering and extraction
Not losing of signal message is ensure that while reducing signal rate) carry out Digital IF Processing, then the base band obtained
Signal is divided into two-way I/Q (inphase/orthogonal component) signal and gives digital signal processor DSP 3.
In fig. 2, in DSP3, linear average cymoscope is realized by 3 steps, first, through digital envelope detection
301, recover the amplitude envelops of the modulated signal on intermediate frequency, i.e. signal；It is simulated fortune subsequently into instrument simulation network 302
Calculate, finally, obtain maximum by peak detection 303, just calculate the output voltage of linear average cymoscope.Last survey
Test result is sent to main frame 4 and carries out display process.
As shown in Figure 2, the instrument simulation network 302 of DSP3 is connected by interface with main frame 4, and main frame 4 is by DSP3 interface
The parameters such as different time constants and sample rate can be set, it is possible to makes DSP3 realize different instrument simulation networks, also
Just can realize different linear average cymoscopes, as long as thus amendment parameter, program need not be revised, be achieved that detection
The reconstruct of device.
The linear average cymoscope that the present invention realizes is completed by program in dsp, and detailed process is as follows:
1) envelope detection: envelope detection is that the inphase/orthogonal component to signal is first the most squared and open radical sign computing again and obtain
's.
2) instrument simulation network: the realization of instrument simulation network is realized by second order iir digital filter, this filter
The transfer function of ripple device is:
Wherein S is transmission signal, TMBeing time constant, time constant is the key parameter of this wave filter, affects wave filter
Filtering performance.In CISPR standard, the time constant of different frequency range is had and following requires to be shown in Table 1:
Time constant requirement to linear average cymoscope in table 1:CISPR standard
For the ease of being realized by program in dsp, this second order iir digital filter is broken down into two identical single orders
The mode that the cascade of iir filter is multiplied:
Described single order I IR filter transfer function is converted into difference equation obtain:
y(n)=b*(x(n)+(x(n-1))-a*y(n-1)
Wherein, x (n) is the input value of wave filter, and y (n) is the output valve of wave filter, b and a is the coefficient of wave filter, by
Time constant TMSampling period T with wave filter input datasDetermining, specific formula for calculation is as follows:
The most in the implementation, time constant T is determined according to table 1 by the wave band of required measurementM, simultaneously according in fig. 2
The sampling period of the I/Q data of FPGA output determines TsFinally obtain the value of coefficient a and b of wave filter.
3) peak detection: the most last process is peak detection, its process was mainly within certain measurement time
Record the maximum of above instrument simulation network, and using this value as the output result of linear average cymoscope.
It should be appreciated that for those of ordinary skills, can be improved according to the above description or be converted,
And all these modifications and variations all should belong to the protection domain of claims of the present invention.
Claims (5)
1. a linear average cymoscope implementation method based on software radio, it is characterised in that comprise the following steps:
Step one: analog intermediate frequency is carried out signal sampling；
Step 2: the signal after sampling is divided into in-phase component signal after Digital IF Processing and orthogonal component signal is sent out
Send；
Step 3: recover intermediate frequency envelope signal from in-phase component signal and orthogonal component signal；
Step 4: intermediate frequency envelope signal is simulated artificial network computing；
Step 5: the output of analog simulation network is obtained maximum detection value, as the output voltage of linear average cymoscope；Institute
Stating step 2: Digital IF Processing, involve extraction including Digital Down Convert, low pass filtered, first, IF input signals is with numeral originally
The signal produced that shakes together enters digital mixer and completes Digital Down Convert；Then, after Digital Down Convert, zero intermediate frequency number is obtained
Word signal, including in-phase component signal and orthogonal component signal two paths of signals altogether, is carried out at digital low-pass filtering and extraction respectively
Reason, to reduce the speed of signal and to ensure that signal message is complete.
2. linear average cymoscope implementation method as claimed in claim 1, it is characterised in that in described step 3, described
Recover concretely comprising the following steps of intermediate frequency envelope signal: by the most squared to in-phase component signal and orthogonal component signal and open radical sign again
Computing obtains.
3. linear average cymoscope implementation method as claimed in claim 1, it is characterised in that in described step 4, described
Analog simulation network operations is realized by second order iir digital filter, and the transfer function of this wave filter is:
Wherein, S is transmission signal, TMFor time constant, H (S) is transfer function.
4. linear average cymoscope implementation method as claimed in claim 3, it is characterised in that described second order IIR numeral is filtered
The transfer function of ripple device is equivalent to the cascade system that is multiplied of two identical first order IIR filtering devices:
Described first order IIR filtering device transfer function is converted into difference equation obtain:
Y (n)=b* (x (n)+(x (n-1))-a*y (n-1)
Wherein, x (n) is the input value of wave filter, and y (n) is the output valve of wave filter, b and a is the coefficient of wave filter, by the time
Constant TMSampling period T with wave filter input datasDetermining, specific formula for calculation is as follows:
Wherein, specify according to CISPR standard, determine time constant T according to the wave band of required measurementM, simultaneously according to in-phase component
The sampling period of signal and orthogonal component signal determines TsFinally obtain the value of coefficient a and b of wave filter.
5. linear average cymoscope implementation method as claimed in claim 4, it is characterised in that in described step 5, described
Obtain maximum detection value to concretely comprise the following steps: be the maximum within the scheduled measurement time in recording step four.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310587171.4A CN103616566B (en) | 2013-11-15 | 2013-11-15 | A kind of linear average cymoscope implementation method based on software radio |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310587171.4A CN103616566B (en) | 2013-11-15 | 2013-11-15 | A kind of linear average cymoscope implementation method based on software radio |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103616566A CN103616566A (en) | 2014-03-05 |
| CN103616566B true CN103616566B (en) | 2016-08-17 |
Family
ID=50167271
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310587171.4A Expired - Fee Related CN103616566B (en) | 2013-11-15 | 2013-11-15 | A kind of linear average cymoscope implementation method based on software radio |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103616566B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105791204A (en) * | 2016-04-29 | 2016-07-20 | 中国人民解放军国防科学技术大学 | Full-digital envelope detection demodulation method of alternating current IRIG-B code and apparatus thereof |
| CN111487476B (en) * | 2020-05-06 | 2022-07-08 | 深圳市鼎阳科技股份有限公司 | Quasi-peak detection method and quasi-peak detector |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1992021984A1 (en) * | 1991-05-27 | 1992-12-10 | Advantest Corporation | Apparatus for measuring average value of impulsive signal |
| JP2010141448A (en) * | 2008-12-10 | 2010-06-24 | Nec Corp | Signal detection device, signal detecting method used for the same, and its program |
-
2013
- 2013-11-15 CN CN201310587171.4A patent/CN103616566B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1992021984A1 (en) * | 1991-05-27 | 1992-12-10 | Advantest Corporation | Apparatus for measuring average value of impulsive signal |
| JP2010141448A (en) * | 2008-12-10 | 2010-06-24 | Nec Corp | Signal detection device, signal detecting method used for the same, and its program |
Non-Patent Citations (2)
| Title |
|---|
| 基于MATLAB的数字下变频器的设计与仿真应用;范文晶等;《电子质量》;20090320(第3期);第37页第2栏至第38页第1栏，图1 * |
| 王铮.CAV检波器及其在GB9254中的使用.《安全与电磁兼容》.2010, * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103616566A (en) | 2014-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106066425B (en) | A kind of impedance measurement device and its method for realizing compensation for calibrating errors | |
| US6316945B1 (en) | Process for harmonic measurement accuracy enhancement | |
| US20100204943A1 (en) | Method and device for the calibration of network analyzers using a comb generator | |
| CN101452067A (en) | Signal generator automatic test and calibration system and method | |
| CN106771709A (en) | A kind of S parameter De- embedding method of multiport network | |
| CN109633262A (en) | Three phase harmonic electric energy gauging method, device based on composite window multiline FFT | |
| CN103457603A (en) | ADC dynamic parameter testing method based on average frequency spectra | |
| CN105487034A (en) | 0.05-level electronic transformer verification method and system | |
| CN104270208B (en) | Method and device for detecting standing-wave ratio of RRU | |
| CN106443548A (en) | Calibration method of using electronic calibration piece to measure frequency conversion loss of mixer by vector network analyzer | |
| CN103217578A (en) | Digital phase meter and method for measuring phase difference between signals and based on phase-locked loop technology | |
| CN105262503A (en) | Group delay calibration based multipath delay generation device and method | |
| CN104483619B (en) | A kind of frequency characteristics measurement system based on virtual instrument | |
| CN106291105A (en) | A kind of sweep generator based on digital zero intermediate frequency | |
| CN103616566B (en) | A kind of linear average cymoscope implementation method based on software radio | |
| CN205880145U (en) | Test panel and testing arrangement are surveyed in TRL calibration of passive link | |
| CN117434426B (en) | Test method, system and device of switched capacitor filter | |
| CN109581062A (en) | High-precision impedance measurement system for oscilloscope calibrator probe | |
| CN107294533A (en) | Analog-digital converter dynamic parameter testing system and method | |
| CN105471525A (en) | Signal processing method for four-channel compressed sensing digital receiver of vector network analyzer | |
| CN101558319A (en) | RF circuit analysis | |
| CN108055091A (en) | A kind of millimeter wave self calibration virtual instrument and its implementation | |
| CN113092892B (en) | Performance testing device and method for electromagnetic interference noise separator | |
| CN204065391U (en) | A kind of wideband current transformer calibrator | |
| CN102201878B (en) | Method and apparatus for separating influences of signal source and analyzer on vector modulation error |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | 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 |
Effective date of registration: 20190314 Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee after: CHINA ELECTRONICS TECHNOLOGY INSTRUMENTS Co.,Ltd. Address before: 266555 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: THE 41ST INSTITUTE OF CHINA ELECTRONICS TECHNOLOGY Group Corp. |
|
| TR01 | Transfer of patent right | ||
| CP03 | Change of name, title or address |
Address after: Huangdao Xiangjiang Road 266555 Shandong city of Qingdao Province, No. 98 Patentee after: CLP kesiyi Technology Co.,Ltd. Address before: 266000 No. 98 Xiangjiang Road, Huangdao District, Shandong, Qingdao Patentee before: CHINA ELECTRONICS TECHNOLOGY INSTRUMENTS Co.,Ltd. |
|
| CP03 | Change of name, title or address | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |