CN116840940A - Channel calibration method and system based on millimeter wave cylindrical calibration body algorithm system - Google Patents
Channel calibration method and system based on millimeter wave cylindrical calibration body algorithm system Download PDFInfo
- Publication number
- CN116840940A CN116840940A CN202310671167.XA CN202310671167A CN116840940A CN 116840940 A CN116840940 A CN 116840940A CN 202310671167 A CN202310671167 A CN 202310671167A CN 116840940 A CN116840940 A CN 116840940A
- Authority
- CN
- China
- Prior art keywords
- channel
- receiving
- millimeter wave
- transmitting
- calibration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 34
- 238000012937 correction Methods 0.000 claims abstract description 32
- 238000007689 inspection Methods 0.000 claims abstract description 15
- 230000006870 function Effects 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 10
- 230000035485 pulse pressure Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 11
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000000605 extraction Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 3
- 239000002184 metal Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007723 transport mechanism Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The application provides a channel calibration method and a system based on a millimeter wave cylindrical calibration algorithm system, wherein the method comprises the following steps: and combining the Newton iteration method and the distance-direction self-adaptive windowing function technology, and iterating out the distance from the position of the receiving and transmitting array element to the cylindrical calibration body by taking the extraction of the effective signal of the echo of the cylindrical calibration body as a criterion, thereby further completing the amplitude phase consistency correction between and in the channels of the receiving and transmitting array element. The technical scheme of the embodiment of the application is suitable for one-dimensional linear array scanning system, two-dimensional electric scanning system and other scanning systems, can be applied to two-dimensional sparse area array configuration, linear array configuration, folding array configuration, arc array and other configurations, has wider application scene, and simultaneously has higher calculation precision and higher efficiency by adopting Newton iteration method, limits clutter interference by self-adaptive windowing technology, has higher subsequent image dynamic by amplitude and phase equalization technology, is convenient to carry in an external field, and is more suitable for millimeter wave human body security inspection application scene.
Description
Technical Field
The application relates to the technical field of millimeter wave imaging, in particular to a channel calibration method and system based on a millimeter wave cylindrical calibration body algorithm system.
Background
The active millimeter wave human body security inspection instrument can be widely applied to security application fields such as subways, airports, customs, public security inspection stations and the like. The millimeter wave human body security inspection instrument is required to have the characteristics of high pass rate, low false alarm rate, no perception and intellectualization in the current human body security inspection occasion. Because the human body security inspection system comprises a frequency source, a frequency converter, a radio frequency front end, a microwave cable and other microwave devices, the amplitude and phase consistency of the devices can drift after long-time working, the accuracy of system correction is reduced, the system correction is finally displayed on a display image, the image quality is deteriorated, further the deep learning intelligent block diagram is affected, the false alarm and false alarm phenomenon is caused, and the correction is needed through a correction body algorithm system. The related art mainly adopts a metal flat plate for calibration, but the portability of the metal flat plate is poor, and the calibration precision is low.
Disclosure of Invention
The present application aims to solve at least one of the technical problems existing in the prior art. Therefore, the application provides a channel calibration method and a system based on a millimeter wave cylindrical calibration body algorithm system, which adopt a cylindrical calibration body with higher portability, accelerate iterative optimization through a Newton iterative method, design a distance-oriented self-adaptive windowing function and have higher calibration precision.
In a first aspect, an embodiment of the present application provides a channel calibration method based on a millimeter wave cylindrical calibration body algorithm system, which is applied to a millimeter wave human body security inspection instrument, where the millimeter wave human body security inspection instrument includes a transceiver array, the transceiver array includes a plurality of transmitting array elements and a plurality of receiving array elements, a cylindrical calibration body is placed in a central position in front of an array face of the transceiver array, and a height of the cylindrical calibration body covers a height range of the transceiver array elements, and the channel calibration method based on the millimeter wave cylindrical calibration body algorithm system includes the following steps:
s1, determining a frequency band to be calibrated, and transmitting radio frequency signals by the transmitting array element according to the frequency band to be calibrated;
s2, receiving background echo data and cylindrical calibration volume echo data by a receiving array element;
s3, determining a channel receiving and transmitting array element distance according to a Newton iteration method, wherein the channel receiving and transmitting array element distance is the shortest distance between each receiving and transmitting array element and the reflecting surface of the cylindrical calibration body;
s4, determining a background cancellation signal, wherein the background cancellation signal is the difference between the cylindrical calibration volume echo data and the background echo data;
s5, performing inverse Fourier transform on the background cancellation signal to obtain a pulse pressure signal;
s6, acquiring the amplitude value and the maximum value index of each group of receiving and transmitting channels in the distance direction of the pulse signal, and carrying out self-adaptive windowing on the amplitude value and the maximum value index of each group of receiving and transmitting channels based on a preset self-adaptive windowing function to obtain pulse pressure signals;
s7, carrying out Fourier transformation on the pulse pressure signal to obtain a frequency domain signal;
s8, acquiring an inter-channel amplitude correction parameter and a preset amplitude correction factor, wherein the inter-channel amplitude correction parameter is the inverse of the frequency dimension power of each group of receiving and transmitting channels;
s9, determining an intra-channel correction parameter, wherein the intra-channel correction parameter is the conjugate of the phase of the receiving-transmitting channel distance and the phase of the channel;
s10, performing conjugate multiplication on the phase of each receiving and transmitting array element distance and the phase of the frequency domain signal, and multiplying the obtained result with the amplitude correction factor to finish correction of the amplitude and the phase in and between channels.
According to some embodiments of the application, the diameter of the cylindrical calibration body is 200mm, and the distance between the center of the cylindrical calibration body and the center of the array surface of the transceiver array is in the range of 0.5 to 0.8m.
According to some embodiments of the application, the expression of the Newton's iterative method isExpression of the channel receiving and transmitting array element distanceThe formula is:
where i is the transmit channel index, j is the receive channel index, (x) t ,y t ,z t ) For transmitting array element coordinates, (x) r ,y r ,z r ) To receive array element coordinates, (x) 0 ,y 0 ,z 0 ) For calibrating the center coordinates of the cylinder, R t For transmitting the distance from the array element coordinates to the reflecting surface of the cylindrical calibration body, R r For receiving the distance from the array element coordinates to the reflecting surface of the cylindrical calibration body, R i,j For the sum of the distance from the transmitting array element coordinates to the reflecting surface of the cylindrical calibration body and the distance from the receiving array element coordinates to the reflecting surface of the cylindrical calibration body, theta n Azimuth angle, θ, for the nth plane iteration n+1 Azimuth angle for n+1th plane iteration, h n For the height direction size of the nth iteration, h n+1 For the height direction size of the n +1 th iteration,representing a jacobian matrix, H (θ) n ,h n ) Representing the hessian matrix, the expression is: /> The expression is as follows: />
According to some embodiments of the application, in the step S6, the windowing function is a rectangular window function, and an expression of the windowing function is:
wherein pSigRangeIndex (i, j) is the index maximum of the transceiving channel.
According to some embodiments of the application, the expression of the background cancellation signal is: pSigCalF (i, j, f) =psigcalf (i, j, f) -pSigCalF (i, j, f), wherein pSigCalF (i, j, f) is the cylinder calibration volume echo data, pSigCalF (i, j, f) is the background echo data, wherein i represents a transmit channel index, j represents a receive channel index, and f represents a frequency point index.
According to some embodiments of the application, the amplitude correction factor is expressed as:
according to some embodiments of the application, the expression of the in-channel correction factor is:
wherein f is E [ f min ，f max ]C is the speed of light.
In a second aspect, an embodiment of the present application provides a channel calibration system based on a millimeter wave cylindrical calibration algorithm system, including at least one control processor and a memory for communication connection with the at least one control processor; the memory stores instructions executable by the at least one control processor to enable the at least one control processor to perform a channel calibration method based on a millimeter wave cylindrical calibration algorithm system as described in the first aspect above.
In a third aspect, an embodiment of the present application provides an electronic device, including a channel calibration system based on the millimeter wave cylindrical calibration algorithm system according to the second aspect.
In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium storing computer-executable instructions for performing the channel calibration method based on the millimeter wave cylindrical calibration algorithm system according to the first aspect.
The channel calibration method based on the millimeter wave cylindrical calibration algorithm system has at least the following beneficial effects: the compensation distance is accelerated and optimized through the Newton iteration method, the calculation accuracy is higher, the efficiency is faster, clutter interference is limited through the self-adaptive windowing technology, the follow-up image is higher in dynamic state through the amplitude and phase equalization technology, the outfield is convenient to carry, and the method is more suitable for millimeter wave human body security inspection application scenes.
Drawings
FIG. 1 is a flow chart of a method for channel calibration based on a millimeter wave cylindrical calibration algorithm system provided by an embodiment of the present application;
FIG. 2 is a schematic illustration of a millimeter wave uniarray arrangement provided in another embodiment of the present application;
FIG. 3 is a layout diagram of transceiver elements according to another embodiment of the present application;
FIG. 4 is an example of a three-dimensional view of a cylindrical calibration system provided in accordance with another embodiment of the present application;
FIG. 5 is a two-dimensional side view of a cylindrical calibration body system provided in accordance with another embodiment of the present application;
fig. 6 is a schematic structural diagram of an algorithm system based on millimeter wave cylindrical calibration according to another embodiment of the present application.
Detailed Description
Embodiments of the present application are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are illustrative only and are not to be construed as limiting the application.
In the description of the present application, it should be understood that the direction or positional relationship indicated with respect to the description of the orientation, such as upper, lower, front, rear, left, right, etc., is based on the direction or positional relationship shown in the drawings, is merely for convenience of describing the present application and simplifying the description, and does not indicate or imply that the terminal or element to be referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present application.
In the description of the present application, a number means one or more, a number means two or more, and greater than, less than, exceeding, etc. are understood to not include the present number, and above, below, within, etc. are understood to include the present number. The description of the first and second is for the purpose of distinguishing between technical features only and should not be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated or implicitly indicating the precedence of the technical features indicated.
In the description of the present application, unless explicitly defined otherwise, terms such as arrangement, installation, connection, etc. should be construed broadly and the specific meaning of the terms in the present application can be reasonably determined by a person skilled in the art in combination with the specific contents of the technical scheme.
According to the technical scheme of the embodiment of the application, the compensation distance is accelerated and optimized by the Newton iteration method, the calculation accuracy is higher, the efficiency is faster, clutter interference is limited by the self-adaptive windowing technology, the subsequent image is higher in dynamic state by the amplitude and phase equalization technology, the outfield is convenient to carry, and the method and the system are more suitable for millimeter wave human body security inspection application scenes.
The control method of the embodiment of the application is further described below with reference to the accompanying drawings.
Referring to fig. 1, fig. 1 is a flowchart of a channel calibration method based on a millimeter wave cylindrical calibration body algorithm system, which is provided by the embodiment of the application, and is applied to a millimeter wave human body security inspection instrument, wherein the millimeter wave human body security inspection instrument comprises a transceiver array, the transceiver array comprises a plurality of transmitting array elements and a plurality of receiving array elements, a cylindrical calibration body is placed in the front center position of an array face of the transceiver array, the height of the cylindrical calibration body covers the height range of the transceiver array elements, and the channel calibration method based on the millimeter wave cylindrical calibration body algorithm system comprises the following steps:
s1, determining a frequency band to be calibrated, and transmitting an RF signal by a transmitting array element according to the frequency band to be calibrated;
s2, receiving background echo data and cylindrical calibration volume echo data by a receiving array element;
s3, determining the channel receiving and transmitting array element distance according to the Newton iteration method, wherein the channel receiving and transmitting array element distance is the shortest distance between each receiving and transmitting array element and the reflecting surface of the cylindrical calibration body;
s4, determining a background cancellation signal, wherein the background cancellation signal is the difference between the cylindrical calibration volume echo data and the background echo data;
s5, performing inverse Fourier transform on the background cancellation signal to obtain a pulse pressure signal;
s6, acquiring the amplitude value and the maximum value index of each group of receiving and transmitting channels in the distance direction of the pulse signal, and carrying out self-adaptive windowing on the amplitude value and the maximum value index of each group of receiving and transmitting channels based on a preset self-adaptive windowing function to obtain a pulse pressure signal;
s7, carrying out Fourier transformation on the pulse pressure signal to obtain a frequency domain signal;
s8, acquiring an inter-channel amplitude correction parameter and a preset amplitude correction factor, wherein the inter-channel amplitude correction parameter is the inverse of the frequency dimension power of each group of receiving and transmitting channels;
s9, determining an intra-channel correction parameter, wherein the intra-channel correction parameter is the conjugate of the phase of the receiving-transmitting channel distance and the phase of the channel;
s10, performing conjugate multiplication on the phase of each receiving and transmitting array element distance and the phase of the frequency domain signal, and multiplying the obtained result by an amplitude correction factor to finish correction of the amplitude and the phase in and between channels.
It should be noted that, in the embodiment of the present application, the millimeter wave sub-array arrangement manner may be shown in fig. 2, the arrangement manner of the transceiver array elements may be shown in fig. 3, in fig. 2 and 3, the dots are the transmitting array element positions, and the stars are the receiving array element positions. Based on the structure shown in fig. 2 and fig. 3, the schematic position diagram obtained after the cylindrical calibration body is placed can be referred to as fig. 4 and fig. 5, and the cylindrical calibration body is placed in the center position in front of the array face of the transceiver array, and the height of the cylindrical calibration body covers the height range of the transceiver array element. As the millimeter wave cylindrical calibration algorithm system does not depend on the frequency band to be changed, the frequency band to be calibrated can be selected to be 30-38 GHz, wherein the frequency is divided into nFreq=128 points, and the center wavelength of the frequency bandThe array element spacing of the transmitting antenna is set to be d less than or equal to 0.75lambda c Setting the distance between the transmitting array element and the receiving array element to be d ant =6mm, the number of mono-array transmit elements is 48, the number of receive elements is 48; the horizontal coverage is 1m, the height direction coverage is 2m, the total number of subarrays is 32, and the total number of transmitting array elements is N t Let=1536, receive array element total N r ＝1536。
In some embodiments, referring to fig. 4 and 5, the diameter of the cylindrical calibration body is 200mm and the center of the cylindrical calibration body is in the range of 0.5 to 0.8m from the center of the array face of the transceiver array.
It should be noted that, the cylinder calibration body is placed in the front center position of the array, the diameter d=200mm of the cylinder calibration body, the center distance range of the cylinder calibration body from the array is 0.5 to 0.8m, the height of the cylinder should cover the height range of the receiving array element, and the effective receiving and transmitting pair of the receiving array element meeting the requirement is ensured.
In some embodiments, the expression of Newton' S iteration in step S3 isThe expression of the channel receiving and transmitting array element distance is as follows: />
Where i is the transmit channel index, j is the receive channel index, (x) t ,y t ,z t ) For transmitting array element coordinates, (x) r ,y r ,z r ) To receive array element coordinates, (x) 0 ,y 0 ,z 0 ) For calibrating the center coordinates of the cylinder, R t For transmitting the distance from the array element coordinates to the reflecting surface of the cylindrical calibration body, R r For receiving the distance from the array element coordinates to the reflecting surface of the cylindrical calibration body, R i,j For the sum of the distance from the transmitting array element coordinates to the reflecting surface of the cylindrical calibration body and the distance from the receiving array element coordinates to the reflecting surface of the cylindrical calibration body, theta n Azimuth angle, θ, for the nth plane iteration n+1 Is n+1st time flatAzimuth angle of plane iteration, h n For the height direction size of the nth iteration, h n+1 For the height direction size of the n +1 th iteration,representing a jacobian matrix, H (θ) n ,h n ) Representing the hessian matrix, the expression is: /> The expression is as follows: />
In order to improve the calculation speed, the millimeter wave human body security inspection instrument can utilize the GPU to perform acceleration calculation, so that the calculation time is saved.
It is noted that in this embodiment, the jacobian matrix takes the first partial derivative and the hessian matrix takes the second partial derivative.
In addition, in an embodiment, the windowing function in step S6 is a rectangular window function, and the expression of the windowing function is:
wherein pSigRangeIndex (i, j) is the index maximum of the transceiving channel.
It should be noted that, by the adaptive windowing technology, clutter interference can be effectively limited, and by the amplitude and phase equalization technology, the subsequent image dynamic is higher.
In addition, in one embodiment, the expression for the background cancellation signal is: pSigCalF (i, j, f) =psigcalf (i, j, f) -pSigCalF (i, h, f), wherein pSigCalF (i, h, f) is the cylinder calibration volume echo data, pSigCalF (i, h, f) is the background echo data, wherein i represents a transmit channel index, j represents a receive channel index, and f represents a frequency point index.
In addition, in one embodiment, the expression of the amplitude correction factor is:
in addition, in one embodiment, the expression of the in-channel correction factor is:
wherein f is E [ f min ，f max ]C is the speed of light.
As shown in fig. 6, fig. 6 is a block diagram of a channel calibration system based on a millimeter wave cylindrical calibration algorithm system according to an embodiment of the present application. The application also provides a channel calibration system based on the millimeter wave cylinder calibration body algorithm system, which comprises:
the processor 601 may be implemented by a general-purpose central processing unit (Central Processing Unit, CPU), a microprocessor, an application specific integrated circuit (Application Specific Integrated Circuit, ASIC), or one or more integrated circuits, etc. for executing related programs, so as to implement the technical scheme provided by the embodiments of the present application;
the Memory 602 may be implemented in the form of a Read Only Memory (ROM), a static storage system, a dynamic storage system, or a random access Memory (Random Access Memory, RAM). The memory 602 may store an operating system and other application programs, and when the technical solution provided in the embodiments of the present disclosure is implemented by software or firmware, relevant program codes are stored in the memory 602, and the processor 601 invokes a channel calibration method based on the millimeter wave cylindrical calibration algorithm system to execute the embodiments of the present disclosure;
an input/output interface 603 for implementing information input and output;
the communication interface 604 is configured to implement communication interaction between the present system and other systems, and may implement communication in a wired manner (e.g. USB, network cable, etc.), or may implement communication in a wireless manner (e.g. mobile network, WIFI, bluetooth, etc.);
a bus 605 for transferring information between the various components of the system (e.g., the processor 601, memory 602, input/output interfaces 603, and communication interfaces 604);
wherein the processor 601, the memory 602, the input/output interface 603 and the communication interface 604 are communicatively coupled to each other within the system via a bus 605.
The embodiment of the application also provides electronic equipment, which comprises the channel calibration system based on the millimeter wave cylindrical calibration body algorithm system.
The embodiment of the application also provides a storage medium, which is a computer readable storage medium, and the storage medium stores a computer program which is executed by a processor to realize the channel calibration method based on the millimeter wave cylindrical calibration algorithm system.
The memory, as a non-transitory computer readable storage medium, may be used to store non-transitory software programs as well as non-transitory computer executable programs. In addition, the memory may include high-speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory optionally includes memory remotely located relative to the processor, the remote memory being connectable to the processor through a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof. The terminal embodiments described above are merely illustrative, in which the elements illustrated as separate components may or may not be physically separate, implemented to reside in one place, or may be distributed over multiple network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
Those of ordinary skill in the art will appreciate that all or some of the steps, systems, and methods disclosed above may be implemented as software, firmware, hardware, and suitable combinations thereof. Some or all of the physical components may be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application specific integrated circuit. Such software may be distributed on computer readable media, which may include computer storage media (or non-transitory media) and communication media (or transitory media). The term computer storage media includes both volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data, as known to those skilled in the art. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage terminals, or any other medium which can be used to store the desired information and which can be accessed by a computer. Furthermore, as is well known to those of ordinary skill in the art, communication media typically include computer readable instructions, data structures, program modules, or other data in a modulated data signal, such as a carrier wave or other transport mechanism, and may include any information delivery media.
While the preferred embodiment of the present application has been described in detail, the present application is not limited to the above embodiments, and those skilled in the art can make various equivalent modifications or substitutions without departing from the spirit and scope of the present application, and these equivalent modifications or substitutions are included in the scope of the present application as defined in the appended claims.
Claims (10)
1. The channel calibration method based on the millimeter wave cylindrical calibration body algorithm system is characterized by being applied to a millimeter wave human body security inspection instrument, wherein the millimeter wave human body security inspection instrument comprises a receiving and transmitting array, the receiving and transmitting array comprises a plurality of transmitting array elements and a plurality of receiving array elements, a cylindrical calibration body is placed in the center position in front of an array face of the receiving and transmitting array, the height of the cylindrical calibration body covers the height range of the receiving and transmitting array elements, and the channel calibration method based on the millimeter wave cylindrical calibration body algorithm system comprises the following steps:
s1, determining a frequency band to be calibrated, and transmitting radio frequency signals by the transmitting array element according to the frequency band to be calibrated;
s2, receiving background echo data and cylindrical calibration volume echo data by a receiving array element;
s3, determining a channel receiving and transmitting array element distance according to a Newton iteration method, wherein the channel receiving and transmitting array element distance is the shortest distance between each receiving and transmitting array element and the reflecting surface of the cylindrical calibration body;
s4, determining a background cancellation signal, wherein the background cancellation signal is the difference between the cylindrical calibration volume echo data and the background echo data;
s5, performing inverse Fourier transform on the background cancellation signal to obtain a pulse pressure signal;
s6, acquiring the amplitude value and the maximum value index of each group of receiving and transmitting channels in the distance direction of the pulse signal, and carrying out self-adaptive windowing on the amplitude value and the maximum value index of each group of receiving and transmitting channels based on a preset self-adaptive windowing function to obtain pulse pressure signals;
s7, carrying out Fourier transformation on the pulse pressure signal to obtain a frequency domain signal;
s8, acquiring an inter-channel amplitude correction parameter and a preset amplitude correction factor, wherein the inter-channel amplitude correction parameter is the inverse of the frequency dimension power of each group of receiving and transmitting channels;
s9, determining an intra-channel correction parameter, wherein the intra-channel correction parameter is the conjugate of the phase of the receiving-transmitting channel distance and the phase of the channel;
s10, performing conjugate multiplication on the phase of each receiving and transmitting array element distance and the phase of the frequency domain signal, and multiplying the obtained result with the amplitude correction factor to finish correction of the amplitude and the phase in and between channels.
2. The method for calibrating a channel based on a millimeter wave cylindrical calibration body algorithm system according to claim 1, wherein the diameter of the cylindrical calibration body is 200mm, and the distance between the center of the cylindrical calibration body and the center of the array surface of the transceiver array is in the range of 0.5 to 0.8m.
3. The method for calibrating a channel based on a millimeter wave cylindrical calibration algorithm system according to claim 1, wherein the expression of the newton's iterative method is thatThe expression of the channel receiving and transmitting array element distance is as follows: r is R i,j (θ,h)＝R t (r c ,θ,h)+R r (r c ,θ,h)，
Where i is the transmit channel index, j is the receive channel index, (x) t ,y t ,z t ) For transmitting array element coordinates, (x) r ,y r ,z r ) To receive array element coordinates, (x) 0 ,y 0 ,z 0 ) For calibrating the center coordinates of the cylinder, R t For transmitting the distance from the array element coordinates to the reflecting surface of the cylindrical calibration body, R r For receiving the distance from the array element coordinates to the reflecting surface of the cylindrical calibration body, R i,j For the sum of the distance from the transmitting array element coordinates to the reflecting surface of the cylindrical calibration body and the distance from the receiving array element coordinates to the reflecting surface of the cylindrical calibration body, theta n Azimuth angle, θ, for the nth plane iteration n+1 Azimuth angle for n+1th plane iteration, h n For the height direction size of the nth iteration, h n+1 For the height direction size of the n +1 th iteration,representing a jacobian matrix, H (θ) n ,h n ]Representing the hessian matrix, the expression is: /> The expression is as follows: />
4. The method for calibrating a channel based on the millimeter wave cylindrical calibration algorithm system according to claim 1, wherein in the step S6, the windowing function is a rectangular window function, and the expression of the windowing function is:
wherein pSigRangeIndex (i, j) is the index maximum of the transceiving channel.
5. The method for calibrating a channel based on a millimeter wave cylindrical calibration algorithm system according to claim 1, wherein the expression of the background cancellation signal is: pSigCalF (i, j, f) =psigcalf (i, j, f) -pSiBackF (i, j, f), wherein pSigCalF (i, j, f) is the cylinder calibration volume echo data, pSigCalF (i, j, f) is the background echo data, wherein i represents a transmit channel index, j represents a receive channel index, and f represents a frequency bin index.
6. The method for calibrating a channel based on a millimeter wave cylindrical calibration algorithm system according to claim 5, wherein the expression of the amplitude correction factor is:
7. the method for calibrating a channel based on a millimeter wave cylindrical calibration algorithm system according to claim 5, wherein the expression of the correction factor in the channel is:
wherein f is E [ f min ，f max ]C is the speed of light.
8. A channel calibration system based on a millimeter wave cylindrical calibration algorithm system, comprising at least one control processor and a memory for communication connection with the at least one control processor; the memory stores instructions executable by the at least one control processor to enable the at least one control processor to perform the millimeter wave cylinder calibrator algorithm-based channel calibration method of any of claims 1 to 7.
9. An electronic device comprising the millimeter wave cylindrical calibration system-based channel calibration system of claim 8.
10. A computer-readable storage medium storing computer-executable instructions for causing a computer to perform the millimeter wave cylinder calibrator algorithm-based channel calibration method according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310671167.XA CN116840940B (en) | 2023-06-07 | 2023-06-07 | Channel calibration method and system based on millimeter wave cylindrical calibration body algorithm system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310671167.XA CN116840940B (en) | 2023-06-07 | 2023-06-07 | Channel calibration method and system based on millimeter wave cylindrical calibration body algorithm system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116840940A true CN116840940A (en) | 2023-10-03 |
| CN116840940B CN116840940B (en) | 2024-04-05 |
Family
ID=88171571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310671167.XA Active CN116840940B (en) | 2023-06-07 | 2023-06-07 | Channel calibration method and system based on millimeter wave cylindrical calibration body algorithm system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116840940B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050128871A1 (en) * | 2003-12-12 | 2005-06-16 | Fernando Garcia-Osuna | Methods and systems for calibrating acoustic receivers |
| CN106597556A (en) * | 2016-12-09 | 2017-04-26 | 北京无线电计量测试研究所 | Method for background elimination of airport runway foreign and debris detection system |
| CN109633630A (en) * | 2018-10-31 | 2019-04-16 | 上海无线电设备研究所 | A kind of E wave band millimeter wave imaging safety check radar system |
| CN109932719A (en) * | 2019-03-18 | 2019-06-25 | 西安电子科技大学 | RCS high-precision measuring method based on SAR imaging |
| CN111025423A (en) * | 2019-12-28 | 2020-04-17 | 北京无线电计量测试研究所 | Calibration system and method of active millimeter wave real-time three-dimensional imaging security inspection system |
| US20200142030A1 (en) * | 2017-04-28 | 2020-05-07 | Shenzhen Institute Of Terahertz Technology And Innovation | Amplitude-phase correction method and system for microwave imaging system |
| CN111679254A (en) * | 2020-05-13 | 2020-09-18 | 苏州理工雷科传感技术有限公司 | Multichannel calibration method for linear frequency modulation signal MIMO radar system |
| CN111708023A (en) * | 2020-08-20 | 2020-09-25 | 湖南雷远电子科技有限公司 | Phase inconsistency correction method and imaging system for millimeter wave multi-transmitting multi-receiving linear array |
| CN111999782A (en) * | 2020-07-30 | 2020-11-27 | 北京遥测技术研究所 | Automatic calibration method of rotary scanning imaging system |
| CN112764116A (en) * | 2020-12-24 | 2021-05-07 | 博微太赫兹信息科技有限公司 | Rapid imaging method of sparse array sparse frequency point planar scanning system |
-
2023
- 2023-06-07 CN CN202310671167.XA patent/CN116840940B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050128871A1 (en) * | 2003-12-12 | 2005-06-16 | Fernando Garcia-Osuna | Methods and systems for calibrating acoustic receivers |
| CN106597556A (en) * | 2016-12-09 | 2017-04-26 | 北京无线电计量测试研究所 | Method for background elimination of airport runway foreign and debris detection system |
| US20200142030A1 (en) * | 2017-04-28 | 2020-05-07 | Shenzhen Institute Of Terahertz Technology And Innovation | Amplitude-phase correction method and system for microwave imaging system |
| CN109633630A (en) * | 2018-10-31 | 2019-04-16 | 上海无线电设备研究所 | A kind of E wave band millimeter wave imaging safety check radar system |
| CN109932719A (en) * | 2019-03-18 | 2019-06-25 | 西安电子科技大学 | RCS high-precision measuring method based on SAR imaging |
| CN111025423A (en) * | 2019-12-28 | 2020-04-17 | 北京无线电计量测试研究所 | Calibration system and method of active millimeter wave real-time three-dimensional imaging security inspection system |
| CN111679254A (en) * | 2020-05-13 | 2020-09-18 | 苏州理工雷科传感技术有限公司 | Multichannel calibration method for linear frequency modulation signal MIMO radar system |
| CN111999782A (en) * | 2020-07-30 | 2020-11-27 | 北京遥测技术研究所 | Automatic calibration method of rotary scanning imaging system |
| CN111708023A (en) * | 2020-08-20 | 2020-09-25 | 湖南雷远电子科技有限公司 | Phase inconsistency correction method and imaging system for millimeter wave multi-transmitting multi-receiving linear array |
| CN112764116A (en) * | 2020-12-24 | 2021-05-07 | 博微太赫兹信息科技有限公司 | Rapid imaging method of sparse array sparse frequency point planar scanning system |
Non-Patent Citations (1)
| Title |
|---|
| 李春林;赵怀坤;吴琳拥;: "利用雷达发射信号对接收阵列实时校正的方法", 电子技术与软件工程, no. 04, 15 February 2020 (2020-02-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116840940B (en) | 2024-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111093267B (en) | IRS-based UE position determination method, communication method and system | |
| CN108139473B (en) | Angle and position sensing using antenna arrays | |
| CN110988862B (en) | Ultra-close millimeter wave radar sensing method and system | |
| US11735828B2 (en) | 24-GHz low-cost continuous beam steering phased array for indoor smart radar and methods relating thereto | |
| US20230065110A1 (en) | Electromagnetic non-line-of-sight imaging method based on time reversal and compressed sensing | |
| US9055523B2 (en) | Apparatus, system and method of calibrating a radio delay of a wireless device | |
| EP3407082B1 (en) | Apparatus and method for determining a distance to an object | |
| EP4266729A1 (en) | Wireless local area network sensing method and apparatus | |
| EP3279687A1 (en) | Beam signal tracking method, device and system | |
| EP3004921B1 (en) | Device, system and method of determining whether a mobile device is located in an indoor location or an outdoor location | |
| CN111679254B (en) | Multichannel calibration method for linear frequency modulation signal MIMO radar system | |
| CN104079307A (en) | Electronic device capable of eliminating interference | |
| US10634780B2 (en) | Electromagnetic wave imaging system and antenna array signal correction method | |
| CN104777467B (en) | Object detection method based on frequency scan antenna | |
| CN113093096A (en) | Method for determining arrival angle and related device | |
| CN112995888B (en) | Positioning method and system based on array antenna, electronic equipment and storage medium | |
| CN114577171A (en) | Three-dimensional angle measuring method, device and equipment and computer storage medium | |
| EP3836454B1 (en) | Method and device for improving phase measurement accuracy | |
| CN116840940B (en) | Channel calibration method and system based on millimeter wave cylindrical calibration body algorithm system | |
| CN109085565B (en) | Miniaturized millimeter wave radio frequency test platform | |
| CN117388810A (en) | One-dimensional millimeter wave sparse linear array correction method, system, equipment and medium | |
| CN111082844B (en) | Side lobe suppression method of microwave direction finding equipment, direction finding method and microwave direction finding equipment | |
| CN114720937A (en) | Method for determining angle of arrival and related device | |
| CN107907880A (en) | One kind is spaceborne to condense polarization SAR fuzziness analysis system, method and memory | |
| Sammartino et al. | Phase compensation and processing in multiple-input–multiple-output radars |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |