Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s mostcited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
Department of Electronic and Communications Systems,Rzeszów University of Technology, Poland
Włodzimierz Kalita
Department of Electronic and Communications Systems,Rzeszów University of Technology, Poland
*Address all correspondence to:
1. Introduction
Current problems that occur in the field of anticollision Radio Frequency IDentification (RFID) prototype systems are solved in experimental way (De Blasi et al., 2010; Lehto et al., 2009; Polivka et al., 2009; Brown, 2007; Clarke et al., 2006; Penttilä et al., 2006; Jones & Chung, 2007). The low efficiency coefficient of identification for the multiple objects localized in the space Ω_{
ID
} doesn't allow to realize practical projects, such as, the identification of Fast Moving Consumer Goods (FMCG) – Fig. 1. In the light of nascent and modified legal communications standards, like for example, Electronic Product Code (EPC) in the area of UHF and HF ISO 180006, ISO 15693, ISO 180003 normalizations, there is a necessity to continue complex theoretical research and experimental investigations in the range of simultaneous analysis of EM field, communication protocols, and electric aspects of operating conditions of efficiency identification in anticollision RFID systems.
To generalise, the operation of passive anticollision inductive (LF, HF), and also propagation (UHF) coupling RFID system is characterized by the interrogation zone (IZ) which is estimated in any direction of 3D space for a group of electronic tags. The elements of algorithm of identification of interrogation zone for anticollision RFID system with the consideration of the energetic (i.e. field and electrical) and communication aspects of operation conditions have been presented in the following chapter. For calculations of the interrogation zone the algorithm based on Monte Carlo (MC) method and a computer program with the use of Mathcad 14 (called JankoRFIDmc’IZ) has been utilized.
2. Determining the interrogation zone using MC method
Unequivocal estimation of the interrogation zone for anticollision RFID system depends on automatic identification process. In accordance with the conditions of the correct operation of any RFID system, different locations of many tags strongly change the functioning of an antenna unit array: read/write device (RWD) and individual tags. The problem of determining the interrogation zone is related to two cases. In the first of them, an assumption is made that the location of the ntags group is determined, whereas in the second case, all possible locations of the group of ntags in a space around the RWD antenna are going to be analyzed. The problem connected with the first case is realizable by the assumption that the process of determining the interrogation zone in RFID system will be carried out in a few feedback cycles which allow to find the proper location of tags. The statement “few feedback cycles” is related to the time which is accepted for determining the interrogation zone for all ntags. The mentioned feedback cycles in the carried out simulation include a modification of the tags location which don't fulfil conditions of the correct RFID system operation. The problem from the second case is almost impossible to solve because the prolonged process of calculations would be ineffective. Seemingly, in that case, a method of "trial and error" during the search of the interrogation zone of RFID system might be easier to apply, however, the presented MC method is a wellfounded alternative.
The presented premises lean towards the necessity of application of the techniques which make use of random numbers (Kalos & Whitlock, 2008). The result of this is the solution of the problem of the ntags group location, and testing the functional efficiency of the antenna unit array: read/write devicetags, that is an estimation of anticollision RFID system interrogation zone for given efficiency of identification η_{
ID
}. The percentage of identification efficiency is given by the equation:
ηID=lIDOKn⋅100%E1
where l_{
IDOK
} is the number of tags for which the desired read/write operations have been properly done.
The problem contained in the MC method has a probabilistic nature, and it's solution is obtained by simulation of the given object (Rubinstein & Kroese, 2007). The simulation object is represented by the antenna unit array: RWDtags with the consideration of a synthesis of this antenna unit array and according to all equations which are going to be determined during the synthesis of its electric model in an anticollision RFID system.
For a laboratory process of automatic objects identification the solution of the problem consists in finding the interrogation zone of given RFID system, with its shape, location and
orientation in 3D space assumed. In the conducted research it was assumed that the demanded area should be square shaped and situated at the z_{
ID
} height, whereas it's location should be axially  symmetrical and parallel to RWD antenna (Fig. 2). Such an assumption results from the orientation of tags that are parallel to symmetrical RWD antenna which has, for example, circular or square shape in inductive coupling RFID systems.
A random layout of ntags at P_{
i
} points of Cartesian space at (x_{
i
},y_{
i
},z_{
ID
}) has been assumed in  considered in a sequence  k steps during the search of the RFID system interrogation zone. The random variables x_{
i
} and y_{
i
}, for i=1..n obtain various values which can’t be predicted, but for which the definite distribution is assumed. The electromagnetic field in any point of communication space is heterogeneous. This effect becomes the clearer, the nearer to the surface of RWD antenna, and the farther from its centre a point is situated. This knowledge allows to make a uniform (rectangular) distribution in intervals: x_{
Ak
},x_{
Bk
}) for the random variable x_{
i
}, and y_{
Ak
},y_{
Bk
} for the random variable y_{
i
}, in kstep for the analyzed area. For uniform distribution of the random variables x_{
i
} and y_{
i
}, and for the definite values of x and y, the distribution functions are given by: (xx_{
Ak
})/(x_{
Bk
}x_{
Ak
}) and (yy_{
Ak
})/(y_{
Bk
}y_{
Ak
}). It should be noticed that the random variables x_{
i
} and y_{
i
} are mutually independent. This means that the random variables x_{
i
} and y_{
i
} are stochastically independent, since the distribution of the x_{
i
} does not depend on the value y_{
i
} and vice versa. In this case, the probability density of a pair of random variables (x_{
i
},y_{
i
}) is equal to the product of the probability density (x_{
i
}) and (y_{
i
}) independently.
In order to determine that the RFID system is functioning correctly for given tags locations it is not enough to achieve the efficiency of identification η_{
ID
}=100% for ntags and fulfill all conditions for a correct operation of anticollision RFID system. It cannot be predicted whether for k area in which all the conditions mentioned above are fulfilled, the coordinates sampling of tags locations on the surface of their arrangement, allows to fulfill the border case of a correct operation of the whole RFID system. In kstep for the analyzed area in which all the conditions of a correct operation of anticollision RFID system for given efficiency of identification are fulfilled, the practical use of the law of large numbers (Kalos & Whitlock, 2008) is the solution to this problem. For the random variables x_{
i
} and y_{
i
} independently, the strong law of large numbers for the analyzed case is given by:
where p denotes the expected value of the random variables x_{
i
} and y_{
i
} (which are equal to zero because the interrogation zone is axially  symmetrical and parallel to RWD antenna), and PP denotes the probability of sampling of variables for m approaching to infinity, but m denotes the number of multiple sampling of tags location (i.e. random variables x_{
i
} and y_{
i
}) for k analyzed area.
What follows from the equations (2) and (3) is that the sequences of random variables S_{
m
}(x_{
i
}) and S_{
m
}(y_{
i
}) converge with probability “1” to the expected value p=0 of the random variables x_{
i
} and y_{
i
}. It can be found that the mtuple increase of the number of the random variables x_{
i
} and y_{
i
} sampling in kstep for the analyzed area lengthens the calculation process during the simulation of an antenna unit array. In accordance with the law of large numbers, the probability of a correct estimation of the interrogation zone for RFID system increases. First of all, this is connected with the examination of a larger number of localized ntags cases. If the conditions of the correct operation of anticollision RFID system are not fulfilled in any of m multiple sampling of tags location for k analyzed area, then the next process of multiple sampling should be stopped, and it becomes necessary to examine the next (k+1)  smaller area of tags location in the xy plane. The MC solution for the analyzed object completes a procedure which confirms the fulfillment of all conditions for the correct operation of anticollision RFID system. The procedure is correct for the given efficiency of identification, and for the area in which all the m multiple sampling of tags location lead to a positive calculation result of the antenna unit array: read/write devicetags.
Correct selection of the m number, which will be satisfactory under the experiment, as well as adequate to calculation time and probability of possible tags locations, is a problem. From equations (2) i (3) which describe the strong law of large numbers, the dependence of probability PP for the random variables x_{
i
} and y_{
i
} (which are stochastically independent, and which have a uniform distribution) in function of the m n has been presented in Fig. 3.
Assuming that the probability PP exceeds the value 0.95 independently for the random variables x_{
i
} and y_{
i
}, the value m n=250 was determined during the calculation of interrogation zone in automatic identification process. These parameters were determined for 10^{6} sampling of 250 independent random variables x_{
i
} and y_{
i
} which have a uniform distribution. For every sampling, the minimum value of probability PP has been searched. The determined value of m n=250 is compatible with a central limit theorem which states that the sum of a sufficiently large number of identically distributed independent random variables, each with finite mean and variance, is going to be approximately normally distributed (Rice, 2006). Uniform distribution of the random variables x_{
i
} and y_{
i
} is in fact different from a normal distribution, but  for this determined value of m n=250  random variables x_{
i
} and y_{
i
} are convergent to a normal distribution. In this case, the obtained compatibility with a central limit theorem confirms the correctness of product m n=250.
The presented idea of ntags analysis at a specifically determined z_{
ID
} height, results from a practical demand for realization of automatic identification process with the anticollision RFID systems. The identification of single products which are located inside a container on a pallet can be the practical example of this process. Identification of single objects separately is impossible in this situation, but their location on a pallet is mostly scheduled  because logistic system has to work satisfactorily (Mo & Lorchirachoonkul, 2010; Shaoping Lu et al., 2007; March, 2005; Jones & Chung, 2007). The development of the presented MC solution on an area Ω_{
ID
} in the xyz space requires investigation of every jsurface independently where tags will be located on the whole area Ω_{
ID
} at points P_{
ij
}(x_{
i
},y_{
i
},z_{
IDj
})  (Fig. 2). If all ntags are in a disordered state in the space Ω_{
ID
}, then the stochastic independence of all the coordinates x_{
i
}, y_{
i
}, z_{
i
} should be assumed. However, this idea is very complicated because it is not enough to assume that all the pairs of random variables are independent. Taking into consideration the practical requirements of different automatic identification processes, the example presented above is marginal, yet very interesting from a scientific point of view.
3. Conditions of correct operation of anticollision RFID system with inductive coupling
Passive RFID systems with inductive coupling are widespread (ID World, 2009; Wolfram et al., 2008; Jones & Chung, 2007; Paret, 2005). These systems can operate in individual and anticollision regime (Finkenzeller, 2003), and the need to design such systems appears more often nowadays. Functioning of RFID systems with inductive coupling is based on the use of energy which is stored in a magnetic field (Chen & Thomas, 2001; Rautio, 2003; Troyke & Edgington, 2000). The kind of executed operation in individual phases of the exchange of data between units of a system is essential during communication in anticollision identification process (JankowskiMihułowicz et al., 2008). The basic condition of effective operation of the system is the proper supply of each tag in a heterogeneous magnetic field created by RWD antenna loop. A minimal value of energy (necessary for proper read/write operation of the tag) is determined by minimal value of magnetic induction B_{
min
} in each point P(x,y,z) of its location (Fig. 4).
Analysis of the general RFID system schema allows for the determination of the complete impedance of RWD antenna Z_{
R
}, taking into account an influence of all coupled tags. Maximal change of the impedance Z_{
R
} under the influence of the tags is expressed by the maximum value of difference in impedance arguments Δφ_{
Rmax
} without the tags and with them, respectively. The value Δφ_{
Rmax
} is limited to assure the correct operation of the system.
A determination of communication conditions for proper operation of RWDtags antenna set is also possible on the basis of the schema analysis, taking into consideration the properties of data transmission process (transmitted frequency band, data flowability and required time relations in selected communication protocol).
The quality factor Q is a measure of tag antenna unit functioning efficiency in areas of energy transfer and communication conditions in RFID systems with inductive coupling (Newman et al., 1975; Redinger et al., 2003). These conditions are expressed by maximal values of quality factors of RWD and tag antennas: Q_{
Rmax
} and Q_{
Tmax
}, respectively. In this case, it is necessary to observe that selection of value Q_{
Tmax
} is compromised by energy (utilization of magnetic field energy) and communication requirements (JankowskiMihułowicz & Kalita, 2009).
During the synthesis of tag antenna, the changes of this parameter can be made by inductance change (indirectly – by changing the effective resistance of antenna loop), adjusted to requirements of the shape and geometrical sizes of an electronic tag. The proper synthesis of interrogation zone of RFID system is closely connected with three aspects which concern the maximum value of Q factor for the operating tag. The first of these aspects is related to the correct operating of the tag supply system, that is, the possibility of radio communication appearing. The second aspect concerns the necessity to obtain the required data transmission bitrate (and also the bandwidth) in direction: tagread/write device. The third aspect concerns the impulse and step response of tag circuit in case of reverse data transmission, to provide a correct identification of commands sent from the RWD. The last two aspects should result directly or indirectly from the electronic tag chip specification for which the antenna is going to be projected. The first aspect should be considered at the stage of antenna synthesis, and the value of Q factor should contain all of the mentioned limitations of operating passive tag. It is essential to ensure the homogeneous proper interrogation zone of RFID system that is a mutual overlay of zones that result from conditions of tag supply (by absorbing the energy of magnetic field) and the radio communication carried in the system (realized with the suitable value of signaltonoise).
Paying attention to the maximum work distance between elements of the RFID system, in particular for systems working in the RFID far field, it is necessary to estimate the simulated and built antenna set RWDtags in relation to the obligatory normalizations of communication and EMC (ETSI EN 300 330, 2010; JankowskiMihułowicz, 2010).
4. Energy transfer in passive anticollision RFID system with inductive coupling – fundamental equations
4.1. EM field aspects
Analysis of the read/write device antenna unit allows to make an assumption that the antenna loop current (I_{
R
}) is constant along the whole flow way (Fig. 5).
Change of the electric charge density in time equals zero, so in that case the electric current density divergence equals zero as well. Making these assumptions permits to apply the magnetostatic laws to magnetic field analysis for any RWD shape. In accordance with vector BiotSavart law, the magnetic induction value of B in any space point P(x,y,z) is given by the equation:
B=μ0IRNR4π∮dl×rb1rb13E4
where: μ_{0}=4 π 10^{7} H/m.
Application of the BiotSavart law for the RWD’s antenna loops is possible by the additional assumptions that the wire diameter of RWD’s antenna is negligible in relation to the geometrical loop sizes, and also that there is full inductive coupling between the individual loop turns (N_{
R
}).
For the circleshaped loop (Fig. 5a), the axial symmetry permits convenient change from Cartesian to cylindrical coordinates. The vector describing the dl location at P point, in which the value of magnetic induction is calculated, is given by the formula:
rb1=rb−r1E5
where the vector describing dl element location that changes in φ angle function, and the vector describing location of point P, are given as follows:
r1=(rR⋅cos(φ)rR⋅sin(φ)0)E6
rb=(xyz)E7
The unit vector connected with dl=rR⋅dφ is given by the formula:
uφ=(−sin(φ)cos(φ)0)E8
Changing the coordinate system leads to final equation, which describes the magnetic vector at any space location with (x,y,z) coordinates for circleshaped RWD loop:
B=μ0IRNR4π∫02πrR⋅uφ×rb1rb13dφE9
In the case of RWD antenna loop, constructed as polygon (Fig. 5b, c), BiotSavart law with principle of superposition permits to add at location P vectors, that descend from individual antenna parts. In this case, the total magnetic induction is calculated from the equation:
B=∑iBiE10
where i denotes analysed side for RWD antenna loop, constructed as polygon.
The obtained vector equations permit numerical calculating of the value of magnetic induction separately for individual components in directions x, y and z (B_{
x
}, B_{
y
}, B_{
z
}). The components of magnetic induction B in any space point P(x,y,z) are given for a circular loop shape (Fig. 5a) by the following equations:
Many practical solutions for identification are characterized by the parallel location of tag antenna and RWD loop (Fig. 6).
This location of individual tags allows the magnetic induction B_{
min
} to reach its minimum value only in relation to zmagnetic induction components (vectors B_{z1}÷B_{zn} for 1 to n tags). This case applies to places in which tags working in anticollision process has been located (P_{
1
}(x_{
1
},y_{
1
},z_{
1
})÷P_{
n
}(x_{
n
},y_{
n
},z_{
n
})).
The presented approach creates numerous limitations connected with the decrease of the interrogation zone in RFID system. This results from too low value of the perpendicular magnetic induction component in relation to tag antenna loop plane. The efficient use of communication space in which anticollision process is going to take place, and also specification of the object marked by passive RFID tag, requires consideration of any tag orientation with regard to the individual components of magnetic induction vector.
The issue of any tag orientation in three dimensions xyz comes down to the tag deviation of α and β angles from parallel location of RWDtag antenna loops (Fig. 7a). In accordance with presented model (Fig. 7b and Fig. 7c), deviation of α angle occurs in zx plane, however deviation of β angle occurs in α y plane. Calculating the value of perpendicular magnetic induction component for tag, which is deviated of α and β angles (B_{
αβ
}) has been divided in two parts. By application of the superposition theorem in the first part, after tag deviation of α angle, the perpendicular magnetic induction component is given by:
Bxzα=Bxα+BzαE17
where the values of vector components are given by:
Bxα=Bx⋅sin(α)E18
Bzα=Bz⋅cos(α)E19
By application of the superposition theorem in the second part again, after tag deviation of β angle, the perpendicular magnetic induction component is given as follows:
Bαβ=Byβ+BxzαβE20
where the values of vector components are given by:
Byβ=By⋅sin(β)E21
Bzα=Bz⋅cos(α)E22
From the equations (1722) comes, that the perpendicular magnetic induction component for passive tag, which is deviated of α and β angles, is given by:
Knowing the magnetic induction separately for individual components in directions x, y and z (B_{
x
}, B_{
y
}, B_{
z
}), the obtained equation (23) permits calculation of the perpendicular magnetic induction component. The aforementioned necessity of changing the tag orientation should be carried out for assurance of correct tag work in the individual space point P(x,y,z). In this way it is possible to calculate the system interrogation zone, which is forced by the specification of identified object, that is the necessity of individual tag location on marked object.
4.2. Electric aspects
The second essential stage of energy transfer from RWD to tags includes an elaboration of electrical model of the whole anticollision RFID system with inductive coupling for the full frequency range (LF – typically 125 kHz, 138 kHz and HF – typically 13.56 MHz). The basic part of the system is circuitry of RWD  tag antenna set in which an adequate operation states should be taken into consideration. The idea of electrical model of analyzed RFID system is based on the circuitry of RWD  tag antenna units (Fig. 8).
The correct operation of tag internal integrated circuit in communication process phases is conditioned by the voltage U_{
T
} which is inducted on the trimming tag antenna. Functioning of the memory and microprocessor of passive tag circuit is connected with the occurrence of the stabilization voltage U_{
TCS
} inside its chip (Villard, 2002; Friedman et al., 1997).
Taking into consideration the fact that the stabilizer is on input to the part of the tag chip and it influences tag’s antenna unit, it is necessary to make the synthesis of a module which was divided into two parts: rectifier and voltage regulator. Rectification of the voltage induced in the tag antenna loop, takes place in halfwave and fullwave rectifier (Jamali, 2006), however, it is controlled by a voltage regulator (Friedman et al., 1997). Functioning of the rectifier and voltage regulator involves the use of other passive elements to the tag's schematic diagram (modulator  index TCM, rectifier  TCR, voltage regulator  TCS). The voltage regulator (represented by the R_{
TCS
} resistance) has the most significant influence on the parameters of the whole electric circuit which is located on the side of connectors of the tag's antenna loop.
The inductive coupling between the antenna loops is expressed by mutual inductance M determined in relation to field conditions at the given system efficiency for arbitrary location and orientation of the tags:
where N_{
T
} is number of tag antenna turns of wire.
In tag antenna loop is induced the voltage U_{
T
} (as effect of the coupling):
UT=jω⋅M⋅IR1+(1(RTCR+RTCS)‖RTCM+jωCT)⋅(jωLT+RT)E26
where the symbol  denotes the parallel connection of resistances representing the generalized circuit of rectifier and voltage regulator (R_{
TCR
}, R_{
TCS
}) and the modulator (R_{
TCM
}) inside the chip structure.
Assuming that the module of U_{
TCS
} is constant, the variable resistance of voltage regulator is given by (27). This equation is true when U_{
T
}>U_{
TCS
} :
The minimum value of U_{
T
} (U_{
Tmin
}) is the base for determining the interrogation zone for the anticollision system. The value of voltage U_{
Tmin
} clearly impacts the analytical relation, which allows for expression of the minimum value of magnetic induction B_{
min
} – the elementary parameter defining the interrogation zone:
The induction B_{
min
} is differentiated on the basis of value U_{
Tmin
}, for the direction of data transmission and the kind of operations in internal tag memory. For experimental verification of calculated value B_{
min
} the special laboratory stand has been made. The stand allows to measure the maximum distance between the RWD and tag antennas, for which the correct operation of the RFID system is ensured.
Operating of the anticollision RFID system with inductive coupling means that working tags, which are located in the system interrogation zone, are going to have an influence on the parameters of RWD antenna unit. The synthesis of the electrical model of an anticollision RFID system includes the replacement of all tag electric circuits by the impedance Z_{
TR
} (Fig. 8). Those circuits are coupled by the mutual inductance M. Location of a specific number of ntags in the interrogation zone of RFID system means that the working tags’ quantitative participation in the influence on the RWD is going to be represented by the sum of all impedances Z_{
TR
} (from Z_{
TR1
} to Z_{
TRn
}). Such an influence of impedances Z_{
TR
} on the read/write device leads to the adverse effect of detuning the RWD antenna. Additionally, the value of antenna loop current I_{
R
} decreases in a read/write devices without output current stabilization. In a read/write device with current stabilization, the antenna voltage U_{
R
} increases to the limit, and after an overflow, the value of antenna loop current I_{
R
} also decreases as in the previous case. The consequence of the influence of operating tags on the RWD in each case is the reduction of the system interrogation zone which is caused by the energetic conditions of devices’ operation. Changes are the stronger the larger is the value of the mutual inductance between tags antenna loop and RWD.
This problem can be illustrated by the total impedance of RWD antenna Z_{
R
}, taking into account the influence of all coupled tags. For an equivalent schematic diagram from Fig. 8, this impedance is given by the following equation:
An influence of working tags causes the RWD antenna unit to detune of the value Δf_{
R
} from the frequency f_{
0
}. This detuning is revealed by the change of impedance’s argument Δφ_{
R
} for the working frequency of RFID system. The border values of these parameters determine the last directive to the efficiency of antenna unit array research for given read/write device. There is a necessity to maintain Δφ_{
R
} within the limits in order to provide correct functioning of the RWD and of the whole anticollision RFID system. Such an assumption is practically used in external antenna tuning devices dedicated to tuning the longrange antennas working in LF and HF inductive coupling RFID systems (Feig, 2005, 2006; Philips, 1996; Texas Instruments, 2002).
The experimental research have been carried out for different RFID elements using the laboratory system which allows to determine single and anticollision identification process for all frequencies in RFID systems with inductive coupling (Fig. 9).
During the search of the interrogation zone of RFID system for a given efficiency of identification η_{
ID
} (1), the appearance of condition Δφ_{
R
}>Δφ_{
Rmax
} makes a correct identification impossible. In MC calculation the parameter Δφ_{
R
} is calculated on the basis of the total impedance Z_{
R
} of RWD antennas arrangement (29), taking into consideration influence of functioning tags on this antenna  Z_{
TR1...n
} calculated from the equation (30). For example, the limit phase value was Δφ_{
Rmax
}=±15º for the Philips HITAG RM 800 read/write device, working at frequency f_{
0
}=125 kHz. The technical documentation available on the basis of an agreement reached between the Department and the Philips Semiconductors has been used in the investigations.
For the correct energy transfer in the anticollision RFID system with inductive coupling, assuming the possibility of using identical tags in automatic identification process, the specified value of minimum magnetic induction B_{
min
} will be the parameter that limits tags’ correct operation area. This parameter should be calculated from the equation (28) for the individual single tag construction and the kind of operations executed in its internal memory (read/write of tags memory). If the value of perpendicular component of magnetic induction vector at point of the location of tag is smaller than his parameter B_{min}, then the correct functioning of this tag in anticollision system is impossible. This denotes that the tag is in the area where communication with RWD is impossible, and efficiency of identification η_{
ID
} is lowered.
Process of determining the interrogation zone using MC method has been preceded by measurement and calculation of B_{
min
} conducted during the process of reading information from the internal memory of tag. The results of these measurements and calculations were presented in the table 1.
In the simulation and measuring part of the experiment respectively, the calculated and measured values B_{min} were the minimum limit of the correct operation of a single tag located in the area of field conditions of functioning of the whole RFID system. In both parts of the experiment locations (in points P_{
i
} of cartesian space at (x_{
i
},y_{
i
},z_{
ID
}) coordinates) of ten tags of a chosen type were selected randomly 25 times (from chapter 2: n∙m=250).
Tag
Measuded z_{max} 1)
Measured Bmin 2)
Calculated Bmin 3)

m
μT
μT
HITAG 1 ISO CARD
0.52
0.74
0.74
HITAG 1 WORLD TAG 50
0.44
1.16
1.16
HITAG 1 WORLD TAG 30
0.27
3.72
3.73
HITAG 1 WORLD TAG 20
0.22
5.63
5.62
1) The measurement of the maximum working distance zmax from the center on axis of symmetry of RWD antenna loop for square read/write device antenna (where a=0.3 m, NR=32, IR=0.213 A) – this is the result of the positive identification of the tag serial number. 2) The measurement by means of analyser Advantest R3132 and Rohde & Schwarz HZ14 near field probe (Rohde & Schwarz, 2003). 3) The values calculated in the JankoRFIDmc’IZ v. 4.08 application (JankowskiMihułowicz, 2007)  on basis of electrical model  equation (28).
Table 1.
Measured and calculated values B_{min} for tags selected to investigations
The example results of the calculated and measured interrogation zone (Fig. 10), were placed on the plane at (x, y, z_{
ID
}) coordinates. The measured interrogation zone is the result of the positive identification of all n=10 tags serial numbers, during conducted experiment, all m=25 multiple sampling of their location. For every multiple sampling of the location of tags in measuring chamber, spatial measurements of z component of magnetic induction B vector were made. On the basis of (Rohde & Schwarz, 2003), the measurement of the component of the vector B perpendicular to the area of the antenna loops of tags was conducted in the 625 points (the resolution of 2 cm on 0.5 m x 0.5 m xy surface – the movable platform in the measuring chamber – Fig. 9b).
All of the calculations and measurements were performed for square antenna of the RWD unit which was tuned in the measuring chamber without the influence of tags, and the achieved value was Δφ_{
R
}=2.5º. In all studied cases, the border value of Δφ_{
Rmax}, wasn't crossed. Thanks to this, the efficiency of identification for the height z_{
ID
} was 100 % in the area of fulfillment of the condition of the magnetic induction minimum value. Difference between the calculated and measured interrogation zone (in the worst case, for the smallest heights z_{
ID
}, on the level ±1.5 cm), is caused mainly by applying an approximate geometrical model of the antenna loop of the RWD. These differences are caused by the fact that the RWD antenna loop was build as loose turns of wire, and that was assumed during synthesis of the geometrical model of the RWD antenna loop.
The measurements in the RWD  tags antennas arrangement required applying many direct and indirect measuring methods. The obtained results always contained certain dispersion of the values, which can always be  in a justified way  ascribed to measured sizes. The multiple results were obtained from many measuring sets.
Generally, the problem of the uncertainty of determining the interrogation zone of the anticollision RFID system with the inductive coupling, has two aspects: simulations and measures. In the process of evaluation of the uncertainty of determining the interrogation zone in the measuring part of the experiment, essential factors are uncertainties of the magnetic induction components u(B) measurements:
u(B)=(∂B∂H⋅u(H))2+(∂B∂μ⋅u(μ))2E31
where:
u(H)=(∂H∂V0⋅u(V0))2+(∂H∂AF⋅u(AF))2E32
where u(V_{
0
})  standard uncertainty of voltage measured by means of Advantest R3132 spectrum analyzer and the R&S HZ14 near magnetic field probe.
This uncertainty includes the systematic influences which cannot be removed during the conducted experiment. They are represented by the set of coefficients read from prepared tables and graphs in the Advantest R3132 spectrum analyzer user manual. u(AF) denotes the uncertainty of antenna coefficient read for measuring frequency (f_{
0
}). For the spatial, multipoint measurements which were made in the measuring chamber of the investigative set, the standard relative uncertainty for the magnetic induction u_{
%
}(B) was on the level 12 %.
In the process of evaluation of uncertainty of the interrogation zone estimation in the simulating part of the experiment, the component factors of the complex uncertainty of the entrance data measurements and output data calculations were considered. They were taken into account in the process of estimating the efficiency of the system antennas arrangement with the MC method, which is made by the JankoRFIDmc’IZ application (JankowskiMihułowicz, 2007).
Explaining this problem, function f which represents the interrogation zone exhibits significant nonlinearity. Therefore, regarding the error propagation, the higher terms in the Taylor's expansion should be taken into account. Their form is as follows:
In indirect measurements every size, calculated or measured directly, brings the different contribution to the uncertainty u(f). The determination of suitable weighting factors resulting from the uncertainty propagation law for the considerably nonlinear function f, according to the higher terms in the Taylor's expansion, is a complicated mathematical question. This is a complicated problem at the present stage of works.
The efficient leading of the automatic identification processes, such as: forwarding mail, materials, articles (in industry); identification of valuable minerals, samples for analysis (in science and medicine), requires the use of a modern radio methods of the simultaneous identification of many objects. The mentioned processes generally belong to the automatic identification group, in which RFID electronic tags are replacing, for example, barcodes. This is caused by the wellknown technical limitations of the objects identification methods used nowadays. The accessibility of electronic tags, the continuous reduction of their production costs and the standardization of work conditions of RFID technology, allows to make a decision about the implementation of quite a new method in the process of automatic identification.
The laboratory research and tests fully confirm the correctness and usefulness of the elaborated (in Department of Electronic and Communication Systems at Rzeszów University of Technology), method of synthesis of anticollision RFID system, where the essential component, based on Monte Carlo method, is the determination of interrogation zone for the system with suitably located tags. It should be noted that the synthesis procedure includes the simultaneous analysis of electromagnetic field, communication protocols and electric aspects of operation conditions in the process of system efficiency identification. Presented part of the problem of interrogation zone synthesis is the base for practical use of projected identification systems, required for specific anticollision RFID applications. The future investigations will be focused on the analysis of efficiency and interrogation zone of the anticollision RFID systems operated in dynamic conditions (speed changes of orientation of suitably located tags). Additionally, the extension of JankoRFIDmc’IZ program on a propagation coupling RFID system is planned. The elements of algorithm of interrogation zone identification for anticollision RFID system taking into consideration the energetic (i.e. field and electrical) and communicational aspects of operation conditions are going to be supplemented by elements of antennas and wave propagation in UHF.
This work was partly supported by the Project "Developing research infrastructure of Rzeszów University of Technology" within the Operational Program Development of Eastern Poland 20072013 of the Priority Axis I Modern Economics of Activity I.3 Supporting Innovation, Contract No. POPW.01.03.0018012/0900.
References
1.De BlasiM.MighaliV.PatronoL.StefanizziM. L.2010 Performance Evaluation of UHF RFID Tags in the Pharmaceutical Supply Chain, In:
2.BrownD.2007
3.ChenS. C. Q.ThomasV.2001 Optimization of inductive RFID technology,
4.ClarkeR. H.TwedeD.TazelaarJ. R.BoyerK. K.2006 Radio frequency identification (RFID) performance: the effect of tag orientation and package contents,
5.ETSI300 300 330 (2010). Short Range Devices (SRD); Technical characteristics and test methods for radio equipment in the frequency range 9 kHz to 25 MHz and inductive loop systems in the frequency range 9 kHz to 30 MHz
6.Feig Electronic GmbH2005 Dynamic Antenna Tuning Board ID ISC.DAT,
7.Feig Electronic GmbH2006 Static Antenna Tuning Controller for longrange antennas ID ISC.SAT.CA,
8.FinkenzellerK.2003
9.FriedmanD.HeinrichH.DuanD. W.1997 A LowPower CMOS Integrated Circuit for FieldPowered Radio Frequency Identification Tags,
10.ID World2009 The Players of the Auto ID Industry, Section RFID,
11.JamaliB.RanasingheD.ColeP.2006 Design and optimisation of power rectifiers for passive RFID systems in monolithic CMOS circuit, In:
12.JankowskiMihułowiczP.2010 Field Conditions of Interrogation Zone in Anticollision Radio Frequency Identification Systems with Inductive Coupling, In:
13.JankowskiMihułowiczP.KalitaW.2009 Efficiency of Tag Antenna Unit in Anticollision Radio Frequency Identification Systems with Inductive Coupling,
14.JankowskiMihułowiczP.KalitaW.PawłowiczB.2008 Problem of dynamic change of tags location in anticollision RFID systems,
15.JankowskiMihułowiczP.2007
16.JonesE. C.ChungC. A.2007
17.KalosM. H.WhitlockP. A.2008
18.LehtoA.NummelaJ.UkkonenL.SydänheimoL.KivikoskiM.2009 Passive UHF RFID in Paper Industry: Challenges, Benefits and the Application Environment.
19.MarchV.2005 Builtin RFID for electronic products,
20.MoJ. P. T.LorchirachoonkulW.2010 RFID Infrastructure for Large Scale Supply Chains Involving Small and Medium Enterprises. In:
21.NewmanE.BohleyP.WalterC.1975 Two Methods for the Measurement of Antenna Efficiency, IEEE Transactions on Antennas and Propagation, 234457461 , 00018926X
22.ParetD.2005
23.PenttiläK.KeskilammiM.SydänheimoL.KivikoskiM.2006 Radio frequency technology for automated manufacturing and logistics control. Part 2: RFID antenna utilisation in industrial applications,
25.PolivkaM.SvandaM.HudecP.ZvanovecS.2009 UHF RF Identification of People in Indoor and Open Areas,
26.Rohde & Schwarz2003 Probe Set for E and H NearField Measurements 9 kHz to 1 GHz HZ14, Operating Manual 1026.7744.03
27.RautioJ. B.2003 Electromagnetic Analysis Speeds RFID Design,
28.RedingerD.FarshchiR.SubramanianV.2003 An allprinted passive component technology for lowcost RFID,
29.RiceJ.2006
30.RubinsteinR. Y.KroeseD. P.2007
31.Shaoping Lu; Yaohua Wu & Yongtao Fu2007 Research and design on palletthroughout system based on RFID.
32.Texas Instruments2002 TI RFID Series 2000 Antenna Tuning Indicator RIACCATI2,
33.TroykeP. R.EdgingtonM.2000 Inductive Links and Drivers for RemotelyPowered Telemetry Systems,
34.VillardP.BourC.DallardE.LattardD.de PontcharraJ.RobertG.RouxS.2002 A lowvoltage mixedmode CMOS/SOI integrated circuit for 13.56 MHz RFID applications,
35.WolframG.GamplB.GabrielP.2008 The RFID Roadmap: The Next Steps for Europe, Springer, 97835407101899783540710189
Written By
Piotr JankowskiMihułowicz and Włodzimierz Kalita
Submitted: November 10th, 2010Published: July 20th, 2011