Properties of Fractional-Order Magnetic Coupling

This paper presents the properties of fractional-order magnetic coupling. The difficulties connected with the analysis of two coils in dynamic states, resulting from the classical approach, provided motivation for studying the properties of fractional-order magnetic coupling. These difficulties arise from failure to comply with the commutation laws, i.e., a sudden power disappearance in the primary winding caused by a switch-mode power supply. Theoretically, under ideal conditions, a sudden power disappearance in the coil is, according to the classical method, manifested by a sudden voltage surge in the form of the Dirac delta function. As is well-known, it is difficult to obtain such ideal conditions in practice; the time of current disappearance does not equal zero due to the circuit breaker’s imperfection (even when electronic circuit breakers are used, the time equals several hundred nanoseconds). Furthermore, it is necessary to take into account phenomena occurring in real inductances, such as the skin effect, the influence of the ferromagnetic core and many other factors. It would be very difficult to model all these phenomena using classical differential calculus. The application of fractional-order differential calculus makes it possible to model them in a simple way by appropriate selection of coefficients and fractional-order derivatives. It should be mentioned that the analysis could be used, for example, in the case of high-voltage generation systems, including spark ignition systems of internal combustion engines. The use of fractional-order differential calculus will allow for more accurate modeling of phenomena occurring in such systems.


Introduction
The classical (standard) approach to the analysis and modeling of electrical systems usually ignores the effects of the non-ideality of elements and allows one to obtain a mathematical model based on integer-order differential equations. The use of this approach does not always lead to precise models. The application of fractional-order derivatives makes it possible to "compensate" for the missed phenomena. Real elements are characterized by loss, dielectric relaxation, electrochemical processes, magnetic heterogeneity, non-linearities, or simply operating conditions, e.g., temperature.
Classical differential and integral calculus is a well-studied area of mathematics, and its application can be found in numerous publications on the analysis of differentiation and integration, solving ordinary and partial differential equations, integral equations, etc. If fractional orders, e.g., 0.5, 1.8, or π, for derivatives and integrals are considered, we deal with differential-integral calculus of non-integer or fractional order. As it turns out, differential-integral calculus of fractional order can also be used with integer orders [1][2][3]. Derivatives of fractional order, also called differintegrals, are an interesting also be used with integer orders [1][2][3]. Derivatives of fractional order, also called differintegrals, are an interesting generalization of well-known dynamic systems, described by ʺclassicalʺ differential equations with derivatives of integer orders.
The use of fractional calculus in theoretical electrical engineering and circuit theory has enjoyed growing interest for several years. The appearance of entirely new electrotechnical components and devices, namely, supercapacitors, memory resistors called ʺmemristorsʺ [4][5][6], inductors with a skin effect [7], and coils with fractional mutual inductance [8], has been a reason for introducing C and L elements, which are described by differential-integral calculus of fractional order as a generalization of classical elements [9][10][11][12][13][14]. Descriptions employing differential calculus of fractional order have been extended to include phenomena occurring in real inductances. This comprises real coils (loss coils) with the skin effect, especially those with soft ferromagnetic cores, or containing electronic active systems, e.g., a generalized impedance converter (GIC) [15] and a coil with a fractional mutual inductance [16]. Other applications of fractional-order elements for materials and solutions described by chemical reactions are presented in [17,18].
The difficulties connected with the analysis of two coils in dynamic states, resulting from the classical approach, provided motivation for studying the properties of fractional-order magnetic coupling. It should be mentioned that the analysis can be used, for example, in the case of highvoltage generation systems, including spark ignition systems of internal combustion engines. The use of fractional-order differential calculus will allow for more accurate modeling of phenomena occurring in such circuits.
The paper is organized as follows. In the Introduction the problem is discussed, and Section 2 deals with the analysis of classical magnetic coupling of two coils and the response of such a system to voltage forcing in the form of a rectangular pulse. The next section contains a synthetic presentation of the mathematical apparatus used in the analysis. The definitions of selected derivatives of fractional order and their Laplace transformations are presented. The basic methods used to determine the inverse Laplace transform of fractional-order systems allow the approximation of the s  factor by the quotient of polynomials with integer powers. The main part of the paper is covered by Section 4, presenting the analysis of fractional-order magnetic coupling of two coils. This section includes the results of theoretical analysis supported by the results of digital simulations. It can be concluded that open-core ferromagnetic coils, which are non-linear elements with ambiguous characteristics, can be more realistically modeled using fractional-order calculus. Experimental tests on a real object (an ignition coil), confirming the purposefulness of the analysis, are presented in Section 5.

The Classic Magnetic Coupling of Two Coils
To compare the magnetic couplings of two coils with matching and opposite magnetic polarity, generally known relationships for classical magnetic coupling are presented. The object studied is the electrical circuit of two coils shown in Figure 1. As we know, by eliminating coupling and using the Laplace transform, it is possible to obtain the diagrams of circuits shown in Figure 2. As we know, by eliminating coupling and using the Laplace transform, it is possible to obtain the diagrams of circuits shown in Figure 2. If the systems in Figure 2 are considered to be two-port networks of T type, the following chain matrix can be determined: Therefore, for the two-port networks in Figures 2a and 2b, the matrix has the form: Assuming that two-port networks are in idle state, the following equation applies: Analyzing the relations, in the context of the elements of the matrices Aa and Ab, it is easy to conclude that the output voltage and the input current of two-port networks a) and b) will differ only by sign.
For the voltage forcing in the form of a rectangular pulse (Figure 3), voltage u10(t) is: If the systems in Figure 2 are considered to be two-port networks of T type, the following chain matrix can be determined: Therefore, for the two-port networks in Figure 2a,b, the matrix has the form: Assuming that two-port networks are in idle state, the following equation applies: Thus, Considering the relations in (2), one obtains: and : Analyzing the relations, in the context of the elements of the matrices A a and A b , it is easy to conclude that the output voltage and the input current of two-port networks a) and b) will differ only by sign.
Energies 2020, 13, x FOR PEER REVIEW 4 of 16 1 1 (7) where U: supply voltage; t0: pulse duration; 1(t): unit function.   The voltage transform is written out as: Thus, the voltage transform at the output and the current transform at the two-terminal input are expressed by the relations: By calculating the inverse transform for time t < t 0 , the system response to rectangular-pulse forcing was obtained. The voltage at the output and current at the input of the two-port network are shown in Figures 4 and 5. Energies 2020, 13, x FOR PEER REVIEW 4 of 16 1 1 (7) where U: supply voltage; t0: pulse duration; 1(t): unit function.   The waveforms were obtained for the real parameters of ignition coils, that is R1 = 0.2 ; L1 = 7 mH; R2 = 1380 ; L2 = 40 mH; U = 13 V [19][20][21][22]. When considering the system response to the forcing assumed for time 0 t t  , the initial condition should be taken into account-i.e., currents passing through coils at 0 t t  . Then, when analyzing the circuits in Figure 2a or 2b, it turns out that only the current passing through the L1 coil should be taken into account, since I2 = 0. Hence, for the case of a), the diagram shown in Figure 6 was obtained. By calculating the voltage component at the two-terminal output, one obtains: By determining the inverse transform of the first component of Equation (11), one obtains the waveforms that are the same as those shown in Figures 4 and 5, but shifted in time by t0 and inverse. The transform of the second component, resulting from the initial condition (it was assumed that , can be written as: As can be seen, one of the components on the right-hand side of Equation (12) is constant, i.e., independent of s, which means that its inverse transform will be the Dirac delta shifted in time by t0.
The result of the digital simulation of the solution is shown in Figure 7. The waveforms were obtained for the real parameters of ignition coils, that is When considering the system response to the forcing assumed for time t ≥ t 0 , the initial condition should be taken into account-i.e., currents passing through coils at t = t 0 . Then, when analyzing the circuits in Figure 2a or Figure 2b, it turns out that only the current passing through the L 1 coil should be taken into account, since I 2 = 0. Hence, for the case of Figure 2a, the diagram shown in Figure 6 was obtained. The waveforms were obtained for the real parameters of ignition coils, that is R1 = 0.2 ; L1 = 7 mH; R2 = 1380 ; L2 = 40 mH; U = 13 V [19][20][21][22]. When considering the system response to the forcing assumed for time 0 t t  , the initial condition should be taken into account-i.e., currents passing through coils at 0 t t  . Then, when analyzing the circuits in Figure 2a or 2b, it turns out that only the current passing through the L1 coil should be taken into account, since I2 = 0. Hence, for the case of a), the diagram shown in Figure 6 was obtained. By calculating the voltage component at the two-terminal output, one obtains: By determining the inverse transform of the first component of Equation (11), one obtains the waveforms that are the same as those shown in Figures 4 and 5, but shifted in time by t0 and inverse. The transform of the second component, resulting from the initial condition (it was assumed that , can be written as: As can be seen, one of the components on the right-hand side of Equation (12) is constant, i.e., independent of s, which means that its inverse transform will be the Dirac delta shifted in time by t0.
The result of the digital simulation of the solution is shown in Figure 7. . By calculating the voltage component at the two-terminal output, one obtains: By determining the inverse transform of the first component of Equation (11), one obtains the waveforms that are the same as those shown in Figures 4 and 5, but shifted in time by t 0 and inverse. The transform of the second component, resulting from the initial condition (it was assumed that I 1 (t 0 ) = U R 1 ), can be written as: As can be seen, one of the components on the right-hand side of Equation (12) is constant, i.e., independent of s, which means that its inverse transform will be the Dirac delta shifted in time by t 0 .
Energies 2020, 13, 1539 The result of the digital simulation of the solution is shown in Figure 7. An electrical circuit with such forcing cannot be considered a classical one because the condition of current continuity in the coil is not met. Therefore, fractional-order calculus should be used to analyze it. Selected elements of the mathematical apparatus used in fractional-order calculus are presented in the next section.

Selected Elements of Fractional-Order Calculus
This section presents definitions of fractional-order derivatives in the Caputo and Riemann-Liouville senses, and the corresponding Laplace transforms. The approximation methods used to determine the inverse Laplace transform of fractional-order systems of equations, including the continued fraction expansion (CFE) and Oustaloup methods, are discussed.

Definitions of Fractional-Order Derivatives
There are many definitions of fractional-order integrals and derivatives. One of them is the Riemann-Liouville differintegral. Only the integral on the left-hand side, mainly used in engineering applications, is applied by the authors. The definition of derivative is obtained from the definition of the Riemann-Liouville fractional integral, which is a generalization of the formula for the n-th-order integral, and is also based on a well-known formula of classical differential calculus, namely: where function f(x) is n-times differentiable in the range [a, b]; k, n ∈N; and k > n. Definition 1. The left-hand side Riemann-Liouville differintegral of the order α∈C, (Re{ α } > 0) of the real function f(x) is defined as: Detailed descriptions of the Riemann-Liouville differintegral can be found, among others, in [3,[23][24][25].
Introduced by Michael Caputo in 1967 in [26], the definition of the Riemann-Liouville fractional integral is also used. Comparison of the Riemann-Liouville derivative and Caputo's (C) derivative shows that, in the latter, the order of integration and differentiation is inverse, namely: An electrical circuit with such forcing cannot be considered a classical one because the condition of current continuity in the coil is not met. Therefore, fractional-order calculus should be used to analyze it. Selected elements of the mathematical apparatus used in fractional-order calculus are presented in the next section.

Selected Elements of Fractional-Order Calculus
This section presents definitions of fractional-order derivatives in the Caputo and Riemann-Liouville senses, and the corresponding Laplace transforms. The approximation methods used to determine the inverse Laplace transform of fractional-order systems of equations, including the continued fraction expansion (CFE) and Oustaloup methods, are discussed.

Definitions of Fractional-Order Derivatives
There are many definitions of fractional-order integrals and derivatives. One of them is the Riemann-Liouville differintegral. Only the integral on the left-hand side, mainly used in engineering applications, is applied by the authors. The definition of derivative is obtained from the definition of the Riemann-Liouville fractional integral, which is a generalization of the formula for the n-th-order integral, and is also based on a well-known formula of classical differential calculus, namely: (14) where function f (x) is n-times differentiable in the range [a, b]; k, n ∈N; and k > n.

Definition 1.
The left-hand side Riemann-Liouville differintegral of the order α ∈ C, (Re{α} > 0) of the real function f(x) is defined as: Detailed descriptions of the Riemann-Liouville differintegral can be found, among others, in [3,[23][24][25]. Introduced by Michael Caputo in 1967 in [26], the definition of the Riemann-Liouville fractional integral is also used. Comparison of the Riemann-Liouville derivative and Caputo's (C) derivative shows that, in the latter, the order of integration and differentiation is inverse, namely: This affects the final structure of the derivative defined as follows: In the Caputo sense, the derivative of the order α (α ∈ R + )is called the operator, The definitions of the fractional-order derivatives are not equivalent in a general case. However, under certain assumptions, there are similarities between them, including the problems with interpretation of initial conditions. An important feature of the Caputo derivative is that initial conditions include derivatives of integer order. Another difference between the Caputo and the Riemann-Liouville definitions is the derivative of the constant function. For the lower limit, α = −∞, the Caputo and Riemann-Liouville definitions are equivalent.

The Laplace Transform of Fractional-Order Derivatives
The Laplace transform is frequently used to solve engineering problems. Its definition, similarly to that of a fractional derivative, is a generalization of concepts from classical calculus. The relations for the Laplace transform and its inverse are well-known and not discussed here. The Laplace transform of fractional-order derivatives can be calculated using the relations for the transform of the n-th-order classical derivative: Furthermore, for the fractional-order derivatives in the Riemann-Liouville sense, The practical application of the transform is limited due to the lack of physical interpretation of the initial values of fractional derivatives for the lower limit t = 0. Likewise, a formula for a fractional derivative in the Caputo sense can be obtained: Its practical application is much wider because the initial conditions for t = 0 refer to derivatives of integer orders whose physical interpretations are well-known.
For derivatives of fractional order, the Laplace transform assuming zero initial conditions is: where F(s) is the transform of the Laplace function f (t).
The detailed properties of the Laplace transform are presented in [3,27].

Selected Approximation Methods in the Analysis of Fractional-Order Electrical Circuits
The methods most commonly used to determine the inverse Laplace transform of fractional order circuits are those that allow the approximation of the s α factor by the quotient of polynomials with integer powers s. Continued fraction expansion (CFE) is such a method. It consists of the expansion of the expression (1 + x), 0 ≤ α ≤ 1 into an infinite continued fraction. By substituting x = s − 1 and taking into account successive terms, the s α approximations of the required order and accuracy are obtained.
In the CFE method, the s α factor can be presented as the quotient of the polynomials of the s variable and the order of the α derivative-here, these variables have integer powers [28,29].
where A: the order of approximation; P Ak (α) and Q Ak (α): α polynomials of A order. It is shown in [30] that approximation of the order of A = 5 gives relatively the best accuracy for fifth-degree polynomials. The higher orders of approximation result in an increased number of components in polynomials and a slight improvement in accuracy. This means that higher orders of approximation cause, apart from an increased number of components in polynomials and, therefore, a greater number of poles and a more complex numerical implementation, a relatively minor change in the solution, which is interpreted as a slight improvement in accuracy. Thus, s α = P 50 s 5 + P 51 s 4 + P 52 s 3 + P 53 s 2 + P 54 s + P 55 Q 50 s 5 + Q 51 s 4 + Q 52 s 3 + Q 53 s 2 + Q 54 s + Q 55 (23) where: The second method used to approximate the s α factor is the expansion into the Oustaloup series, where for frequency range (ω l ,ω h ), it is shown as transfer function H(s) [31][32][33][34]: where N: the order of approximation; C 0 : the gain described by the following relation: where ω ' k : zeros and ω k : the poles of the H(s) transfer function, which are: It should be stressed that since the polynomials in the above approximations contain an integer of s, the inverse transform can be determined by standard methods [35,36]. The following equations can be written out for the model:

Fractional-Order Magnetic Coupling of Two Coils
after the Laplace transformation with zero initial conditions, one obtains:

( 3 0 )
After transformation, the equations assume the form:

( 3 1 )
After ordering, the following relation is obtained:

( 3 2 )
This system of equations corresponds to the electrical circuit in Figure 9. It should be noted that in [8], the same system was obtained with the use of a different algorithm. If we assume that the circuit in Figure 9 is a T-type two-port network, then on the basis of The following equations can be written out for the model: after the Laplace transformation with zero initial conditions, one obtains: After transformation, the equations assume the form: After ordering, the following relation is obtained: This system of equations corresponds to the electrical circuit in Figure 9. The following equations can be written out for the model: after the Laplace transformation with zero initial conditions, one obtains:

( 3 0 )
After transformation, the equations assume the form: After ordering, the following relation is obtained:

( 3 2 )
This system of equations corresponds to the electrical circuit in Figure 9. It should be noted that in [8], the same system was obtained with the use of a different algorithm. If we assume that the circuit in Figure 9 is a T-type two-port network, then on the basis of Equation (1), its chain matrix has the form: Figure 9. Circuit with magnetic coupling eliminated.
Energies 2020, 13,1539 It should be noted that in [8], the same system was obtained with the use of a different algorithm. If we assume that the circuit in Figure 9 is a T-type two-port network, then on the basis of Equation (1), its chain matrix has the form: Following the same procedure as in the case of classical coupling, the relations determining output voltage and input current transforms at the input of the two-port network in Figure 9 take the form: Next, if the rectangular-pulse forcing shown in Figure 3 is considered, the relations are: As can be seen, the expression for the current at the input of the two-port network has the same form as classical coupling.
After substituting the s α factor by the quotient of polynomials according to the CFE method (22) or the Oustaloup series (24), the relation has the form: By calculating the inverse transform for the time t < t 0 , voltage curves at the output of a two-port network with rectangular-pulse forcing were obtained for α = 0.95, 0.9, 0.8, 0.75, as shown in Figure 10. Following the same procedure as in the case of classical coupling, the relations determining output voltage and input current transforms at the input of the two-port network in Figure 9 take the form: Next, if the rectangular-pulse forcing shown in Figure 3 is considered, the relations are: As can be seen, the expression for the current at the input of the two-port network has the same form as classical coupling.
After substituting the s α factor by the quotient of polynomials according to the CFE method (22) or the Oustaloup series (24), the relation has the form: By calculating the inverse transform for the time t < 0 t , voltage curves at the output of a twoport network with rectangular-pulse forcing were obtained for  = 0.95, 0.9, 0.8, 0.75, as shown in Figure 10.    After considering the circuit response to rectangular-pulse forcing in time 0 t t  and taking into account the initial condition, i.e., the diagram presented in Figure 6, the following relations are written out: The expression determining the current at the input of the two-port network does not contain the s  factor, so the reverse transform can be determined directly: ( 3 9 ) Figure 12 shows the current curve at the input of the two-port network for the whole time interval.  Figure 11. Voltage curves at the output of a two-port network with rectangular-pulse forcing for α = 0.5, 0.4, 0.3.
After considering the circuit response to rectangular-pulse forcing in time t ≥ t 0 and taking into account the initial condition, i.e., the diagram presented in Figure 6, the following relations are written out: The expression determining the current at the input of the two-port network does not contain the s α factor, so the reverse transform can be determined directly: Figure 12 shows the current curve at the input of the two-port network for the whole time interval.
Energies 2020, 13, x FOR PEER REVIEW 11 of 16 Figure 11. Voltage curves at the output of a two-port network with rectangular-pulse forcing for 0.5, 0.4, 0.3.
After considering the circuit response to rectangular-pulse forcing in time 0 t t  and taking into account the initial condition, i.e., the diagram presented in Figure 6, the following relations are written out: The expression determining the current at the input of the two-port network does not contain the s  factor, so the reverse transform can be determined directly: ( 3 9 ) Figure 12 shows the current curve at the input of the two-port network for the whole time interval. The voltage transform at the output of the two-port network is: By calculating the inverse transform for the time t ≥ t 0 , voltage waveforms for the assumed forcing at the output of the two-port network are obtained. Figure 13 shows voltage curves for α = 0.95, 0.9, 0.8, 0.75, and Figure 14 for α = 0.5, 0.4, 0.3.
Energies 2020, 13, x FOR PEER REVIEW 12 of 16 By calculating the inverse transform for the time t ≥t 0 , voltage waveforms for the assumed forcing at the output of the two-port network are obtained. Figure 13 shows voltage curves for = α0.95, 0.9, 0.8, 0.7  The figures presented in the paper show that the simulation results are consistent with those given in the references.

Experimental Research on a Real Object
The digital simulations were verified experimentally by means of an ignition coil used to generate high voltages in typical ignition systems of spark-ignition engines. The tested coil was a circuit of two induction coils with an open ferromagnetic core of matching magnetization, as shown in Figure 8. The experimental tests consisted in recording the u2(t) voltage at the system output for the input voltage in the form of a rectangular pulse. The recorded curves of u2(t) voltage and i1(t) current are shown in Figure 15. By calculating the inverse transform for the time t ≥t 0 , voltage waveforms for the assumed forcing at the output of the two-port network are obtained. Figure 13 shows voltage curves for = α0.95, 0.9, 0.8, 0.7  The figures presented in the paper show that the simulation results are consistent with those given in the references.

Experimental Research on a Real Object
The digital simulations were verified experimentally by means of an ignition coil used to generate high voltages in typical ignition systems of spark-ignition engines. The tested coil was a circuit of two induction coils with an open ferromagnetic core of matching magnetization, as shown in Figure 8. The experimental tests consisted in recording the u2(t) voltage at the system output for the input voltage in the form of a rectangular pulse. The recorded curves of u2(t) voltage and i1(t) current are shown in Figure 15. The figures presented in the paper show that the simulation results are consistent with those given in the references.

Experimental Research on a Real Object
The digital simulations were verified experimentally by means of an ignition coil used to generate high voltages in typical ignition systems of spark-ignition engines. The tested coil was a circuit of two induction coils with an open ferromagnetic core of matching magnetization, as shown in Figure 8. The experimental tests consisted in recording the u 2 (t) voltage at the system output for the input voltage in the form of a rectangular pulse. The recorded curves of u 2 (t) voltage and i 1 (t) current are shown in Figure 15. Energies 2020, 13, x FOR PEER REVIEW 13 of 16 Figure 15. Oscillograms of the voltage at the output of a two-port network and the current at the input with rectangular-pulse forcing.
As can be seen, the curves of i1(t) current and u2(t) voltage obtained as a result of the experiments contain oscillations and are insignificantly different from the results of the numerical experiments shown in Figures 13 and 14. In order to compare the experimental results with the simulations, a detailed curve of u2(t) voltage at the moment of its induction for the fading impulse is presented in Figure 16. The comparison of the obtained curve with the results of numerical simulations shows that the curves are similar in shape. The differences between the real and simulated curves result from the existence of micro-capacities in the tested circuit (e.g., between coils) and the capacity of the voltage probe used. Complementing the model of the circuit in Figure 8 with capacitive and resistive elements and then selecting α to make the model as close as possible to the real curves will allow for an accurate mathematical description of the analyzed circuit. As can be seen, the curves of i 1 (t) current and u 2 (t) voltage obtained as a result of the experiments contain oscillations and are insignificantly different from the results of the numerical experiments shown in Figures 13 and 14. In order to compare the experimental results with the simulations, a detailed curve of u 2 (t) voltage at the moment of its induction for the fading impulse is presented in Figure 16. As can be seen, the curves of i1(t) current and u2(t) voltage obtained as a result of the experiments contain oscillations and are insignificantly different from the results of the numerical experiments shown in Figures 13 and 14. In order to compare the experimental results with the simulations, a detailed curve of u2(t) voltage at the moment of its induction for the fading impulse is presented in Figure 16. The comparison of the obtained curve with the results of numerical simulations shows that the curves are similar in shape. The differences between the real and simulated curves result from the existence of micro-capacities in the tested circuit (e.g., between coils) and the capacity of the voltage probe used. Complementing the model of the circuit in Figure 8 with capacitive and resistive elements and then selecting α to make the model as close as possible to the real curves will allow for an accurate mathematical description of the analyzed circuit. The comparison of the obtained curve with the results of numerical simulations shows that the curves are similar in shape. The differences between the real and simulated curves result from the existence of micro-capacities in the tested circuit (e.g., between coils) and the capacity of the voltage probe used. Complementing the model of the circuit in Figure 8 with capacitive and resistive elements and then selecting α to make the model as close as possible to the real curves will allow for an accurate mathematical description of the analyzed circuit.

Conclusions
To sum up, it can be stated that the use of fractional-order magnetic coupling between two coils makes it possible to avoid the Dirac impulse that is used in the classical approach. In the case of first-order magnetic coupling (classical approach), the immediate appearance and immediate disappearance of voltage are assumed (Figure 3). In practice, no mechanical or even electronic switch is able to realize this state. The use of fractional-order magnetic coupling makes it possible to consider the non-ideality of coils and switches. Thus, results closer to those obtained experimentally can be achieved. In future studies, fractional-order magnetic coupling will be used to model the spark-ignition system of an internal combustion engine, and the results will be experimentally verified in a system with a standard ignition coil. As can be seen, the real curves presented in Figure 15 contain oscillations and are insignificantly different from the results of numerical experiments. The differences between the real and simulated curves result from the existence of micro-capacities (e.g., between coils) in the tested circuit. The next stage of research will include the development of a model containing capacitive and resistive elements and the selection of α to make the developed model as close as possible to the real curves.
An additional statement resulting from the study may be formulated as follows: The current waveform at the input of the electrical circuit is the same as with classical coupling, and the coupling value α considerably affects the value of voltage induced in the second coil. This means that the higher the α value, the "stronger" the coupling.

Conflicts of Interest:
The authors declare no conflict of interest. Nomenclature L 1 and L 2 self-inductance of coils M mutual inductance of coils R 1 and R 2 coil resistance U supply voltage t 0 pulse duration 1(t) unit function F(s) the transform of the Laplace function f (t) A, N the order of approximation P Ak (α), Q Ak (α) α polynomials of A order C 0 the gain described by the relation ω k ' zeros ω k the poles H(s) transfer function