# Minimal Solution of Singular LR Fuzzy Linear Systems

Authors: M. Nikuie; M. Z. Ahmad
Journal: The Scientific World Journal
 (2014)

---

## Abstract

In this paper, the singular LR fuzzy linear system is introduced. Such systems are divided into two parts: singular consistent LR fuzzy linear systems and singular inconsistent LR fuzzy linear systems. The capability of the generalized inverses such as Drazin inverse, pseudoinverse, and {1}-inverse in finding minimal solution of singular consistent LR fuzzy linear systems is investigated.


---

## Body

## 1. Introduction

The concept of fuzzy numbers and fuzzy arithmetic operations was first introduced and investigated by Zadeh [1, 2], Dubois and Prade [3], and Nahmias [4]. Some different approaches to fuzzy numbers and the structure of fuzzy number spaces were given by Puri and Ralescu [5], Goetschel and Voxman [6, 7], and Cong-Xin and Ming [8, 9]. Fuzzy systems are used to a variety of problems ranging from control chaotic systems [10] to fuzzy metric spaces [11], fuzzy linear systems, fuzzy differential equations [12], particle physics [13], and so on. System of simulations linear equations plays major role in various areas such as mathematics, physics, statistics, engineering, and social sciences [14]. Linear systems of equations plays major role in various areas of sciences and engineering such as fuzzy mathematics, and so forth [15]. The linear system of equations
(1)
A
x
~
=
b
~
,
where A is a n
×
n crisp matrix and the right-hand side is a vector of fuzzy numbers in parametric form, is called a fuzzy linear system (FLS). Friedman et al. [14] use an embedding method for solving a special case of (1) when A is a nonsingular matrix. They introduce new notions for solution of a FLS (1) and transform the FLS to a crisp linear system. In their method, the ordinary inverse matrix is used to solve the FLS. The fuzzy linear system (1) where A is a n
×
n singular crisp matrix and the right-hand side is a vector of fuzzy numbers in parametric form is called a singular fuzzy linear system (SFLS) [16]. The capability of pseudoinverse [17] and Drazin inverse [16] in solving fuzzy linear systems and {
1
}-inverse [18] and {
1,3
}-inverse [19] in finding fuzzy least-square solutions of inconsistent fuzzy linear systems has been studied.Ordinary inverse for anyn
×
n nonsingular matrix exists and is unique. Unlike the case of the nonsingular matrix, which has a single unique inverse for all purposes, there are different generalized inverses for different purposes. For some purposes, as in the examples of solutions of linear systems, there is not a unique inverse, but any matrix of a certain class will do. Generalized inverses have various applications in the theory of finite Markov chains, the study of singular differentia and difference equations, the investigation of Cesaro-Neumann iterations, cryptograph, iterative methods in numerical analysis, multibody system dynamics, and others [20–23].The fuzzy linear system (1) with LR fuzzy variables and LR fuzzy right-hand side is called a fuzzy LR linear system (FLRLS). Ghanbari et al. [24, 25] give the exact and approximated solutions of fuzzy LR linear systems. They give new algorithms using a least-squares model and the ABS approach. LR-type fuzzy systems of linear algebraic equations are discussed by Abramovich et al. in [26], and based on their work in [27] Allahviranloo et al. study LR fuzzy linear systems. They prove existence of fuzzy number vector solution of LR fuzzy linear system by solving optimization problem and obtain the approximated and exact solutions of mentioned systems. This approximated solution in general is different from the weak fuzzy solution considered by Friedman et al. [14, 25].Here, the singular LR fuzzy linear systems are introduced. Such systems are divided into two parts: singular consistent LR fuzzy linear systems and singular inconsistent LR fuzzy linear systems. In this paper, the capability of the generalized inverses in solving singular consistent LR fuzzy linear systems is investigated. In Section2, we recall some preliminaries. The purpose of this paper is to propose a method for solving singular LR fuzzy linear systems on the LR fuzzy linear systems discussed in Section 3. The proposed method is given in Section 4. Numerical examples are given in Section 5 followed by suggestions and concluding remarks in Section 6.


## 2. Preliminaries

This section gives a brief summary of index of matrix, Drazin inverse, and LR fuzzy numbers. We refer the reader to [28–31] for more information.Definition 1.
LetA
∈
C
n
×
n. The index of A is the size of the largest Jordan block corresponding to the zero eigenvalue of A and is denoted by ind
(
A
).Definition 2.
LetA
∈
C
n
×
n, with ind
(
A
)
=
k. The matrix X of order n is the Drazin inverse of A, denoted by A
D, if X satisfies the following conditions:
(2)
A
X
=
X
A
,
X
A
X
=
X
,
A
k
X
A
=
A
k
.Definition 3.
LetA
∈
C
m
×
n. The matrix X of order n
×
m is the pseudoinverse of A, denoted by A
+, if X satisfies the following conditions:
(3)
A
X
A
=
A
,
X
A
X
=
X
,
(
A
X
)
*
=
A
X
,
(
X
A
)
*
=
X
A
.Definition 4.
LetA
∈
C
m
×
n. The matrix X of order n
×
m is the {
1
}-inverse of A, denoted by A
(
1
), if X satisfies the following condition:
(4)
A
X
A
=
A
.Definition 5.
The fuzzy numberA
~ is of LR-type if there exist shape functions L (for left), R (for right), and scalers α
>
0, β
>
0 with
(5)
μ
A
~
(
x
)
=
{
L
(
a
-
x
α
)
,
x
≤
a
,
R
(
x
-
a
β
)
,
x
≥
a
.
The mean value of A
~, a, is a real number and α, β are called the left and right spreads, respectively. A
~ is denoted by (
a
,
α
,
β
)
LR. A crisp number a is specified to be (
a
,
0,0
)
LR.Definition 6.
The fuzzy numberA
~
=
(
m
,
α
,
α
)
L
L is a symmetric fuzzy number and is denoted by A
~
=
(
m
,
α
)
L.Definition 7.
The set of LR fuzzy numbers is denoted byF
(
R
1
)
LR. For arbitrary LR fuzzy numbers M
~
=
(
m
,
α
,
β
)
LR, N
~
=
(
n
,
γ
,
δ
)
LR, and λ
∈
R
+, we have (1)
λ
M
~
=
(
λ
m
,
λ
α
,
λ
β
)
LR,
(2)
-
M
~
=
(
-
m
,
β
,
α
)
LR,
(3)
M
~
⊕
N
~
=
(
m
+
n
,
α
+
γ
,
β
+
δ
)
LR.


## 3. Singular LR Fuzzy Linear Systems

In this section, singular LR fuzzy linear systems are introduced and the existing methods for solving LR fuzzy linear systems are investigated.


### 3.1. Shortcomings of the Existing Methods

In this subsection, the shortcomings of the existing methods for solving LR fuzzy linear systems are pointed out.(
1
) Ghanbari and Mahdavi-Amiri [24] considered the LR fuzzy linear system (FLRLS) of the form
(6)
A
x
~
=
b
~
,
where A is a m
×
n crisp matrix, the unknown vector x
~
=
(
x
~
1
,
…
,
x
~
n
)
T consists of n fuzzy numbers, and the constant b
~
=
(
b
~
1
,
…
,
b
~
m
)
T is a vector consisting of n LR fuzzy numbers. They convert the LR fuzzy linear system (6) into the corresponding crisp linear system
(7)
A
x
=
b
,
and the constrained least-square problem
(8)
(
B
+
-
B
-
-
B
-
B
+
)
(
α
β
)
=
(
b
l
b
r
)
,
(
α
β
)
>
0
,
where [
B
+
] and [
B
-
] are determined as follows:(9)
[
B
+
]
i
j
=
{
a
i
j
,
a
i
j
≥
0
,
0
,
a
i
j
<
0
,
[
B
-
]
i
j
=
{
a
i
j
,
a
i
j
<
0
,
0
,
a
i
j
≥
0
,
for all i
=
1
,
…
,
m and j
=
1
,
…
,
n. Also, let x
=
(
x
1
,
…
,
x
n
)
T, b
=
(
b
1
,
…
,
b
m
)
T, α
=
(
α
1
,
…
,
α
n
)
T, β
=
(
β
1
,
…
,
β
n
)
T, b
l
=
(
b
1
l
,
…
,
b
m
l
)
T, and b
r
=
(
b
1
r
,
…
,
b
m
r
)
T. They define(10)
r
x
~
=
∥
B
+
α
-
B
-
β
-
b
l
∥
2
+
∥
-
B
-
α
+
B
+
β
-
b
r
∥
2
.For solving FLRLS (6), they must solve systems (7) and (8).(i)
If the crisp system (7) lacks solution, then the FLRLS (6) does not have any solution.
(ii)
Ifx is a solution of the system A
x
=
b and α, β are the solutions of the constrained least-square problem (8), then the corresponding x
~
=
(
x
,
α
,
β
)
LR is said to be an exact solution of (6).
(iii)
Ifr
x
~
>
0 for all x
~
∈
F
(
R
n
)
LR, then the system (10) lacks a solution, and hence FLRLS (6) lacks an exact solution. In this case, if x
~
=
(
x
,
α
,
β
)
LR is a solution of the system A
x
=
b and α, β are the solutions of the constrained least-square problem (8).This method [24] is applicable when in the crisp linear system (7) A
∈
R
m
×
n and rank
⁡
(
A
)
=
m.Example 8.
Consider the LR fuzzy linear system(11)
x
~
1
-
x
~
2
=
(
1,1
,
1
)
LR
,
x
~
1
+
3
x
~
2
=
(
5,1
,
2
)
LR
.By Ghanbari’s method, we obtain an exact solutionx
~
1
=
(
2,0.625,0.875
) using ABS algorithm and x
~
2
=
(
1,0.125,0.375
) usinglsqlin function of MATLAB software.(2) Allahviranloo et al. [27] transform the LR fuzzy linear system
(12)
A
(
x
,
α
x
,
β
x
)
=
(
b
,
α
b
,
β
b
)
into the 1-level crisp linear system (7) and the following minimization problem:
(13)
min
⁡
z
1
+
z
2
s
.
t
.
-
z
1
≤
∑
a
i
j
>
0
‍
a
i
j
α
x
j
-
∑
a
i
j
<
0
‍
a
i
j
β
x
j
-
α
b
i
≤
z
1
,
i
=
1
,
…
,
n
,
-
z
2
≤
∑
a
i
j
>
0
‍
a
i
j
β
x
j
-
∑
a
i
j
<
0
‍
a
i
j
α
x
j
-
β
b
i
≤
z
2
,
i
=
1
,
…
,
n
,
α
x
j
,
β
x
j
,
z
1
,
z
2
≥
0
,
j
=
1
,
…
,
n
.
The spreads of fuzzy number vector solution are obtained from the model (13). Finally, by allocating these spreads to vector solution of 1-level system, the fuzzy number vector solution results and if the computed solution satisfies the FLS, it can be as exact or approximated solution. (i)
If the crisp system (7) lacks solution, then the FLRLS (12) does not have any solution.
(ii)
If model (13) has more than two optimal solutions then it has a set of solutions.Example 9.
Consider the LR fuzzy linear system(14)
x
~
1
-
x
~
2
=
(
1,2
,
3
)
LR
,
x
~
1
+
x
~
2
=
(
3,3
,
2
)
LR
.Using Allahviranloo’s method, the vector solution of the 1-level system is(
2,1
)
t. This fuzzy linear system has infinite number of fuzzy number vector solutions as follows:
(15)
(
x
~
1
x
~
2
)
=
(
λ
(
2,1
,
1
)
+
(
1
-
λ
)
(
2,2
,
2
)
λ
(
1,2
,
1
)
+
(
1
-
λ
)
(
1,1
,
0
)
)
,
λ
∈
[
0,1
]
.
The existing methods are incapable of finding the minimal solution for the following LR fuzzy linear system:
(16)
x
~
1
⊕
(
-
1
)
x
~
2
=
(
1,1
,
1
)
LR
,
2
x
~
1
⊕
(
-
2
)
x
~
2
=
(
2,2
,
2
)
LR
.


### 3.2. Singular LR Fuzzy Linear Systems

In this subsection, the singular consistent or inconsistent LR fuzzy linear systems is introduced.Definition 10.
Then
×
n linear system
(17)
(
a
11
⋯
a
1
n
⋮
⋱
⋮
a
n
1
⋯
a
n
n
)
(
x
~
1
⋮
x
~
n
)
=
(
b
~
1
⋮
b
~
n
)
,
where the n
×
n crisp singular matrix A
=
(
a
i
j
) and LR fuzzy vector b
~ are given and x
~
=
F
(
R
n
)
LR is an unknown LR fuzzy vector to be found, is called a singular LR fuzzy system of equations (SFLRLS).In Definition10, let x
~
i
=
(
x
i
,
α
x
i
,
β
x
i
) and b
~
i
=
(
b
i
,
α
b
i
,
β
b
i
),1
≤
i
≤
n. The n
×
n singular LR fuzzy system of (17) can be extended into the 3
n
×
3
n crisp linear system
(18)
(
A
0
0
0
B
C
0
C
B
)
(
x
x
α
x
β
)
=
(
b
b
α
b
β
)
,
abbreviated
as
S
X
=
Y
,
wherein
(19)
[
B
]
i
j
=
{
a
i
j
,
a
i
j
≥
0
,
0
,
a
i
j
<
0
,
[
C
]
i
j
=
{
a
i
j
,
a
i
j
<
0
,
0
,
a
i
j
≥
0
,
x
α
=
(
α
x
1
⋮
α
x
n
)
,
x
β
=
(
β
x
1
⋮
β
x
n
)
,
b
α
=
(
α
b
1
⋮
α
b
n
)
,
b
β
=
(
β
b
1
⋮
β
b
n
)
.Definition 11.
The singular LR fuzzy linear system (17) is called a singular consistent LR fuzzy linear system while the systems (18) are consistent crisp linear systems, that is, rank
[
S
]
=
rank
[
S
∣
Y
], inconsistent otherwise.


## 3.1. Shortcomings of the Existing Methods

In this subsection, the shortcomings of the existing methods for solving LR fuzzy linear systems are pointed out.(
1
) Ghanbari and Mahdavi-Amiri [24] considered the LR fuzzy linear system (FLRLS) of the form
(6)
A
x
~
=
b
~
,
where A is a m
×
n crisp matrix, the unknown vector x
~
=
(
x
~
1
,
…
,
x
~
n
)
T consists of n fuzzy numbers, and the constant b
~
=
(
b
~
1
,
…
,
b
~
m
)
T is a vector consisting of n LR fuzzy numbers. They convert the LR fuzzy linear system (6) into the corresponding crisp linear system
(7)
A
x
=
b
,
and the constrained least-square problem
(8)
(
B
+
-
B
-
-
B
-
B
+
)
(
α
β
)
=
(
b
l
b
r
)
,
(
α
β
)
>
0
,
where [
B
+
] and [
B
-
] are determined as follows:(9)
[
B
+
]
i
j
=
{
a
i
j
,
a
i
j
≥
0
,
0
,
a
i
j
<
0
,
[
B
-
]
i
j
=
{
a
i
j
,
a
i
j
<
0
,
0
,
a
i
j
≥
0
,
for all i
=
1
,
…
,
m and j
=
1
,
…
,
n. Also, let x
=
(
x
1
,
…
,
x
n
)
T, b
=
(
b
1
,
…
,
b
m
)
T, α
=
(
α
1
,
…
,
α
n
)
T, β
=
(
β
1
,
…
,
β
n
)
T, b
l
=
(
b
1
l
,
…
,
b
m
l
)
T, and b
r
=
(
b
1
r
,
…
,
b
m
r
)
T. They define(10)
r
x
~
=
∥
B
+
α
-
B
-
β
-
b
l
∥
2
+
∥
-
B
-
α
+
B
+
β
-
b
r
∥
2
.For solving FLRLS (6), they must solve systems (7) and (8).(i)
If the crisp system (7) lacks solution, then the FLRLS (6) does not have any solution.
(ii)
Ifx is a solution of the system A
x
=
b and α, β are the solutions of the constrained least-square problem (8), then the corresponding x
~
=
(
x
,
α
,
β
)
LR is said to be an exact solution of (6).
(iii)
Ifr
x
~
>
0 for all x
~
∈
F
(
R
n
)
LR, then the system (10) lacks a solution, and hence FLRLS (6) lacks an exact solution. In this case, if x
~
=
(
x
,
α
,
β
)
LR is a solution of the system A
x
=
b and α, β are the solutions of the constrained least-square problem (8).This method [24] is applicable when in the crisp linear system (7) A
∈
R
m
×
n and rank
⁡
(
A
)
=
m.Example 8.
Consider the LR fuzzy linear system(11)
x
~
1
-
x
~
2
=
(
1,1
,
1
)
LR
,
x
~
1
+
3
x
~
2
=
(
5,1
,
2
)
LR
.By Ghanbari’s method, we obtain an exact solutionx
~
1
=
(
2,0.625,0.875
) using ABS algorithm and x
~
2
=
(
1,0.125,0.375
) usinglsqlin function of MATLAB software.(2) Allahviranloo et al. [27] transform the LR fuzzy linear system
(12)
A
(
x
,
α
x
,
β
x
)
=
(
b
,
α
b
,
β
b
)
into the 1-level crisp linear system (7) and the following minimization problem:
(13)
min
⁡
z
1
+
z
2
s
.
t
.
-
z
1
≤
∑
a
i
j
>
0
‍
a
i
j
α
x
j
-
∑
a
i
j
<
0
‍
a
i
j
β
x
j
-
α
b
i
≤
z
1
,
i
=
1
,
…
,
n
,
-
z
2
≤
∑
a
i
j
>
0
‍
a
i
j
β
x
j
-
∑
a
i
j
<
0
‍
a
i
j
α
x
j
-
β
b
i
≤
z
2
,
i
=
1
,
…
,
n
,
α
x
j
,
β
x
j
,
z
1
,
z
2
≥
0
,
j
=
1
,
…
,
n
.
The spreads of fuzzy number vector solution are obtained from the model (13). Finally, by allocating these spreads to vector solution of 1-level system, the fuzzy number vector solution results and if the computed solution satisfies the FLS, it can be as exact or approximated solution. (i)
If the crisp system (7) lacks solution, then the FLRLS (12) does not have any solution.
(ii)
If model (13) has more than two optimal solutions then it has a set of solutions.Example 9.
Consider the LR fuzzy linear system(14)
x
~
1
-
x
~
2
=
(
1,2
,
3
)
LR
,
x
~
1
+
x
~
2
=
(
3,3
,
2
)
LR
.Using Allahviranloo’s method, the vector solution of the 1-level system is(
2,1
)
t. This fuzzy linear system has infinite number of fuzzy number vector solutions as follows:
(15)
(
x
~
1
x
~
2
)
=
(
λ
(
2,1
,
1
)
+
(
1
-
λ
)
(
2,2
,
2
)
λ
(
1,2
,
1
)
+
(
1
-
λ
)
(
1,1
,
0
)
)
,
λ
∈
[
0,1
]
.
The existing methods are incapable of finding the minimal solution for the following LR fuzzy linear system:
(16)
x
~
1
⊕
(
-
1
)
x
~
2
=
(
1,1
,
1
)
LR
,
2
x
~
1
⊕
(
-
2
)
x
~
2
=
(
2,2
,
2
)
LR
.


## 3.2. Singular LR Fuzzy Linear Systems

In this subsection, the singular consistent or inconsistent LR fuzzy linear systems is introduced.Definition 10.
Then
×
n linear system
(17)
(
a
11
⋯
a
1
n
⋮
⋱
⋮
a
n
1
⋯
a
n
n
)
(
x
~
1
⋮
x
~
n
)
=
(
b
~
1
⋮
b
~
n
)
,
where the n
×
n crisp singular matrix A
=
(
a
i
j
) and LR fuzzy vector b
~ are given and x
~
=
F
(
R
n
)
LR is an unknown LR fuzzy vector to be found, is called a singular LR fuzzy system of equations (SFLRLS).In Definition10, let x
~
i
=
(
x
i
,
α
x
i
,
β
x
i
) and b
~
i
=
(
b
i
,
α
b
i
,
β
b
i
),1
≤
i
≤
n. The n
×
n singular LR fuzzy system of (17) can be extended into the 3
n
×
3
n crisp linear system
(18)
(
A
0
0
0
B
C
0
C
B
)
(
x
x
α
x
β
)
=
(
b
b
α
b
β
)
,
abbreviated
as
S
X
=
Y
,
wherein
(19)
[
B
]
i
j
=
{
a
i
j
,
a
i
j
≥
0
,
0
,
a
i
j
<
0
,
[
C
]
i
j
=
{
a
i
j
,
a
i
j
<
0
,
0
,
a
i
j
≥
0
,
x
α
=
(
α
x
1
⋮
α
x
n
)
,
x
β
=
(
β
x
1
⋮
β
x
n
)
,
b
α
=
(
α
b
1
⋮
α
b
n
)
,
b
β
=
(
β
b
1
⋮
β
b
n
)
.Definition 11.
The singular LR fuzzy linear system (17) is called a singular consistent LR fuzzy linear system while the systems (18) are consistent crisp linear systems, that is, rank
[
S
]
=
rank
[
S
∣
Y
], inconsistent otherwise.


## 4.  Minimal Solution of Singular LR Fuzzy Linear Systems

Let(20)
A
x
~
=
b
~
be a singular consistent LR fuzzy linear system and let (18) be the extended crisp linear system of it. The coefficient matrix system (18) is a singular matrix [16, 32]. In this section, the capability of generalized inverses in finding minimal solution of the singular consistent LR fuzzy linear system (20) is investigated. Equivalent systems of linear equations have exactly the same solutions [33]. Therefore, the fundamental theorem for the singular consistent LR fuzzy linear system is given.Theorem 12 (fundamental theorem).
The system (20) is equivalent to the singular consistent crisp linear system
(21)
A
x
=
b
,
k
=
ind

(
A
)
,
and the singular consistent constrained linear system
(22)
(
B
C
C
B
)
(
x
α
x
β
)
=
(
b
α
b
β
)
,
(
x
α
x
β
)
>
0
.


### 4.1. Notes

Note 1.
If model (22) has a solution, then (20) has an exact LR fuzzy solution and (22) gives exact solutions; otherwise we get a set of approximated solutions. For computing approximated solution instead of the singular constrained linear system (22), we solve the following singular constrained linear system:
(23)
(
B
C
C
B
)
(
x
α
x
β
)
=
(
b
α
b
β
)
,
(
x
α
x
β
)
≥
ɛ
e
,
where e
=
(
1
,
…
,
1
)
T
∈
R
n and ɛ
>
0 is a user-defined parameter.Note 2.
If model (22) has a set of solutions and (21) is a consistent system, then (20) has a set of LR fuzzy exact solutions.Note 3.
If the crisp linear system (21) is a singular inconsistent linear system, the singular LR fuzzy linear system (20) has no solution.Note 4.
We observe that(
x
,
x
α
,
x
β
) is the exact minimal LR fuzzy solution of the singular consistent LR fuzzy linear system (20) if and only if x is a minimal LR fuzzy solution of the singular consistent crisp linear system (21) and the following constrained least-square problem:
(24)
z
=
min
⁡
{
∥
B
x
α
-
C
x
β
-
b
α
∥
2
+
∥
-
B
x
α
+
C
x
β
-
b
β
∥
2
}
s
.
t
.
x
α
>
0
,
x
β
>
0
,
has a solution with Z
*
=
0; otherwise, the singular consistent LR fuzzy linear system has a set of approximated minimal LR fuzzy solutions.


### 4.2. Solution Strategy

According to Theorem12, for solving (20) we must solve the system (21) and singular constrained linear system (22). Finding positive solutions for linear systems of equations has been investigated [34]. Nonnegativity of matrices and their group or Drazin [35] inverses and pseudoinverse [36] were investigated.


#### 4.2.1. Solving Singular Crisp Linear Systems (21)

Theorem 13 (see [20]).
A
D
b is a solution of (21) if and only if b
∈
range
(
A
k
), and A
D
b is an unique solution of (21) provided that x
∈
range
(
A
k
).Theorem 14 (see [37]).
The minimal solution of the system (21) is x
=
A
+
b.Theorem 15 (see [28]).
LetA
∈
C
n
×
n, b
∈
C
n. Then (21) is consistent if and only if, for some A
(
1
),
(25)
A
A
(
1
)
b
=
b
,
in which case the general solution of (21) is
(26)
x
=
A
(
1
)
b
+
(
I
-
A
(
1
)
A
)
y
,
for arbitrary y
∈
C
n.


#### 4.2.2. Solving Singular Constrained Linear System (22)

Illustrative examples for finding positive exact solution for singular linear systems,(27)
(
B
C
C
B
)
(
x
α
x
β
)
=
(
b
α
b
β
)
,
using generalized inverses are given in the next section. It is possible that the generalized inverses do not obtain the positive solution for the system (27). In this case, we use thelsqlin function of MATLAB software. The capability of thelsqlin function of MATLAB software in finding exact or approximated solutions is illustrated in next the section.


## 4.1. Notes

Note 1.
If model (22) has a solution, then (20) has an exact LR fuzzy solution and (22) gives exact solutions; otherwise we get a set of approximated solutions. For computing approximated solution instead of the singular constrained linear system (22), we solve the following singular constrained linear system:
(23)
(
B
C
C
B
)
(
x
α
x
β
)
=
(
b
α
b
β
)
,
(
x
α
x
β
)
≥
ɛ
e
,
where e
=
(
1
,
…
,
1
)
T
∈
R
n and ɛ
>
0 is a user-defined parameter.Note 2.
If model (22) has a set of solutions and (21) is a consistent system, then (20) has a set of LR fuzzy exact solutions.Note 3.
If the crisp linear system (21) is a singular inconsistent linear system, the singular LR fuzzy linear system (20) has no solution.Note 4.
We observe that(
x
,
x
α
,
x
β
) is the exact minimal LR fuzzy solution of the singular consistent LR fuzzy linear system (20) if and only if x is a minimal LR fuzzy solution of the singular consistent crisp linear system (21) and the following constrained least-square problem:
(24)
z
=
min
⁡
{
∥
B
x
α
-
C
x
β
-
b
α
∥
2
+
∥
-
B
x
α
+
C
x
β
-
b
β
∥
2
}
s
.
t
.
x
α
>
0
,
x
β
>
0
,
has a solution with Z
*
=
0; otherwise, the singular consistent LR fuzzy linear system has a set of approximated minimal LR fuzzy solutions.


## 4.2. Solution Strategy

According to Theorem12, for solving (20) we must solve the system (21) and singular constrained linear system (22). Finding positive solutions for linear systems of equations has been investigated [34]. Nonnegativity of matrices and their group or Drazin [35] inverses and pseudoinverse [36] were investigated.


### 4.2.1. Solving Singular Crisp Linear Systems (21)

Theorem 13 (see [20]).
A
D
b is a solution of (21) if and only if b
∈
range
(
A
k
), and A
D
b is an unique solution of (21) provided that x
∈
range
(
A
k
).Theorem 14 (see [37]).
The minimal solution of the system (21) is x
=
A
+
b.Theorem 15 (see [28]).
LetA
∈
C
n
×
n, b
∈
C
n. Then (21) is consistent if and only if, for some A
(
1
),
(25)
A
A
(
1
)
b
=
b
,
in which case the general solution of (21) is
(26)
x
=
A
(
1
)
b
+
(
I
-
A
(
1
)
A
)
y
,
for arbitrary y
∈
C
n.


### 4.2.2. Solving Singular Constrained Linear System (22)

Illustrative examples for finding positive exact solution for singular linear systems,(27)
(
B
C
C
B
)
(
x
α
x
β
)
=
(
b
α
b
β
)
,
using generalized inverses are given in the next section. It is possible that the generalized inverses do not obtain the positive solution for the system (27). In this case, we use thelsqlin function of MATLAB software. The capability of thelsqlin function of MATLAB software in finding exact or approximated solutions is illustrated in next the section.


## 4.2.1. Solving Singular Crisp Linear Systems (21)

Theorem 13 (see [20]).
A
D
b is a solution of (21) if and only if b
∈
range
(
A
k
), and A
D
b is an unique solution of (21) provided that x
∈
range
(
A
k
).Theorem 14 (see [37]).
The minimal solution of the system (21) is x
=
A
+
b.Theorem 15 (see [28]).
LetA
∈
C
n
×
n, b
∈
C
n. Then (21) is consistent if and only if, for some A
(
1
),
(25)
A
A
(
1
)
b
=
b
,
in which case the general solution of (21) is
(26)
x
=
A
(
1
)
b
+
(
I
-
A
(
1
)
A
)
y
,
for arbitrary y
∈
C
n.


## 4.2.2. Solving Singular Constrained Linear System (22)

Illustrative examples for finding positive exact solution for singular linear systems,(27)
(
B
C
C
B
)
(
x
α
x
β
)
=
(
b
α
b
β
)
,
using generalized inverses are given in the next section. It is possible that the generalized inverses do not obtain the positive solution for the system (27). In this case, we use thelsqlin function of MATLAB software. The capability of thelsqlin function of MATLAB software in finding exact or approximated solutions is illustrated in next the section.


## 5. Numerical Results

In this section, solving singular LR fuzzy linear systems is illustrated.Example 16.
Consider the following singular consistent LR fuzzy linear system:(28)
x
~
1
⊕
(
-
1
)
x
~
2
=
(
1,1
,
1
)
LR
,
2
x
~
1
⊕
(
-
2
)
x
~
2
=
(
2,2
,
2
)
LR
.
According to Tables1 and 2 Drazin inverse, pseudoinverse, a {
1
}-inverse and thelsqlin function gives a set of exact LR fuzzy solution for (28).
The exact minimal LR fuzzy solution of the system (28) is x
~
1
,
x
~
2
=
(
0.5000
,
-
0.5000,0.5000
).Table 1
Finding the left and right spreads.
Using
α
x
1
α
x
2
β
x
1
β
x
2
∥
(
x
α
,
x
β
)
t
∥
2
Drazin inverse
0.3333
0.6667
0.3333
0.6667
1.0541
Pseudoinverse
0.5000
0.5000
0.5000
0.5000
1.0000
{
1
}-inverse
0.5000
0.5000
0.5000
0.5000
1.0000
Other{
1
}-inverse
2.0000
2.0000
−1.0000
−1.0000
3.1623
lsqlin function
0.5000
0.5000
0.5000
0.5000
1.0000Table 2
Finding the mean values.
Using
x
1
x
2
∥
x
∥
2
Drazin inverse
−1.0000
−2.0000
2.2361
Pseudoinverse
0.5000
−0.5000
0.7071
{
1
}-inverse
0
−1.0000
1.0000Example 17.
Consider the following singular consistent LR fuzzy linear system:(29)
(
1
-
1
-
1
1
-
4
4
3
-
3
-
3
)
(
x
~
1
x
~
2
x
~
3
)
=
(
(
5,1
,
2
)
LR
(
-
2,10,1
)
LR
(
15,3
,
6
)
LR
)
.
The solutions obtained by generalized inverses for the system (29) are listed in Tables 3 and 4.
The generalized inverses are incapable of finding an exact LR fuzzy solution for (29). Therefore, we use thelsqlin function of MATLAB software for finding a set of approximated solutions (see Table 5).
The approximated minimal LR fuzzy solution withɛ
=
1
0
-
3 is
(30)
x
~
1
=
(
1.5918,0.1631,0.1703
)
,
x
~
2
=
(
-
1.2551,0.0941,0.7236
)
,
x
~
3
=
(
-
2.1531,1.7356,0.1133
)
.Table 3
Finding the mean values.
Using
x
1
x
2
x
3
∥
x
∥
2
Drazin inverse
−1.0476
−2.9048
−3.1429
4.4060
Pseudoinverse
1.5918
−1.2551
−2.1531
2.9572
{
1
}-inverse
4.6667
0.6667
−1.0000
4.8189Table 4
Finding the left and right spreads.
Using
α
x
1
α
x
2
α
x
3
β
x
1
β
x
2
β
x
3
∥
(
x
α
,
x
β
)
t
∥
2
Drazin inverse
0.4762
0.7143
1.4286
−0.1429
0.9524
−0.4286
1.9720
Pseudoinverse
0.0952
0.3095
1.4524
0.2381
1.0238
−0.1190
1.8257
{
1
}-inverse
−1.0000
−1.0833
1.7500
1.3333
1.0000
1.0000
3.0023Table 5
Finding the left and right spreads.
Using
α
x
1
α
x
2
α
x
3
β
x
1
β
x
2
β
x
3
ε
lsqlin function
0.1631
0.0941
1.7356
0.1703
0.7236
0.1133
10−3
lsqlin function
0.1793
0.1000
1.7460
0.1540
0.7092
0.1115
10−2
lsqlin function
0.1867
0.1006
1.7528
0.1466
0.7005
0.1128
10−1Example 18.
Now consider a system where the coefficient matrix is the same as in the previous system, but with the LR fuzzy vector on the right-hand side leads to an exact minimal LR fuzzy solution(31)
(
1
-
1
-
1
1
-
4
4
3
-
3
-
3
)
(
x
~
1
x
~
2
x
~
3
)
=
(
(
1,1
,
1
)
LR
(
-
1,3
,
3.25
)
LR
(
3,3
,
3
)
LR
)
.
According to Tables 6 and 7, the system (31) has a set of the exact minimal LR fuzzy solutions.
Hence forλ
∈
[
0,1
], the singular LR fuzzy linear system (31) has a set of exact minimal LR fuzzy solution as follows:(32)
(
x
~
1
x
~
2
x
~
3
)
=
(
λ
(
0.3061,0.2885,0.2948
)
LR
+
(
1
-
λ
)
(
0.3061,0.2891,0.2942
)
LR
λ
(
-
0.1837,0.3638,0.3365
)
LR
+
(
1
-
λ
)
(
-
0.1837,0.3682,0.3401
)
LR
λ
(
-
0.5102,0.3414,0.3750
)
LR
+
(
1
-
λ
)
(
-
0.5102,0.3376,0.3707
)
LR
)
.Table 6
Finding the left and right spreads.
Using
α
x
1
α
x
2
α
x
3
β
x
1
β
x
2
β
x
3
∥
(
x
α
,
x
β
)
t
∥
2
Drazin inverse
0.2857
−0.1548
0.8571
0.2976
−0.1786
0.8929
1.3259
Pseudoinverse
0.2891
0.3682
0.3376
0.2942
0.3401
0.3707
0.8203
{
1
}-inverse
−1.0000
−1.0833
1.7500
1.3333
1.0000
1.0000
3.0023
lsqlin function
0.2885
0.3638
0.3414
0.2948
0.3365
0.3750
0.8203Table 7
Finding the mean values.
Using
x
1
x
2
x
3
∥
x
∥
2
Drazin inverse
−0.2381
−0.5238
−0.7143
0.9172
Pseudoinverse
0.3061
−0.1837
−0.5102
0.6227
{
1
}-inverse
−1.0000
−1.0000
−1.0000
1.7321Example 19 (singular dual LR fuzzy linear systems).
Usually, there is no inverse element for an arbitrary LR fuzzy number [38]. Also, if u
~
=
F
(
R
1
)
LR is an arbitrary LR fuzzy number, then we have
(33)
u
~
⊕
-
u
~
≠
0
.
Therefore, the following two fuzzy linear systems are not equivalent
(34)
A
1
x
~
+
b
1
=
A
2
x
~
+
b
2
,
(
A
1
-
A
2
)
x
~
=
(
b
2
-
b
1
)
.
In [38] the necessary and sufficient conditions for A
1
x
~
=
A
2
x
~
+
b
~
2 were given, where A
1 and A
2 are square matrices and in [25] on the DFLRLS (17) was discussed. Now, we consider an economic application of minimal solution of general dual fuzzy linear systems. The market price of goods and the quantity produced are determined by the equality between supply and demand. Suppose that demand and supply are linear functions of the price
(35)
q
d
+
a
=
b
*
p
+
c
,
q
s
+
d
=
e
*
p
+
f
,
where q
s is the quantity supplied, which is required to be equal to q
d, the quantity requested, p is the price, a, b, c, d, e, and f are coefficients to be estimated, where the coefficients a, c, d, and f are represented by LR fuzzy numbers, and b and e are crisp numbers. By imposing the equality between quantity supplied and requested, the following general dual fuzzy linear system should be solved:
(36)
x
~
1
⊕
(
5,5
,
4.5
)
LR
=
x
~
2
⊕
(
1,1
,
1
)
LR
,
2
x
~
1
⊕
(
10,10,9
)
LR
=
2
x
~
2
⊕
(
2,2
,
2
)
LR
.
Using Drazin inverse we get the exact solution x
~
1
=
(
4,4
,
3.5
)
LR and x
~
2
=
(
8,8
,
7
)
LR. For economic analysis, we refer to [17, 39].


## 6. Conclusions and Suggestions

In this paper, then
×
n singular LR fuzzy linear system A
x
~
=
b
~ is transformed into the 3
n
×
3
n singular crisp linear system S
X
=
Y. Such systems are divided into two parts: singular consistent LR fuzzy linear systems and singular inconsistent LR fuzzy linear systems. The singular consistent LR fuzzy linear system A
x
~
=
b
~ has an exact minimal LR fuzzy solution or a set of the approximated minimal LR fuzzy solutions. It is possible that the generalized inverses be incapable of finding the exact or approximated minimal LR fuzzy solution of the system A
x
~
=
b
~. In this case, we use thelsqlin function of MATLAB software. We hope that this paper will be a new path to improve solving singular consistent LR fuzzy linear systems using preconditioning techniques and iterative methods.



---

*Source: 517218-2014-03-11.xml*