Applied Computational Intelligence and Soft Computing

Volume 2018, Article ID 2797038, 13 pages

https://doi.org/10.1155/2018/2797038

## Two Iterative Methods for Solving Linear Interval Systems

^{1}Department of Computer Science, Yazd University, Yazd, Iran^{2}Parallel Processing Laboratory, Yazd University, Yazd, Iran

Correspondence should be addressed to Seyed Abolfazl Shahzadeh Fazeli; ri.ca.dzay@ilezaf

Received 28 May 2018; Revised 28 August 2018; Accepted 24 September 2018; Published 8 October 2018

Academic Editor: Shyi-Ming Chen

Copyright © 2018 Esmaeil Siahlooei and Seyed Abolfazl Shahzadeh Fazeli. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Conjugate gradient is an iterative method that solves a linear system , where is a positive definite matrix. We present this new iterative method for solving linear interval systems , where is a diagonally dominant interval matrix, as defined in this paper. Our method is based on conjugate gradient algorithm in the context view of interval numbers. Numerical experiments show that the new interval modified conjugate gradient method minimizes the norm of the difference of and at every step while the norm is sufficiently small. In addition, we present another iterative method that solves , where is a diagonally dominant interval matrix. This method, using the idea of steepest descent, finds exact solution for linear interval systems, where ; we present a proof that indicates that this iterative method is convergent. Also, our numerical experiments illustrate the efficiency of the proposed methods.

#### 1. Introduction

Solving system of linear equations is a well-known problem in linear algebra. Many practical problems are modeled as system of linear equations. These problems have been studied by many scientists and several methods have been proposed to solve them [1, 2]. But, in everyday life, measuring instruments just estimate values, and usually measured values are not accurate. Sometimes for more accuracy, intervals are used to represent the actual values. For example, the length of a metal rod is estimated between and centimeters, or temperature measured by a thermometer is between and Celsius.

There are many types of uncertainty and there are many different mathematical systems that calculate uncertainties [3], e.g., rough sets theory [4], fuzzy numbers [5], probability theory [6], interval valued numbers [7], and dual fuzzy numbers [8]. In this paper, uncertainty is considered as interval numbers. Some articles used this model to solve their problems [9, 10].

Using interval numbers in algebra was initially developed in the mid-1960s. In 1966, Moore presented his book on interval analysis [11]. Then Hansen offered a solution on interval linear algebraic equations [12]. Then many authors published their methods for solving linear interval systems, such as Neumaier, Abolmasoumi and Alavi, Nirmala, and Ganesan [13–15]. Nowadays, interval analysis methods have been applied to engineering problems, such as dynamic response analysis [16, 17], geotechnical structures [18], and control systems [19].

In this paper, we introduce two new iterative methods for solving a linear interval system of equations that is a linear system involving uncertain coefficients appearing as interval numbers. The solution of this system is an interval vector.

The present paper is organized as follows. The basic definitions related to interval numbers and linear interval systems are discussed in Section 3. In Section 3.1, we introduce two additional interval operations. These are inverse operations of “−” and “+.” Section 3.2 recalls some properties of interval numbers and interval arithmetic. Also, it introduces some new properties on introduced operations. Next, we review conjugate gradient method and then present a new method for solving linear interval systems based on conjugate gradient method. Another new iterative method using steepest descent idea is proposed in Section 3.4. Section 4 shows the experimental results of the proposed methods and discusses the accuracy and efficiency of the new methods.

#### 2. Interval Arithmetic

Interval numbers and arithmetic are explained in [7, 11, 13, 20]. We review main definition here. Given , where , the real bounded set is called a* proper interval*. The set of all proper intervals on is denoted by . In this paper, all elements in are shown with a hat, i.e., .

The* infimum* and* supremum* of are and , respectively.

The* magnitude* is defined to be the maximum value of for all . ThusThe width and midpoint of are denoted by and , respectively, and defined asAn interval number is defined as a subset of interval number when and and denoted by . Equality held when .

Each real number can be viewed as a special interval number . The interval simply can be denoted by without confusion.

For each binary operation which is defined on , a binary operation is defined on asfor . Note that we will use operator for both interval and real numbers when there is no confusion.

Therefore, basic interval arithmetic is defined in the following. Let ; thenwhereNote that there is a simple alternative way to formulate the interval multiplication.

An* interval vector* is a vector whose elements are interval numbers. Similarly, an* interval matrix* is a matrix whose elements are interval numbers.

Consider , , , and . We say if each element of vector belongs to corresponding element of and if each element of matrix belongs to the corresponding element of .

Norm of interval vector is defined in [21] asWe define a square interval matrix as* diagonally magnitude dominant* or for simplicity* diagonally dominant* if

Proposition 1. *Interval matrix is diagonally dominant if and only if*

*Proof. *Equation (13) directly follows (12); just set and when .

If (13) holds then andFinally,General matrix form of linear interval systems can be written as follows:where is an interval matrix, is an interval vector solution, and is an interval vector.

To solve a linear interval system, there exist some different approaches [22–24]. Four of the well-known approaches are as follows:(1)Find an interval vector such that, for all , and exist such that (*united solution set*).(2)Find an interval vector such that, for all , for all , exists such that (*tolerable solution set*). Equivalently,(3)Find an interval vector such that, for all , for all , exists, such that (*controllable solution set*). Equivalently,(4)Find an interval vector such that multiplication is equal to (*exact solution set*). United solutions are studied in [13, 25, 26]. In [22, 27] methods have been published to find tolerable solutions. Also, some methods have been developed to obtain controllable solutions [14, 28] and exact solutions [12, 14, 29–31] of interval linear systems.

Our proposed method calculates a solution for in the sense of the exact solution, where is a diagonally dominant interval matrix.

#### 3. Materials and Methods

##### 3.1. Inverse Operations for “+” and “”

Suppose we are given an equation , where is a known interval number, and is an unknown interval number, how can one find ?

Obviously is the solution to this equation. However, which is not equal to unless . If the operation “−” is inverse of the operation “+,” then must equal . But and unless .

As it can be observed, the operator “−” is not inverse of the operation “+.” Similarly, it can be observed that the operation “” is not inverse of the operation “.” Therefore, we propose the inverse of operations “+” and “.”

Before introducing inverse operations for “+” and “,” we need another general type of intervals that is defined in [32, 33]. A* general closed interval * is identified by two real numbers and is defined asThe* general interval number space over real numbers* is denoted by . Obviously, . All definitions (supremum, infimum, and magnitude) and operations in are defined the same as but computation over is more complicated than . The arithmetic over is well discussed in [20, 32, 33]. For example, as in (7), or as in (8). The width of a general interval number can be negative; for example, . Also, according to definition of subset, .

In this paper, by* general interval,* we mean general closed intervals.

For , the operation “” is defined in [20, 33] asFor , we define the operation “” asNote that this definition can handle many cases where belongs to , and just in three cases compute undefined or unbounded solution. This property lets us design algorithms with high stability and consistent when working with interval values that contain .

In the following two theorems, we show that operations “” and “” are inverse of operations “+” and “,” respectively [34].

Theorem 2. *Given , then the interval equation has a unique solution , where . The equation has a unique solution in if ; otherwise the equation has no solution in but has exactly one in .*

*Proof. * That is, , which solves . Similarly, ; therefore, we have , where .

Therefore if and , then can be found with this inverse operation uniquely.

Theorem 3. *Given , then the interval equation has a solution , where . The equation has a solution in if ; otherwise the equation has no solution in but has in .*

*Remark*. Solution is not usually unique; e.g., the interval equation has solutions . Actually, this equation has unlimited solutions. All intervals and where and can be solution of this equation.

In cases that is not unique, is considered as the longest interval as possible. In the above example, must be .

*Proof. *Suppose . Before starting the proof, note that multiplication result, , depends on the sign of and . Table 1 shows sign of related to sign of and and also shows formulas of multiplication result.

Suppose that, given , , and and . Now, from Table 1, it is observed that the following must occur: and . ThenThis means in this case is uniquely calculated by inverse operations of “.” For other cases except cases with condition , one can provide similar proofs.

For cases with condition , there is an ambiguity. From Table 1, there are three cases in column . They have the same conditions . The solution can be in forms , , and , respectively. The solution should be in form of , because it is the longest possible interval solution and includes both of the others. Now, find and . If with condition is the solution of , then . It follows thatandThen,If , then , and . This contradicts with . HenceSimilarly,and therefore, we have .

In addition, special case should be considered separately. If and , no exists such that . If , then is free and expression is consistent for all . At this point is considered as .