The disclosure relates to a method for modeling the structure of an item of currency. In particular, the disclosure relates to a method for predicting the behavior of a currency sensing system as related to the structure of a tested item of currency. The disclosure also relates to a sensing apparatus used for sensing characteristics of an item of currency.
Automated transaction machines (e.g. vending machines, gaming machines, ATMs, etc.) typically accept items of currency in exchange for goods and/or services. Items of currency are typically inserted into an automated transaction machine, and are evaluated by an authentication unit to determine if they are genuine or non-genuine. In some forms of currency (e.g. banknotes) there can be inks used for printing images and other features deemed necessary by a respective banking authority. It is known that some inks used for printing can exhibit electromagnetic properties such that a sensing system can be used to verify its presence or characteristics. Banknotes are sometimes constructed using multiple layers of different materials to form a substrate. In some cases one or more of these layers exhibit electromagnetic properties such that a sensing system can be used to verify its presence or characteristics.
Other items of currency (e.g. coins or tokens) can be constructed using at least one component or material that exhibits electromagnetic properties. Some currently circulating coins are constructed using more that one material (e.g. cladded coins, platted coins, or bi-color coins), and in some cases at least one of the materials used exhibit electromagnetic properties. In automated transaction machines, there can be provided a sensing unit that is capable of verifying the presence or characteristics of a given material in an item of currency. For the purposes of the disclosure the term “item of currency” includes, but is not limited to, banknotes, bills, coupons, security papers, checks, valuable documents, coins, tokens, and gaming chips.
The authentication of items of currency can also occur in processing equipment used by central banking institutions for sorting and evalaution. This equipment can include an authentication unit configured to sense at least one electromagnetic property of an item of currency for the purpose of recognition and/or authentication.
Various aspects of the invention are set forth in the claims.
In some implementations there can be provided a method for predicting the response from an item of currency when using a specified currency sensing system. There can be provided a mathematical model of an item of currency and a mathematical model of a given sensing system such that for a specified set of parameters of an item of currency, the response of the specified currency sensing system can be predicted.
In some implementations, there can be provided a method for determining a particular construction of an item of currency based on theoretical responses from such an item of currency being tested with a theoretical sensing system. In some implementations, there can be provided a method and system for determining the structure of an item of currency based on theoretical responses of such an item of currency being evaluated by a theoretical sensing system and further based on a set of known items of currency.
In some implementations there can be provided a method of predicting the response of a specified currency sensing system for a given inputted set of parameters of an item of currency. More particularly, an item of currency can be constructed using at least one component (e.g. material layer) exhibiting electromagnetic properties. In some implementations, there can be a mathematical model of an item of currency such that at least one component of an item of currency can be described relative to its respective electromagnetic properties. It is possible that for a specific item of currency there can be a plurality of components (e.g. 3 layers) exhibiting electromagnetic properties. With an item of currency having a plurality of layers, each layer can be inspected to determine the material thickness and type. In some implementations, the inductance relative to frequency can be used to characterize at least one electromagnetic component present in an item of currency. In some implementations, an item of currency can be characterize by a complex impedance measurement (or estimation) relative to frequency when being evaluated (i.e. sensed) by a Pulse Eddy Current (PEC) sensing system.
In some implementations, currency handling machine 10 includes a sensing device 25 including a PEC sensor 60. PEC sensor 60 can be arranged to include input 61, a coil 63, core 65, and output 68 as is commonly known in the arts. In some implementations input 61 can be configured to use broad band techniques for driving PEC sensor 60. In other implementations, input 61 can be configured to use other techniques (e.g. spread spectrum, frequency hopping) for driving PEC sensor 60.
In some implementations, the input 61 and output 68 of coil 63 can be used as inputs to a model (e.g. Equation (A)) to obtain electromagnetic properties of at least one material (i.e. component) of currency item 50. The material properties obtained from the model can then be used as inputs to a classification technique (e.g. Malahanobis Distance, Linear Discriminant Analysis, Feature Vector Selection) to obtain statistical information on item of currency 50 relative to at least one known other item of currency (e.g. other classes, forgeries, other denominations). In other implementations, the sennsing system 25 is arranged such that a numerical solution of the Maxwell equations are required in order to obtain the material properties of currency item 50. In such implementations, the material properties can be used as inputs to a classification technique or algorithm (e.g. Malahanobis Distance, Linear Discriminant Analysis, Feature Vector Selection).
In some implementations there can be provided a sensing system 25 configured to discriminate and/or classify an item of currency 50. Sensing system 25 can be arranged to include a processing unit 80 for driving the input 61 and receiving signals at output 68. In some implementations, sensing system 25 includes a memory unit 90 electrically coupled to processing unit 80. In some implementations, processing unit 80 is arranged as a component of authentication device 20 and electrically coupled to sensing system 25. In other implementations, processing unit 80 is integrated as a component of sensing system 25. Either arrangement is not intended to be a limitation of the scope of the disclosure.
Processing unit 80 uses the signals of input 61 and output 68 and a specified model (e.g. Equation (A)) to compute material properties of currency item 50. Processing unit 80 can be further configured to use the computed material properties of currency item 50 as inputs to a classification algorithm in order to discriminate or classify item of currency 50 from at least one other known item of currency. For example, authentication device 20, can be arranged to accept $1, $5, $10, and $20 US banknotes. In such an implementation, currency item 50 is evaluated by authentication device 20 and processing unit 80 can be arranged to determine if currency item 50 belongs to one of the aforementioned US denominations (i.e. classes). In some implementations, other classes can be used including, but not limited to, genuine, non-genuine, fit for circulation, not fit for circulation or any other class as required for the given application for authentication unit 20.
Using the example of a sensing device 25 configured for employing a PEC sensor 60, the structure of an item of currency can be estimated.
Solving Maxwell equations, a particular model can be created for a specified sensing system. For example, the complex impedance Z(ω), represented by equation (A), can be used.
It should be understood that for other types of sensing system (i.e. sensor configurations), a different model can be established by solving the Maxwell equations (as commonly known is the arts) given the particular constraints of such a sensing system (e.g. different coil geometries). Alternatively, where a model from solving the Maxwell equations is not practical, the Maxwell equations themselves can be numerically solved. In equation (A) J1(x) is the Bessel function of the first kind, first order. U12 is the first line, second column of the matrix U, and U22 is its the second line, second column and f is the frequency. In addition, μn is the nth material layer permeability [H/m], and σn its associated conductivity [S/m]. Finally N is the amount of turn for the coil wire.
If the first layer of equation (A) is assumed to have an infinite thickness, it can be thought of as acting as a half space. Choosing μ1=1 and σ1=0, the first layer of currency item 50 becomes air like. Equation (A) is an exact mathematical solution for an air-core coil for sensor 60. If the coil is inside of a ferrite pot, equation (A) still can be used as a good approximation, assuming μ0 and coil 65 geometrical dimensions are changed accordingly to fit the actual coil impedance. For example, this can be accomplished by trial and error in a known situation until a good fit has been reached.
Using an example sensing system 25 having a sensor 60 (as shown in
In some implementations, differential impedance ΔZ(ω), rather than the absolute one Z(ω) can be used. Such an approach can be used to remove the effect of the sire resistance and other common factors (e.g. temperature drift). The differential impedance can be represented by equation (B).
ΔZ(ω)=Zcoin(ω)−Zair(ω) Equation (B):
In an example of an implementation, item of currency 50 is a multi-layer coin. In equation (B), Zair(ω) corresponds to the situation where there is no coin 50 near sensor 60, while Zcoin(ω) corresponds to the situation having coin 50 present. Zair(ω) is computed just before processing coin 50, for example as an idle background processor task of sensing system 25. In such an example, Zair(ω) is an estimation at the current system temperature and set up of sensing system 25.
To illustrate the technique of the disclosure, and example will now be described. It should be understood however, that the following example is an example implementation, and in no way is intended to be limiting on the scope of the disclosure or claims.
The differential impedance of equation (B) is a complex function and therefore can be split into two terms. In an implementation of the disclosure, the differential impedance can be investigated using an inductive part ΔL(ω) and a resistive part ΔR(ω). Each can be represented by equations (C) and (D) respectively.
For example, when there is not item of currency 50 (e.g. a two layer coin) in the presence of sensing system 25, ΔL(ω)=ΔR(ω)=0∀ω≧0. When an item of currency 50 is in the presence of sensing system 25, the system becomes non-linear in which ΔL(ω) and ΔR(ω) evolve with the pulsation of ω. Such a situation is a result of eddy currents developing inside of each material of currency item 50. For example, in
In some implementations, it can be necessary to account for the distance between sensor 60 and an item of currency 50. For example, in the example described and shown in
where: where ΔLO could be chosen among different definitions, such as:
ΔL0=ΔL(ω0)−ΔL(2πfθ)ω0≠2πθ E(a)
Or the simplification, for small θ:
From the above set of equations and definitions, a simplified version of Equation (E) can be represented by equation (F).
In some implementations, the structure of an item of currency 50 can be further estimated using Model inversion techniques as commonly known in the arts. Applying such techniques to equation (A) and/or equation (B) allows for the estimation of the structure of an item of currency 50 from experimental data. As an example, inversion of Z(ω) will now be described, although it is not intended as a limitation of the disclosure. For example, a similar process can be used for ΔZ(ω) without varying in scope from the present disclosure.
In the present example, experimental data is gathered from an item of currency 50 (e.g. coin or banknote) using swept frequency techniques, direct signal spread spectrum, or any suitable signals. Furthermore, in the current example the frequency domain will be focused on, however the same procedure can be used for the time domain using the inverse Fourier Transform. Once experimental data is obtained for an item of currency 50, an estimation of the coil impedance {circumflex over (Z)}(ω) can be obtained. This can be accomplished using a non-parametric approach such as Fast Fourier Transform (FFT) or by a parametric approach such as ARMAX. In the present example, the inversion can be viewed as a non-linear regression. In order to accomplish this, the empirical risk (equation (F)) associated with a pointwise loss function (equation (G)) need to be minimized.
In the above equations (F) and (G), M is the amount of samples and {right arrow over (θ)} is the parameter vector, where {right arrow over (θ)} regroups all the unknown values, which can each be layer characteristics μn, σn, zn, as well as the lift off and the geometry of coil 65 if no prior knowledge is available. Therefore the inversion solution can be represented by equation (H).
Equation (H) is a classical unconstrained least mean square (LMS) optimization, however other optimization techniques known in the arts can be used. For example, inversion techniques can include constraints and regularization since inversion problems are often ill posed, especially in a noisy condition.
In some implementations, classification of an item of currency can be made using a simpler approximation of Z(ω) (or ΔZ(ω)) which avoids the inversion of equation (A). For example, ARMAX or OE error or any other known model for retrieving {circumflex over (Z)}(ω) can be used. Although the aforementioned models are linear, by increasing their orders (i.e. poles and zeros) they can fit more complex functions and therefore give a reasonable approximation of Z(ω).
In some implementations, the coefficients of the model can be used as inputs for recognition and/or classification. In other implementations, a spectral version of {circumflex over (Z)}(ω), either from the above models or an FFT, can be used. In such implementations, it may be important to select the most relevant frequencies, to reduce the amount of computation based on the performance requirements of authentication device 20 (e.g. processing time or acceptance/rejection rates).
In some implementations of the disclosure, there can be provided a method of using currency item 50 modeling techniques as described above in conjunction with modeling of an electromagnetic sensing system 25 for predicting the complex impedance Z(ω) and/or differential impedance ΔZ(ω). In such implementations, by specifying (or proposing) a structure for a given (or theoretical) item of currency 50, a theoretical differential impedance ΔZ(ω) can be estimated. By using known Finite Element modeling techniques, the differential impedance ΔZ(ω) (or Z(ω)) can be estimated for any hypothetical item of currency. In other implementations where sensing system 25 is arranged such that model inversion is impractical, the estimated material properties can be obtained from directing solving the Maxwell equations given the constraints of the sensing system 25.
In some implementations there can be provided a method of estimating (i.e. predicting) how a proposed or new item of currency 50 structure would be sensed by a specific sensing system 25. More particularly, it is contemplated that using the methods of the disclosure one could estimate if a proposed structure (e.g. 5 layer coin of varying materials) would be sensed, and thus classified) as an already known (and possibly circulating) item of currency or any other class of item as relevant to the specific application of an authentication device 20. Such an analysis can provide a useful tool in developing new items of currency such that the probability of a newly designed item of currency being classified as another item of currency (or as a known forgery) is minimized.
Having the parameters (μn, σn, zn) of a hypothetical currency item 500 identified, along with a specific sensing system 250, an approximation of a differential impedance ΔZ(ω) (or any other relevant model of the disclosure) can be computed in step 140. In some implementations, the outcome of the method of the disclosure results in a comparison of the hypothetical item of currency 500 with known items of currency in circulation (or any subset thereof) in step 150. In some implementations the output results in a set of coefficients from the associated model that can be used with a classification technique to determine how well the hypothetical currency item 500 can be discriminated from known items of currency in optional step 160. In some implementations, the set of coefficients from the associated model can be used with a classification algorithm or fitness function (e.g. Malahanobis Distance, Feature Vector Selection, Linear Disciminant Analysis, Support Vector Machine).
In other implementations, there is provided a method for determining the structure of a new item of currency 5000 based on a differential impedance ΔZ(ω) (or any other relevant model of the disclosure) with a given sensing system 1250.
Having the parameters (μn, σn, zn) of a hypothetical currency item 5000 identified, along with a specific sensing system 2500, at least one complex impedance can be computed for the possible configurations of item of currency 5000 by varying each parameter. A proposed solution can be output for a each material layer 3000 based a comparison of the at least one complex impedance of the hypothetical item of currency 5000 and of known items of currency in circulation (or any subset thereof) in step 1400. In other implementations of the disclosure the outcome of the method of the disclosure results is a suggestion (or guidance) of other currency item parameters including, but not limited to number of material layers, type of material, and thickness of material. Such an output is based on the given constraints used (e.g. only 3 layers, or only copper and steel, or any combination of structural characteristics).
In some implementations, the theoretical material properties of currency item (e.g. currency item 5000) obtained from a model inversion are used as inputs to a classification method or algorithm. For example, when the material properties are used as inputs to a classification technique such as Linear Discriminant Analysis (LDA), a statistical separation is obtained from at least one other class of currency items. Other classification techniques can be used including, but not limited to, Malahanobis Distance, Support Vector Machine, Feature Vector Selection. In order to determine the optimal structure for an item of currency, an optimization technique (gradient distance, or a genetic algorithm) can be used to find the optimal statistical separation of currency item 5000 from all other known items of currency (or any subset thereof). For example, if LDA is used a disciminant axis and the distance between an item of currency 5000 and an at least one known item of currency (classes), at least one material property (e.g. material thickness) of currency item 5000 can be varied to determine which value of the material property maximizes the staticical separation of currency item 5000 from the respective known currency items. In such an implementation, a solution can be obtained for the establishment of a new currency item 5000 having at least one material property (e.g. material layer thickness) having been optimized and identified based on finding the maximum statistical separation of currency item 5000 from the known class used. A process flow implementations of the disclosure is shown in
In some implementations the output from step 830 can be used as one of the inputs to an optimization step 840 for example, employing a gradient distance algorithm. The optimization step 840 uses as inputs the design constraints from step 800 and how they can be varied, the Maxwell model being used in step 810, and the fitness factor from step 830. The optimization step 840 finds the optimal design parameter that result in the best fitness factor based on the constraints of all the inputs to step 840. For example, when using gradient distance, the algorithm uses the gradient to converge on a solution that optimizes the fitness factor from step 830. In some configurations there may be local maximum found using the optimization step 840 and thus other optimization techniques can be further included to determine if the local maximum found is in fact the true maximum as is commonly known in the arts.
It is contemplated that any combination of design parameters can be fixed and/or varied to establish a new item of currency 5000 as required for a given application. For example, there can be certain design constraints that are known such as manufacturing tolerances, processing of certain materials, and/or manufacturing costs.
In other implementations, the optimization step 840 from
In some implementations, the Maxwell model from step 810 requires a direct numerical solution of the Maxwell equations to determine the simulated sensor 60 signals. Such a need arises when the model deduced from the Maxwell equations is open form and/or depending on the particular sensor arrangement.
In some implementations, there can be provided a method and apparatus for classifying items of currency as shown in
In some implementations, an item of currency is inserted into currency handling machine 10 and transported to authentication device 20. Sensing system 25 obtains response information from currency item 50 and corresponding signals are obtained from sensor 60. Authentication unit 20, using processor 80, selects an initial set of design parameters in step 900. The initial set of parameters can be selected at random or in a perdetermined manner. The design parameters from step 900 are used in step 910 to produce simulated signals for an item of currency having such design parameters. The simulated signals from step 910 and the actual signals from sensor 60 are provided as inputs to step 915 for comparison. For example, the error between the two signals can be computed. The output from step 915 (e.g. computed error) is provided as an input to step 916 where by an optimization (e.g. minimization through gradient distance) is made in order to select new design parameters (or modify the initial ones) to be inputs to step 910. Since in some implementations there is not an existing mathematical inversion of the Maxwell model from step 910, an annealing technique can be implemented to iteratively cycle from between steps 900, 910, 915, and 915 until a desired minimum error (for example) is reached. The design parameters from step 900 that are selected (or identified) by the optimization technique, are then used to produce simulated signals to be provided to step 920 as inputs. Step 920 uses feature extraction to select predetermined features from the signals from step 910 and provide them as inputs to step 930. Step 930 is a classification step whereby the inserted currency item 50 is compared with at least one known currency item to determine if it is a member of that class. In some implementations, the step 930 employs a classification technique including, but not limited to, Malahanobis distance, Linear Discriminat Analysis, Support Vector Machine, and Feature Vector Selection. In some implementations, step 930 is a fitness filter. The output of step 930 provides a fitness value for use in discriminating between at least one known currency item and an inserted item of currency 50. For example, when Malahanobis Distance is used, inserted currency item can be evaluated for belonging to a certain class if the fitness value obtained from step 930 falls within a predetermined threshold.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/US2009/059354 | 10/2/2009 | WO | 00 | 6/3/2011 |
Number | Date | Country | |
---|---|---|---|
61102685 | Oct 2008 | US |