Behavioral Modeling of Memristor-Based Rectiﬁer Bridge

: In electrical engineering, radio engineering, robotics, computing, control systems, etc., a lot of nonlinear devices are synthesized on the basis of a nanoelement named memristor that possesses a number of useful properties, such as passivity, nonlinearity, high variability of parameters, nonvolatility, compactness. The efﬁciency of this electric element has led to the emergence of many memristor technologies based on different physical principles and, as a result, to the occurrence of different mathematical models describing these principles. A general approach to the modeling of memristive devices is represented. The essence is to construct a behavioral model that approximates nonlinear mapping of the input signal set into the output signal set. The polynomials of split signals, which are adaptive to the class of input signals, are used. This adaptation leads to the model’s simpliﬁcation important in practice. Multi-dimensional polynomials of split signals are built for the rectiﬁer bridge at harmonic input signals. The modeling error is estimated in the mean-square norm. It is shown that the accuracy of the modeling is increased in the case of using the piecewise polynomial with split signals.


Introduction
Recently, the element base for the designing of nonlinear dynamic devices is changed due to the emergence of elements with new properties and technologies for their creation. More and more devices are synthesized based on memristors or memristive structures [1][2][3][4]. This process results from the effective characteristics of a memristor as a passive electrotechnical element. The memristor possesses the following useful properties [1][2][3][4]: • Being an analog element, its resistance can take any values. This is a positive property in comparison with any binary element whose value can be either 0 or 1. Such variability of resistance is realized in one element, the memristor size is reduced to several nanometers and the response rate is reduced to nanoseconds.

•
The memristor does not store a charge. This means that it is not prone to charge leaks, which must be dealt with when going to nanometer-scale microcircuits.

•
The memristor is a non-volatile element, and data can be stored as long as the materials from which it is made exist. • Memristors placed on crossing conductors (crossbars) can be used to form densely packed memory. • Many memristor materials are compatible with complementary metal-oxide-semiconductor (CMOS) technology.
A memristor was theoretically described first by L. Chua in 1971 [5]. There were a lot of investigations devoted to elements with resistive switching [6][7][8][9][10][11] and only in 2008 the link between resistive switching and the theory of memristors was established in the Hewlett Packard (HP) laboratory [12]. In [12] the authors represented a variant of the ideal memristor model, which doesn't completely describe any resistive switching element [13].
of these cases the current flow on the AC side is non-linear which causes higher order current harmonics. For mains powered devices the magnitude of these harmonics is to be limited [42]. An additional power factor correction (PFC) circuit is needed. The utilization of memristors instead of diodes for the rectifier bridge is a new investigation trend. Due to the nature of a physically realized memristor it is not yet usable for power electronic applications. In order to get an unidirectional valve effect for the current, as it is the case for the rectifier diode, a bipolar switching memristor is required and the operating frequency needs to be in the lower Hertz region [43,44]. For low power applications the memristorbased rectifier bridge can be used in synapse-based neural networks with programmable synaptic weights [45,46]. Another application is image processing; by adding a capacitor at the output of the memristive bridge embodies a low-pass filter with time-varying cut-off frequency and transfer function [47]. Thus, the use of memristors in the well-known Graetz circuit essentially expands its application. For investigation of memristor-based devices, both empirical models adapted from the ideal memristor model and models based on physical considerations [48][49][50][51][52][53] play an important role.

Multi-Dimensional Split Polynomial as a Behavioral Nonlinear Model
A split polynomial approximates the operator of a nonlinear system under an assigned error. This error is estimated in accordance with a norm. We take the well-known meansquare norm as the most popular in practice. The model maps the domain of the input signals into the corresponding domain of the output signals. The model parameters result from solving the following approximation task where s o out (t), s out (t) are the output signals of a nonlinear system and a model correspondingly; C is a vector of the model parameters. The solution of the optimization task (1) is determined on using the given sets of the input and output signals. This solution is global optimum in the case of linear-in-the-parameter model [31][32][33][34]. It encourages the popularity of polynomial models.
A designed split polynomial is to be unique, dynamic and nonlinear. The unique property means there is no mapping of one input signal into some output signals. The model block-scheme depicted in Figure 1 ensures the mentioned properties. The first block named Splitter transforms every scalar input signal s in (A, t) into the vector signal S sp (A, t), named split signal: where A, A ∈ G a is the parameters vector of the input signal s in (A, t) from the domain G a of parameters; t, t ∈ G t is the time variable existing in the time domain G t .
Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 17 many others [42]. A rectifier bridge can operate as a passive circuit using diodes or as an active circuit using bipolar junction transistors (BJTs) or field effect transistors (FETs). In any of these cases the current flow on the AC side is non-linear which causes higher order current harmonics. For mains powered devices the magnitude of these harmonics is to be limited [42]. An additional power factor correction (PFC) circuit is needed. The utilization of memristors instead of diodes for the rectifier bridge is a new investigation trend. Due to the nature of a physically realized memristor it is not yet usable for power electronic applications. In order to get an unidirectional valve effect for the current, as it is the case for the rectifier diode, a bipolar switching memristor is required and the operating frequency needs to be in the lower Hertz region [43,44]. For low power applications the memristor-based rectifier bridge can be used in synapse-based neural networks with programmable synaptic weights [45,46]. Another application is image processing; by adding a capacitor at the output of the memristive bridge embodies a low-pass filter with timevarying cut-off frequency and transfer function [47]. Thus, the use of memristors in the well-known Graetz circuit essentially expands its application. For investigation of memristor-based devices, both empirical models adapted from the ideal memristor model and models based on physical considerations [48][49][50][51][52][53] play an important role.

Multi-Dimensional Split Polynomial as a Behavioral Nonlinear Model
A split polynomial approximates the operator of a nonlinear system under an assigned error. This error is estimated in accordance with a norm. We take the well-known mean-square norm as the most popular in practice. The model maps the domain of the input signals into the corresponding domain of the output signals. The model parameters result from solving the following approximation task where s out o (t), s out (t) are the output signals of a nonlinear system and a model correspondingly; C is a vector of the model parameters. The solution of the optimization task (1) is determined on using the given sets of the input and output signals. This solution is global optimum in the case of linear-in-the-parameter model [31][32][33][34]. It encourages the popularity of polynomial models.
A designed split polynomial is to be unique, dynamic and nonlinear. The unique property means there is no mapping of one input signal into some output signals. The model block-scheme depicted in Figure 1 ensures the mentioned properties. The first block named Splitter transforms every scalar input signal s in (A,t) into the vector signal S sp (A,t), named split signal: where A, A ∈ G a is the parameters vector of the input signal s in (A,t) from the domain G a of parameters; , t ∈ G t is the time variable existing in the time domain G t . The mentioned transformation is carried out according to the following conditions [35][36][37][38]: • for all A ∈ G a , t ∈ G t , there is an inequality The mentioned transformation is carried out according to the following conditions [35][36][37][38]: • for all A ∈ G a , t ∈ G t , there is an inequality 4 of 15 Equations (3) and (4) mathematically describe the situation, that the phase portraits, which are built on the axes of the elements of split vector (2), don't cross and touch each other as well as don't pass through zero. The splitter is a dynamic circuit that provides unique and dynamic properties of a model. Splitting signals can be implemented by linear and nonlinear, stationary and non-stationary dynamic circuits [35][36][37][38]. The main implementation criterion is, the lower the length of a split vector (2), the lower the dimensionality and consequently complexity of a model. The splitter depends on the input signal, since it splits the particular input signal. Therefore, one can say, the split model is adaptive to the input signal, thus it can be simpler in comparison with general polynomial, regression or neural models.
The splitter is followed by the block named Nonlinear Memoryless Transformer (NMT), that creates the model nonlinearity. The NMT transforms the split vector signal into the scalar output signal of a model. This nonlinearity can be described by a polynomial, a fraction and a regression, as well as neural networks. The most popular form is a polynomial, that is used for modeling in this paper. The polynomial is written The power of polynomial (5) equals J = ∑ m r=1 J r . The elements of the split vector signal (2) form basis functions in (5), that is why the polynomial is multi-dimensional. The dimension of the model is the length of the split vector (2).
Further, we use the described split polynomial for the mathematical modeling of the memristor-based rectifier bridge. By simulating the rectifier bridge in the LTspice program, we obtain the sets of the input and output signals required for building this model. This program is famous for simulating a majority of analog circuits and it includes some of the well-known memristor models [43].

Forming the Sets of Input and Output Signals of the Memristor-Based Rectifier Bridge
The circuit of the memristor-based rectifier bridge simulated in LTspice is shown in Figure 2. A sinusoidal voltage source with variable amplitude A and frequency f acts as the excitation (the input signal). A simple 1 kΩ resistor acts as a load. The voltage of the resistive load is the output signal. At a positive half period of the sinusoidal input voltage, the memristance of the two P-xsv-N memristors (see Figure 2) is decreasing, while the memristance of the two N-xsv-P memristors is increasing until saturation is reached. At a negative half period of the input voltage the process is reversed. During one period there is a change of processes between both memristances. Thus, in the ideal case the output voltage is only positive, and this pulsed voltage has double the frequency of the input signal. In the real case there are the temporary short negative voltage-time areas in the output signal, due to the time shift in the state variable change. The circuit comprises memristors described by the Yakopcic model. The parameters of the Yakopcic model are represented in LTspice and have the following values: [54,55]. The equations of the memristor model are implemented in a sub-circuit called memris-tor_yakopcic.sub. The Yakopcic model will be explained in detail below. The directives (.param and .tran) seen in Figure 2 are only used as an example. called memristor_yakopcic.sub. The Yakopcic model will be explained in detail below. The directives (.param and .tran) seen in Figure 2 are only used as an example. For the investigation of the rectifier, we have set the input signal parameters to the amplitude/frequency pairs (12 V/1 Hz), (13 V/2 Hz), (14 V/3 Hz), (15 V/4 Hz) and (16 V/5 Hz). The time window for simulation was set depending on the frequency and the number of signal periods N p . We have chosen N p = 3 in order to let startup transients settle. The LTspice output signal does not have equally spaced timestamps, as the program adjusts the time steps automatically dependent on the temporal derivative of the signal. Only the final period of the output voltage was kept for post-processing. For interpolation we used N s = 2001 equally spaced samples. As a final step, all waveforms were normalized to the same angular frequency ω = 6π rad. By this the different maxima and minima are aligned.

The Yakopcic Model of a Memristor in LTspice
The memristor model used in the circuit depicted in Figure 2 plays the important role for signal processing. In the past, several different mathematical memristor models were published; each of them aiming for a certain application. The first model was created by HP Laboratory in 2008 [12]. The basic current (I)-voltage (V) relationship can be found there. In 2009, the original HP model was extended by introducing the factor η , which describes the state variable motion with respect to the polarity of the applied voltage. Furthermore, a window function for the state variable motion was added [56]. The window function was modified by Biolek for a more realistic behavior of the state variable motion near the boundaries [57]. The first hardware related model was introduced by the Air Force Research Lab [58]. It does not use a state variable; the resistance change rate was defined directly. Later, models were implemented that describe a hyperbolic sine I-V relationship, like the one by the University of Michigan [59]. Previous models could only emulate sinusoidal input signal correctly, not general signal shapes. The hyperbolic sine model can approximate the I-V relationship of the MIM junction better.
For the implementation of a memristor model in LTspice we use the generalized hyperbolic sine model by Yakopcic et al. [55]. In comparison with the above-mentioned memristor models, the Yakopcic model uses more general equations for the memristor dynamics. Via fitting parameters, the model can be adapted to a variety of published memristor devices. The LTspice schematic of the sub-circuit can be seen in Figure 3. Only the final period of the output voltage was kept for post-processing. For interpolation we used N s = 2001 equally spaced samples. As a final step, all waveforms were normalized to the same angular frequency ω = 6π rad. By this the different maxima and minima are aligned.

The Yakopcic Model of a Memristor in LTspice
The memristor model used in the circuit depicted in Figure 2 plays the important role for signal processing. In the past, several different mathematical memristor models were published; each of them aiming for a certain application. The first model was created by HP Laboratory in 2008 [12]. The basic current (I)-voltage (V) relationship can be found there. In 2009, the original HP model was extended by introducing the factor η, which describes the state variable motion with respect to the polarity of the applied voltage. Furthermore, a window function for the state variable motion was added [56]. The window function was modified by Biolek for a more realistic behavior of the state variable motion near the boundaries [57]. The first hardware related model was introduced by the Air Force Research Lab [58]. It does not use a state variable; the resistance change rate was defined directly. Later, models were implemented that describe a hyperbolic sine I-V relationship, like the one by the University of Michigan [59]. Previous models could only emulate sinusoidal input signal correctly, not general signal shapes. The hyperbolic sine model can approximate the I-V relationship of the MIM junction better.
For the implementation of a memristor model in LTspice we use the generalized hyperbolic sine model by Yakopcic et al. [55]. In comparison with the above-mentioned memristor models, the Yakopcic model uses more general equations for the memristor dynamics. Via fitting parameters, the model can be adapted to a variety of published memristor devices. The LTspice schematic of the sub-circuit can be seen in Figure 3. The I-V relationship of the generalized memristor model depends on the state variable x(t). Due to the MIM sandwich structure the memristive device shows an increase in conductivity beyond a certain voltage threshold. This behavior is considered by using a hyperbolic sine function. The I-V relationship of the generalized memristor model depends on the state variable x(t). Due to the MIM sandwich structure the memristive device shows an increase in conductivity beyond a certain voltage threshold. This behavior is considered by using a hyperbolic sine function.
where a 1 , a 2 and b are the parameters used to fit the model function to real memristive devices. The state variable x(t) is related to the resistance of the device, i.e., the relation of the length of the doped region to the overall length l d /l mem . In the LTspice sub-circuit, the state variable is named XSV. For the LTspice netlist Equation (6) is written as [55] .func IVRel(V1,V2) = IF(V1 >= 0, a1*V2*sinh(b*V1), a2*V2*sinh(b*V1)) Here V1 stands for the voltage across the memristor v(t) and V2 stands for the state variable x(t). The state variable change is described in (7) by two functions f (x(t)) and g(v(t)) : where η is a factor which can be +1 or −1. It represents the direction of the motion of the state variable relative to the polarity of v(t).
In order to get the derivative of the state variable, one can use a 1F capacitor in LTspice (Figure 3). Then the derivative of the node voltage V(XSV) is proportional to the current i x (t) in the capacitor which is equal to the current of the voltage dependent source Gx. Eventually, we can write: The initial state x 0 can be thought of as the voltage across the capacitor v C (t = 0) [55]: In order to change the state variable, a certain voltage threshold must be exceeded. The function g(v(t)) accounts for this behavior: where V p and V n are the positive and negative voltage thresholds; A p and A n are multipliers for state variable motion intensity in positive and negative direction. The netlist implementation is as follows [55]: It is harder to change the state variable as it approaches the boundaries. This behavior is modeled by the second function f (x(t)), where x p and x n are the points at which the state variable motion is reduced; α p and α n are the rate in positive and negative direction, at which the state variable motion decays. To assure that f (x(t)) approaches zero when x(t) = 1, this function includes two window functions: These functions are written into the sub-circuit netlist as [55] . To complete the sub-circuit netlist, the voltage dependent current sources Gx and Gm get the values [55] Gx 0 XSV value = {eta*F(V(P,N),V(XSV,0))*G(V(P,N))} Gm P N value = {IVRel(V(P,N),V(XSV,0))} Next, we describe the generation of input and output signals for the behavioral modeling of the memristor-based rectifier bridge.

Smoothing Output Signals of the Memristor-Based Rectifier Bridge
To construct polynomial models, the set of the test input signals of the rectifier and the corresponding set of the output signals are formed. The harmonic input signal is written where A is an amplitude; ω 0 = 2π · f 0 is an angular frequency; f 0 is a frequency; t ∈ [0, T) is time interval on a signal period; T is a signal period. The set of the test input signals is generated on specifying the following vectors of amplitude and frequency: A = [12; 13; 14; 15; 16] V, F 0 = [1; 2; 3; 4; 5] Hz, where A and F 0 are the amplitude and frequency vectors of the harmonic test exciters. The set of the corresponding output signals results from simulating the circuit shown in Figure 4. To build the mathematical models (a split polynomial), the number N p = 2001·5 = 10, 005 of the input-output pairs is used. We use high amplitudes (9) of the input signal (8), since the rectifier bridge model based on the Yakopcic model of a memristor is severely nonlinear and the Yakopcic model requires high nonlinearity. With the given set of parameters this nonlinearity is achieved at high input amplitudes. One of our targets is to compare the split multi-dimensional polynomial of the rectifier bridge with the model based on the Yakopcic memristor model. It is well known that the memristor resistance can change for smaller voltages, leading to the so-called voltage-time dilemma and the resistance may degrade even at zero voltage [60]. It should be pointed out that the investigated "black box" approach is independent of the mentioned problem, since a behavioral model does not take the memristor model into account.
The input and output signals are sampled and normalized according to equation time-normalized range t n ∈ [0, 2π). The rectifier outputs v out (t) are also converted to the time-normalized signals v out t n .
=     0 1; 2; 3; 4; 5 F Hz, where A and F 0 are the amplitude and frequency vectors of the harmonic test exciters. The set of the corresponding output signals results from simulating the circuit shown in Figure 4. To build the mathematical models (a split polynomial), the number N p = 2001 • 5 = 10,005 of the input-output pairs is used.  We use high amplitudes (9) of the input signal (8), since the rectifier bridge model based on the Yakopcic model of a memristor is severely nonlinear and the Yakopcic model requires high nonlinearity. With the given set of parameters this nonlinearity is achieved at high input amplitudes. One of our targets is to compare the split multi-dimensional polynomial of the rectifier bridge with the model based on the Yakopcic memristor model. It is well known that the memristor resistance can change for smaller voltages, leading to the so-called voltage-time dilemma and the resistance may degrade even at zero voltage [60]. It should be pointed out that the investigated "black box" approach is independent of the mentioned problem, since a behavioral model does not take the memristor model into account.
The input and output signals are sampled and normalized according to equation The time-normalized output signals of the rectifier at the excitation (8), in which the amplitude and frequency are assigned by the combination of elements from vectors (9) and (10), are shown in Figure 4a. As seen from Figure 4a, the output signals of the rectifier simulated with the Yakopcic model of memristors are not smooth functions. Every output signal has a point similar to discontinuity. This discontinuity is known to be generated by high-frequency harmonic components. To smooth signals, a low-pass filter, connected at the output of the rectifier, is used. Let's apply an ideal low-pass filter. We simulate the filter action as follows. Each output signal of the rectifier is processed with the discrete Fourier transform of length N = 2001. We keep the first 15 harmonics and zero the remaining ones. Then, we carry out the inverse discrete Fourier transform. As a result, the rectifier outputs are smoothed out. Figure 4b shows signals smoothed in comparison with those which are depicted in Figure 4a. It should be pointed out that the spectrum of the rectifier output signals contains even harmonics, since the nonlinear operator of a full-wave rectifier is even. This nonlinear operator transforms the harmonic input signal (8) into a polyharmonic one comprising only even harmonics.
Besides normalizing the time variable, the amplitude of the input signals is normalized where v(t n ) is the discrete input signal with the normalized variables of time and amplitude, v(t n ) ∈ [0, 1]; t n ∈ [0, 2π). Moreover, the amplitude of the output signals is normalized where A is the amplitude of the harmonic input signal (8). The amplitude normalization is effective to compare the results obtained by different models. Figure 4c shows the smoothed signals with the amplitude-and time-normalized variables.

Results of the Rectifier Modeling Based on Multi-Dimensional Polynomials
The rectifier polynomial describes the mapping between the input signals set and the set of the smoothed rectifier output signals. The rectifier models are two-dimensional polynomials of different powers. Note there are the polynomials of the even-only degree, since the full-wave rectifier is described by an even operator. The two-dimensional polynomial of the 6th degree (the number of parameters is 15) is used as an example. The polynomial is of the form where C i 1 i 2 ... i m are the model parameters; L is the even power of the polynomial; m is an even degree; I 1 = I 2 = · · · = I m = 1 are constants that determine the dimension of the polynomial in the corresponding coordinates. As the analysis of expression (11) shows, the presented polynomial is two-dimensional for any degree, its basis functions are the samples of the input signal at the current and previous time instants. The model parameters result from solving the optimization problem: where v out,sm (t n ) is discrete amplitude and time-normalized smoothed rectifier output signal, C is a vector containing the parameters of the model (8).
For the 3rd trial input signal, as shown in Figure 5a, the blue curve is the response y 3 (t n ) of polynomial (11) of the 6th power, the red curve is the smoothed output signal v out,sm,3 (t n ), the black curve is the uniform error ∆ 3 (t n ) from Equation (12) in the range t n ∈ [0, 2π). Figure 5b shows correspondingly the response of the split polynomial, the non-smoothed load voltage v out,3 (t n ), the uniform error obtained without smoothing by the low-pass filter. By comparing Figure 5a with Figure 5b, it can be seen that the non-smoothed output signal is described with a significantly greater error (black curve) than the smoothed output signal.
The maximum absolute and root-mean-square errors of modeling are represented in Table 1. For the 3rd trial input signal, as shown in Figure 5a, the blue curve is the response y 3 (t ̅ n ) of polynomial (11) of the 6th power, the red curve is the smoothed output signal v out, sm,3 (t ̅ n ), the black curve is the uniform error ∆ 3 (t ̅ n ) from Equation (12) in the range t ̅ n ∈[0,2π). Figure 5b shows correspondingly the response of the split polynomial, the nonsmoothed load voltage v out,3 (t ̅ n ), the uniform error obtained without smoothing by the low-pass filter. By comparing Figure 5a with Figure 5b, it can be seen that the nonsmoothed output signal is described with a significantly greater error (black curve) than the smoothed output signal.
(a) (b) Figure 5. For the 3rd trial input signal, the normalized output signals (the red curve is the load voltage of the circuit, the blue curve is the response of the split polynomial) and the uniform error (the black curve) in the case of smoothing by the low-pass filter (a) and without smoothing (b).
As follows from the analysis of Table 1 and Figure 5a,b, the higher the polynomial power, the higher the modeling accuracy. However, high-power polynomials are known to lead to ill-conditioned solutions when approximating a nonlinear operator; therefore, high-power polynomials are not constructive in practice. In this regard, we turn to piecewise polynomial models to achieve high modeling accuracy using the set of the lowpower polynomials.

Results of the Rectifier Modeling Based on Piecewise Multi-Dimensional Polynomials
To construct a piecewise polynomial model of the rectifier, we divide the range of the normalized modulus of the input signal amplitude into several sub-ranges. The thresholds of the sub-ranges we denote as λ j−1 ; j = 1, 2, . . . , J; 0 < λ j−1 < 1; λ 0 = 0; λ J = 1. The input sub-signals correspond to the mentioned amplitude sub-ranges, and every input signal is described by the sum where v (j) (t n ) is the input sub-signal with the normalized amplitude and time variables, As a result, we form the sub-sets from the input sub-signals v (j) (t n ), t n ∈ [0, 2π), j = 1, 2, . . . , J. The sub-sets of the output sub-signals v (j) out,sm t n , t n ∈ [0, 2π), j = 1, 2, . . . , J correspond to the sub-sets of the input sub-signals.
The block-scheme of dividing the input signal into sub-signals according to the amplitudes and the formation of the response of the split piecewise-polynomial model is shown in Figure 6. In the root-mean-square norm, using the sub-sets of the input sub-signals v j (t ̅ n ), t ̅ n ∈[0, 2π), j = 1, 2, …, J and the output sub-signals v out,sm j (t ̅ n ), we build models: In the root-mean-square norm, using the sub-sets of the input sub-signals v (j) (t n ), t n ∈ [0, 2π), j = 1, 2, . . . , J. and the output sub-signals v (j) out,sm t n , we build models: (13) where y (j) pw (t n ) and C (j) i 1 i 2 ... i m ; j = 1, 2, . . . , J are the output sub-signals and parameters of the piecewise model on the sub-sets of the input signal respectively; L is the even power of the polynomial; m is an even degree.
The output signal of the piecewise model results from the following summation where y (j) pw (t n ) is described in (13). Table 2 demonstrates the errors of modeling the memristive rectifier with the split piecewise-polynomial of the 6th power at three trial excitations. We used three subranges of the normalized amplitude of the input signal, the thresholds of sub-ranges are λ 0 = 0, λ 1 = 0.8, λ 2 = 0.9, λ 3 = 1. The total number of the model parameters is 45. As shown in Figure 7, the blue curve is the response of the polynomial (14) of the 6th power, the red curve is the smoothed output signal v out,sm,3 (t n ) of the circuit depicted in Figure 2, the black curve is the uniform error ∆ 3 (t n ). As shown in Figure 7, the blue curve is the response of the polynomial (14) of the 6th power, the red curve is the smoothed output signal v out,sm,3 (t ̅ n ) of the circuit depicted in Figure 2, the black curve is the uniform error ∆ 3 (t n ).

Figure 7.
For the 3rd trial input signal, the normalized output signals (the red curve is the load voltage of the circuit, the blue curve is the response of the split piecewise-polynomial) and the uniform error (the black curve). Tables 1 and 2, as well as in Figures 5a and 7 results in the fact, that the piecewise-polynomial model ensures lower errors than a general polynomial model does.

Conclusions
Following previous research, we demonstrate the capabilities of the behavioral model in the form of the multi-dimensional split polynomial for the memristor-based rec- Figure 7. For the 3rd trial input signal, the normalized output signals (the red curve is the load voltage of the circuit, the blue curve is the response of the split piecewise-polynomial) and the uniform error (the black curve). Tables 1 and 2, as well as in Figures 5a and 7 results in the fact, that the piecewise-polynomial model ensures lower errors than a general polynomial model does.

Conclusions
Following previous research, we demonstrate the capabilities of the behavioral model in the form of the multi-dimensional split polynomial for the memristor-based rectifier bridge. This model describes the nonlinear input-output mapping of the device. The behavioral model is built with the help of the sets of the input and output signals which are calculated by the LTspice software. The rectifier bridge is represented as an electrical circuit with memristive elements described by the Yakopcic model at a harmonic excitation. For smoothing the output signals of the rectifier simulated in LTspice, an ideal low-pass filter is used. Then, both amplitude and time samples of the signals from the obtained input and output sets are normalized.
The behavioral model is represented as a multi-dimensional polynomial of split signals. The scalar input signal is converted to a vector by splitting. Further, the vector signal elements are used as base functions for constructing the multi-dimensional polynomial. The model parameters result from solving the approximation problem of a nonlinear rectifier operator in the mean-square norm. The advantages of the multi-dimensional split polynomial are as follows: 1.
The form of the model and the method of its construction are universal, because they do not depend on the technology of memristor implementation.

2.
Since a polynomial is linear-in-parameters, these parameters are defined as globally optimal by solving the approximation problem.

3.
The splitting property allows the polynomial to be adapted to the assigned signal class, hence to construct a simpler model than other behavioral models.
To improve the accuracy of the modeling, the power of the split multi-dimensional polynomial should be increased. However, this operation leads to the problem of poor condition of the approximation problem when finding the model parameters. Therefore, the split piecewise-polynomial model free of this problem is proposed. Moreover, the piecewise-polynomial model is useful under the switching variability in memristors. As a result, we observe a significant reduction in the error of the rectifier bridge modeling.
The memristor-based rectifier bridge relates to, so-called, imperfect uncertain dynamic systems [61,62], since its working conditions affect the device parameters, the threshold and switching effects, robustness, etc. and, eventually, the device behavior. Thus, we face the problem of control taking into account the device limits. This problem can be solved in the framework of the control paradigm of imperfect and uncertain dynamic systems to gain a satisfactory trade-off between behavior and performance of the device [61,62].