Next Article in Journal
A Novel Hybrid Method for KPI Anomaly Detection Based on VAE and SVDD
Next Article in Special Issue
Multi-Objective UAV Trajectory Planning in Uncertain Environment
Previous Article in Journal
Classification of Metaplectic Fusion Categories
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling RL Electrical Circuit by Multifactor Uncertain Differential Equation

1
Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China
2
School of Information and Statistics, Guangxi University of Finance and Economics, Nanning 530003, China
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(11), 2103; https://doi.org/10.3390/sym13112103
Submission received: 2 October 2021 / Revised: 28 October 2021 / Accepted: 3 November 2021 / Published: 5 November 2021
(This article belongs to the Special Issue Uncertainty Theory: Symmetry and Applications)

Abstract

:
The symmetry principle of circuit system shows that we can equate a complex structure in the circuit network to a simple circuit. Hence, this paper only considers a simple series RL circuit and first presents an uncertain RL circuit model based on multifactor uncertain differential equation by considering the external noise and internal noise in an actual electrical circuit system. Then, the solution of uncertain RL circuit equation and the inverse uncertainty distribution of solution are derived. Some applications of solution for uncertain RL circuit equation are also investigated. Finally, the method of moments is used to estimate the unknown parameters in uncertain RL circuit equation.

1. Introduction

For the sake of rationally handling the likelihood that something will happen, uncertainty theory was established by Liu [1] and then perfected by Liu [2]. To this day, uncertainty theory has been extensively studied and has spawned numerous theoretical branches.
Among the theoretical branches of uncertainty theory, an uncertain differential equation is a type of differential equation driven by Liu process. It was first presented by Liu [3] to model the dynamic systems with continuous-time noises. For the purpose of investigating the solution of uncertain differential equation, Chen and Liu [4] first verified the existence and uniqueness theorem of solution of uncertain differential equation under linear growth condition and Lipschitz condition. Later, the theorem was extended to local linear growth condition and local Lipschitz condition by Gao [5]. After that, the first concept of stability of uncertain differential equation was proposed by Liu [2]. Afterwards, Yao and Gao [6] proved some theorems to develop the stability analysis of uncertain differential equation, and then other types of stability were presented by Sheng and Wang [7], Yao et al. [8], Yang et al. [9], etc. As the most significant contribution, the Yao–Chen formula was shown by Yao and Chen [10] to associate uncertain differential equation with ordinary differential equations, which greatly promoted the development of uncertain differential equation. Based on the Yao–Chen formula, the numerical method of solving uncertain differential equation was first presented by Yao and Chen [10], which was then developed by Yang and Shen [11], Yang and Ralescu [12], Gao [13], etc. As a topic that has attracted much attention in recent years, how to estimate the unknown parameters in an uncertain differential equation that fits the observed data as much as possible is a core problem. In order to deal with this problem, Yao and Liu [14] first presented the method of moments to estimate the unknown parameters. After that, several methods were proposed such as generalized moment estimation (Liu [15]), minimum cover estimation (Yang et al. [16]), least squares estimation (Sheng et al. [17]) and maximum likelihood estimation (Liu and Liu [18]). Based on that, uncertain differential equation was used to model the cumulative number of COVID-19 infections in China by scholars such as Lio and Liu [19], Chen et al. [20] and Jia and Chen [21].
As an important research field of physics systems, electrical circuit is nearly everywhere in our lives. To name only a few, electrical circuit is an essential component of power supply network, car, television, and mobile phone. In the actual electrical circuit system, due to the electromagnetic interference of external environment and the dimension error of potential source device, there is always noise in the circuit system (Bennett [22]). In order to study the effect of noise on circuit system, white noise was first used to describe the classification and numerical simulation of electrical circuit by Kampowsky et al. [23]. Later, the noise was modeled as Wiener process and stochastic differential equation was applied in the modeling of electrical circuit by Penski [24], Kolarova [25], Kolarova [26], etc. However, Liu [27] proved that stochastic circuit differential equation would lead to a paradox and suggested that the noise should be modeled by uncertain differential equation.
In order to model the series RL circuit system more reasonably, this paper first proposes the uncertain RL circuit equation based on multifactor uncertain differential equation by considering the external noise and internal noise. The main contributions of this paper are to derive the solution of uncertain RL circuit equation and the inverse uncertainty distributions of solution, supremum of solution and time integral of solution. The overall structure of this paper takes the form of five sections, including this introductory section. Section 2 begins by deriving the uncertain RL circuit equation, and then gives the solution of uncertain RL circuit equation and the inverse uncertainty distribution of solution. Section 3 investigates two applications of solution of uncertain RL circuit equation. In Section 4, parameter estimation for uncertain RL circuit equation is considered and a numerical experiment is provided. Finally, a concise conclusion is given in Section 5, and the future scope of the study is provided in Section 6.

2. Uncertain RL Circuit Equation

Based on the symmetry principle of circuit system, we only consider a simple series RL circuit, which is an electrical circuit composed of a resistor and inductor and is driven by a potential source. The resistance of the resistor R is denoted by r and measured in ohms ( Ω ), and the inductance of the inductor L is denoted by l and measured in henrys (H). The current through the resistor and inductor at time t is denoted by I t and measured in amperes (A), the potential source provided to the circuit system is denoted by F ( t ) and measured in volts (V).

2.1. Uncertain RL Circuit Model

Next, we will derive the uncertain RL circuit model based on multifactor uncertain differential equation. Let the time interval [ 0 , t ] be divided into n subintervals such that
0 = t 0 < t 1 < < t n 1 < t n = t .
Then, we write
Δ = max 1 j n ( t j t j 1 ) .
For each index j with 1 j n , we consider the behavior of the circuit system in time interval ( t j 1 , t j ] . Due to the electromagnetic interference of external environment to the potential source and the dimension error of potential source device, the potential source provided to the circuit system is often affected by some noise, which is modeled by
σ 1 C 1 t j C 1 t j 1 t j t j 1 .
Here,
C 1 t j C 1 t j 1 t j t j 1
is a normal uncertain variable of mean zero and variance one, and σ 1 > 0 is known as the intensity of noise. Thus, the actual potential source provided to the circuit system is
F ( t j 1 ) + σ 1 C 1 t j C 1 t j 1 t j t j 1
when Δ is sufficiently small. When the current through the resistor, the resistor will generate heat, which will then disturb the resistance. Thus, the resistance of resistor can be regarded as
r + σ 2 C 2 t j C 2 t j 1 t j t j 1
where σ 2 > 0 is also known as the intensity of noise. By using Ohm’s law, the resistor voltage is
I t j 1 r + σ 2 C 2 t j C 2 t j 1 t j t j 1 .
According to the definition of inductance, the voltage across the inductor is
l I t j I t j 1 t j t j 1 .
It follows from Kirchhoff’s second law that
I t j 1 r + σ 2 C 2 t j C 2 t j 1 t j t j 1 + l I t j I t j 1 t j t j 1 = F ( t j 1 ) + σ 1 C 1 t j C 1 t j 1 t j t j 1 ,
which can be rewritten as
I t j I t j 1 = F ( t j 1 ) r I t j 1 l t j t j 1 + σ 1 l C 1 t j C 1 t j 1 σ 2 I t j l C 1 t j C 1 t j 1 .
Thus,
I t I 0 = j = 1 n I t j I t j 1 = j = 1 n F ( t j 1 ) r I t j 1 l t j t j 1 + σ 1 l C 1 t j C 1 t j 1 σ 2 I t j l C 1 t j C 1 t j 1 = j = 1 n F ( t j 1 ) r I t j 1 l t j t j 1 + j = 1 n σ 1 l C 1 t j C 1 t j 1 j = 1 n σ 2 I t j l C 1 t j C 1 t j 1 .
Letting Δ 0 , we have
I t I 0 = 0 t F ( s ) r I s l d s + 0 t σ 1 l d C 1 s 0 t σ 2 I s l d C 2 s .
This is equivalent to the following multifactor uncertain differential equation
d I t = F ( t ) l r l I t d t + σ 1 l d C 1 t σ 2 l I t d C 2 t
which is called the uncertain RL circuit equation.

2.2. Solution of Uncertain RL Circuit Equation

Next, we will give the solution for the uncertain RL circuit Equation (1).
Theorem 1.
Let F ( t ) be a continuous functions on ℜ. Then the uncertain RL circuit equation,
d I t = F ( t ) l r l I t d t + σ 1 l d C 1 t σ 2 l I t d C 2 t ,
has a solution
I t = I 0 exp r l t σ 2 l C 2 t + 1 l 0 t F ( s ) exp r l s t + σ 2 l C 2 s C 2 t d s + σ 1 l 0 t exp r l s t + σ 2 l C 2 s C 2 t d C 1 s
where I 0 is the initial current through the circuit system.
Proof of Theorem 1.
We first define two uncertain processes U t and V t via uncertain differential equations,
d U t = r l U t d t σ 2 l U t d C 2 t , d V t = F ( t ) l U t d t + σ 1 l U t d C 1 t .
Since
d ( U t V t ) = V t d U t + U t d V t = r l V t U t d t σ 2 l V t U t d C 2 t + F ( t ) l d t + σ 1 l d C 1 t = F ( t ) l r l V t U t d t + σ 1 l d C 1 t σ 2 l V t U t d C 2 t ,
the uncertain process U t V t is a solution of the uncertain RL circuit Equation (1). Note that
U t = U 0 exp 0 t r l d s 0 t σ 2 l d C 2 s = U 0 exp r l t σ 2 l C 2 t
and
V t = V 0 + 0 t F ( s ) l U s d s + 0 t σ 1 l U s d C 1 s = V 0 + 1 l U 0 0 t F ( s ) exp r l s + σ 2 l C 2 s d s + σ 1 l U 0 0 t exp r l s + σ 2 l C 2 s d C 1 s .
Hence, we have
U t V t = V 0 U 0 exp r l t σ 2 l C 2 t + 1 l 0 t F ( s ) exp r l s t + σ 2 l C 2 s C 2 t d s + σ 1 l 0 t exp r l s t + σ 2 l C 2 s C 2 t d C 1 s .
Taking U 0 = 1 and V 0 = I 0 , we can get the solution (2). The theorem is proved.

2.3. Inverse Uncertainty Distribution of Solution

We first introduce an important lemma before giving the inverse uncertainty distribution of solution (2).
Lemma 1.
Let I t be the solution of the uncertain RL circuit equation
d I t = F ( t ) l r l I t d t + σ 1 l d C 1 t σ 2 l I t d C 2 t .
Then, we have
M { I t Ψ t 1 ( α ) , t } = α
and
M { I t > Ψ t 1 ( α ) , t } = 1 α ,
where
Ψ t 1 ( α ) = I 0 + σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r exp σ 2 Φ 1 ( α ) r l t σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r + 1 l 0 t F ( s ) exp σ 2 Φ 1 ( α ) r l t s d s
and
Φ 1 ( α ) = 3 π ln α 1 α
is the inverse standard normal uncertainty distribution.
Proof of Lemma 1.
For any α ( 0 , 1 ) , we define some events as follows,
Λ 1 + = γ | d C 1 t ( γ ) d t Φ 1 ( α ) , t ,
Λ 1 = γ | d C 1 t ( γ ) d t > Φ 1 ( α ) , t ,
Λ 2 + = γ | d C 2 t ( γ ) d t < Φ 1 ( 1 α ) , t ,
Λ 2 = γ | d C 2 t ( γ ) d t Φ 1 ( 1 α ) , t .
Since C 1 t and C 2 t are two uncertain independent increment processes, we have
M { Λ 1 + } = M { Λ 2 } = α
and
M { Λ 1 } = M { Λ 2 + } = 1 α .
On the one hand, we have
M { I t Ψ t 1 ( α ) , t } = M I 0 exp r l t σ 2 l C 2 t + 1 l 0 t F ( s ) exp r l s t + σ 2 l C 2 s C 2 t d s + σ 1 l 0 t exp r l s t + σ 2 l C 2 s C 2 t d C 1 s I 0 + σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r exp σ 2 Φ 1 ( α ) r l t σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r
+ 1 l 0 t F ( s ) exp σ 2 Φ 1 ( α ) r l t s d s , t M C 2 t t Φ 1 ( 1 α ) , C 2 t C 2 s ( t s ) Φ 1 ( 1 α ) , d C 1 s ( γ ) d s Φ 1 ( α ) , 0 s t , t M γ | d C 2 s ( γ ) d s Φ 1 ( 1 α ) , d C 1 s ( γ ) d s Φ 1 ( α ) , 0 s t , t = M γ | d C 2 t ( γ ) d t Φ 1 ( 1 α ) , d C 1 t ( γ ) d t Φ 1 ( α ) , t = M γ | d C 1 t ( γ ) d t Φ 1 ( α ) , t   M   γ | d C 2 t ( γ ) d t Φ 1 ( 1 α ) , t = M { Λ 1 + }   M   { Λ 2 } = α .
On the other hand, we have
M { I t > Ψ t 1 ( α ) , t } = M I 0 exp r l t σ 2 l C 2 t + 1 l 0 t F ( s ) exp r l s t + σ 2 l C 2 s C 2 t d s + σ 1 l 0 t exp r l s t + σ 2 l C 2 s C 2 t d C 1 s > I 0 + σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r exp σ 2 Φ 1 ( α ) r l t σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r + 1 l 0 t F ( s ) exp σ 2 Φ 1 ( α ) r l t s d s , t M C 2 t < t Φ 1 ( 1 α ) , C 2 t C 2 s < ( t s ) Φ 1 ( 1 α ) , d C 1 s ( γ ) d s > Φ 1 ( α ) , 0 s t , t M γ | d C 2 s ( γ ) d s < Φ 1 ( 1 α ) , d C 1 s ( γ ) d s > Φ 1 ( α ) , 0 s t , t = M γ | d C 2 t ( γ ) d t < Φ 1 ( 1 α ) , d C 1 t ( γ ) d t > Φ 1 ( α ) , t = M γ | d C 1 t ( γ ) d t > Φ 1 ( α ) , t M γ | d C 2 t ( γ ) d t < Φ 1 ( 1 α ) , t = M { Λ 1 } M { Λ 2 + } = 1 α .
Thus, we can get
M { I t Ψ t 1 ( α ) , t } α
and
M { I t > Ψ t 1 ( α ) , t } 1 α .
Note that
{ I t Ψ t 1 ( α ) , t }
and
{ I t > Ψ t 1 ( α ) , t }
are two disjoint events with each other. Thus, we obtain
{ I t > Ψ t 1 ( α ) , t } { I t Ψ t 1 ( α ) , t } c .
According to the duality axiom and monotonicity theorem, we have
M { I t Ψ t 1 ( α ) , t } + M { I t > Ψ t 1 ( α ) , t } 1 .
It follows from (3), (4) and (5) that
M { I t Ψ t 1 ( α ) , t } = α
and
M { I t > Ψ t 1 ( α ) , t } = 1 α .
The lemma is proved. □
Theorem 2.
Let I t be the solution of the RL uncertain circuit equation
d I t = F ( t ) l r l I t d t + σ 1 l d C 1 t σ 2 l I t d C 2 t .
Then I t has an inverse uncertainty distribution
Ψ t 1 ( α ) = I 0 + σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r exp σ 2 Φ 1 ( α ) r l t σ 1 Φ 1 ( α ) σ 2 Φ 1 ( α ) r + 1 l 0 t F ( s ) exp σ 2 Φ 1 ( α ) r l t s d s .
Proof of Theorem 2.
For each α ( 0 , 1 ) and t > 0 , it is easy to infer that
{ I s Ψ s 1 ( α ) , s } { I t Ψ t 1 ( α ) }
and
{ I s > Ψ s 1 ( α ) , s } { I t > Ψ t 1 ( α ) } .
It follows from the monotonicity theorem and Lemma 1 that
M { I t Ψ t 1 ( α ) } { I s Ψ s 1 ( α ) , s } = α
and
M { I t > Ψ t 1 ( α ) } { I s > Ψ s 1 ( α ) , s } = 1 α .
Since { I t Ψ t 1 ( α ) } and { I t > Ψ t 1 ( α ) } are two opposite events, we have
M { I t Ψ t 1 ( α ) } + M { I t > Ψ t 1 ( α ) } = 1 .
Thus,
M { I t Ψ t 1 ( α ) } = α
and Ψ t 1 ( α ) is the inverse uncertainty distribution of I t immediately. □

2.4. Examples of Uncertain RL Circuit Equation

Next we provide some examples to illustrate the results of the solution of uncertain RL circuit Equation (1) and the inverse uncertainty distribution of solution.
Example 1.
Suppose that F ( t ) = 5 and σ 1 = σ 2 = 1 and then the potential source is 5 + Δ C 1 t Δ t . Then, the differential equation
d I t = 5 l r l I t d t + 1 l d C 1 t 1 l I t d C 2 t
is an uncertain RL circuit equation. It follows from Theorem 1 that uncertain RL circuit Equation (6) has a solution
I t = I 0 exp r l t 1 l C 2 t + 5 l 0 t exp r l s t + 1 l C 2 s C 2 t d s + 1 l 0 t exp r l s t + 1 l C 2 s C 2 t d C 1 s .
By using Theorem 2, I t has an inverse uncertainty distribution
Ψ t 1 ( α ) = I 0 + Φ 1 ( α ) + 5 Φ 1 ( α ) r exp Φ 1 ( α ) r l t Φ 1 ( α ) + 5 Φ 1 ( α ) r
which is shown in Figure 1.
Example 2.
Suppose that F ( t ) = exp ( t ) and σ 1 = σ 2 = 0.1 and then the potential source is
exp ( t ) + 0.1 Δ C 1 t Δ t .
Then, the differential equation
d I t = exp ( t ) l r l I t d t + 1 10 l d C 1 t 1 10 l I t d C 2 t
is an uncertain RL circuit equation. It follows from Theorem 1 that uncertain RL circuit Equation (7) has a solution
I t = I 0 exp r l t 1 10 l C 2 t + 1 l 0 t exp r l l s r l t + 1 10 l C 2 s C 2 t d s + 1 10 l 0 t exp r l s t + 1 10 l C 2 s C 2 t d C 1 s .
By using Theorem 2, I t has an inverse uncertainty distribution
Ψ t 1 ( α ) = I 0 + Φ 1 ( α ) Φ 1 ( α ) 10 r + 10 Φ 1 ( α ) 10 r + 10 l exp Φ 1 ( α ) 10 r 10 l t 10 Φ 1 ( α ) 10 r + 10 l exp ( t ) Φ 1 ( α ) Φ 1 ( α ) 10 r
which is shown in Figure 2.

3. Applications of Solution

In this section, two applications of solution are introduced. First is the supremum of solution and the other is the time integral of solution.

3.1. Supremum of Solution

When current flows through a circuit system, we usually care about the maximum value of the current in the time interval [0, t] to prevent it from exceeding the maximum fault current that the circuit system can withstand, which is showed in the following theorem.
Theorem 3.
Let I t be the solution of the uncertain RL circuit equation
d I t = F ( t ) l r l I t d t + σ 1 l d C 1 t σ 2 l I t d C 2 t .
Then, the maximum current through the circuit over time [ 0 , t ] ,
sup 0 s t I s ,
has an inverse uncertainty distribution
Υ t 1 ( α ) = sup 0 s t Ψ s 1 ( α )
where Ψ t 1 ( α ) is the inverse uncertainty distribution of I t .
Proof of Theorem 3.
For any given time t > 0 and α ( 0 , 1 ) , it follows from the basic property of extreme value that
I s Ψ s 1 ( α ) , s sup 0 s t I s sup 0 s t Ψ s 1 ( α )
and
I s > Ψ s 1 ( α ) , s sup 0 s t I s > sup 0 s t Ψ s 1 ( α ) .
By using the monotonicity theorem and Lemma 1, we can get
M sup 0 s t I s sup 0 s t Ψ s 1 ( α ) M I s Ψ s 1 ( α ) , s = α
and
M sup 0 s t I s > sup 0 s t Ψ s 1 ( α ) M I s > Ψ s 1 ( α ) , s = 1 α .
Since
sup 0 s t I s sup 0 s t Ψ s 1 ( α )
and
sup 0 s t I s > sup 0 s t Ψ s 1 ( α )
are two opposite events, we have
M sup 0 s t I s sup 0 s t Ψ s 1 ( α ) + M sup 0 s t I s > sup 0 s t Ψ s 1 ( α ) = 1 .
Hence
M sup 0 s t I s sup 0 s t Ψ s 1 ( α ) = α
and Υ t 1 ( α ) is the inverse uncertainty distribution of
sup 0 s t I s
immediately. □

3.2. Time Integral of Solution

Sometimes, we also care about the amount of charge passing through an electrical component such as resistor in the time period [ 0 , t ] , which can be calculated by the time integral of I t .
Theorem 4.
Let I t be the solution of the uncertain RL circuit equation
d I t = F ( t ) l r l I t d t + σ 1 l d C 1 t σ 2 l I t d C 2 t .
Then, the amount of charge passing through the resistor over [ 0 , t ] ,
0 t I s d s ,
has an inverse uncertainty distribution
G t 1 ( α ) = 0 t Ψ s 1 ( α ) d s
where Ψ t 1 ( α ) is the inverse uncertainty distribution of I t .
Proof of Theorem 4.
For any given time t > 0 and α ( 0 , 1 ) , it follows from the basic property of time integral that
I s Ψ s 1 ( α ) , s 0 t I s d s 0 t Ψ s 1 ( α ) d s
and
I s > Ψ s 1 ( α ) , s 0 t I s d s > 0 t Ψ s 1 ( α ) d s .
By using the monotonicity theorem and Lemma 1, we can get
M 0 t I s d s 0 t Ψ s 1 ( α ) d s M I s Ψ s 1 ( α ) , s = α
and
M 0 t I s d s > 0 t Ψ s 1 ( α ) d s M I s > Ψ s 1 ( α ) , s = 1 α .
Since
0 t I s d s 0 t Ψ s 1 ( α ) d s
and
0 t I s d s > 0 t Ψ s 1 ( α ) d s
are two opposite events, we have
M 0 t I s d s 0 t Ψ s 1 ( α ) d s + M 0 t I s d s > 0 t Ψ s 1 ( α ) d s = 1 .
Hence
M 0 t I s d s 0 t Ψ s 1 ( α ) d s = α
and G t 1 ( α ) is the inverse uncertainty distribution of
0 t I s d s
immediately. □

3.3. Example of Applications

Next we will provide an example to illustrate the results of the applications of solution.
Example 3.
Consider the uncertain RL circuit equation in Example 1. According to Theorem 3, the inverse uncertainty distribution of the maximum current through the circuit over time [ 0 , t ] is
Υ t 1 ( α ) = I 0 , if α 1 + exp π 5 r I 0 3 I 0 + 1 1 I 0 + Φ 1 ( α ) + 5 Φ 1 ( α ) r exp Φ 1 ( α ) r l t Φ 1 ( α ) + 5 Φ 1 ( α ) r , if α > 1 + exp π 5 r I 0 3 I 0 + 1 1
and is shown in Figure 3. By using Theorem 4, we can get the inverse uncertainty distribution of the amount of charge passing through the resistor over [ 0 , t ] as
G t 1 ( α ) = l Φ 1 ( α ) r I 0 + Φ 1 ( α ) + 5 Φ 1 ( α ) r exp Φ 1 ( α ) r l t 1 Φ 1 ( α ) + 5 Φ 1 ( α ) r t
which is shown in Figure 4.

4. Parameter Estimation

How to estimate the unknown parameters of uncertain RL circuit equation based on observed data is a core problem in practice. In this section, we will apply the method of moments to estimating the unknown parameters.
Consider the uncertain RL circuit Equation (1)
d I t = F ( t ) l r l I t d t + σ 1 l d C 1 t σ 2 l I t d C 2 t .
Suppose that F ( t ) is a continuous functions on , r, l, σ 1 > 0 and σ 2 > 0 are unknown parameters to be estimated. Ye [28] suggested that we should consider the difference form of uncertain RL circuit Equation (1), which can be written as
I t j + 1 I t j = F ( t j ) l r l I t j ( t j + 1 t j ) + σ 1 l ( C 1 t j + 1 C 1 t j ) σ 2 l I t j ( C 2 t j + 1 C 2 t j ) ,
i.e.,
l I t j + 1 l I t j + r I t j F ( t j ) ( t j + 1 t j ) = σ 1 ( C 1 t j + 1 C 1 t j ) σ 2 I t j ( C 2 t j + 1 C 2 t j ) .
Dividing both sides of the above equation by
σ 1 + σ 2 | I t j | ( t j + 1 t j ) ,
we obtain
l I t j + 1 l I t j + r I t j F ( t j ) ( t j + 1 t j ) σ 1 + σ 2 | I t j | ( t j + 1 t j ) = σ 1 ( C 1 t j + 1 C 1 t j ) σ 2 I t j ( C 2 t j + 1 C 2 t j ) σ 1 + σ 2 | I t j | ( t j + 1 t j ) .
Suppose that there are n observed data i t 1 , i t 2 , , i t n of the solution I t at time-points t 1 < t 2 < < t n . By substituting the observed data into Equation (8), we write
h j ( r , l , σ 1 , σ 2 ) = l i t j + 1 l i t j + r i t j F ( t j ) ( t j + 1 t j ) σ 1 + σ 2 | i t j | ( t j + 1 t j )
and
ξ j = σ 1 ( C 1 t j + 1 C 1 t j ) σ 2 i t j ( C 2 t j + 1 C 2 t j ) σ 1 + σ 2 | i t j | ( t j + 1 t j )
for j = 1 , 2 , , n 1 . It is easy to infer that h j ( r , l , σ 1 , σ 2 ) can be regarded as a sample of ξ j with 1 j n 1 . On the other hand, it follows from the definition of Liu process and C 1 t and C 2 t are two independent Liu processes that
ξ j N ( 0 , σ 1 + σ 2 | i t j | ( t j + 1 t j ) ) σ 1 + σ 2 | i t j | ( t j + 1 t j ) = N ( 0 , 1 )
for j = 1 , 2 , , n 1 . Thus, h j ( r , l , σ 1 , σ 2 ) , j = 1 , 2 , , n 1 can be regarded as n 1 samples of the standard normal uncertainty distribution N ( 0 , 1 ) . For the purpose of estimating the unknown parameters, Yao and Liu [14] suggested that we should equate the first 4 sample moments to the corresponding first 4 population moments. Then, we can get the following system of equations:
1 n 1 j = 1 n 1 h j ( r , l , σ 1 , σ 2 ) = 0 1 n 1 j = 1 n 1 ( h j ( r , l , σ 1 , σ 2 ) ) 2 = 1 1 n 1 j = 1 n 1 ( h j ( r , l , σ 1 , σ 2 ) ) 3 = 0 1 n 1 j = 1 n 1 ( h j ( r , l , σ 1 , σ 2 ) ) 4 = 21 5 .
The solution ( r * , l * , σ 1 * , σ 2 * ) of the above system of equations is the estimate of unknown parameters, which is called the moment estimate of uncertain RL circuit Equation (1).
Next, we will provide an example to illustrate the moment estimate of uncertain RL circuit equation.
Example 4.
Suppose the potential source provided to the series RL circuit is 5 V, and Simulink (MATLAB R2021a, 9.10.0.1602886, maci64, Simulink, 10.3.) has provided a credible collection of current through the circuit as shown in Table 1. Since the potential source is constant at 5 V, we have
h j ( r , l , σ 1 , σ 2 ) = l i t j + 1 l i t j + r i t j 5 ( t j + 1 t j ) σ 1 + σ 2 | i t j | ( t j + 1 t j )
for j = 1 , 2 , , 29 . By solving the system of Equation (11), we have
r * = 0.9888 , l * = 1.1621 , σ 1 * = 0.0015 , σ 2 * = 0.0318 .
Therefore, the uncertain RL circuit equation is
d I t = 4.3024 0.8508 I t d t + 0.0013 d C 1 t 0.0274 I t d C 2 t .
Next, we test whether the estimates of unknown parameters are appropriate. By substituting the estimates into (12), we can get 29 residuals ε 1 , ε 2 , , ε 29 of uncertain RL circuit Equation (13) as shown in Table 2. Since we have regarded h 1 ( r , l , σ 1 , σ 2 ) , h 2 ( r , l , σ 1 , σ 2 ) ,   , h 29 ( r , l , σ 1 , σ 2 ) as samples of the population N ( 0 , 1 ) , then if the estimates are appropriate, the residuals should be samples of N ( 0 , 1 ) .
Suppose that the residuals are samples of N ( e , σ ) . Then, we consider the following hypotheses:
H 0 : e = 0   and   σ = 1   versus   H 1 : e 0   or   σ 1 .
For the level of significance α = 0.05 , it follows from uncertain hypothesis test (Ye and Liu [29]) that the test for the hypotheses (14) is
W = { ( z 1 , z 2 , , z 29 ) : there   is   at   least 2 of   indexes j s   with 1 j 29 such   that z j < 2.0198 or z j > 2.0198 } .
As shown in Figure 5, we can see that only
ε 1 [ 2.0198 , 2.0198 ] .
Thus, ( ε 1 , ε 2 , , ε 29 ) W and then we think the estimated parameters ( r * , l * , σ 1 * , σ 2 * ) is appropriate.
In addition, if we divide the residuals into two groups to test whether ε 1 , ε 2 , , ε 9 and ε 10 , ε 11 , , ε 29 are from the same population, then the result of the two-sample Kolmogorov–Smirnov test is to reject the hypothesis that they are from the same population with a significance level α = 0.05 , and the p-value is 0.0038 . That is, the residuals are not identically distributed in the sense of probability theory. Therefore, we should regard the residuals as uncertain variables and use uncertain differential equation to model the circuit system.

5. Conclusions

By considering the external noise and internal noise, this paper first applied a multifactor uncertain differential equation to modeling the time evolution of a simple series RL electrical circuit system. Based on the uncertain RL circuit equation, some results were obtained, such as the solution of uncertain RL circuit equation, the inverse uncertainty distribution of solution and the applications of solution. Furthermore, parameter estimation for the uncertain RL circuit equation was also considered.

6. Future Work

Further research may consider modeling more complex circuit systems such as RLC circuit systems, delay circuit systems, and circuit networks.

Author Contributions

Conceptualization, Y.L.; methodology, Y.L.; software, Y.L.; validation, Y.L. and L.Z.; formal analysis, Y.L.; investigation, Y.L. and L.Z.; resources, Y.L. and L.Z.; data curation, Y.L.; writing—original draft preparation, Y.L.; writing—review and editing, Y.L. and L.Z.; supervision, L.Z.; funding acquisition, L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China Grant No. 61873329.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this paper is generated by Simulink.

Acknowledgments

The authors especially thank the editors and anonymous referees for their kindly review and helpful comments. In addition, the authors would like to acknowledge the gracious support of this work by the National Natural Science Foundation of China. Any remaining errors are ours.

Conflicts of Interest

We declare that we have no relevant or material financial interests that relate to the research described in this paper. Neither the entire paper nor any part of its content has been published or has been accepted elsewhere. It is also not being submitted to any other journal.

References

  1. Liu, B. Uncertainty Theory, 2nd ed.; Springer: Berlin, Germany, 2007. [Google Scholar]
  2. Liu, B. Some research problems in uncertainty theory. J. Uncertain Syst. 2009, 3, 3–10. [Google Scholar]
  3. Liu, B. Fuzzy process, hybrid process and uncertain process. J. Uncertain Syst. 2008, 2, 3–16. [Google Scholar]
  4. Chen, X.; Liu, B. Existence and uniqueness theorem for uncertain differential equations. Fuzzy Optim. Decis. Ma. 2010, 9, 69–81. [Google Scholar]
  5. Gao, Y. Existence and uniqueness theorem on uncertain differential equations with local Lipschitz condition. J. Uncertain Syst. 2012, 6, 223–232. [Google Scholar]
  6. Yao, K.; Gao, J.; Gao, Y. Some stability theorems of uncertain differential equation. Fuzzy Optim. Decis. Ma. 2013, 12, 3–13. [Google Scholar] [CrossRef]
  7. Sheng, Y.; Wang, C. Stability in the p-th moment for uncertain differential equation. J. Intell. Fuzzy Syst. 2014, 26, 1263–1271. [Google Scholar] [CrossRef]
  8. Yao, K.; Ke, H.; Sheng, Y. Stability in mean for uncertain differential equation. Fuzzy Optim. Decis. Ma. 2015, 14, 365–379. [Google Scholar] [CrossRef]
  9. Yang, X.; Ni, Y.; Zhang, Y. Stability in inverse distribution for uncertain differential equations. J. Intell. Fuzzy Syst. 2017, 32, 2051–2059. [Google Scholar] [CrossRef] [Green Version]
  10. Yao, K.; Chen, X. A numerical method for solving uncertain differential equations. J. Intell. Fuzzy Syst. 2013, 25, 825–832. [Google Scholar] [CrossRef]
  11. Yang, X.; Shen, Y. Runge-Kutta method for solving uncertain differential equations. J. Uncertain. Anal. Appl. 2015, 3, 17. [Google Scholar] [CrossRef] [Green Version]
  12. Yang, X.; Ralescu, D. Adams method for solving uncertain differential equations. Appl. Math. Comput. 2015, 270, 993–1003. [Google Scholar] [CrossRef]
  13. Gao, R. Milne method for solving uncertain differential equations. Appl. Math. Comput. 2016, 274, 774–785. [Google Scholar]
  14. Yao, K.; Liu, B. Parameter estimation in uncertain differential equations. Fuzzy Optim. Decis. Ma. 2020, 19, 1–12. [Google Scholar] [CrossRef]
  15. Liu, Z. Generalized moment estimation for uncertain differential equations. Appl. Math. Comput. 2021, 392, 125724. [Google Scholar] [CrossRef]
  16. Yang, X.; Liu, Y.; Park, G. Parameter estimation of uncertain differential equation with application to financial market. Chaos Solitons Fract. 2020, 139, 110026. [Google Scholar] [CrossRef]
  17. Sheng, Y.; Yao, K.; Chen, X. Least squares estimation in uncertain differential equations. IEEE Trans. Fuzzy Syst. 2020, 28, 2651–2655. [Google Scholar] [CrossRef]
  18. Liu, Y.; Liu, B. Estimating Unknown Parameters in Uncertain Differential Equation by Maximum Likelihood Estimation. 2020. submitted. [Google Scholar]
  19. Lio, W.; Liu, B. Initial value estimation of uncertain differential equations and zero-day of COVID-19 spread in China. Fuzzy Optim. Decis. Ma. 2021, 20, 177–188. [Google Scholar]
  20. Chen, X.; Li, J.; Xiao, C.; Yang, P. Numerical solution and parameter estimation for uncertain SIR model with application to COVID-19. Fuzzy Optim. Decis. Ma. 2021, 20, 189–208. [Google Scholar]
  21. Jia, L.; Chen, W. Uncertain SEIAR model for COVID-19 cases in China. Fuzzy Optim. Decis. Mak. 2021, 20, 243–259. [Google Scholar] [CrossRef]
  22. Bennett, W. Electrical Noise; McGraw-Hill: New York, NY, USA, 1960. [Google Scholar]
  23. Kampowsky, W.; Rentrop, P.; Schmidt, W. Classification and numerical simulation of electric circuits. Sur. Math. Ind. 1992, 2, 23–65. [Google Scholar]
  24. Penski, C. A new numerical method for SDEs and its application in circuit simulation. J. Comput. Appl. Math. 2000, 115, 461–470. [Google Scholar] [CrossRef] [Green Version]
  25. Kolarova, E. Modelling RL electrical circuits by stochastic differential equations. In Proceedings of the The International Conference on “Computer as a Tool”, Belgrade, Serbia, 21–24 November 2005; pp. 1236–1238. [Google Scholar]
  26. Kolarova, E. An application of stochastic integral equations to electrical networks. Acta Electrotech. Inform. 2008, 8, 14–17. [Google Scholar]
  27. Liu, Y. Uncertain circuit equation. J. Uncertain Syst. 2021, 14, 2150018. [Google Scholar] [CrossRef]
  28. Ye, T. Parameter Estimation in Multifactor Uncertain Differential Equation with Application to Stock Market; Technical Report; 2021. [Google Scholar]
  29. Ye, T.; Liu, B. Uncertain hypothesis test with application to uncertain regression analysis. Fuzzy Optim. Decis. Mak. 2021. [Google Scholar] [CrossRef]
Figure 1. Inverse uncertainty distribution of solution with I 0 = 0 , r = 5 and l = 10 in Example 1.
Figure 1. Inverse uncertainty distribution of solution with I 0 = 0 , r = 5 and l = 10 in Example 1.
Symmetry 13 02103 g001
Figure 2. Inverse uncertainty distribution of solution with I 0 = 0 , r = 5 and l = 10 in Example 2.
Figure 2. Inverse uncertainty distribution of solution with I 0 = 0 , r = 5 and l = 10 in Example 2.
Symmetry 13 02103 g002
Figure 3. Inverse uncertainty distribution of the maximum current with α ( 0 , 0.8 ] , I 0 = 0 , r = 5 and l = 10 in Example 3.
Figure 3. Inverse uncertainty distribution of the maximum current with α ( 0 , 0.8 ] , I 0 = 0 , r = 5 and l = 10 in Example 3.
Symmetry 13 02103 g003
Figure 4. Inverse uncertainty distribution of the amount of charge passing through the resistor with I 0 = 0 , r = 5 and l = 10 in Example 3.
Figure 4. Inverse uncertainty distribution of the amount of charge passing through the resistor with I 0 = 0 , r = 5 and l = 10 in Example 3.
Symmetry 13 02103 g004
Figure 5. Residual plot in Example 4.
Figure 5. Residual plot in Example 4.
Symmetry 13 02103 g005
Table 1. Observations of a simple series RL circuit in Example 4.
Table 1. Observations of a simple series RL circuit in Example 4.
t j (s)0.16150.36150.56150.76151.00001.19241.39241.5924
i t j (A)0.72481.47432.08792.59033.07183.36963.62413.8324
t j (s)1.79242.03812.22832.62832.82833.20003.60003.8000
i t j (A)4.00304.18274.33984.58694.67804.80584.89944.9338
t j (s)4.00004.40004.80005.00005.40005.80006.00006.4000
i t j (A)4.96215.04625.10275.12345.06355.02345.00865.0090
t j (s)6.80007.00007.20007.40008.00009.2000
i t j (A)5.00925.00935.03765.06085.10804.8716
Table 2. Residual plot in Example 4.
Table 2. Residual plot in Example 4.
j12345678
ε j 2.92050.4789−0.2392−1.1031−1.6516−1.7423−1.7655−1.7768
j910111213141516
ε j −1.48990.71020.06400.43880.16620.15360.28090.2700
j1718192021222324
ε j 0.94520.94791.0107−0.6577−0.6759−0.7379−0.2892−0.2904
j2526272829
ε j −0.28910.72990.71560.5866−1.0873
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Y.; Zhou, L. Modeling RL Electrical Circuit by Multifactor Uncertain Differential Equation. Symmetry 2021, 13, 2103. https://doi.org/10.3390/sym13112103

AMA Style

Liu Y, Zhou L. Modeling RL Electrical Circuit by Multifactor Uncertain Differential Equation. Symmetry. 2021; 13(11):2103. https://doi.org/10.3390/sym13112103

Chicago/Turabian Style

Liu, Yang, and Lujun Zhou. 2021. "Modeling RL Electrical Circuit by Multifactor Uncertain Differential Equation" Symmetry 13, no. 11: 2103. https://doi.org/10.3390/sym13112103

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