Next Article in Journal
An Overdispersed Black-Box Variational Bayesian–Kalman Filter with Inaccurate Noise Second-Order Statistics
Previous Article in Journal
Estimating Resting HRV during fMRI: A Comparison between Laboratory and Scanner Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research of Probability-Based Tunneling Magnetoresistive Sensor Static Hysteresis Model

1
School of Instrumentation and Optoelectronic Engineering, Beihang University, Beijing 100191, China
2
High Voltage Research Institute, China Electric Power Research Institute, Beijing 100192, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(22), 7672; https://doi.org/10.3390/s21227672
Submission received: 16 September 2021 / Revised: 15 November 2021 / Accepted: 16 November 2021 / Published: 18 November 2021
(This article belongs to the Section Electronic Sensors)

Abstract

:
Tunneling magnetoresistive (TMR) sensors have broad application prospects because of their high sensitivity and small volume. However, the inherent hysteresis characteristics of TMR affect its applications in high accuracy scenarios. It is essential to build a model to describe the attributes of hysteresis of TMR accurately. Preisach model is one of the popular models to describe the behavior of inherent hysteresis for TMR, whereas it presents low accuracy in high-order hysteresis reversal curves. Furthermore, the traditional Preisach model has strict congruence constraints, and the amount of data seriously affects the accuracy. This paper proposes a hysteresis model from a probability perspective. This model has the same computational complexity as the classic Preisach model while presenting higher accuracy, especially in high-order hysteresis reversal curves. When measuring a small amount of data, the error of this method is significantly reduced compared with the classical Preisach model. Besides, the proposed model’s congruence in this paper only needs equal vertical chords.

1. Introduction

TMR is an emerging commercial magnetoresistive sensor with good performance. It has a broad application prospect in automobiles, electronic storage, biomedicine, and aerospace. Hysteresis is a characteristic widely present in magnetic materials, which has a significant impact on the features of the sensor. On one side, the hysteresis characteristic affects the measurement accuracy in specific applications; on the other side, the hysteresis could be utilized for data storage.
Although TMR has been extensively studied, the research on the hysteresis of TMR is still on the way. Current investigations on TMR mainly focuses on its performance, while there is relatively little research on its magnetic hysteresis. There are two main types of research on hysteresis. Portugal INESC [1], Mark Tondra and James M. Daughton [2] improved the hysteresis characteristics of GMR by adding a bias magnetic field; Jean Lamour et al. [3] improved the hysteresis characteristics of GMR by pinning the direction of the layer. Shu-Hsien Liao et al. [4] reduced the hysteresis of TMR with an AC modulation magnetic field. Fei Xie, Roland Weiss, and Robert Weigel [5] used controlled current pulses to compensate for hysteresis.
To reduce the characteristics of hysteresis of TMR in practical application, the essential work should be the analysis and mathematical description of hysteresis characteristics. The original research on hysteresis characteristics focused on the physical model, such as the Stoner–Wohlfarth model [6]. The physical model needs deep expert knowledge on the inherent aspects of TMR and is difficult for regular researchers. The other kind of model is the phenomenon model. The phenomenon model constructs modeling hysteresis by studying the relationship between the input and output. The phenomenon model can be further classified into two types.
One type of phenomenon model is the black-box model based on a statistical model or machine learning model. This method is relatively simple to implement and does not need to pay attention to the explanation. However, this method has little apparent explanation, requires a large amount of data for training, and has stability problems and high hardware requirements to implement [7,8]. The other type of phenomenon model is the white-box method, which tries to describe the characteristics of the hysteresis based on a particular equation. This method has two primary kinds. One is based on the differential equation. It uses differential equations to describe the relationship between the input and output. Such as the Bouc-Wen model [9,10,11,12]. The other is operator-based equations. Examples include the Preisach model, KP (Krasnosel’skii-Pokrovskii) model, etc. The Preisach model is a simple, effective, and widely used model. Some researchers intend to extern its scope of application. Refs. [13,14,15,16] focused on the dynamic hysteresis model, and [17] focused on the vector hysteresis model. Some researchers try to apply it in practice. For example, ref. [18] used First Order Reversal Curve (FORC) to analyze the features of magnetic tunnel junction (MTJ) thin-film stacks. Ref. [19] used Preisach theory in reproducing the magnetization processes of the transformer core. Refs. [20,21,22] focus on the inverse model to compensate for the hysteresis. Some researchers try to improve its accuracy. Ref. [23] proposed a polynomial model that Prandtl-Ishlinskii (PI) operator can describe. Ref. [24] used a modified reduced Preisach model based on discrete empirical interpolation method (DEIM) to resolve the contradiction between the model accuracy and the number of the hysteresis kernels. Ref. [25] proposed an input-dependent Preisach model. This model is also called the nonlinear Preisach model. [26] identified the parameters of the Preisach model by the analysis method. Refs. [27,28] identified the parameters of the Preisach model by the machine learning method.
Although the Preisch model-based methods have been widely adopted to describe and analyze the characteristics of the hysteresis of TMR, these methods are either complicated or require more data. Therefore, based on the Preisach model, we explores the input and output relationship of the TMR sensor from the perspective of probability, and proposes a static hysteresis model based on probability; this model has higher accuracy than the classical Preisach model, especially in the case of a high-order rotation curve. Two kinds of sensors are tested.
The rest of this paper is organized as follows: Section 2 introduces the basic principle of the Preisach model. Section 3 details the proposed probability model. Section 4 presents experimental results and discussion. Finally, Section 5 concludes this work.

2. Preisach Introduction

In this section, the Preisach model is introduced based on the literature [6]. The magnet is considered to be composed of finite simple hysteresis operators γαβ, shown in Figure 1. The values of the operator can only be −1 or +1. The simple hysteresis operator has two switch values, the up-switch value α, and down-switch value β, with αβ. The value of the operator is 1 when the input Hα and 0 when the input Hβ; otherwise, it shows no change. In mathematical terms, it can be represented by Equation (1).
γ α β = { 1                                 H > α m a i n t a i n   v a l u e       O t h e r 1                                 H < β .
The Preisach model can be represented as
f ( t ) = α β u ( α , β ) γ α β U ( t ) d α d β .
where U(t) is the input, γαβU(t) is the value of the operators, and u ( α , β ) is the Preisach function.
The geometric interpretation of the Preisach model is shown in Figure 2. The triangle T is composed of α = β, α = α0, β = β0. The value of u(α,β) is non-zero inside the triangle T and zero outside the triangle T. A line parallel to the β-axis moves from α = β0 to α = α1 when the input monotonically increases from β0 to α1. Correspondingly, the value of all operators swept by the line becomes 1. A vertical line moves from β = α1 to β = β1 when the input monotonically decreases from α1 to β1. Correspondingly, the value of all operators swept by the line becomes −1. After a serial change of β0,α1,β1,α2,β2,α3,β3, the staircase line in the figure is formed. The points (α1,β1), (α2,β2) and (α3,β3) are reversal points. The staircase line separates the triangle T into two parts, S+ and S. The values of the points in S+ are 1 and in S are −1.
The limiting ascending branch f+ (α) is formed by the monotonically increased input from negative saturation to some value α. The limiting descending branch f (β) is formed by the monotonically decreased input from positive saturation to some value β.
The first-order descending curve f+ (α,β) is formed by a monotonically decreased input following the limiting ascending branch f+ (α). The first-order ascending curve f (α,β) is formed by a monotonically increased input following the limiting descending branch. The high-order curve can then be defined accordingly.
In the Preisach model, has
f + ( α ) = f ( α ) .
f + ( α ,   β ) = f ( β , α ) .
Equation (3) shows that the limiting ascending branch and the descending branch are odd-symmetric. Equation (4) shows that the first-order ascending curve and the descending branch are symmetric about α = −β.
There are two essential properties of the Preisach model, the wipe-out property and the congruency property. By the wipe-out property, only the alternating series of dominant input extrema are stored by the Preisach model, while all other input extrema are wiped out. By the congruency property, all minor hysteresis loops corresponding to back-and-forth variations of inputs between the same two consecutive extremum values are congruent. In the non-linear model, the congruency property becomes the property of equal vertical chords. All minor loops resulting from back-and-forth input variations between the same two consecutive extrema have equal vertical chords (output increments) for the same input values.
Equations (5) and (6) are used to calculate the output in the Preisach model. Equation (5) is used when the last stage is monotonically decreased, and Equation (6) is used when the last stage is monotonically increased.
f ( t ) = f + + k = 1 n 1 ( f M k m k f M k m k 1 ) + f M k u ( t ) f M k m n 1 ,
f ( t ) = f + + k = 1 n 1 ( f M k m k f M k m k 1 ) + f m n 1 f m n 1 , u ( t ) ,
where f(t) is the output, f+ is the output of the positive saturation, f M k m k is the first-order curve values of reversal points and f m n 1 is the limiting curve value.

3. Probability Model

Suppose that the magnet consists of a finite simple hysteresis operator; then, the hysteresis operator is shown in Figure 3.
The operator has two values, 0 and 1, while the up- and down-switch values are α and β . The operator can be represented by Equation (7):
γ α β = { 1   H α m a i n t a i n   v a l u e   o t h e r 0   H β .
The output can be written as
h ( t ) = β α u ( α , β ) γ α β U ( t ) d α d β ,
where γαβU(t) is the value of the operator, which can be 0 or 1, and u(α,β) is the power function; therefore, h(t) is the number of operators for which the value is 1. h(t) needs to be obtained using Equation (9) because the output of the sensor has negative values.
h ( t ) = f ( t ) f ( H s a t ) ,
where f(t) is the experiment datum, −Hsat is the negative saturation input and +Hsat is the positive saturation input. When the input is smaller than the negative saturation, all of the operators are 0 and h(t) = 0. When the input is bigger than the positive saturation, all of the operators are 1 and h(t) = M, where M is the number of operators.
Let
g ( t ) = h ( t ) M .
We know that h(t) is the number of the operators for which the value is 1, and M is the number of operators. Therefore, g(t) is the probability of the operators having a value of 1. Suppose the up-switch value and down-switch value of the operator are random variables, A and B, then we obtain 3.
A [ H s a t , + H s a t ] .
B [ H s a t , + H s a t ] .
According to the probability theory,
P ( A = a ) = P { α | α   a } .
P ( B = b ) = P { β | β b } .
where α and β are the up-switch value and down-switch value of the operators. When the input monotonically increases from some value below the negative saturation to some value a, the limiting ascending curve g(a) is formed. From now on, if g(t) is not specified, it starts from negative saturation. All of the operators are 0 when the input is smaller than the negative saturation. Then, the operators with up-switch values α smaller than a become 1, when the input monotonically increases to a. At this time, all of the operators with a value of 1 have up-switch values smaller than a. Therefore, g(a) is the probability of the operators with up-switch values smaller than a. This can be written as Equation (15):
g ( a ) = P ( A = a ) .
When the input is monotonically decreased from some value above the positive saturation to some value b, the limiting descending curve g(b) is formed. All of the operators are 1 when the input is more significant than the positive saturation. Then, the operators with down-switch values bigger than b become 0, as the input monotonically decreases to b. This means that the operators with a value of 1 are those with a down-switch value smaller than b. Therefore, g(b) is the probability of the operators having down-switch values smaller than b. We can then write
g ( b ) = P ( B = b ) .
The first-order decreased curve g(a,b) is formed when g(a) is followed by a subsequent monotonically decreasing to b. Therefore, g(b) can be considered g(Hmax,b). All operators with a value of 1 at g(a) have up-switch values smaller than a. Among these operators, those with down-switch values bigger than b become 0, as in g(a,b). Therefore, the operators with a value of 1 have up-switch values smaller than a and down-switch values smaller than b. We then obtain
g ( a ,   b ) = P ( A = a , B = b ) .
The first-order increased curve g(b,a) is formed when the input monotonically increases to b after g(b). Therefore, g(a) can be thought of as g (−Hmax,a). All of the operators with a value of 1 at g(b) have down-switch values smaller than b. This means that the operators with down-switch values bigger than b should be 0. As g(b,a), the operators with up-switch values smaller than a in the 0-value operators turn to 1. Therefore,
g ( b , a ) = P ( B = b ) + P ( A = a , B > b ) = P ( B = b ) + P ( A = a ) P ( A = a , B = b ) = P ( B + A ) .              
Equation (17) can be written as
g ( a , b ) = P ( A = a , B = b ) = P ( A = a ) ( P ( A = a ) P ( A = a , B = b ) .
We can see that the increase in g(b,a) is equal to the decrease in g(a,b). In addition, we obtain Equation (20),
g ( a , b ) + g ( b , a ) = P ( A = a ) + P ( B = b ) = g ( a ) + g ( b ) .
Let a = −b; then,
g ( b , b ) + g ( b , b ) = g ( b ) + g ( b ) = h ( b ) h ( H s a t ) + h ( b ) h ( H s a t ) = f ( b ) f ( H s a t ) f ( H s a t ) f ( H s a t ) + f ( b ) f ( H s a t ) f ( H s a t ) f ( H s a t ) .
If
f ( b ) = f ( b )
then, Equation (21) is equal to 1.
The second-order curve g(a1,b1,a2) is formed by monotonically increasing the input from b1 to a2 in the first-order curve g(a1,b1). As previously discussed, in the circumstance of first-order curve g(a1,b1), the values of operators with up-switch values smaller than a1 and down-switch values smaller than b1 are 1. As the input increases, the curve is changed to the seconde-order g(a1,b1,a2), and the values of operators with up-switch values smaller than a2 used to be 0 are changed to 1.
Let A1 represent the event {α < a1}, B1 represent the event {β < b1} and A2 represent the event {α < a2}. Then, we obtain
g ( a 1 , b 1 , a 2 ) = P ( A 1 B 1 + A 2 ) = P ( A 1 B 1 ) + P ( A 2 ) P ( A 1 B 1 A 2 )
If a2 > a1, then
g ( a 1 , b 1 , a 2 ) = P ( A 1 B 1 ) + P ( A 2 ) P ( A 1 B 1 ) = P ( A 2 ) .
This is in accordance with the wipe-out property.
If a2 < a1, then
g ( a 1 , b 1 , a 2 ) = P ( A 1 B 1 ) + P ( A 2 ) P ( A 2 B 1 ) = g ( a 1 , b 1 ) + g ( a 2 ) g ( a 2 , b 1 ) .
This means that the output increase from g(a1,b1) to g(a1,b1,a2) equals the output decrease from g(a2) to g(a2,b1).
The third-order curve g(a1,b1,a2,b2) is formed by the input monotonically decreasing from a2 to b2 after g(a1,b1,a2) (obviously, a2 < a1). Let B2 represent the event {β < b2}.
g ( a 1 , b 1 , a 2 , b 2 ) = P ( ( A 1 B 1 + A 2 )   B 2 ) = P ( A 1 B 1 B 2 ) + P ( A 2 B 2 ) P ( A 1 B 1 A 2 B 2 ) .
If b1 > b2, then
g ( a 1 , b 1 , a 2 , b 2 ) = P ( A 1 B 2 ) + P ( A 2 B 2 ) P ( A 2 B 2 ) = P ( A 1 B 2 ) = g ( a , d ) .
This satisfies the wipe-out propriety.
If b1 < b2, then
      g ( a 1 , b 1 , a 2 , b 2 ) = P ( A 1 B 1 ) + P ( A 2 B 2 ) P ( A 2 B 1 ) = P ( A 1 B 1 ) + P ( A 2 ) P ( A 2 B 1 ) ( P ( A 2 ) P ( A 2 B 2 ) ) = g ( a 1 , b 1 , a 2 ) ( P ( A 2 ) P ( A 2 B 2 ) ) .
This means that the change in the third-order curve equals that of the first-order curve. This is equivalent to the property of equal vertical chords in the non-linear model. As mentioned before, monotonically increasing to an means adding the probability of the new event An = {α < an} to the original event. When the input monotonically decreases to bn, the operators with values of 1 should meet the event Bn = {β < bn}. Thus, any order curve can express as
g ( a 1 , b 1 , , a n , b n ) = P ( ( ( ( A 1 B 1 + A 2 ) B 2 + ) B n 1 + A n ) B n )
= P ( { ( B n 1 + A n ) B n } ) .
In Equation (30), {(Bn−1 + An)Bn} is short for (((A1B1 + A2)B2 + …)Bn−1 + An)Bn (from Equation (29)). Next, we prove that the wipe-out property and equal vertical chords property are also true for higher-order curves.
Suppose the 2n-order curve satisfies the two properties. Then, Equation (31) can be obtained:
g ( a 1 , b 1 , , a n , b n , u ( t ) ) = P ( { ( B n 1 + A n ) B n } + U )
Then, the (2n + 1)-order curve can be written as Equation (32):
g ( a 1 , b 1 , , a n , b n , a n + 1 , u ( t ) ) = P ( { ( B n + A n + 1 ) } U ) = P ( { ( B n U + A n + 1 U ) } ) = P { B n U } + P ( A n + 1 U ) P ( { B n U A n + 1 U } ) .
If u(t) > bn, then
g ( a 1 , b 1 , , a n , b n , a n + 1 , u ( t ) ) = P ( { B n } ) + P ( A n + 1 U ) P ( { B n A n + 1 } ) = P ( { B n } ) + P ( A n + 1 ) P ( { B n } A n + 1 ) + P ( A n + 1 U ) P ( A n + 1 )   = P ( { B n } + A n + 1 ) ( P ( A n + 1 ) P ( A n + 1 U ) )                                 = g ( a 1 , b 1 , , a n , b n , a n + 1 ) ( g ( a n + 1 ) g ( a n + 1 , u ( t ) ) ) .        
Equation (33) shows that the (2n + 1)-order curve meets the equal vertical chords property.
If u(t) < bn, then
g ( a 1 , b 1 , , a n , b n , a n + 1 , u ( t ) ) = P ( { U } ) + P ( A n + 1 U ) P ( { U A n + 1 } ) = P ( { U + A n + 1 U } ) = P ( { A n + A n + 1 } U )   = P ( { A n } U )                                 = g ( a 1 , b 1 , , a n , u ( t ) ) .        
This is in accordance with the wipe-out property.
The (2n + 2)-order is a monotonically increased curve and can be written as follows:
g ( a 1 , b 1 , , a n , b n , a n + 1 , b n + 1 , u ( t ) ) = P ( { ( B n + A n + 1 ) B n + 1 } + U ) = P ( { ( B n + A n + 1 ) B n + 1 } ) + P ( U ) P ( { ( B n + A n + 1 ) B n + 1 } U ) = P ( { ( B n + A n + 1 ) B n + 1 } ) + P ( U ) P ( { ( B n + A n + 1 ) B n + 1 U } ) = P ( { ( B n + A n + 1 ) B n + 1 } ) + P ( U ) P ( { ( B n U + A n + 1 B n + 1 U ) } ) = P ( { ( B n + A n + 1 ) B n + 1 } ) + P ( U ) P ( { ( B n U ) } ) P ( A n + 1 B n + 1 U )   + P ( { ( B n U A n + 1 ) } ) .
If u(t) < an+1,
g ( a 1 , b 1 , , a n , b n , a n + 1 , b n + 1 , u ( t ) ) = P ( { ( B n + A n + 1 ) B n + 1 } ) + P ( U ) P ( B n + 1 U ) = g ( a 1 , b 1 , , b n + 1 ) + g ( u ( t ) ) g ( u ( t ) , b n + 1 ) .
The equal vertical chords property is thus proved.
If u(t) > an+1, then
g ( a 1 , b 1 , , a n , b n , a n + 1 , b n + 1 , u ( t ) ) = P ( { ( B n + A n + 1 ) B n + 1 } ) + P ( U ) P ( { ( B n U ) } ) P ( B n + 1 A n + 1 )   + P ( { ( B n A n + 1 ) } )         = P ( { ( B n B n + 1 + A n + 1 B n + 1 ) } ) + P ( U ) P ( { ( B n U ) } ) P ( B n + 1 A n + 1 )   + P ( { ( B n A n + 1 ) } )   = P ( { ( B n B n + 1 ) } ) + P ( A n + 1 B n + 1 ) P ( { ( B n A n + 1 B n + 1 ) } ) + P ( U )   P ( { ( B n U ) } ) P ( B n + 1 A n + 1 ) + P ( { ( B n A n + 1 ) } ) = P ( { ( B n ) } ) + P ( U ) P ( { ( B n U ) } ) = P ( { ( B n ) } + U ) = g ( a 1 , b 1 , b n , u ( t ) ) .
This is in accordance with the wipe-out property. Therefore, the two properties are valid for any order curve. We can use these two properties to calculate any curve starting from the saturation input.

4. Experiment and Results

4.1. Experiment Introduction

The experimental scheme is shown in Figure 4. The high-precision power supply drives the Helmholtz coil through the power amplifier to generate a magnetic field. Sensor 1, sensor two, and Tesla meter probes are fixed in the uniform magnetic field area of the Helmholtz coil. Measure the magnetic field with a Tesla meter. Two high-precision voltmeters are used to measure the outputs of the sensors. The computer controls the whole acquisition process.
The flow chart of acquisition is shown in Figure 5. First, the computer reads the required magnetic field, calculates an approximate voltage value, controls the high-precision power supply to output the voltage value, and collects the Tesla meter reading. Adjust the output voltage value according to the difference between the Tesla meter reading and the target value when the reading is stable. This operation continues until the target value is obtained. Then read the data of the two sensors, wait until the sensor data is sound, record the data, and then read the next target value.

4.2. Results and Discussion

The limiting curves f(u(t)) acquired following the experiment are shown in Figure 6. fifty-one sample points were chosen, denoted as f(a1), f(a2), , f(a21), and the odd symmetry of the limiting descending branch is also plotted in Figure 6 to study the symmetry of the limiting curves.
In Figure 6, the y-axis represents the sensor’s output, and the x-axis represents the magnetic field. We can see that the odd symmetry of the descending branch is not identical to that of the ascending branch. This means that the limiting ascending curves and limiting descending curves are not symmetric. This might be due to the material and manufacturing process of the sensor. There are 50 first-order curves corresponding to 51 acquisition points. The first-order curves are shown in Figure 7. Every color curve represents a first-order curve.
According to the discussion in Section 3, we obtain Equation (38):
g ( u ( t ) ) = h ( u ( t ) ) h ( + Hsat ) = f ( u ( t ) ) f ( H s a t ) f ( H s a t ) f ( H s a t ) .
The limiting curve g(ai) and several first-order curves g(ai,bj) can be obtained by Equation (38). The limiting curve g(ai) is shown in Figure 8 and the first-order curves g(ai,bj) are shown in Figure 9. Every color curve represents a first-order curve.
Any limiting ascending curve g(a) can be obtained by linearly interpolating g(ai). Then, any first-order curve g(a,b) needs to be calculated. One method is to directly interpolate the first-order curves g(ai,bj). However, because the distance between the first-order revolution curves does not have apparent characteristics, as shown in Figure 10, the accuracy of this method is not high.
Another method to calculate any first-order curve g(a,b) is as follows. Any first-order curve g(a,b) represents the probability of up-switch values smaller than a and down-switch values smaller than b. This can be written as Equation (39):
g ( a , b ) = P ( A = a , B = b ) .
According to the property of probability, Equation (40) can be written as
g ( a , b ) = P ( A = a , B = b ) = P ( B = b | A = a ) P ( A = a ) .
In Equation (40), P(A) can be obtained by interpolating g(ai). P(A,B) can be calculated when P(B|A) is determined. Then, experimental P(Bj|Ai) is needed and P(Ai,Bj) are obtained. Therefore, we can obtain the conditional probability using Equation (41). The conditional probability curves are shown in Figure 11.
P ( B j |   A i ) = P ( A i = a i , B j = b j ) P ( A i = a i ) = g ( a i , b j ) g ( a i ) .
In Figure 11, the x axis represents the magnetic field, and the y axis represents the probability. Every curve in the figure represents one conditional probability. The conditional probability curves intersect the line at y = 100% at the reversal points, such as A1(a1,1), A2(a2,1), A3(a3,1), and A4(a4,1). Some features can be seen in the figure. The probability is 0 at negative saturation and 1 when the magnetic saturation is more significant than the reversal points. The probability slope with more minor reversal points is more significant than that with more prominent reversal points—the greater the reversal points, the denser the curves. The spaces between the curves decrease with the increase in the reversal points. The ratio of the spacing between adjacent curves remains roughly unchanged. For example, the ratio of A2A3 and A3A4 is roughly equal to the ratio of CF and FG; this can be written as Equation (42):
C F F G = A 2 A 3 A 3 A 4 .
Next, we illustrate how to calculate P(B|A). As shown in Figure 10, the straight line x = b intersects curves A2C, A3D and A4E at the points C (b,P(b|a2)), D (b,P(b|a3)) and E (b,P(b|a4)), respectively. Suppose the probability of point D is needed. This means that curves A2C and A4E are known; however, curve A3D is not. The probabilities of points C and E can be calculated by interpolating the curves A2C and A4E, respectively. Then, according to point C’ and E’s probability, the abscissa g and m of point G and point M are obtained by interpolating curve A4E and curve A2M respectively. According to the previous assumptions, the lengths of EK and CF can be obtained. Moreover, the probability of P(B|A3) = P(b|a3) needs to be calculated. It can be supposed that curve KF is straight line. Then, the triangle CFD is similar to triangle DEK. Therefore, the following equation can be obtained:
C D D E = C F K E = C G * A 2 A 3 A 2 A 4 M E * A 3 A 4 A 2 A 4 = C G * A 2 A 3 M E * A 3 A 4 .
Substitute the coordinate value to Equation (43); then, we obtain
P ( b |   a 2 ) P ( b |   a 3 ) P ( b |   a 3 ) P ( b |   a 4 ) = b g * ( a 2 a 3 ) ( m b ) * ( a 3 a 4 ) .
Then, P(b| a 3 ) can be calculated by Equation (45):
P ( b |   a 3 ) = ( m b ) * ( a 3 a 4 ) P ( b |   a 2 ) + ( m b ) * ( a 3 a 4 ) P ( b |   a 4 ) ( b g ) * ( a 2 a 3 ) + ( m b ) * ( a 3 a 4 ) .
If the point C (b,P(b|a2)) is needed, then curves A1H and A3D are known and curve A2C is not. The triangle BCA2 is approximately similar to triangle BDA3. Therefore,
B C C D = B A 2 H D = B A 2 I D * A 2 A 3 A 1 A 3 .
Substitute the coordinate value to Equation (46); then,
1 P ( b |   a 2 ) P ( b |   a 2 ) P ( b |   a 3 ) = b a 2 ( i b ) * a 2 a 3 a 1 a 3 .
Then,
P ( b |   a 2 ) = ( i b ) ( a 2 a 3 ) + ( b a 2 ) ( a 1 a 3 ) P ( b |   a 3 ) ( b a 2 ) ( a 1 a 3 ) + ( i b ) ( a 2 a 3 ) .
According to the above, any conditional probability P(B|A) can be calculated. Moreover, any P(A) can be obtained by interpolating the limiting curve. Then, the probability P(AB) of any first-order reversal curve can be obtained and, according to Section 3, any order reversal curves can be predicted. In our experiment, two types of sensors are used. The test curves were acquired, and the probability model and the classical Preisach model were used to predict the curve. The results of the two models were compared with the experiment data; then, the experimental data and the results of the classical Preisach model were turned into probability data by Equation (38). The results of the first kind of sensor are shown in Figure 12. The first sensor’s error of the probability model and that of the classic Preisach model are shown in Figure 13, in which the y axis represents the probability, and its unit is 1. The results of the second sensor are shown in Figure 14. The error is shown in Figure 15, in which the y axis represents the probability, and its unit is 1.
The results show that the error of the probability model is clearly smaller than that of the classic Preisach model. At the beginning stage, the errors of the two models are similar. That is because the beginning stage is the limiting curve. Both models use the limiting curve linear interpolation method to predict the limiting curve. However, in the subse-quent data, the error of the probability model is less than that of the classical Preisach model, especially in the high-order reversal curve. It can be seen that in the last section of data, the error fluctuation range of classical Preisach tends to enlarge. In contrast, the error fluctuation of the probability model tends to narrow. Both models can predict the sensor’s output; however, the probability model is more accurate than the classic Preisach model.
Various indicators of error are shown in Table 1. show that the mean error of one sensor between the Preisach model and probability model are 0.0245 and 0.0197; the variance is 4.5650 × 10−5 and 4.1695 × 10−5, and the maximum error is 0.0347 and 0.0209. And, the mean error of the other sensor between the Preisach model and probability model are 0.0156 and 0.0146, respectively, the variance is 8.3547 × 10−6 and 8.2343 × 10−6, and the maximum error is 0.0209 and 0.0187. The reduced ratio is listed as the relevant index of the Preisach model minus the relevant index of the probability model and then divided by the appropriate index of the Preisach model. It can be seen that the lift of sensor 1 is greater than that of sensor 2. This is because the magnetic hysteresis of sensor 1 is relatively large relative to its range, while that of sensor 2 is relatively small relative to its range.
The amount of data of the first-order reversal curve, that is, the number of reversal points, contradicts the the model’s accuracy. The model’s accuracy is higher with more reversal points, but the time required is significantly increased. With fewer turning points, the time required is less, but the accuracy may be reduced. So we also verified the model by using sparse data. This paper uses the first-order reversal curve data with an interval of 5 OE between reversal points to test two models and two sensors. The results and errors of sensor one are shown in Figure 16 and Figure 17, and the results and mistakes of sensor two are shown in Figure 18 and Figure 19.
The error indicators are shown in Table 2. The mean error of the two models on one sensor is 0.0135 and 0.0058, the variance is 1.7869 × 10−4 and 6.3849 × 10−5, and the maximum error is 0.0667 and 0.0204. On the other sensor, the mean error of the Preisach model and probability model are 0.0082 and 0.0070, respectively, the variance is 2.7725 × 10−5 and 1.1834 × 10−5, and the maximum error is 0.0342 and 0.0132. It can be seen that under the condition of using sparse data, the improvement of the probability model is more evident than that of the Preisach model. Especially in the case of a high-order reversal curve, the error of the Preisach model increases. In contrast, the error of the probability model hardly changes, and even the fluctuation range tends to narrow. Moreover, the complexity of the model and the required data of the probability model are approximately the same as the classic Preisach model.

5. Conclusions

This paper analyzes hysteresis from the perspective of probability and proposes a probability-based hysteresis modeling method. This method has the same complexity and required data level as the classical model but achieves higher accuracy. The test results show that the error of the probability model is smaller than that of the Preisach model, especially under the condition of a high-order reversal curve.
It can be seen that the error of the probability model relative to the Preisach model is reduced, especially in the case of a high-order reversal curve. The amount of the first-order reversal points is a primary reason affecting the accuracy of the model. This paper uses two steps, e.g., 1 OE and 5 OE to evaluate the two models, which shows that the probability model has higher accuracy than the classical model in sparse data, which shows that this method can effectively improve test efficiency while improving accuracy.
Furthermore, the reduction of the number of first-order reversal points will significantly shorten the acquisition time. The results show that the improvement of the probability model is more evident than that with a break of 1 OE. This helps to shorten the acquisition time. In addition, the congruence property of the Preisach model requires the consistency of small loops, while the probability model only requires the equality of vertical chords.

Author Contributions

Y.L. proposed the method, performed the experiments, analyzed the data, wrote the paper; L.W., H.Y. and Z.Q. review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ku, W.; Silva, F.; Bernardo, J.; Freitas, P.P. Integrated giant magnetoresistance bridge sensors with transverse permanent magnet biasing. J. Appl. Phys. 2000, 87, 5352–5355. [Google Scholar] [CrossRef]
  2. Tondra, M.; Daughton, J.M.; Nordman, C.; Wang, D.; Taylor, J. Micromagnetic design of spin dependent tunnel junctions for optimized sensingperformance. J. Appl. Phys. 2000, 87, 4678–4681. [Google Scholar] [CrossRef] [Green Version]
  3. Negulescu, B.; Lacour, D.; Montaigne, F.; Gerken, A.; Paul, J.; Spetter, V.; Marien, J.; Duret, C.; Hehn, C. Wide range and tunable linear magnetic tunnel junction sensor using two exchange pinned electrodes. Appl. Phys. Lett. 2009, 95, 112502. [Google Scholar] [CrossRef] [Green Version]
  4. Liao, S.-H.; Huang, H.-S.; Sokolov, A.; Yang, Y.; Liu, Y.-F.; Yin, X.; Wilson, A.; Trujilo, A.; Ewing, D.; Liou, S.-H. Hysteresis Reduction in Tunneling Magnetoresistive Sensor with AC Modulation Magnetic Field. IEEE Trans. Magn. 2021, 57, 1–4. [Google Scholar] [CrossRef]
  5. Xie, F.; Weiss, R.; Weigel, R. Hysteresis Compensation Based on Controlled Current Pulses for Magnetoresistive Sensors. IEEE Trans. Ind. Electron. 2015, 62, 7804–7809. [Google Scholar] [CrossRef]
  6. Mayergoyz, D. Mathematical Models of Hysteresis and Their Applications; Academic Press: Cambridge, MA, USA, 2003. [Google Scholar]
  7. Xu, Q.; Wong, P.K. Hysteresis modeling and compensation of a piezostage using least squares support vector machines. Mechatronics 2011, 21, 1239–1251. [Google Scholar] [CrossRef]
  8. Branko, K.; Alenka, M.; Nebojsa, M. Mathematical modelling of frequency-dependent hysteresis and energy loss of FeBSiC amorphous alloy. J. Magn. Magn. Mater. 2017, 422, 37–42. [Google Scholar]
  9. Yu, G.R.; You, C.S.; Hong, R.J. Self-tuning fuzzy control of a piezoelectric actuator system. In Proceedings of the 2006 IEEE International Conference on Systems, Man and Cybernetics, Taipei, Taiwan, 8–11 October 2006; Volume 2, pp. 1108–1113. [Google Scholar]
  10. Payam, A.F.; Fathipour, M.; Yazdanpanah, M.J. A backstepping controller for piezoelectric actuators with hysteresis in nanopositioning. In Proceedings of the 2009 4th IEEE International Conference on Nano/Micro Engineered and Molecular Systems, Shenzhen, China, 5–8 January 2009; pp. 711–716. [Google Scholar]
  11. Chang, K.M. Adaptive control for piezoelectric positioning system. In Proceedings of the 2006 IEEE International Symposium on Industrial Electronics, Montreal, QC, Canada, 9–13 July 2006; pp. 37–42. [Google Scholar]
  12. Rakotondrabe, M. Bouc-Wen modeling and inverse multiplicative structure to compensate hysteresis nonlinearity in piezoelectric actuators. IEEE Trans. Autom. Sci. Eng. 2011, 8, 428–431. [Google Scholar] [CrossRef] [Green Version]
  13. Grech, C.; Buzio, M.; Pentella, M.; Sammut, N. Dynamic Ferromagnetic Hysteresis Modelling Using a Preisach-Recurrent Neural Network Model. Materials 2020, 13, 2561. [Google Scholar] [CrossRef]
  14. Dang, X.; Tan, Y. Neural networks dynamic hysteresis model for piezoceramic actuator based on hysteresis operator of first-order differential equation. Physica B 2005, 365, 173–184. [Google Scholar] [CrossRef]
  15. Mayergoyz, I.D. Dynamic preisach models of hysteresis. IEEE Trans. Magn. 1988, 24, 6. [Google Scholar] [CrossRef]
  16. Dupré, L.; Bertotti, G.; Basso, V.; Fiorillo, F.; Melkebeek, J. Generalisation of the dynamic Preisach model toward grain oriented Fe-Si alloys. Physica B 2000, 275, 202–206. [Google Scholar] [CrossRef]
  17. Zhu, L.; Jiang, J.; Wu, W.; Li, W.; Lu, K.; Koh, C.S. An Improved Anisotropic Vector Preisach Model for Nonoriented Electrical Steel Sheet Based on Iron Loss Separation Theory. Math. Probl. Eng. 2020, 2020, 7348648. [Google Scholar] [CrossRef]
  18. Abugri, J.B.; Visscher, P.B.; Gupta, S.; Chen, P.J.; Shull, R.D. FORC+ analysis of perpendicular magnetic tunnel junctions. J. Appl. Phys. 2018, 124, 043901. [Google Scholar] [CrossRef] [Green Version]
  19. Andreev, M.; Askarov, A.; Suvorov, A. Design of the magnetic hysteresis mathematical model based on Preisach theory. Electr. Eng. 2019, 101, 3–9. [Google Scholar] [CrossRef]
  20. Davino, D.; Giustiniani, A.; Visone, C. Fast Inverse Preisach Models in Algorithms for Static and Quasistatic Magnetic-Field Computations. IEEE Trans. Magn. 2008, 44, 6. [Google Scholar] [CrossRef]
  21. Sarkera, P.C.; Guo, Y.; Lu, H.; Zhu, J.G. A generalized inverse Preisach dynamic hysteresis model of Fe-based amorphous magnetic materials. J. Magn. Magn. Mater. 2020, 514, 167290. [Google Scholar] [CrossRef]
  22. Reimers, A.; Della Torre, E. Fast Preisach-Based Magnetization Model and Fast Inverse Hysteresis Model. IEEE Trans. Magn. 1998, 34, 6. [Google Scholar] [CrossRef]
  23. Ru, C.; Chen, L.; Shao, B.; Rong, W.; Sun, L. A hysteresis compensation method of piezoelectric actuator: Model, identification and control. Control Eng. Pract. 2009, 17, 1107–1114. [Google Scholar] [CrossRef]
  24. Wang, D.; Wang, Y.; Lu, Y.; Chen, B.; Wang, L.; Wu, H. Reduced Preisach Model: Beyond Discrete Empirical Interpolation Method. IEEE Access 2019, 7, 155563–155573. [Google Scholar] [CrossRef]
  25. Friedman, G.; Mayergoyz, I.D. Input-dependent Preisach model and hysteretic energy losses. J. Appl. Phys. 1991, 8, 4611. [Google Scholar] [CrossRef]
  26. Ktenaa, A.; Fotiadisa, D.I.; Spanosb, P.D.; Massalasc, C.V. A Preisach model identification procedure and simulation of hysteresis in ferromagnets and shape-memory alloys. Physica B 2001, 306, 84–90. [Google Scholar] [CrossRef]
  27. Hergli, K.; Marouani, H.; Zidi, M.; Fouad, Y.; Elshazly, M. Identification of Preisach hysteresis model parameters using genetic algorithms. J. King Saud Univ. Sci. 2019, 31, 746–752. [Google Scholar] [CrossRef]
  28. E, J.; Qian, C.; Zhu, H.; Peng, Q.; Zuo, W.; Liu, G. Parameter-identification investigations on the hysteretic Preisach model improved by the fuzzy least square support vector machine based on adaptive variable chaos immune algorithm. J. Low Freq. Noise Vib. Act. Control 2017, 36, 227–242. [Google Scholar]
Figure 1. Preisach operator.
Figure 1. Preisach operator.
Sensors 21 07672 g001
Figure 2. The geometric interpretation of the Preisach model.
Figure 2. The geometric interpretation of the Preisach model.
Sensors 21 07672 g002
Figure 3. The hysteresis operator of the probability model.
Figure 3. The hysteresis operator of the probability model.
Sensors 21 07672 g003
Figure 4. Experimental scheme.
Figure 4. Experimental scheme.
Sensors 21 07672 g004
Figure 5. Flow chart.
Figure 5. Flow chart.
Sensors 21 07672 g005
Figure 6. The limiting curves.
Figure 6. The limiting curves.
Sensors 21 07672 g006
Figure 7. The first-order curves.
Figure 7. The first-order curves.
Sensors 21 07672 g007
Figure 8. The limiting curve of the probability model.
Figure 8. The limiting curve of the probability model.
Sensors 21 07672 g008
Figure 9. The first-order curves of the probability model.
Figure 9. The first-order curves of the probability model.
Sensors 21 07672 g009
Figure 10. A partial enlarged view of the first-order curves.
Figure 10. A partial enlarged view of the first-order curves.
Sensors 21 07672 g010
Figure 11. The conditional probability curves.
Figure 11. The conditional probability curves.
Sensors 21 07672 g011
Figure 12. Test results of the first kind sensor.
Figure 12. Test results of the first kind sensor.
Sensors 21 07672 g012
Figure 13. Errors of the first kind sensor.
Figure 13. Errors of the first kind sensor.
Sensors 21 07672 g013
Figure 14. Test results of the second kind sensor.
Figure 14. Test results of the second kind sensor.
Sensors 21 07672 g014
Figure 15. Errors of the second kind sensor.
Figure 15. Errors of the second kind sensor.
Sensors 21 07672 g015
Figure 16. Test results of the first kind sensor for sparse data.
Figure 16. Test results of the first kind sensor for sparse data.
Sensors 21 07672 g016
Figure 17. Errors of the first kind sensor for the sparse data.
Figure 17. Errors of the first kind sensor for the sparse data.
Sensors 21 07672 g017
Figure 18. Test results of the second kind sensor for the sparse data.
Figure 18. Test results of the second kind sensor for the sparse data.
Sensors 21 07672 g018
Figure 19. Errors of the second kind sensor for the sparse data.
Figure 19. Errors of the second kind sensor for the sparse data.
Sensors 21 07672 g019
Table 1. Errors.
Table 1. Errors.
SensorsModelsAverage ValueReduced ProportionAbsolute MaximumReduced ProportionVarianceReduced Proportion
Sensor 1Preisach model0.024519.6%0.034739.8%4.5650 × 10−58.6%
Probability model0.01970.02094.169 × 10−5
Sensor 2Preisach model0.01566.5%0.020910.5%8.3547 × 10−614.4%
Probability model0.01460.01878.2343e × 10−6
Table 2. Errors.
Table 2. Errors.
SensorsModelsAverage ValueReduced ProportionAbsolute MaximumReduced ProportionVarianceReduced Proportion
Sensor 1Preisach model0.013557.0%0.066769.4%1.7869 × 10−464.3%
Probability model0.00580.02046.3849 × 10−5
Sensor 2Preisach model0.008214.6%0.034261.4%2.7725 × 10−557.3%
Probability model0.00700.01321.1834 × 10−5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Wang, L.; Yu, H.; Qian, Z. Research of Probability-Based Tunneling Magnetoresistive Sensor Static Hysteresis Model. Sensors 2021, 21, 7672. https://doi.org/10.3390/s21227672

AMA Style

Li Y, Wang L, Yu H, Qian Z. Research of Probability-Based Tunneling Magnetoresistive Sensor Static Hysteresis Model. Sensors. 2021; 21(22):7672. https://doi.org/10.3390/s21227672

Chicago/Turabian Style

Li, Yutao, Liliang Wang, Hao Yu, and Zheng Qian. 2021. "Research of Probability-Based Tunneling Magnetoresistive Sensor Static Hysteresis Model" Sensors 21, no. 22: 7672. https://doi.org/10.3390/s21227672

APA Style

Li, Y., Wang, L., Yu, H., & Qian, Z. (2021). Research of Probability-Based Tunneling Magnetoresistive Sensor Static Hysteresis Model. Sensors, 21(22), 7672. https://doi.org/10.3390/s21227672

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop