3D target location based on RFID polarization phase model

The three-dimensional target position is vital to automatic driving, which can efficiently detect the geometrical shape of obstacles and achieve automatic environment recognition for both indoor and outdoor scenarios. In this paper, we propose an RFID 3D localization prototype based on the polarization phase model. To eliminate the impacts of tag angle, we leverage the polarization effect to obtain polarized phases. After that, we use polarized phase decomposition to get pure phases, i.e., the phase corresponding to the line of sight path. Hence, we can calculate the raw phases using the HMFCW algorithm. Finally, we can estimate the target’s 3D position according to unwrapped phase information. In this way, we implement a low-cost, high-precision, easy-for-deployment 3D position solution using commercial UHF RFID hardware. The experiments show that the method can obtain effective positioning results in various environments and achieve positioning accuracy with an average error of about 9 cm.

make a relatively good positioning effect at a low cost. It can transform the positioning of objects into the positioning of tags, simplify the theoretical complexity of positioning and improve the positioning efficiency.
In view of the wide application prospect and implementation value of indoor positioning, especially three-dimensional positioning technology, this paper proposes a low-cost, convenient and efficient three-dimensional positioning method. Based on the existing commercial RFID equipment, by constructing the phase model based on polarization effect and combined with the signal propagation model in space, this paper proposes a three-dimensional positioning method based on pure signal extraction to achieve low cost. Highly reliable Indoor 3D positioning method.
The main work of this paper is based on RFID in three-dimensional space. The main problems and solutions are as follows: First, for the problems of high cost and difficult popularization in current positioning, it is necessary to use low-cost and large-scale deployable RFID tags as a substitute for target positioning to bind with the target object. Considering that the information that RFID tag can receive includes received signal strength, phase, timestamp and Doppler frequency shift, and the phase data have the highest accuracy, phase information is used as the basic parameter of positioning. Secondly, the problem of phase folding should be solved first when using phase for positioning. Because the phase changes periodically, the phase in a cycle can only be used to locate objects within more than ten centimeters, and the positioning range is too small. Therefore, HMFCW algorithm needs to be used to unwrap to improve the range of phase positioning. Thirdly, the phase information from the tag received by the reader receiving antenna contains not only the tag information, but the result of reflection superposition of multiple paths in space. Therefore, it is necessary to separate the pure phase information from the mixed information. Finally, because the existing phase models cannot describe the impact of label angle on phase, and label angle does affect phase in actual experiments, it is necessary to reconstruct the phase model related to label angle to ensure the accuracy of the phase model. The method proposed in this paper is mainly divided into the following two steps: (1) Theoretical phase modeling: obtain the phase information of the tag to be tested from the existing commercial reader equipment, establish a polarization-based phase model for the tag phase and preprocess the phase with deconvolution and discrete wavelet transform denoising; (2) Using phase positioning: decompose the phase on the basis of (1), extract the pure phase of the decomposed phase, eliminate the complex multipath information from the measured phase, use HMFCW algorithm to unwrap the phase, and finally, complete the positioning.
The main part of this paper is the process of reshaping the theoretical model according to the experiment. The first part is the modeling part, that is, the phase model based on the polarization effect. According to the experimental phenomena obtained in the experiment, it is found that the label angle has a great impact on the phase, but this is not mentioned in the existing phase model, therefore, based on the existing phase model, a phase model based on polarization effect is proposed, and the influence of label angle on phase is explained in detail. The second part is the positioning part. According to the theoretical model proposed above, it is found in the actual experiment that the phase information provided by the existing commercial RFID equipment cannot directly match the theoretical model. The existence of multiple reflectors in the environment has a great impact on the experimental results, that is, the multipath effect. Therefore, this paper proposes a positioning algorithm based on pure phase extraction. Aiming at the multipath information in complex environments, this paper theoretically analyzes the relationship between pure phase and multipath phase, uses the collected phase information of different frequencies to fit the theoretical phase and corrects the error to ensure the accuracy of the pure phase. Finally, the pure phase at different frequencies is unwrapped to locate the label.
This paper designs a large number of experiments, modifies the theoretical model through the actual collected experimental data and designs a series of algorithms to reduce noise and error as much as possible. Finally, the feasibility of the algorithm is verified by experiments. Compared with the existing work, the research innovations and contributions of this paper are as follows: (1) By analyzing the phenomenon of signal propagation and polarization in RFID system, a phase model based on the polarization effect is proposed, which accurately reflects the influence of tag angle on phase. (2) By analyzing the theoretical relationship between actual phase information, pure phase and multipath phase, a positioning algorithm based on pure phase extraction is proposed, which eliminates the multipath phase from the phase obtained by commercial RFID readers, and solves the problem of the significant error increase of traditional positioning methods in multipath environment. (3) There is no need to refer to labels or anchor points, and there is no need to collect a large number of training data, which improves the convenience of system deployment. (4) Using commercial RFID equipment, a three-dimensional positioning algorithm with easy deployment, high precision and low cost is realized, and the average positioning error is about 9 cm.

System model based on polarization effect
In the previous section, we briefly introduce our basic idea and contribution. And then, we will describe the polarization effect and the phase model based on this effect. In order to solve the influence of label angle on phase ignored in the existing phase model, a phase model based on polarization effect is proposed in this paper.
Because different label directions lead to different phases, the existing phase model cannot accurately quantify the relationship between phase and distance, and the fundamental reason for this phenomenon is polarization effect.
In order to explain the influence of tag direction on phase, this paper theoretically analyzes the signal propagation model of RFID. There can be a circular polarization between two antennas π/2 phase difference. If the directions of two orthogonal linear polarization antennas are marked as sum and the polarization direction of the label is marked as, the signal sent by the transmitting antenna is u, v and w: When the transmitted signal reaches the tag through propagation, the signal received at the tag end S ′ (t) is where t d is the signal propagation time/s. After the transmitted signal is printed on the tag, the tag collects energy, activates the working circuit and modulates its own information and polarization direction on the received signal, and the modulated received signal S ′ (t) is where θ t is the phase offset/rad caused by hardware defect of tag.
After modulating the information, the tag will backscatter the signal into space. After the receiving antenna captures the signal of the tag in space, it can be decomposed into linear polarization components in two orthogonal directions. The signal decomposed by the receiving antenna in each direction vector is: The signal of the receiving antenna received by the reader is the superposition of the signals of the two linearly polarized antennas. Since the circularly polarized antenna is abstracted into two phases, the phase difference is π/2 . The line polarized antenna, therefore, actually receives the signal can be R(t) written as: where θ t is the phase offset/rad due to tag antenna hardware and θ r is the phase offset/ rad due to reader antenna hardware.
Equation (5) describes the influence of the label angle on the signal in the actual received signal, and the phase model based on polarization effect can be obtained where is the wavelength of electromagnetic wave/M; d is distance from antenna to tag/M; u, v are unit direction vectors of reader antenna; w is the unit direction vector of the label.
Note that due to the hardware limitation, e.g., ImpinJ R420 has a phase resolution about 0.0015 radians, the theoretical ranging resolution for COTS RFID readers is
Different from the existing phase model, this phase model based on polarization effect takes the angle of the label in three dimensions into account, which effectively reduces the error caused by the label angle. As shown in Fig. 1, it is a schematic diagram of antenna placement and coordinate selection in the system. The x-axis parallel to the ground in the X-Y plane of the top antenna is defined as the x-axis of the space, the positive half axis is the horizontal right direction. The y-axis perpendicular to the ground is defined as the Y-axis of the space, the positive half axis is the vertically upward direction. And the axis perpendicular to the X-Y plane is defined as the z-axis, the positive half axis is the direction in which the antenna faces.
In this way, the polarization decomposition vector coordinates of three antennas can be obtained where v vector of all three antennas is (0,1,0), u vector of top antenna . The antenna deployment mode is not unique, as long as the polarization decomposition vector of three antennas can cover three axes in the space. However, in order to reduce the computational complexity and improve the execution efficiency of the system, the above deployment mode is selected in this paper. The number of antennas can also continue to increase. The increase in the number can continue to improve the positioning accuracy, but it will increase the cost. In this paper, three antennas are the minimum requirements for positioning, so three antennas are used in the process of the experiment.

Methods
In the previous section, we state the polarization effect and the phase model based on this effect. In this section, we will detail the method about how to calculate the pure phase and how to obtain the measured phase corresponding to the distance.

Pure phase calculation
In the previous section, theoretical analysis demonstrates that there is no direct linear relationship between the measured phase β , pure phase θ and multipath phase α . In order to simplify the computational complexity, two new variables, multipath variable α and mirror phase θ , are selected to establish a linear relationship with the measured Fig. 1 Schematic diagram of antenna triaxial coordinates phase β . In this section, we will demonstrate how to utilize the internal relationship under multichannel of the measured phase β to calculate the actual value of the mirror phase θ .
First, the measured phase in multichannel β has the following relationship: Assuming that a reader has N channels, the measurement phase β n under each channel n can be expressed with: where θ is the mirror phase/rad of the first channel; α n is the multipath variable/rad of the nth channel.
The phase change of the signal will change with the change of frequency. Even under the same environmental reflection, the phase at different frequencies is still different. Moreover, path loss and reflection attenuation are also related to frequency. The theoretical path loss in free space is: where d is the path length/M; f is the frequency/Hz.
The reflection attenuation in free space is: where Z in is impedance/ and it is frequency dependent. According to Eqs. (8) and (9), both factors affecting multipath effect will change with frequency. Therefore, multipath variable α n is also different in different channels and varies nonlinearly. However, the phase offset between two adjacent channels �θ varies linearly. It consists of three parts, including the phase difference affected by distance �θ d , phase difference affected by reader antenna hardware �θ A and phase difference affected by tag antenna hardware �θ T . When the channel changes from n to n + 1 , the wavelength will also change, so the received signal will also produce a phase difference affected by the wavelength at the same distance �θ D . According to the above theoretical equation, the three factors that can impact �θ are all affected by the frequency change. Hence, when the frequency changes linearly in different channels, the phase offset �θ will also vary linearly.
In order to prove the above theory, the following experiments are done in this paper. A label is placed in two environments with open space and dense multipath, and we measure the phase of each of the 16 channels under the conditions of the same distance, the same hardware equipment and the same label angle β n . As shown in Fig. 2, the phase information measured in the open area is roughly on a straight line, because the multipath impact of the environment in the open area is minor and can be almost ignored, that is α n ≪ �θ . Therefore, the following expression can be obtained: where y ′ is the ideal fitting line.
The ideal fitting line can be obtained only when there is no multipath interference, but the phase value measured in the complex actual experimental environment does not necessarily meet the fitting line due to the interference of the multipath effect. In the complex multipath environment, due to the influence of the multipath effect, the measured phase value and frequency do not meet the direct linear relationship. A shake can be observed between channel 3 and channel 8 as well as some small inconsistency along other channels. This is caused by frequency selective fading. It is because, under such complex multipath environments, signals of different central frequencies may travel through different paths causing their propagation distance to vary from channel to channel. Therefore, the measured phase obtained in the actual experiment is not necessarily accurate, and even there will be large errors. In order to obtain an accurate mirror phase from the measured phase θ , quantitative analysis of the measured phase and mirror phase is required, and we can obtain: Obviously, Eq. (11) is a system of non-homogeneous linear equations, in which there are n + 2 unknowns and n equations, so there are countless solutions of unknowns X. In order to find the effective solution of X, two additional conditions need to be established to fit this set of measurements β n to a straight line. Hence, an actual fitting line is defined in this paper y n = κ · n + d , the purpose of this line is to be as close as possible to the fitting line in the ideal case. Therefore, the difference between the actual fitting line and the ideal fitting line is also defined. The equation is as follows: where ω n is the weight of the nth channel.
Because the frequency of each channel is different, the multipath effect of each channel is also different. Therefore, it is necessary to define the channel weight function to reduce the influence of the possible outliers of some channels on the fitting line. The standard of designing channel weight function is very simple, that is, to weaken the influence of the multipath effect as much as possible. The generation of the multipath effect is often composed of two parts, dynamic reflector (such as moving objects and people) and static reflector (such as walls, furniture and ceiling in a room). Compared with the static reflector, the dynamic reflector will produce more uncontrollable and unpredictable errors on the measurement phase, while the static reflector will have a relatively stable error distribution, and the more serious multipath effect will produce more discrete points. Therefore, the purpose of designing the channel weight function is to eliminate the reflection of dynamic multipath as much as possible and estimate the influence of static multipath signal. Therefore, multipath variables can be α n It is divided into two parts, the static multipath component composed of the reflection superposition of the static reflector α s n And a dynamic multipath component composed of the reflection superposition of the dynamic reflectorα d n , i.e., where β 0 n is the measurement phase/rad of the nth channel; T is the number of samples; and overlineβ n is the average value/rad of measured phase under t samples.
After the above analysis, the weight function can be obtained as where p n is the degree of dispersion of the nth channel for all channels.
Degree of dispersion p n The measurement equation is: In order to make the result of the weight function easy to solve and improve the calculation efficiency, the sum of all weight functions ω n is equal to the number of channels N. Under this calculation method, the influence of the channel affected by the unpredictable and uncontrollable dynamic multipath on the final result is effectively reduced. And (15) the dynamic multipath effect generally obeys Gaussian distribution, so it can be considered as N n=1 ω n ·α d n when the number of samples is large enough. This weight function can be used to obtain more accurate fitting results.
As Eq. (15) suggests, when the error φ gets its minimum, the slope κ and intercept d closest to the ideal fitting line can be obtained. Therefore, by demanding φ on the slope κ and intercept d partial derivatives, we can obtain: where e 1 is the difference of slope between the fitting line and the ideal line; e 2 is the difference of intercept between the fitted line and the ideal line. After that, we leverage linear regression (LR) to estimate the original phase. Although the phase error does not follow the Gaussian distribution, we can feel comfortable using ordinary least squares (OLS) since Gaussian-Markov theorem (GMT) [14] guarantees it is the best linear unbiased estimator (BLUE) in this case. The GMT states that the OLS estimator has the lowest sampling variance within the class of linear unbiased estimators, if the errors in the linear regression model are uncorrelated and have equal variances and expectation value of zero. We know that the measurement errors are independent. The mean residual of phase error is about 10 −7 . The residual can pass the Breusch-Pagan test, which means equal variances. As shown in Fig. 3, the LR line, which uses OLS, and the LAR line, which uses the least absolute residual solution for robust estimation, almost perfectly coincide. Hence, we leverage OLS to estimate the phase for convenience.
In this system, the calculated tag angle is used to obtain the phase model based on the polarization effect on the three antennas. Then the measured phase actually received is phase decomposed and pure phase calculated. The multipath information in the environment is eliminated, and finally, the distance is calculated by using the phase difference at different frequencies. (19)

Phase unwrapping
In the RFID system, the measured phase value received by the reader is actually a roll phase, that is: where θ true is the true phase/rad; θ is the measured phase/rad. The difference between the true phase and the measured phase is going to be k × 2π . Therefore, if the measured phase is directly used for positioning, the effective positioning distance does not exceed half wavelength, about 16 cm, which is far less than the distance requirements of general positioning. So, the calculation of k is required. This process is called phase unwrapping, also known as integer ambiguity resolution. In this paper, HMFCW algorithm is used to solve the integer ambiguity. [7] As shown in Fig. 4, the schematic diagram of HMFCW is provided. Under the n different frequency, we sample the phase of the same tag. Ideally, the theoretical distance between the tag and the reader antenna is: where i is the first i frequency; d h is the integer wavelength distance calculated from integer cycles/M; d m is the distance calculated from measured phase/M; k i is the first i integer ambiguity at two frequencies; i is the first i wavelength at two frequencies/M; and θ i is the first i measured phase/rad at two frequencies.
The overall idea of HMFCW is getting n measurement of phase under n different frequency and calculating the distance corresponding to each phase measurement as θ i 2π × i . Then adding 0, 1, 2...N wavelengths, respectively, to obtain n + 1 distance values. Because the real distance is fixed, the values with similar distances obtained at different frequencies can be grouped into a cluster, and the group with the best clustering result can be selected as the ranging result, so as to realize positioning. However, the phase obtained in the actual scene often has a certain error, so it is necessary to have a certain error tolerance for the phase, and the phase tolerance is proposed �θ i ; the theoretical distance in the real scene can be expressed as:

Fig. 4 Schematic diagram of HMFCW
where d is the true distance/M and �θ i is the frequency Phase tolerance under f i . It is assumed that the absolute value of phase tolerance at all frequencies is less than an upper limit tolerance �θ max ; then, subtract the theoretical distance at the two frequencies to obtain: where n i is the number of full cycles at ith frequency and n j is the number of full cycles at jth frequency. According to Eq. (23), if �θ max is too large, there will be multiple clusters of whole weeks that meet the filtering criteria. Therefore, if you want to have a unique solution to Eq. (23), you need to limit �θ max as follows: where f max is the maximum frequency of frequency conversion signal/MHz and f min is the minimum frequency of frequency conversion signal/MHz. Equation (24) shows the relationship between phase tolerance and frequency. The distance between two frequencies can be obtained by making a difference: where R max is the maximum ranging range/m.
Equation (25) shows that when the measurement phase error at all frequencies is less than �θ max , there must be a unique whole cycle solution, and then, the positioning result is calculated. However, in the actual measurement, the phase error may still be greater than the phase tolerance. At this time, it is necessary to slowly expand the phase tolerance to ensure that there is a unique whole cycle cluster to meet the conditions.
The size of phase error is an important index affecting the accuracy of HMFCW algorithm. Once some frequencies are greatly offset under the influence of complex multipath environments, it will directly have a great impact on the positioning results. Therefore, it is necessary to study the impact of different sizes of phase error on ranging accuracy.
Suppose there are two frequencies f i and f j , the wavelengths corresponding to these two frequencies are i and j , and f i < f j , i > j . Denote the corresponding phase errors as −�θ max and �θ max . Therefore, the distances calculated by these two frequencies can be subtracted and get: Equation (26) shows that when the phase error is �θ max . It is also the phase tolerance, and there must be a unique solution for the whole cycle. So when the phase error is less than �θ max , a unique solution must exist. When the whole cycle number cluster needs to be solved for positioning, the whole cycle number of the similar distance at each other frequency needs to be calculated based on the distance corresponding to the maximum frequency. Let the maximum frequency be f max . The wavelength at this frequency is min . The corresponding phase error at the maximum frequency is �θ max . The corresponding number of whole cycles at frequency f i is: where r is the whole cycle difference count between f i and f max . Because the number of whole cycles at different frequencies will be different, r needs to be calculated. The calculation equation is: where round is the rounding function.
The real distance between the antenna and the tag can be obtained after the above integer cycle processing. The whole steps of the positioning algorithm based on pure phase extraction are summarized as shown in Algorithm 1:

Results and discussion
In this section, the above methods are tested in different scenarios, and the experimental results are used to evaluate the system.

Environmental configuration
This paper mainly studies the three-dimensional positioning of RFID tags in indoor scenes. The system is mainly composed of an Impinj RFID reader and three corresponding matched circular polarization antennas Laird s9028pcr. As shown in Fig. 5, we conduct our experiments in three different scenarios: (a) outdoor case, (b) indoor case with weak multipath effect and (c) indoor case with strong multipath effect. For each experiment, we utilize three antennas and one tag. The tag is positioned on a test stand with rotation functionality. Three antennas are positioned as one on the top and two at the bottom. The bottom two antennas are horizontally placed with an angle of 90 degrees. And the entire antenna setup is facing toward the tag. We locate 10 kinds of RFID tags in the space of 2 × 2 × 2 M 3 . The reader Impinj r420 follows EPC global C1G2 protocol and operates at 920.625-924.375 mhz with 16 channels in total. The bandwidth of the each channel is 231.25 kHz, and the transmission power of the directional antenna is set to 30 dBm by default.

Impact of different environments
In the practical application of the RFID system, the real scene must be far from the laboratory environment. In the actual environment, the system will face a variety of environmental conditions, but the difference in environmental conditions is actually attributed to the difference of multipath. Therefore, in order to evaluate the positioning effect under different multipath conditions, this paper selects three representative environments with typical real application conditions as the test, namely outdoor open environment, indoor weak multipath environment and indoor strong multipath environment. In this paper, the system is arranged in these three environments, respectively, and the control variable method is adopted to only change the environment to ensure that all other parameters in the system except distance are exactly the same. The aln-9640 tag carrying m730 chip was used as the positioning target in the experiment ×2× Select several positions within 1m for positioning test, collect data for 3 times at each position, and calculate the distance error of each collected data, respectively, �d(i) Finally, the distance error of each position is recorded as: � = 3 i=1 �d(i)/3. As shown in Fig. 6, it can be seen from experience that the positioning errors are often different in different environments. The more complex the environment is, the richer the multipath information is, and the greater the interference to the RF signal is. The same is the experimental result. In the outdoor open environment, only the tag itself and the ground can reflect the RF signal, and the multipath information is the least, so the positioning error is the smallest, with an error of 6.79 cm. The indoor weak multipath environment contains labels, ground, walls and a small number of other objects, which has more multipath information than the outdoor open environment. Therefore, the positioning error is higher than the outdoor open environment, with an error of 8.90 cm. In the indoor strong multipath environment, many reflectors are added on the basis of the indoor weak multipath environment to simulate the warehouse placement in the real scene. A large number of items are often stacked in the real warehouse, and the number of reflectors is countless, which often has strong multipath information. Therefore, the error is the largest among the three environments, with an error of 9.90 cm. The errors in the three environments are less than half the wavelength of the RF signal. Within the expected positioning error range, it is proved that the system is feasible in a variety of environments.

Impact of dynamic multipath
In the real scene, there cannot be only static multipath information in the environment near the label to be located, and the uncertainty and unpredictability of dynamic multipath information are often the biggest factors causing errors.
In order to simulate the real environment, in this experiment, a person walks around the label to be tested irregularly to simulate the dynamic multipath information in the real environment, and tests are carried out in the three environments in the previous section. The results are shown in Fig. 7. Under the three environments, the positioning error with dynamic multipath information in the environment to be tested is greater than that without dynamic multipath information, in which in the outdoor open environment, the positioning error with dynamic multipath information is 8.39 cm, which is 23.41% higher than that without dynamic multipath information. The positioning error with dynamic multipath information in an indoor weak multipath environment is 12.28 cm, which is 46.72% higher than that without dynamic multipath information. The positioning error with dynamic multipath information in indoor strong multipath Fig. 6 Positioning error under different environments environment is 12.86 cm, which is 29.22% higher than that without dynamic multipath information.
The above experimental results show that in an open environment, even if there is dynamic multipath information, the impact of dynamic multipath on the positioning will not increase significantly due to the small number of reflectors. In the more complex indoor environment, due to the large number of reflectors, once the dynamic multipath information is added, the effect of affecting the positioning accuracy will increase significantly. However, when the number of reflectors reaches a certain extent, a small amount of dynamic multipath will affect the positioning accuracy, but the impact will not be particularly large.

Impacts of different distances
In order to verify the influence of the positioning distance on the algorithm and ensure a good positioning effect in the measurement range, we conducted experiments using the control variable method.
This experiment only changes the distance between the tag and the antenna under the experimental conditions of the fixed experimental environment and tag type. And the distance is distributed between 20 cm and 200 cm. In the experiment, every 20 cm is divided as an interval, 20 samples are measured in each interval, and the positioning error of each sample is calculated. As shown in Fig. 8, the positioning error at all distances is below half wavelength, the minimum average error is 6.96 cm, and the maximum average error is 12.41 cm. With the increase in the distance, the positioning error has a certain positive correlation, but more than 70% of the positioning error is concentrated in 8-12 cm, and the maximum positioning error is only 16.15 cm. Within the range of expected positioning error, it shows that the work has high positioning accuracy within the existing agreed range.

Comparison with existing algorithms
In order to illustrate the availability and reliability of our solution, as shown in Table 1, we compare some existing works, including 3DinSAR, 3DTagArray, Phase Relock, RF-MVO and 3DLRA. These works are only tested in a single environment. In order to reduce the comparison error, this paper selects the indoor weak multipath environment scene as the work comparison. Among them, 3DinSAR uses the principle of interferometric synthetic aperture radar (InSAR) to simulate a multiantenna array by moving antenna in the plane to achieve the purpose of positioning. Compared with the method in this paper, 3DinSAR has the problems of large error and inconvenient deployment. 3DTagArray adds energy collection circuit and motion sensor on the commercial label and binds 2-4 label arrays on the object to locate the target. Although it has good positioning effect, the transformation of the label undoubtedly increases the cost of deployment and the difficulty of popularization. Phase Relock binds a single antenna to the robot and simulates the antenna array through the movement of the robot. Compared with 3DinSAR algorithm, it has lower computational complexity and higher real-time performance, but at the same time, the error also increases significantly, which is not a suitable positioning choice. RF-MVO combines RF signal positioning with monocular camera. By installing monocular camera on the mobile antenna, the phase and received signal strength of the tag are  measured through the antenna in the process of moving, and the visual information is captured by the monocular camera to improve the positioning accuracy. Although it can have high positioning accuracy, the positioning cost is also the highest. 3DLRA uses deep learning to input the phase, received signal strength and time stamp as features into the convolutional neural network (CNN) for positioning. Although it has high positioning accuracy, the method based on deep learning often requires a long training time and is strongly dependent on the environment. Once the environment is changed, it needs to be retrained, and the portability is poor. The three-dimensional positioning algorithm proposed in this paper has a good positioning effect in various environments and does not need additional hardware. All system deployment is completely carried out on the basis of existing commercial equipment and has good environmental adaptability and portability. However, we should point out that our solution can achieve 3D localization for the stable object. Nevertheless, due to the Doppler frequency shift, the proposed solution cannot track moving objects.

Conclusion
In this paper, we propose a 3D localization system for RFID. The main work and innovations of this paper are summarized as follows: (1) A phase model based on the polarization effect is proposed to make up for the influence of the phase error caused by ignoring the label angle in the existing phase model, and the phase accuracy is improved by phase unwrapping and phase filtering based on discrete wavelet transform. (2) By analyzing the vector relationship between the received phase, the pure phase and the multipath variable, a positioning algorithm based on the pure phase extraction is proposed, which removes the multipath variable from the received phase, solves the problem that the error of the traditional positioning method increases significantly in the multipath environment, and uses the HMFCW algorithm to unwrap the phase and recover the real phase information. (3) There is no need to refer to labels or anchors, nor to collect training data, which greatly enhances the environmental portability of the method. Using commercial RFID equipment, a three-dimensional positioning algorithm with easy deployment, high precision and low cost is realized, and the average positioning error is about 9 cm.

Abbreviations
The technical terms we used in this manuscript are list in the glossary, as shown in Table 2.
Author's contributions W.S. proposed the basic idea and analyzed the polarized features of the phase. K.Z. and Z.C. leveraged HMFCW to unwrap the polarized phase. Z.G. analyzed and verified that the OLS was appreciated to fit the phase data. Y.Q., P.W. and et al. performed the experiments. All authors read and approved the final manuscript.

Funding
This work is supported by the Natural Science Foundation of China under Grants. 61802299, 62002284, and the Fundamental Research Funds for the Central Universities xzy012020019.