Solving Arbitrary Coupled Trapezoidal Fully Fuzzy Sylvester Matrix Equation with Necessary Arithmetic Multiplication Operations

Abstract

A couple of Sylvester matrix equations (CSME) are required to be solved simultaneously in many applications, especially in analysing the stability of control systems. However, there are some situations in which the crisp CSME are not well equipped to deal with the uncertainty problem during the stability analysis of control systems. Thus, this paper proposes a new method for solving Arbitrary Coupled Trapezoidal Fully Fuzzy Sylvester Matrix Equation (ACTrFFSME). New arithmetic fuzzy multiplication operations are developed and applied to convert the ACTrFFSME to a reduced system of non-linear equations. Then it is converted to a system of absolute equations where the arbitrary fuzzy solution is obtained by solving that system. The proposed method can solve many arbitrary fuzzy equations, such as fully fuzzy matrix equations, Sylvester and Lyapunov fully fuzzy matrix equations with triangular and trapezoidal fuzzy numbers without any restrictions. We illustrate the proposed method by solving two numerical examples.

KEYWORDS:

• Arbitrary fuzzy systems
• couple Sylvester matrix equations
• near-zero fuzzy numbers
• trapezoidal fuzzy numbers
• trapezoidal fuzzy multiplication

Subject classifications codes:

• 03B52
• 94D05
• 08A72
• 15B15

1. Introduction

Sylvester matrix equation (SME) in the form $AX+XB=C$ has many important applications in the design and analysis of linear control systems [1], observer design [2], reduction of large-scale dynamical systems [3], restoration of noisy images [4,5], medical imaging data acquisition, model reduction [6] and stochastic control, image processing and filtering [5].

However, there are many applications where CSME in the form $\{\begin{array}{c}AX+YB=C\\ DX+YE=F\end{array}$ are required to be solved simultaneously. The CSME is important in making the computational process less complicated, especially in analysing the stability of control systems so that the control system always performs well according to its specifications [7]. Researchers for many years have proposed many analytical and numerical methods for solving the CSME with crisp numbers [8–13]. However, in many applications, some system parameters are represented by fuzzy numbers rather than crisp numbers due to uncertainty problems such as conflicting requirements during the system process, the distraction of any elements and noise.

Fuzzy logic has been studied since the 1920s, as infinite valued logic by Lukasiewicz and Tarski [14], and the fuzzy set theory was introduced by Lotfi Zadeh [15] in 1965 while the set theory was developed by Georg Cantor [16]. Fuzzy Relation Equations (FREs) with the max–min composition was first studied by Sanchez [17]. The theory and applications of FREs can be found in Di Nola et al. [18], which indicated that if the solvability of max-continuous t-norm FREs is assumed, then the solution set for the FREs can be fully determined from a unique greatest solution and all minimal solutions, and the number of minimal solutions is always finite. Since then, FREs based on various compositions have been investigated. Some common compositions include max–min [19–24], max-product [25–28], max-Archimedean t-norm [29, 30], u-norm [31], max t-norm [32] and max arithmetic mean [33,34]. The conditions for the existence of a solution to the inverse problem concerned with FRE are investigated in [35], a finite system of FREs with sup-T composition was studied in [36], and a system of FREs was investigated in [37, 38].

When all parameters of the CSME are in the fuzzy form, it is called coupled fully fuzzy Sylvester matrix equation (CFFSME).

Definition 1.1:

A matrix equation that can be written as (1.1) $\{\begin{array}{c}\tilde{A}\tilde{X}+\tilde{Y}\tilde{B}=\tilde{C}\\ \tilde{D}\tilde{X}+\tilde{Y}\tilde{E}=\tilde{F}\end{array}$(1.1) is called an arbitrary coupled trapezoidal fully fuzzy Sylvester matrix equation (ACTrFFSME) if $\tilde{A}=({\tilde{a}}_{ij}{)}_{m\times n}$, $\tilde{B}=({\tilde{b}}_{ij}{)}_{r\times p}$, $\tilde{C}=({\tilde{c}}_{ij}{)}_{m\times p}$, $\tilde{D}=({\tilde{d}}_{ij}{)}_{m\times n}$, $\tilde{E}=({\tilde{e}}_{ij}{)}_{r\times p}$, $\tilde{F}=({\tilde{f}}_{ij}{)}_{m\times p}$, $\tilde{X}=({\tilde{x}}_{ij}{)}_{n\times p}$ and $\tilde{Y}=({\tilde{y}}_{ij}{)}_{m\times r}$ are arbitrary trapezoidal fuzzy matrices.

Equation (1.1) is of interest in many different applications. But until now, there is less study for the arbitrary solution of this equation. Most of the methods are proposed for its special cases, such as fully fuzzy Sylvester matrix equations, fully fuzzy continuous-time Lyapunov matrix equations and fully fuzzy matrix equations.

Definition 1.2:

The fully fuzzy matrix equation that can be written as (1.2) $\tilde{A}\tilde{X}+\tilde{X}\tilde{B}=\tilde{E}$(1.2) where, $\tilde{A}=({\tilde{a}}_{ij}{)}_{p\times p}$, $\tilde{B}=({\tilde{b}}_{ij}{)}_{n\times n}$, $\tilde{X}=({\tilde{x}}_{ij}{)}_{p\times n}$ and $\tilde{E}=({\tilde{e}}_{ij}{)}_{p\times n}$. is called Fully Fuzzy Sylvester Matrix Equation (FFSME).

Definition 1.3:

The fully fuzzy matrix equation that can be written as (1.3) $\tilde{A}\tilde{X}+\tilde{X}{\tilde{A}}^T=\tilde{E}$(1.3) where, $\tilde{A}=({\tilde{a}}_{ij}{)}_{p\times p}$, ${\tilde{A}}^T=({\tilde{a}}_{ij}^T{)}_{p\times p}$, $\tilde{X}=({\tilde{x}}_{ij}{)}_{p\times p}$ and $\tilde{E}=({\tilde{e}}_{ij}{)}_{p\times p}$. is called a Fully Fuzzy Continuous-Time Lyapunov Matrix Equation (FFCTLME).

Definition 1.4:

The fully fuzzy matrix equation that can be written as (1.4) $\tilde{A}\tilde{X}=\tilde{E}$(1.4) where, $\tilde{A}=({\tilde{a}}_{ij}{)}_{m\times p}$, $\tilde{X}=({\tilde{x}}_{ij}{)}_{p\times n}$ and $\tilde{E}=({\tilde{e}}_{ij}{)}_{m\times n}$ is called a Fully Fuzzy Matrix Equation (FFME). To the best of our knowledge, the ACFFSME in Equation (1.1) is not solved by any known method. Thus, researchers restrict themselves to some of its special cases. For example, the authors in [39–42] considered the solvability for a Fully Fuzzy Linear System (FFLS) and authors in [43,44] considered the FFME in Equation (1.4) with triangular fuzzy numbers (TFNs). However, these methods cannot be extended to FFME with trapezoidal fuzzy numbers (TrFNs). In addition, several analytical methods have been proposed for triangular FFSME (TFFSME). The most important of which is the Kronecker product and matrix inversion [45], Kronecker product and associated linear system [46–48]. However, the methods were restricted only to positive triangular fuzzy numbers (TFNs) only. Positive and negative solutions for trapezoidal FFSME (TrFFSME) are considered in [49]. However, the method cannot be extended to arbitrary TrFFSME. A few studies have been conducted for solving a pair of fuzzy matrix equations. [50] proposed a method for solving a pair of fuzzy matrix equations in the form: $\{\begin{array}{c}A\tilde{X}+\tilde{X}B=\tilde{C}\\ D\tilde{X}E=\tilde{F}\end{array}$

In addition, [51, 52] proposed analytical methods for solving a pair of fully fuzzy matrix equations (PFFME) in the form: $\{\begin{array}{c}\tilde{A}\tilde{X}+\tilde{X}\tilde{B}=\tilde{C}\\ \tilde{D}\tilde{X}\tilde{E}=\tilde{F}\end{array}$

In that study, a direct method was proposed to solve the positive PFFME with arbitrary coefficients by applying the Kronecker product and Vec-operator. However, both methods obtained a positive fuzzy solution only. Recently, The positive fuzzy solution to the CTrFFSME was obtained by El Sayed, Malkawi and Ahmad [53]. However, this method cannot be applied to ACTrFFSME. Therefore, limitations of the existing methods are pointed out.

1. Existing methods for solving TFFSME and TrFFSME can only obtain positive or negative fuzzy solutions and cannot be extended to fuzzy systems with arbitrary solutions. Consequently, the existing methods are limited to some real-life situations. However, there exist several instances in real scenarios where the coefficients may not be entirely positive or nonnegative.

2. Analytical methods proposed for solving different arbitrary systems are based on Vec-operator and Kronecker product. It is worth mentioning that the Vec-operator and Kronecker product approach is applicable for fuzzy systems with positive or negative fuzzy numbers only and cannot be applied to fuzzy systems with near-zero fuzzy numbers. In addition, the Kronecker product method for $m\times n$ fuzzy system requires getting the inverse of $mn\times mn$ matrix, which is impossible for large systems.

3. Using the existing methods, it is not possible to check whether the obtained solution is unique or not. There may be several cases where the fuzzy systems generate unique or infinitely many solutions.

4. The existing methods can only be applied to a certain fuzzy number and cannot be extended to other types of fuzzy numbers.

Consequently, the developed methods in the literature limit the size and sign of fuzzy systems. Therefore, it is important to develop a new method without considering Vec-operator and Kronecker product for solving the CTrFFSME and its special cases. Therefore, to deal with this shortcoming, in this paper, reduced fuzzy multiplication operations are applied to convert the arbitrary CTrFFSME into a system of non-linear equations. The reduced system is converted to an equivalent system of absolute equations where the fuzzy solutions are obtained by solving that system.

The proposed method is applicable for solving ACTrFFSME with large-size matrices. In addition, it can also be applied to its special cases in Equations (1.2)–(1.4) with both TFNs and TrFNs without any restriction. This paper is organised as follows: Section 2 introduces preliminary arithmetic operations of trapezoidal fuzzy numbers. In Section 3, a proposed method for solving ACTrFFSME is developed, along with a presentation of its algorithm. In Section 4, a numerical example is presented to illustrate the proposed method. Section 5 is dedicated to the conclusion.

2. Preliminaries

The following are basic definitions and results related to TrFNs in fuzzy theory [54–56] and fuzzy matrix [57,58].

Definition 2.1:

Let $X$ be a universal set. Then, the fuzzy subset $\tilde{A}$ of $X$ is defined by its membership function ${\mu }_{\tilde{A}}:X\to [0,1]$ which assigns to each element $x\in X$ a real number ${\mu }_{\tilde{A}}(x)$ in the interval $[0,1]$, where the function value of ${\mu }_{\tilde{A}}(x)$ represents the grade of membership of $x$ in $\tilde{A}.$ A fuzzy set $\tilde{A}$ is written as $\tilde{A}=\{(x,{\mu }_{\tilde{A}}(x)),x\in X,{\mu }_{\tilde{A}}(x)\in [0,1]\}$.

Definition 2.2:

A fuzzy set $\tilde{A}$, defined on the universal set of real number $R$, is said to be a fuzzy number if its membership function has the following characteristics:

1. $\tilde{A}$ is convex, i.e. ${\mu }_{\tilde{A}}(\lambda x_1+(1-\lambda )x_2\ge min({\mu }_{\tilde{A}}(x),{\mu }_{\tilde{A}}(x))\mathrm{\forall }{x}_1,x_2\in R,\mathrm{\forall }\lambda \in [0,1].$

2. $\tilde{A}$ is normal, i.e. $\mathrm{\exists }{x}_0\in R\text{such}\text{that}{\mu }_{\tilde{A}}(x_0)=1.$

3. ${\mu }_{\tilde{A}}$ is piecewise continuous.

Definition 2.3:

A fuzzy number $\tilde{A}=(a_1,a_2,a_3,a_4)$ is a TrFN if its membership function is: ${\mu }_{\tilde{A}}(x)=\{\begin{array}{ll}0& x<a_1\\ \frac{x-a_1}{a_2-a_1}& a_1\le x\le a_2\\ 1& a_2\le x\le a_3\\ \frac{a_4-x}{a_4-a_3}& a_3\le x\le a_4\\ 0& x>a_4\end{array}$ Figure 1 representation of TrFN in the form ($a_1,a_2,a_3,a_4$).

Figure 1. Representation of TrFN ($a_1,a_2,a_3,a_4$).

Definition 2.4:

The sign of the TrFN $\tilde{A}=(a_1,a_2,a_3,a_4)$ can be classified as:

• $\tilde{A}$ is positive (negative) iff $a_1>0,(a_4<0).$

• $\tilde{A}$ is zero iff $(a_1,a_2,a_{3}and{a}_4=0).$

• $\tilde{A}$ is near zero iff $a_1\le 0\le a_4.$

Definition 2.5:

Operations of TrFNs.

The arithmetic operations of TrFNs are presented as follows. Let $\tilde{A}=(a_1,a_2,a_3,a_4)$, $\tilde{B}=(b_1,b_2,b_3,b_4)$ be two TrFNs then:

1. Addition (2.1) $\tilde{A}+\tilde{B}=(a_1+b_1,a_2+b_2,a_3+b_3,a_4+b_4).$(2.1)

2. Subtraction $\tilde{A}-\tilde{B}=(a_1-b_4,a_2-b_3,a_3-b_2,a_4-b_1).$

3. Symmetric image $-\tilde{A}=(-a_4,-a_3,-a_2,-a_1).$

4. Scalar multiplication: Let $\lambda \in \mathbb{R}$ then, $\lambda \otimes (a_1,a_2,a_3,a_4)=\{\begin{array}{cc}(\lambda a_1,\lambda a_2,\lambda a_3,\lambda a_4),& \lambda \ge 0\\ (\lambda a_4,\lambda a_3,\lambda a_2,\lambda a_1).& \lambda <0\end{array}$

5. Multiplication: The multiplication between fuzzy numbers is neither commutative nor associative. Thus, Arithmetic Multiplication Operations (AMO) [59] for TrFNs can be classified as follows:

Let $\tilde{A}=(a_1,a_2,a_3,a_4)$, $\tilde{B}=(b_1,b_2,b_3,b_4)$ be TrFNs,

Case (I) If $\tilde{A}$, $\tilde{B}$ both arbitrary TrFNs then: (2.2) $\tilde{A}\tilde{B}=(a,h,m,d)$(2.2) where $\begin{array}{rl}a& =min(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4),\end{array}$ $\begin{array}{rl}h& =min(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3),\end{array}$ $\begin{array}{rl}m& =\text{max}(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3),\end{array}$ $\begin{array}{rl}d& =max(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4).\end{array}$ Case (II) If $\tilde{A}$, $\tilde{B}$ both restricted TrFNs then: (2.3a) $\begin{array}{rl}\tilde{A}\tilde{B}& =(a_1{b}_1,a_2{b}_2,a_3{b}_3,a_4{b}_4)\tilde{A},\tilde{B}>0,\end{array}$(2.3a) (2.3b) $\begin{array}{rl}\tilde{A}\tilde{B}& =(a_4{b}_4,a_3{b}_3,a_2{b}_2,a_1{b}_1)\tilde{A},\tilde{B}<0,\end{array}$(2.3b) (2.3c) $\begin{array}{rl}\tilde{A}\tilde{B}& =(a_4{b}_1,a_3{b}_2,a_2{b}_3,a_1{b}_4)\tilde{A}>0,\tilde{B}<0,\end{array}$(2.3c) (2.3d) $\begin{array}{rl}\tilde{A}\tilde{B}& =(a_1{b}_4,a_2{b}_3,a_3{b}_2,a_4{b}_1)\tilde{A}<0,\tilde{B}>0.\end{array}$(2.3d)

• (VI) Equality: The fuzzy numbers $\tilde{A}=(a_1,a_2,a_3,a_4)$ and $\tilde{B}=(b_1,b_2,b_3,b_4)$ are equal iff (2.4) $a_1=b_1,a_2=b_2,a_3=b_{3}and{a}_4=b_4$(2.4)

Definition 2.6:

A matrix $\tilde{A}=({\tilde{a}}_{ij}{)}_{m\times n}$ is called a trapezoidal fuzzy matrix, if each element of $\tilde{A}$ is a TrFN.

Definition 2.7:

A fuzzy matrix $\tilde{A}$ will be:

1. Positive (negative) and denoted by $\tilde{A}>0$, $(\tilde{A}<0)$ if each element of $\tilde{A}$ is positive (negative) TrFN.

2. Nonnegative (non-positive) and denoted by $\tilde{A}\ge 0$, $(\tilde{A}\le 0)$ if each element of $\tilde{A}$ is non-negative (non-positive) TrFNs.

3. Arbitrary, if at least one element of $\tilde{A}$ is near zero TrFNs.

Definition 2.8

[60]: A TrFN $\tilde{A}=(a_1,a_2,a_3,a_4)$ is said to be near-zero TrFN if $a_1<0<a_4$, and can be classified as follows: (2.5a) $\begin{array}{rl}\text{I})\tilde{A}& =(a_1,a_2,a_3,a_4)\text{is called}{N}_1-zero\text{TrFN iff}{a}_1\le a_2\le a_3<0<a_4.\end{array}$(2.5a) (2.5b) $\begin{array}{rl}\text{II})\tilde{A}& =(a_1,a_2,a_3,a_4)\text{is called}{N}_2-zero\text{TrFN iff}{a}_1\le a_2\le a_3<0<a_4.\end{array}$(2.5b) (2.5c) $\begin{array}{rl}\text{III})\tilde{A}& =(a_1,a_2,a_3,a_4)\text{is called}{N}_3-zero\text{TrFN iff}{a}_1\le a_2\le a_3<0<a_4.\end{array}$(2.5c)

Definition 2.9:

The Vec-operator generates a column vector from a matrix $A$ by stacking the column vectors of $A=\left(\begin{array}{ccc}{a}_{11}& \cdots & a_{1n}\\ ⋮& \ddots & ⋮\\ a_{n1}& \cdots & a_{nn}\end{array}\right)$ as, $Vec(A)=\left(\begin{array}{c}{a}_{11}\\ a_{21}\\ ⋮\\ a_{nn}\end{array}\right).$ In addition, If,hen $A=\left(\begin{array}{ccc}{a}_{11}& \cdots & a_{1n}\\ ⋮& \ddots & ⋮\\ a_{n1}& \cdots & a_{nn}\end{array}\right)$.

Definition 2.10:

For any integers $x$ and $y$, $\text{min}(x,y)$ and $\text{max}(x,y)$ denote the minimum and maximum of $x$ and $y$, respectively as follows [41], $min(a,b)=\left(\frac{a+b}{2}-\frac{|a-b|}{2}\right),max(a,b)=\left(\frac{a+b}{2}+\frac{|a-b|}{2}\right).$

Definition 2.11:

An absolute system (or A system of absolute equations) is a collection of equations such that one of them at least is an absolute equation.

Definition 2.12:

The solution set of a system is called a finite solution or alternative solution, whereby the number of solutions is more than one and not infinite solutions.

Remark 2.1:

Two TrFNs are said to be restricted if their signs are either positive or negative.

Remark 2.2:

Two TrFNs are said to be semi-restricted if the sign of one of them is either positive or negative and the other one is arbitrary.

The AMO in Definition 2.5, Equation (2.2) can find the product of arbitrary TrFNs. However, when applying these operations to solve arbitrary fuzzy systems, the fuzzy systems are converted to a min–max non-linear system which is very challenging to be solved. Therefore, reducing the min–max non-linear system from four terms to two terms makes the solution to the system much easier in terms of computational timing and memory usage. Thus, a further modification to the AMO in Definition 2.5, Equation (2.2) need to be done to make it more practical in solving arbitrary fuzzy systems. Therefore, in the next Section 3, AMO is reduced to what so-called Reduced Arithmetic Multiplication Operations (RAMO) based on the sign of the TrFNs that are positive, negative or near-zero.

3. Reduced Multiplication Operations for Semi-Restricted TrFNs

In this section, the fuzzy arithmetic multiplication operations for arbitrary TrFNs in Equation (2.2) are reduced for semi-restricted TrFNs. In the following corollaries, the multiplication between semi-restricted TrFNs is reduced from four terms into two terms only. This reduction contributes significantly to the solution of the arbitrary CTrFFSME in Equation (1.1) in Section 4.

Corollary 3.1:

Suppose $\tilde{A}=(a_1,a_2,a_3,a_4)$ is positive TrFN and $\tilde{B}=(b_1,b_2,b_3,b_4)$ is arbitrary TrFN then: (3.1) $\tilde{A}\tilde{B}=(min(a_1{b}_1,a_4{b}_1),min(a_2{b}_2,a_3{b}_2),max(a_2{b}_3,a_3{b}_3),max(a_1{b}_4,a_4{b}_4))$(3.1)

Proof:

Given $\tilde{B}$ is arbitrary TrFN, consequently, its sign could be positive, negative or near-zero respectively. Thus, the proof for this corollary is done in three steps as follows:

1. Both $\tilde{A}$ and $\tilde{B}$ are positive. Then by Equation (2.3a), the following is obtained, $\tilde{A}\tilde{B}=(a_1{b}_1,a_2{b}_2,a_3{b}_3,a_4{b}_4).$

2. If $\tilde{A}$ is positive and $\tilde{B}$ is negative. Then by Equation (2.3c), the following is obtained, $\tilde{A}\tilde{B}=(a_4{b}_1,a_3{b}_2,a_2{b}_3,a_1{b}_4).$

3. If $\tilde{A}$ is positive and $\tilde{B}$ is near-zero. Then by Definition 2.7 the product $\tilde{A}\tilde{B}$ is classified as follows:

Case (I) Let $\tilde{A}$ be positive TrFN and $\tilde{B}$ be $N_1-zero$ TrFN then, $0<a_1<a_4\text{and}{b}_1<0<b_4.$

Therefore $a_4{b}_1<a_1{b}_1<0$ and $0<a_1{b}_4<a_4{b}_4$.

Consequently, (3.2a) $a_4{b}_1<a_1{b}_1<0<a_1{b}_4<a_4{b}_4.$(3.2a)

In addition, since $0<a_2\le a_3$ and $b_2<b_3<0$, then, (3.2b) $a_3{b}_2<a_2{b}_2<a_2{b}_3<0\text{and}{a}_3{b}_2<a_3{b}_3<a_2{b}_3<0.$(3.2b)

By Equations (3.2a) and (3.2b), AMO in Equation(2.2) is reduced as follows: $\begin{array}{rl}a& =min(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4)=a_4{b}_1,\end{array}$ $\begin{array}{rl}h& =min(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3)=a_3{b}_2,\end{array}$ $\begin{array}{rl}m& =max(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3)=a_2{b}_3,\end{array}$ $\begin{array}{rl}d& =max(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4)=a_4{b}_4.\end{array}$

Therefore, if $\tilde{A}$ is positive TrFN and $\tilde{B}$ is $N_1-zero$ then: (3.2c) $\tilde{A}\tilde{B}=(a_4{b}_1,a_3{b}_2,a_2{b}_3,a_4{b}_4).$(3.2c)

Case (II) Let $\tilde{A}$ be positive TrFN and $\tilde{B}$ be $N_2-zero$ TrFN then, $0<a_1<a_4\text{and}{b}_1<0<b_4$

Therefore $a_4{b}_1<a_1{b}_1<0$ and $0<a_1{b}_4<a_4{b}_4$ Consequently, (3.3a) $a_4{b}_1<a_1{b}_1<0<a_1{b}_4<a_4{b}_4.$(3.3a)

In addition, since $0<a_2\le a_3$ and $b_2<0<b_3$. $a_3{b}_2\le a_2{b}_2<0$ and $0<a_2{b}_3<a_3{b}_3$. Consequently, (3.3b) $a_3{b}_2\le a_2{b}_2<0<a_2{b}_3<a_3{b}_3$(3.3b)

By Equations (3.3a) and (3.3b), AMO in Equation (3.1) is reduced as follows: $\begin{array}{rl}a& =min(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4)=a_4{b}_1,\end{array}$ $\begin{array}{rl}h& =min(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3)=a_3{b}_2,\end{array}$ $\begin{array}{rl}m& =\text{max}(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3)=a_3{b}_3,\end{array}$ $\begin{array}{rl}d& =\text{max}(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4)=a_4{b}_4.\end{array}$

Therefore, if $\tilde{A}$ is positive TrFN and $\tilde{B}$ is $N_2-zero$ TrFN then: (3.3c) $\tilde{A}\tilde{B}=(a_4{b}_1,a_3{b}_2,a_3{b}_3,a_4{b}_4)$(3.3c)

Case (III) Let $\tilde{A}$ be positive TrFN and $\tilde{B}$ is $N_3-zero$ TrFN then, $0<a_1<a_4\text{and}{b}_1<0<b_4$

Therefore $a_4{b}_1<a_1{b}_1<0$ and $0<a_1{b}_4<a_4{b}_4$ Consequently (3.4a) $a_4{b}_1<a_1{b}_1<0<a_1{b}_4<a_4{b}_4.$(3.4a)

In addition, since $0<a_2\le a_3$ and $0<b_2<b_3$, then, (3.4b) $0<a_2{b}_2<a_2{b}_3<a_3{b}_3,\text{and}0<a_2{b}_2<a_3{b}_2<a_3{b}_3$(3.4b)

By Equations (3.4a) and (3.4b), AMO in Equation (2.2) is reduced as follows: $\begin{array}{rl}a& =min(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4)=a_4{b}_1,\end{array}$ $\begin{array}{rl}h& =min(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3)=a_2{b}_2,\end{array}$ $\begin{array}{rl}m& =\text{max}(a_2{b}_2,a_2{b}_3,a_3{b}_2,a_3{b}_3)=a_3{b}_3,\end{array}$ $\begin{array}{rl}d& =\text{max}(a_1{b}_1,a_1{b}_4,a_4{b}_1,a_4{b}_4)=a_4{b}_4.\end{array}$

Therefore, if $\tilde{A}$ is positive TrFN and $\tilde{B}$ is $N_3-zero$ then: (3.4c) $\tilde{A}\tilde{B}=(a_4{b}_1,a_2{b}_2,a_3{b}_3,a_4{b}_4).$(3.4c)

By combining Equations (2.3a), (2.3c), (3.2c), (3.3c) and (3.4c), If $\tilde{A}$ is positive TrFN and $\tilde{B}$ is arbitrary TrFN, then the product $\tilde{A}\tilde{B}$ is: $\tilde{A}\tilde{B}=(min(a_1{b}_1,a_4{b}_1),min(a_2{b}_2,a_3{b}_2),max(a_2{b}_3,a_3{b}_3),max(a_1{b}_4,a_4{b}_4)).$

Example 3.1:

Let $\tilde{\mathrm{A}}=(1,2,4,5)$ and $\tilde{\mathrm{B}}=(-6,-4,6,7)$ two TrFNs respectively, then $\widetilde{\mathrm{AB}}$ is found using the RAMO in Equation (3.1) as follows, $\begin{array}{rl}\widetilde{\mathrm{AB}}& =(\text{min}(1\times -6,5\times -6),\text{min}(2\times -4,4\times -4),max(2\times 6,4\times 6),\\ & \text{ max}(1\times 7,5\times 7)).\end{array}$ Therefore, $\tilde{A}\tilde{B}=(-30,-16,24,35).$

Corollary 3.2:

If $\tilde{A}=(a_1,a_2,a_3,a_4)<0$ and $\tilde{B}=(b_1,b_2,b_3,b_4)$ arbitrary fuzzy number, then: (3.5) $\tilde{A}\tilde{B}=(min(a_1{b}_4,a_4{b}_4),min(a_2{b}_3,a_3{b}_3),max(a_3{b}_2,a_2{b}_2),max(a_4{b}_1,a_1{b}_1))$(3.5)

Proof: Given $\tilde{B}$ is arbitrary TrFN, it could be positive, negative or near-zero TrFN respectively. Thus, the proof for this corollary is done in three as follows:

1. Both $\tilde{A}$ and $\tilde{B}$ are negative. Then by Equation (2.3b), the following is obtained, $\tilde{A}\tilde{B}=(a_4{b}_4,a_3{b}_3,a_2{b}_2,a_1{b}_1).$

2. If $\tilde{A}$ is negative and $\tilde{B}$ is positive. Then by Equation (2.3d), the following is obtained, $\tilde{A}\tilde{B}=(a_1{b}_4,a_2{b}_3,a_3{b}_2,a_4{b}_1).$

3. If $\tilde{A}$ is negative and $\tilde{B}$ is near-zero. Then by the definition of near-zero TrFNs in Definition 2.7 the product of $\tilde{A}\tilde{B}$ is classified as follows:

Case (I) If $\tilde{A}$ is negative and $\tilde{B}$ is $N_1-zero$ then: (3.6a) $\tilde{A}\tilde{B}=(a_1{b}_4,a_3{b}_3,a_2{b}_2,a_4{b}_1).$(3.6a)

Proof: Since $\tilde{A}$ is negative TrFN and $\tilde{B}$ is $N_1-zero$ TrFN then, $a_1<a_4<0\text{and}{b}_1<0<b_4.$

Thus, $0<a_1{b}_1$, $0<a_4{b}_1$, $a_1{b}_4<0$ and $a_4{b}_4<0$.

In addition, $a_4{b}_1>a_1{b}_1>0$ and $a_1{b}_4<a_4{b}_4<0$. Therefore, (3.6b) $a_1{b}_4<a_4{b}_4<0<a_1{b}_1<a_4{b}_1$(3.6b)

In addition, since $a_2\le a_3<0$ and $b_2<b_3<0$. Therefore, (3.6c) $0<a_3{b}_3<a_3{b}_2<a_2{b}_2\text{and}0<a_3{b}_3<a_2{b}_3<a_2{b}_2$(3.6c)

By Equations (3.6b) and (3.6c), the following is obtained $\tilde{A}\tilde{B}=(a_4{b}_1,a_3{b}_3,a_2{b}_2,a_1{b}_4).$

Case (II) If $\tilde{A}$ is negative and $\tilde{B}$ is $N_2-zero$ then: (3.7) $\tilde{A}\tilde{B}=(a_1{b}_4,a_2{b}_3,a_2{b}_2,a_1{b}_1).$(3.7)

Proof: Since $\tilde{A}$ is negative TrFN and $\tilde{B}$ is $N_2-zero$ TrFN then,

$a_1<a_4<0$ and $b_1<0<b_4$. Therefore $0<a_4{b}_1<a_1{b}_1$ and $a_1{b}_4<a_4{b}_4<0$ Consequently, (3.7a) $a_1{b}_4<a_4{b}_4<0<a_4{b}_1<a_1{b}_1$(3.7a)

In addition, since $a_2\le a_3<0$ and $b_2<0<b_3$. $0<a_3{b}_2\le a_2{b}_2$ and $a_2{b}_3<a_3{b}_3<0$. Consequently, (3.7b) $a_2{b}_3\le a_3{b}_3<0<a_3{b}_2<a_2{b}_2$(3.7b)

By Equations (3.7a) and (3.7b), the following is obtained $\tilde{A}\tilde{B}=(a_1{b}_4,a_1{b}_1,a_2{b}_3,a_2{b}_2).$

Case (III) If $\tilde{A}$ is positive and $\tilde{B}$ is $N_3-zero$ then: (3.8) $\tilde{A}\tilde{B}=(a_1{b}_4,a_2{b}_3,a_3{b}_2,a_1{b}_1).$(3.8)

Proof: Since $\tilde{A}$ is positive TrFN and $\tilde{B}$ is $N_1-zero$ TrFN then,

$a_1<a_4<0$ and $b_1<0<b_4$. Therefore $0<a_4{b}_1<a_1{b}_1$ and $a_1{b}_4<a_4{b}_4<0$ Consequently, (3.8a) $a_1{b}_4<a_4{b}_4<0<a_4{b}_1<a_1{b}_1$(3.8a)

In addition, since $a_2\le a_3<0$ and $0<b_2<b_3$. Then, $a_2{b}_3<a_2{b}_2<a_3{b}_2$ and $a_2{b}_3<a_2{b}_2<a_3{b}_3$. In addition, $a_3{b}_3<a_3{b}_2$. Consequently, (3.8b) $a_2{b}_3\le a_2{b}_2<a_3{b}_3<a_3{b}_2$(3.8b)

By Equations (3.8a) and (3.8b), the following is obtained $\tilde{A}\tilde{B}=(a_1{b}_4,a_2{b}_3,a_3{b}_2,a_1{b}_1).$

By Equations (2.3b), (2.3d), (3.6a), (3.7) and (3.8) the following is obtained. $\tilde{A}\tilde{B}=(min(a_1{b}_4,a_4{b}_4),min(a_2{b}_3,a_3{b}_3),max(a_3{b}_2,a_2{b}_2),max(a_4{b}_1,a_1{b}_1)).$

Example 3.2:

Let $\tilde{\mathrm{A}}=(-11,-7,-4,-2)$ and $\tilde{\mathrm{B}}=(-1,4,5,7)$ two TrFNs respectively, then $\widetilde{\mathrm{AB}}$ is found using the RAMO in Equation (3.5) as follows, $\begin{array}{rl}\widetilde{\mathrm{AB}}& =(\text{min}(-11\times 7,-2\times 7),\text{min}(-7\times 5,-4\times 5),max(-4\times 4,-7\times 4),\\ & \text{ max}(-2\times -1,-11\times -1)).\end{array}$

Therefore, $\tilde{A}\tilde{B}=(-77,-35,-16,11).$

Corollary 3.3:

If $\tilde{A}$ is near-zero TrFN and $\tilde{B}$ arbitrary TrFN and based on the definition of near-zero TrFNs in Definition 2.7 the product $\tilde{A}\tilde{B}$ can be classified as follow:

Case (I) If $\tilde{A}$ is $N_1-zero$ TrFN and $\tilde{B}$ arbitrary TrFN then: (3.9) $\tilde{A}\tilde{B}=(min(a_1{b}_4,a_4{b}_1),min(a_2{b}_3,a_3{b}_3),max(a_3{b}_2,a_2{b}_2),max(a_4{b}_4,a_1{b}_1))$(3.9) Case (II) If $\tilde{A}$ is $N_2-zero$ TrFN and $\tilde{B}$ arbitrary TrFN then: (3.9a) $\tilde{A}\tilde{B}=(min(a_1{b}_4,a_4{b}_1),min(a_2{b}_3,a_3{b}_2),max(a_2{b}_2,a_3{b}_3),max(a_4{b}_4,a_1{b}_1))$(3.9a) Case (III) If $\tilde{A}$ is $N_3-zero$ TrFN and $\tilde{B}$ arbitrary TrFN then: (3.9b) $\tilde{A}\tilde{B}=(min(a_1{b}_4,a_4{b}_1),min(a_2{b}_2,a_3{b}_2),max(a_2{b}_3,a_3{b}_3),max(a_4{b}_4,a_1{b}_1))$(3.9b)

Proof: Straightforward similar to Corollary 3.1 and 3.2.

Corollary 3.4:

If $\tilde{A}=(a_1,a_2,a_3,a_4)$ arbitrary TrFN and $\tilde{B}>0$ then: (3.10) $\tilde{A}\tilde{B}=(min(a_1{b}_1,a_1{b}_4),min(a_2{b}_2,a_2{b}_3),max(a_3{b}_2,a_3{b}_3),max(a_4{b}_1,a_4{b}_4)$(3.10)

Proof: Straightforward similar to Corollary 3.1 and 3.2.

Corollary 3.5:

If $\tilde{A}=(a_1,a_2,a_3,a_4)$ arbitrary TrFN and $\tilde{B}<0$ then: (3.11) $\tilde{A}\tilde{B}=(min(a_4{b}_1,a_4{b}_4),min(a_3{b}_2,a_3{b}_3),max(a_2{b}_3,a_2{b}_2),max(a_1{b}_4,a_1{b}_1)$(3.11) Proof: Straightforward similar to Corollary 3.1 and 3.2.

Corollary 3.6:

If $\tilde{A}$ is arbitrary TrFNs and $\tilde{B}$ near-zero TrFN and based on the definition of near-zero TrFN in Definition 2.7 the product $\tilde{A}\tilde{B}$ is classified as follows:

Case (I) If $\tilde{A}$ is arbitrary TrFN and $\tilde{B}$ is $N_1-zero$ TrFN then: (3.12) $\tilde{A}\tilde{B}=(min(a_4{b}_1,a_1{b}_4),min(a_3{b}_2,a_3{b}_3),max(a_2{b}_3,a_2{b}_2),max(a_1{b}_1,a_4{b}_4)$(3.12) Case (II) If $\tilde{A}$ is arbitrary TrFN and $\tilde{B}$ is $N_2-zero$ TrFN then: (3.12a) $\tilde{A}\tilde{B}=(min(a_4{b}_1,a_1{b}_4),min(a_3{b}_2,a_2{b}_3),max(a_2{b}_2,a_3{b}_3),max(a_1{b}_1,a_4{b}_4)$(3.12a) Case (III) If $\tilde{A}$ is arbitrary TrFN and $\tilde{B}$ is $N_3-zero$ TrFN then: (3.12b) $\tilde{A}\tilde{B}=(min(a_4{b}_1,a_1{b}_4),min(a_2{b}_2,a_2{b}_3),max(a_3{b}_2,a_3{b}_3),max(a_1{b}_1,a_4{b}_4)$(3.12b)

Proof: Straightforward similar to Corollary 3.1 and 3.2.

The following Figure 2 summarises the reduced fuzzy multiplication operations for TrFNs.

Figure 2. Summary of the reduced fuzzy multiplication operations for TrFNs.

4. Proposed Method

In the following section, a new method is proposed for solving the CTrFFSME in Equation (1.1).

4.1. Arbitrary Coupled Trapezoidal Fully Fuzzy Sylvester Matrix Equation

In this section, the arbitrary CTrFFSME in Equation (1.1) is converted to an equivalent reduced system of non-linear equations based on the RAMO in Section 3. In addition, the obtained non-linear system is converted to a system of absolute equations. Solving the system of absolute equations will give the fuzzy solution to the arbitrary CTrFFSME in Equation (1.1).

In the following Definition 6.1.2, the system of non-linear equations is introduced.

Definition 4.1:

The system of equations in the form, $\{\begin{array}{c}min(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})+min(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)})=e_{ij}^{(1)},\\ min\left(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)}\right)+min\left(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)}\right)=e_{ij}^{(2)},\\ max\left(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)}\right)+max\left(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)}\right)=e_{ij}^{(3)},\\ max\left(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)}\right)+max\left(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)}\right)=e_{ij}^{(4)},\\ min(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)})+min(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)})=f_{ij}^{(1)},\\ min\left(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)}\right)+min\left(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)}\right)=f_{ij}^{(2)},\\ max\left(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)}\right)+max\left(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)}\right)=f_{ij}^{(3)},\\ max\left(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)}\right)+max\left(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)}\right)=f_{ij}^{(4)}.\end{array}$ is called a system of non-linear equations.

In the following Theorem 4.1, the arbitrary CTrFFSME in Equation (1.1) is converted to an equivalent system of non-linear equations.

Theorem 4.1:

Fundamental theorem of arbitrary coupled trapezoidal fully fuzzy Sylvester matrix equation.

Suppose that $\tilde{A},\tilde{B},\tilde{C},\tilde{D},\tilde{X},\tilde{Y},\tilde{E}and\tilde{F}$are arbitrary trapezoidal fuzzy matrices respectively, then the arbitrary CTrFFSME $\{\begin{array}{c}\tilde{A}\tilde{X}+\tilde{Y}\tilde{B}=\tilde{C}\\ \tilde{D}\tilde{X}+\tilde{Y}\tilde{E}=\tilde{F}\end{array}$ is equivalent to the following system of non-linear equations: (4.1) $\{\begin{array}{c}min(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})+min(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)})=e_{ij}^{(1)},\\ min(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})+min(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)})=e_{ij}^{(2)},\\ max(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})+max(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)})=e_{ij}^{(3)},\\ max(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})+max(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)})=e_{ij}^{(4)},\\ min(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)})+min(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)})=f_{ij}^{(1)},\\ min(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)})+min(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)})=f_{ij}^{(2)},\\ max(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)})+max(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)})=f_{ij}^{(3)},\\ max(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)})+max(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)})=f_{ij}^{(4)}.\end{array}$(4.1)

Proof:

Let $\tilde{A},\tilde{B},\tilde{C},\tilde{D},\tilde{X},\tilde{Y},\tilde{E}\text{and}\tilde{F}$ in the ACTrFFSME $\{\begin{array}{c}\tilde{A}\tilde{X}+\tilde{Y}\tilde{B}=\tilde{E}\\ \tilde{C}\tilde{X}+\tilde{Y}\tilde{D}=\tilde{F}\end{array}$ be arbitrary trapezoidal fuzzy matrices respectively, then the RAMO and EAMO in Sections 3.1.2, 3.1.3 and 3.2 respectively can be applied to obtain ${\tilde{a}}_{ij}{\tilde{x}}_{ij},$ ${\tilde{y}}_{ij}{\tilde{b}}_{ij}$, ${\tilde{c}}_{ij}{\tilde{x}}_{ij}$ and ${\tilde{y}}_{ij}{\tilde{d}}_{ij}$ as follows: $\tilde{A}\tilde{X}=\sum _{k=1}^n{\tilde{a}}_{ik}{\tilde{x}}_{kj}\mathrm{\forall }1<i<m,1<j<p.$ which can be written as $\tilde{A}\tilde{X}=(M_{ij},N_{ij},P_{ij},Q_{ij})$ where, $\begin{array}{rl}{M}_{ij}& =min(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)}),\end{array}$ $\begin{array}{rl}{N}_{ij}& =min(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}{P}_{ij}& =max(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}{Q}_{ij}& =max(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)}).\end{array}$

and, $\tilde{Y}\tilde{B}=\sum _{k=1}^n{\tilde{y}}_{ik}{\tilde{b}}_{kj}\mathrm{\forall }1<i<m,1<j<p.$ which can be written as $\tilde{Y}\tilde{B}=(F_{ij},L_{ij},H_{ij},R_{ij})$ where, $\begin{array}{rl}{F}_{ij}& =min(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)}),\end{array}$ $\begin{array}{rl}{L}_{ij}& =min(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}{H}_{ij}& =max(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}{R}_{ij}& =max(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)}).\end{array}$

Similarly, $\tilde{D}\tilde{X}=\sum _{k=1}^n{\tilde{d}}_{ik}{\tilde{x}}_{kj}\mathrm{\forall }1<i<m,1<j<p.$ which can be written as $\tilde{D}\tilde{X}=(R_{ij},S_{ij},T_{ij},V_{ij})$ where, $\begin{array}{rl}{R}_{ij}& =min(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)}),\end{array}$ $\begin{array}{rl}{S}_{ij}& =min(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}{T}_{ij}& =max(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}{V}_{ij}& =max(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)}).\end{array}$

and, $\tilde{Y}\tilde{E}=\sum _{k=1}^n{\tilde{y}}_{ik}{\tilde{e}}_{kj}=\sum _{k=1}^n(U_{ij},W_{ij},Y_{ij},Z_{ij})\mathrm{\forall }1<i<m,1<j<p.$ which can be written as $\tilde{Y}\tilde{E}=(U_{ij},W_{ij},Y_{ij},Z_{ij})$ where, $\begin{array}{rl}{U}_{ij}& =min(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)})\end{array}$ $\begin{array}{rl}{W}_{ij}& =min(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)})\end{array}$ $\begin{array}{rl}{Y}_{ij}& =max(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)})\end{array}$ $\begin{array}{rl}{Z}_{ij}& =max(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)})\end{array}$

Combining $\tilde{A}\tilde{X}$ and $\tilde{Y}\tilde{B}$, $\tilde{D}\tilde{X}$ and $\tilde{Y}\tilde{E}$ we get: (4.2) $\{\begin{array}{c}\sum _{k=1}^n{\tilde{a}}_{ik}{\tilde{x}}_{kj}+\sum _{k=1}^n{\tilde{y}}_{ik}{\tilde{b}}_{kj}={\tilde{c}}_{ij}\mathrm{\forall }1\le i\le m,1\le j\le p,\\ \sum _{k=1}^n{\tilde{d}}_{ik}{\tilde{x}}_{kj}+\sum _{k=1}^n{\tilde{y}}_{ik}{\tilde{e}}_{kj}={\tilde{f}}_{ij}\mathrm{\forall }1\le i\le m,1\le j\le p.\end{array}$(4.2)

By Definition 2.7 and Definition 1.1, the arbitrary CTrFFSME in Eq. (1.1) is equivalent to the following system: $\{\begin{array}{c}min(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})+min(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)})=e_{ij}^{(1)},\\ min(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})+min(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)})=e_{ij}^{(2)},\\ max(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})+max(y_{ij}^{(2)}{b}_{ij}^{(2)},y_{ij}^{(2)}{b}_{ij}^{(3)},y_{ij}^{(3)}{b}_{ij}^{(2)},y_{ij}^{(3)}{b}_{ij}^{(3)})=e_{ij}^{(3)},\\ max(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})+max(y_{ij}^{(1)}{b}_{ij}^{(1)},y_{ij}^{(1)}{b}_{ij}^{(4)},y_{ij}^{(4)}{b}_{ij}^{(1)},y_{ij}^{(4)}{b}_{ij}^{(4)})=e_{ij}^{(4)},\\ min(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)})+min(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)})=f_{ij}^{(1)},\\ min(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)})+min(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)})=f_{ij}^{(2)},\\ max(d_{ij}^{(2)}{x}_{ij}^{(2)},d_{ij}^{(2)}{x}_{ij}^{(3)},d_{ij}^{(3)}{x}_{ij}^{(2)},d_{ij}^{(3)}{x}_{ij}^{(3)})+max(y_{ij}^{(2)}{e}_{ij}^{(2)},y_{ij}^{(2)}{e}_{ij}^{(3)},y_{ij}^{(3)}{e}_{ij}^{(2)},y_{ij}^{(3)}{e}_{ij}^{(3)})=f_{ij}^{(3)},\\ max(d_{ij}^{(1)}{x}_{ij}^{(1)},d_{ij}^{(1)}{x}_{ij}^{(4)},d_{ij}^{(4)}{x}_{ij}^{(1)},d_{ij}^{(4)}{x}_{ij}^{(4)})+max(y_{ij}^{(1)}{e}_{ij}^{(1)},y_{ij}^{(1)}{e}_{ij}^{(4)},y_{ij}^{(4)}{e}_{ij}^{(1)},y_{ij}^{(4)}{e}_{ij}^{(4)})=f_{ij}^{(4)}.\end{array}$

To solve the ACTrFFSME in Equation (1.1), we consider the corresponding system of non-linear equations in Equation (4.2).

Remark 4.1:

The number of non-linear equations obtained from $m\times n$ ACTrFFSME is equal to $8\times m\times n$ equations. For example, $2\times 2$ CTrFFSME can be converted to a system of 32 non-linear equations.

Authors in [41] showed that the reduced non-linear system of equations can be converted to an equivalent system of absolute equations. Therefore, the obtained non-linear system in Equation (3.2) can be converted to an equivalent system of absolute equations where the nature of the fuzzy solutions of the ACTrFFSME depends upon the nature of the solutions of the absolute system which may be no solution, unique, many definite or infinitely many solutions.

The following algorithm can be used to obtain the fuzzy solution to the ACTrFFSME in Equation (1.1) as follows:

Algorithm 1 Solving arbitrary CTrFFSME

Step 1: Convert the CTrFFSME in Equation (1.1) to the reduced non-linear system in Equation (4.2) using Theorem 4.1.

Step 2: Convert the reduced non-linear system to an absolute system using Equation (2.6) in Definition 2.10.

Step 3: Solve the system of absolute equations.

Step 4: By solving the system of absolute equations and eliminating the non-fuzzy solutions, the following arbitrary fuzzy solution is obtained: (4.3) $\{\begin{array}{c}\tilde{X}=\left(\begin{array}{ccc}(x_{11}^{(1)},x_{11}^{(2)},x_{11}^{(3)},x_{11}^{(4)})& \cdots & (x_{1p}^{(1)},x_{1p}^{(2)},x_{1p}^{(3)},x_{1p}^{(4)})\\ ⋮& \ddots & ⋮\\ (x_{n1}^{(1)},x_{n1}^{(2)},x_{n1}^{(3)},x_{n1}^{(4)})& \cdots & (x_{np}^{(1)},x_{np}^{(2)},x_{np}^{(3)},x_{np}^{(4)})\end{array}\right)\\ \tilde{Y}=\left(\begin{array}{ccc}(y_{11}^{(1)},y_{11}^{(2)},y_{11}^{(3)},y_{11}^{(4)})& \cdots & (y_{1p}^{(1)},y_{1p}^{(2)},y_{1p}^{(3)},y_{1p}^{(4)})\\ ⋮& \ddots & ⋮\\ (y_{m1}^{(1)},y_{m1}^{(2)},y_{m1}^{(3)},y_{m1}^{(4)})& \cdots & (y_{mr}^{(1)},y_{mr}^{(2)},y_{mr}^{(3)},y_{mr}^{(4)})\end{array}\right)\end{array}$(4.3)

Feasibility of the arbitrary fuzzy solution to the CTrFFSME

The obtained arbitrary fuzzy solution in Equation (4.3) to the CTrFFSME in Equation (1.1) is feasible (strong arbitrary fuzzy solution) if the following conditions are satisfied:

1. $x_{ij}^{(1)}\le x_{ij}^{(2)}\le x_{ij}^{(3)}\le x_{ij}^{(4)}\mathrm{\forall }1\le i\le n,1\le j\le p.$

2. $y_{ij}^{(1)}\le y_{ij}^{(2)}\le y_{ij}^{(3)}\le y_{ij}^{(4)}\mathrm{\forall }1\le i\le m,1\le j\le r.$

3. At least one element of $\tilde{X}$ is near-zero TrFN.

4. At least one element of $\tilde{Y}$ is near-zero TrFN.

Remark 4.2:

If the solution fails to satisfy the feasibility conditions, then it is an infeasible fuzzy solution (weak and also known as a non-fuzzy solution).

Algorithm 1 can be applied to ACTrFFSME with different sizes. In the following Section 4.2, the proposed method is explicitly applied to a $2\times 2$ arbitrary CTrFFSME.

4.2. Applications of the Proposed Method to Other Fuzzy Systems and Fuzzy Numbers.

The proposed method is able to solve different arbitrary fuzzy systems with triangular and trapezoidal fuzzy numbers, without any restrictions. It can be directly applied to the TrFFME in Equation (1.4). In the following Theorem 4.2 the TrFFME in Equation (1.4) is converted to an equivalent system of non-linear matrix equations.

Theorem 4.2:

If $\tilde{A}=({\tilde{a}}_{ij}{)}_{m\times n},\mathrm{\forall }1\le i,j\le m,n\text{ and }\tilde{X}=({\tilde{x}}_{ij}{)}_{n\times r}$, $\mathrm{\forall }1\le i,j\le n,r$. Then the arbitrary TrFFME in Equation (1.4) can be written as follows: (4.5) $\{\begin{array}{c}\sum _{i,j=1}^{n}min(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})=e_{ij}^{(1)}\\ \sum _{i,j=1}^{n}min(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})=e_{ij}^{(2)}\\ \sum _{i,j=1}^{n}max(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})=e_{ij}^{(3)}\\ \sum _{i,j=1}^{n}max(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})=e_{ij}^{(4)}\end{array}$(4.5)

Proof:

Assuming $\tilde{A}=({\tilde{a}}_{ij}{)}_{m\times n}=(a_{ij}^{(1)},a_{ij}^{(2)},a_{ij}^{(3)},a_{ij}^{(4)})$ and $\tilde{X}=({\tilde{x}}_{ij}{)}_{n\times r}=(x_{ij}^{(1)},x_{ij}^{(2)},x_{ij}^{(3)},x_{ij}^{(4)})$ and $\tilde{E}=({\tilde{e}}_{ij}{)}_{m\times r}=(e_{ij}^{(1)},e_{ij}^{(2)},e_{ij}^{(3)},e_{ij}^{(4)})$ arbitrary TrFNs. By Equation (2.2) in Definition 2.5, we have, $\tilde{A}\tilde{X}=(M,N,P,Q)$ where, $\begin{array}{rl}M& =min(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)}),\end{array}$ $\begin{array}{rl}N& =min(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}P& =max(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)}),\end{array}$ $\begin{array}{rl}Q& =max(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)}).\end{array}$

By the definition of arbitrary TrFFME in Eq. (1.4), it can be written as: $\{\begin{array}{c}\sum _{i,j=1}^{n}min(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})=e_{ij}^{(1)}\\ \sum _{i,j=1}^{n}min(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})=e_{ij}^{(2)}\\ \sum _{i,j=1}^{n}max(a_{ij}^{(2)}{x}_{ij}^{(2)},a_{ij}^{(2)}{x}_{ij}^{(3)},a_{ij}^{(3)}{x}_{ij}^{(2)},a_{ij}^{(3)}{x}_{ij}^{(3)})=e_{ij}^{(3)}\\ \sum _{i,j=1}^{n}max(a_{ij}^{(1)}{x}_{ij}^{(1)},a_{ij}^{(1)}{x}_{ij}^{(4)},a_{ij}^{(4)}{x}_{ij}^{(1)},a_{ij}^{(4)}{x}_{ij}^{(4)})=e_{ij}^{(4)}\end{array}$ and the proof is completed.

Based on algorithm 1, the following algorithm is used to obtain the arbitrary fuzzy solution to the TrFFME in Equation (1.4) as follows:

Algorithm 2 Solving arbitrary TrFFME.

Step 1: Convert the arbitrary TrFFME in Equation (1.4) to a reduced non-linear system in Equation (4.5) using Equation (2.2).

Step 2: Convert the reduced non-linear system to an absolute system using Definition 2.6.

Step 3: Solve the system of absolute equations and check which solution(s) satisfy the following.

1. $x_{ij}^{(1)}\le x_{ij}^{(2)}\le x_{ij}^{(3)}\le x_{ij}^{(4)}\mathrm{\forall }1\le i\le n,1\le j\le m.$

2. At least one element of $\tilde{X}$ is near-zero TrFN.

Step 4: By solving the system of absolute equations and eliminating the non-fuzzy solutions, the following arbitrary fuzzy solution(s) is obtained: $\tilde{X}=\left(\begin{array}{ccc}(x_{11}^{(1)},x_{11}^{(2)},x_{11}^{(3)},x_{11}^{(4)})& \cdots & (x_{1p}^{(1)},x_{1p}^{(2)},x_{1p}^{(3)},x_{1p}^{(4)})\\ ⋮& \ddots & ⋮\\ (x_{n1}^{(1)},x_{n1}^{(2)},x_{n1}^{(3)},x_{n1}^{(4)})& \dots & (x_{np}^{(1)},x_{np}^{(2)},x_{np}^{(3)},x_{np}^{(4)})\end{array}\right).$

Similarly, the proposed method is able to solve different arbitrary fuzzy systems with triangular and trapezoidal fuzzy numbers without any restrictions. It can be directly applied to the following:

Fully fuzzy Sylvester matrix equation $\tilde{A}\tilde{X}+\tilde{X}\tilde{D}=\tilde{E}$.

Fully fuzzy continuous-time Lyapunov matrix equation $\tilde{A}\tilde{X}+\tilde{X}{\tilde{A}}^T=\tilde{E}$.

To illustrate the proposed methods, two examples are solved in the following Section 5. The first example is for ACTrFFSME, and the second example is for FFME with TFNs.

5. Numerical Example

In this section, the proposed method in Section 3 is illustrated by solving two examples.

Example 5.1:

Consider the following arbitrary CTrFFSME and solve it by the proposed method: $\{\begin{array}{c}\tilde{A}\tilde{X}+\tilde{Y}\tilde{B}=\tilde{C}\\ \tilde{D}\tilde{X}+\tilde{Y}\tilde{E}=\tilde{F}\end{array}$ where, $\begin{array}{rl}\tilde{A}& =\left(\begin{array}{cc}(5,7,9,11)& (2,4,7,10)\\ (-11,-9,-5,-2)& (-4,-3,2,4)\end{array}\right),\tilde{B}=\left(\begin{array}{cc}(2,4,7,9)& (2,3,5,7)\\ (-5,-4,-3,-2)& (-4,-2,4,6)\end{array}\right),\\ \tilde{E}& =\left(\begin{array}{cc}(1,4,5,7)& (3,4,5,7)\\ (-5,-4,-3,-2)& (-6,-3,2,4)\end{array}\right),\tilde{D}=\left(\begin{array}{cc}(5,7,10,12)& (2,4,5,6)\\ (-10,-7,-4,-2)& (-3,-2,3,4)\end{array}\right),\\ \tilde{F}& =\left(\begin{array}{cc}\left(-106,-47,43,102\right)& \left(-79,-22,80,142\right)\\ \left(-108,-50,23,86\right)& \left(-138,-68,-10,50\right)\end{array}\right),\\ \tilde{C}& =\left(\begin{array}{cc}\left(-119,-51,50,126\right)& \left(-87,-25,89,160\right)\\ \left(-125,-65,24,95\right)& \left(-134,-76,-6,65\right)\end{array}\right).\end{array}$

Solution

The solution to the given CTrFFSME is obtained by Algorithm 1 as follows:

Step 1: Converting the given arbitrary $2\times 2$ CTrFFSME to a reduced system of non-linear equations using RAMO in Section 3 as follows: $\begin{array}{rl}& min(5x_{11}^{(1)},11x_{11}^{(1)})+min(2x_{21}^{(1)},10x_{21}^{(1)})+min(2y_{11}^{(1)},9y_{11}^{(1)})+min(-5y_{12}^{(4)},-2y_{12}^{(4)})=-119,\end{array}$ $\begin{array}{rl}& min(7x_{11}^{(2)},9x_{11}^{(2)})+min(4x_{21}^{(2)},7x_{21}^{(2)})+min(7y_{11}^{(2)},4y_{11}^{(2)})+min(-4y_{12}^{(3)},-3y_{12}^{(3)})=-51,\end{array}$ $\begin{array}{rl}& max(7x_{11}^{(3)},9x_{11}^{(3)})+max(4x_{21}^{(3)},7x_{21}^{(3)})+max(4y_{11}^{(3)},7y_{11}^{(3)})+max(-4y_{12}^{(2)},-3y_{12}^{(2)})=50,\end{array}$ $\begin{array}{rl}& max(5x_{11}^{(4)},11x_{11}^{(4)})+max(2x_{21}^{(4)},10x_{21}^{(4)})+max(2y_{11}^{(4)},9y_{11}^{(4)})+max(-5y_{12}^{(1)},-2y_{12}^{(1)})=126,\end{array}$ $\begin{array}{rl}& min(5x_{12}^{(1)},11x_{12}^{(1)})+min(2x_{22}^{(1)},10x_{22}^{(1)})+min(2y_{11}^{(1)},7y_{11}^{(1)})+min(6y_{12}^{(1)},-4y_{12}^{(4)})=-87,\end{array}$ $\begin{array}{rl}& min(7x_{12}^{(2)},9x_{12}^{(2)})+min(4x_{22}^{(2)},7x_{22}^{(2)})+min(5y_{11}^{(2)},3y_{11}^{(2)})+min(4y_{12}^{(2)},-2y_{12}^{(3)})=-25,\end{array}$ $\begin{array}{rl}& max(7x_{12}^{(3)},9x_{12}^{(3)})+max(4x_{22}^{(3)},7x_{22}^{(3)})+max(3y_{11}^{(3)},5y_{11}^{(3)})+max(-2y_{12}^{(2)},4y_{12}^{(3)})=89,\end{array}$ $\begin{array}{rl}& max(5x_{12}^{(4)},11x_{12}^{(4)})+max(2x_{22}^{(4)},10x_{22}^{(4)})+max(2y_{11}^{(4)},7y_{11}^{(4)})+max(-4y_{12}^{(1)},6y_{12}^{(4)})=160,\end{array}$ $\begin{array}{rl}& min(-2x_{11}^{(1)},-11x_{11}^{(4)})+min(4x_{21}^{(1)},-4x_{21}^{(4)})+min(2y_{21}^{(1)},9y_{21}^{(1)})\\ & +min(-5y_{22}^{(4)},-2y_{22}^{(4)})=-125\end{array}$ $\begin{array}{rl}& min(-5x_{11}^{(3)},-9x_{11}^{(3)})+min(2x_{21}^{(2)},-3x_{21}^{(3)})+min(4y_{21}^{(2)},7y_{21}^{(2)})\\ & +min(-4y_{22}^{(3)},-3y_{22}^{(3)})=-65,\end{array}$ $\begin{array}{rl}& max(-9x_{11}^{(2)},-5x_{11}^{(2)})+max(-3x_{21}^{(2)},2x_{21}^{(3)})+max(4y_{21}^{(3)},7y_{21}^{(3)})\\ & +max(-4y_{22}^{(2)},-3y_{22}^{(2)})=24,\end{array}$ $\begin{array}{rl}& max(-11x_{11}^{(1)},-2x_{11}^{(1)})+max(-4x_{21}^{(1)},4x_{21}^{(4)})+max(2y_{21}^{(4)},9y_{21}^{(4)})\\ & +max(-5y_{22}^{(1)},-2y_{22}^{(1)})=95,\end{array}$ $\begin{array}{rl}& min(-2x_{12}^{(1)},-11x_{12}^{(4)})+min(4x_{22}^{(1)},-4x_{22}^{(4)})+min(2y_{21}^{(1)},7y_{21}^{(1)})\\ & +min(6y_{22}^{(1)},-4y_{22}^{(4)})=-134,\end{array}$ $\begin{array}{rl}& min(-9x_{12}^{(3)},-5x_{12}^{(3)})+min(2x_{22}^{(2)},-3x_{22}^{(3)})+min(3y_{21}^{(2)},5y_{21}^{(2)})+min(4y_{22}^{(2)},-2y_{22}^{(3)})=-76,\end{array}$ $\begin{array}{rl}& max(-9x_{12}^{(2)},-5x_{12}^{(2)})+max(-3x_{22}^{(2)},2x_{22}^{(3)})+max(3y_{21}^{(3)},5y_{21}^{(3)})\\ & +max(-2y_{22}^{(2)},4y_{22}^{(3)})=-6,\end{array}$ $\begin{array}{rl}& max(-11x_{12}^{(1)},-2x_{12}^{(1)})+max(-4x_{22}^{(1)},4x_{22}^{(4)})+max(2y_{21}^{(4)},7y_{21}^{(4)})\\ & +max(-4y_{22}^{(1)},6y_{22}^{(4)})=65,\end{array}$ $\begin{array}{rl}& min(5x_{11}^{(1)},12x_{11}^{(1)})+min(2x_{21}^{(1)},6x_{21}^{(1)})+min(y_{11}^{(1)},7y_{11}^{(1)})+min(-5y_{12}^{(4)},-2y_{12}^{(4)})=-106,\end{array}$ $\begin{array}{rl}& min(7x_{11}^{(2)},10x_{11}^{(2)})+min(4x_{21}^{(2)},5x_{21}^{(2)})+min(5y_{11}^{(2)},4y_{11}^{(2)})+min(-4y_{12}^{(3)},-3y_{12}^{(3)})=-47,\end{array}$ $\begin{array}{rl}& max(7x_{11}^{(3)},10x_{11}^{(3)})+max(4x_{21}^{(3)},5x_{21}^{(3)})+max(4y_{11}^{(3)},5y_{11}^{(3)})+max(-4y_{12}^{(2)},-3y_{12}^{(2)})=43,\end{array}$ $\begin{array}{rl}& max(5x_{11}^{(4)},12x_{11}^{(4)})+max(2x_{21}^{(4)},6x_{21}^{(4)})+max(y_{11}^{(4)},7y_{11}^{(4)})+max(-5y_{12}^{(1)},-2y_{12}^{(1)})=102,\end{array}$ $\begin{array}{rl}& min(5x_{12}^{(1)},12x_{12}^{(1)})+min(2x_{22}^{(1)},6x_{22}^{(1)})+min(3y_{11}^{(1)},7y_{11}^{(1)})+min(4y_{12}^{(1)},-6y_{12}^{(4)})=-79,\end{array}$ $\begin{array}{rl}& min(7x_{12}^{(2)},10x_{12}^{(2)})+min(4x_{22}^{(2)},5x_{22}^{(2)})+min(5y_{11}^{(2)},4y_{11}^{(2)})+min(2y_{12}^{(2)},-3y_{12}^{(3)})=-22\end{array}$ $\begin{array}{rl}& max(7x_{12}^{(3)},10x_{12}^{(3)})+max(4x_{22}^{(3)},5x_{22}^{(3)})+max(4y_{11}^{(3)},5y_{11}^{(3)})+max(-3y_{12}^{(2)},2y_{12}^{(3)})=80,\end{array}$ $\begin{array}{rl}& max(5x_{12}^{(4)},12x_{12}^{(4)})+max(2x_{22}^{(4)},6x_{22}^{(4)})+max(3y_{11}^{(4)},7y_{11}^{(4)})+max(-6y_{12}^{(1)},4y_{12}^{(4)})=142,\end{array}$ $\begin{array}{rl}& min(-2x_{11}^{(1)},-10x_{11}^{(4)})+min(4x_{21}^{(1)},-3x_{21}^{(4)})+min(y_{21}^{(1)},7y_{21}^{(1)})\\ & +min(-5y_{22}^{(4)},-2y_{22}^{(4)})=-108,\end{array}$ $\begin{array}{rl}& min(-4x_{11}^{(3)},-7x_{11}^{(3)})+min(3x_{21}^{(2)},-2x_{21}^{(3)})+min(4y_{21}^{(2)},5y_{21}^{(2)})\\ & +min(-4y_{22}^{(3)},-3y_{22}^{(3)})=-50,\end{array}$ $\begin{array}{rl}& max(-7x_{11}^{(2)},-4x_{11}^{(2)})+max(-2x_{21}^{(2)},3x_{21}^{(3)})+max(4y_{21}^{(3)},5y_{21}^{(3)})\\ & +max(-4y_{22}^{(2)},-3y_{22}^{(2)})=23,\end{array}$ $\begin{array}{rl}& max(-10x_{11}^{(1)},-2x_{11}^{(1)})+max(-3x_{21}^{(1)},4x_{21}^{(4)})+max(y_{21}^{(4)},7y_{21}^{(4)})\\ & +max(-5y_{22}^{(1)},-2y_{22}^{(1)})=86,\end{array}$ $\begin{array}{rl}& min(-2x_{12}^{(1)},-10x_{12}^{(4)})+min(4x_{22}^{(1)},-3x_{22}^{(4)})+min(3y_{21}^{(1)},7y_{21}^{(1)})\\ & +min(4y_{22}^{(1)},-6y_{22}^{(4)})=-138,\end{array}$ $\begin{array}{rl}& min(-7x_{12}^{(3)},-4x_{12}^{(3)})+min(3x_{22}^{(2)},-2x_{22}^{(3)})+min(4y_{21}^{(2)},5y_{21}^{(2)})+min(2y_{22}^{(2)},-3y_{22}^{(3)})=-68,\end{array}$ $\begin{array}{rl}& max(-7x_{12}^{(2)},-4x_{12}^{(2)})+max(-2x_{22}^{(2)},3x_{22}^{(3)})+max(4y_{21}^{(3)},5y_{21}^{(3)})\\ & +max(-3y_{22}^{(2)},2y_{22}^{(3)})=-10,\end{array}$ $\begin{array}{rl}& max(-10x_{12}^{(1)},-2x_{12}^{(1)})+max(-3x_{22}^{(1)},4x_{22}^{(4)})+max(3y_{21}^{(4)},7y_{21}^{(4)})\\ & +max(-6y_{22}^{(1)},4y_{22}^{(4)})=50,\end{array}$

Subject to:

$x_{11}^{(1)}-x_{11}^{(2)}\le 0$,

$x_{11}^{(2)}-x_{11}^{(3)}\le 0$,

$x_{11}^{(3)}-x_{11}^{(4)}\le 0$,

$x_{12}^{(1)}-x_{12}^{(2)}\le 0$,

$x_{12}^{(2)}-x_{12}^{(3)}\le 0$,

$x_{12}^{(3)}-x_{12}^{(4)}\le 0$,

$x_{21}^{(1)}-x_{21}^{(2)}\le 0$,

$x_{21}^{(2)}-x_{21}^{(3)}\le 0$,

$x_{21}^{(3)}-x_{21}^{(4)}\le 0$,

$x_{22}^{(1)}-x_{22}^{(2)}\le 0$,

$x_{22}^{(2)}-x_{22}^{(3)}\le 0$,

$x_{22}^{(3)}-x_{22}^{(4)}\le 0$,

$y_{11}^{(1)}-y_{11}^{(2)}\le 0$,

$y_{11}^{(2)}-y_{11}^{(3)}\le 0$,

$y_{11}^{(3)}-y_{11}^{(4)}\le 0$,

$y_{12}^{(1)}-y_{12}^{(2)}\le 0$,

$y_{12}^{(2)}-y_{12}^{(3)}\le 0$,

$y_{12}^{(3)}-y_{12}^{(4)}\le 0$,

$y_{21}^{(1)}-y_{21}^{(2)}\le 0$,

$y_{21}^{(2)}-y_{21}^{(3)}\le 0$,

$y_{21}^{(3)}-y_{21}^{(4)}\le 0$,

$y_{22}^{(1)}-y_{22}^{(2)}\le 0$,

$y_{22}^{(2)}-y_{22}^{(3)}\le 0$,

$y_{22}^{(3)}-y_{22}^{(4)}\le 0$.

Step 2: The reduced non-linear system in step 1 is reduced to an absolute system using Equation (2.6) in Definition 2.10; see Appendix A for the obtained absolute system of equations.

Step 3: Solve the system of absolute equations in step 2 and choose the solution(s) that satisfy the following.

1. $x_{ij}^{(1)}\le x_{ij}^{(2)}\le x_{ij}^{(3)}\le x_{ij}^{(4)}\mathrm{\forall }1\le i\le 2,1\le j\le 2.$

2. $y_{ij}^{(1)}\le y_{ij}^{(2)}\le y_{ij}^{(3)}\le y_{ij}^{(4)}\mathrm{\forall }1\le i\le 2,1\le j\le 2.$

3. At least one element of $\tilde{X}$ is near zero TrFN.

4. At least one element of $\tilde{Y}$ is near zero TrFN.

Step 4: By solving the system of absolute equations and eliminating the non-fuzzy solutions, the following arbitrary fuzzy solution is obtained: $\{\begin{array}{c}\tilde{X}=\left(\begin{array}{cc}(-3,-2,3,4)& (3,4,5,6)\\ (-2,2,3,5)& (-5,-4,2,3)\end{array}\right),\\ \tilde{Y}=\left(\begin{array}{cc}(-4,-3,2,4)& (2,4,5,6)\\ (-4,-3,-2,3)& (-3,-2,2,5)\end{array}\right).\end{array}$

Figure 3 shows the arbitrary trapezoidal fuzzy solution $\tilde{X}$ for example 4.1.

Figure 3. The arbitrary fuzzy solution $\tilde{X}$ of Example 5.1

Verification of the solution:

To verify the obtained arbitrary fuzzy solution, we first multiply $\tilde{A}\tilde{X}$, $\tilde{Y}\tilde{B}$, $\tilde{D}\tilde{X}$ and $\tilde{Y}\tilde{E}$ as follows: $\begin{array}{rl}\tilde{A}\tilde{X}& =\left(\begin{array}{cc}(5,7,9,11)& (2,4,7,10)\\ (-11,-9,-5,-2)& (-4,-3,2,4)\end{array}\right)\left(\begin{array}{cc}(-3,-2,3,4)& (3,4,5,6)\\ (-2,2,3,5)& (-5,-4,2,3)\end{array}\right)\\ & =\left(\begin{array}{cc}(-53,-10,48,94)& (-35,0,59,96)\\ (-64,-36,24,53)& (-86,-53,-8,14)\end{array}\right)\end{array}$ $\begin{array}{rl}\tilde{Y}\tilde{B}& =\left(\begin{array}{cc}(-4,-3,2,4)& (2,4,5,6)\\ (-4,-3,-2,3)& (-3,-2,2,5)\end{array}\right)\left(\begin{array}{cc}(2,4,7,9)& (2,3,5,7)\\ (-5,-4,-3,-2)& (-4,-2,4,6)\end{array}\right)\\ & =\left(\begin{array}{cc}\left(\begin{array}{c}-66,-41,2,32\end{array}\right)& \left(\begin{array}{c}-52,-25,30,64\end{array}\right)\\ \left(\begin{array}{c}-61,-29,0,42\end{array}\right)& \left(\begin{array}{c}-48,-23,2,51\end{array}\right)\end{array}\right)\end{array}$ $\begin{array}{rl}\tilde{D}\tilde{X}& =\left(\begin{array}{cc}(5,7,10,12)& (2,4,5,6)\\ (-10,-7,-4,-2)& (-3,-2,3,4)\end{array}\right)\left(\begin{array}{cc}(-3,-2,3,4)& (3,4,5,6)\\ (-2,2,3,5)& (-5,-4,2,3)\end{array}\right)\\ & =\left(\begin{array}{cc}\left(\begin{array}{c}-48,-12,45,78\end{array}\right)& \left(\begin{array}{c}-15,8,60,90\end{array}\right)\\ \left(\begin{array}{c}-55,-27,23,50\end{array}\right)& \left(\begin{array}{c}-80,-47,-8,9\end{array}\right)\end{array}\right).\end{array}$

Adding, $\tilde{A}\tilde{X}+\tilde{Y}\tilde{B}$ and $\tilde{D}\tilde{X}+\tilde{Y}\tilde{E}$ gives: $\begin{array}{rl}\tilde{A}\tilde{X}+\tilde{Y}\tilde{B}& =\left(\begin{array}{cc}\left(\begin{array}{c}-53,-10,48,94\end{array}\right)& \left(\begin{array}{c}-35,0,59,96\end{array}\right)\\ \left(\begin{array}{c}-64,-36,24,53\end{array}\right)& \left(\begin{array}{c}-86,-53,-8,14\end{array}\right)\end{array}\right)\\ & +\left(\begin{array}{cc}\left(\begin{array}{c}-66,-41,2,32\end{array}\right)& \left(\begin{array}{c}-52,-25,30,64\end{array}\right)\\ \left(\begin{array}{c}-61,-29,0,42\end{array}\right)& \left(\begin{array}{c}-48,-23,2,51\end{array}\right)\end{array}\right)\\ & =\left(\begin{array}{cc}\left(\begin{array}{c}-119,-51,50,126\end{array}\right)& \left(\begin{array}{c}-87,-25,89,160\end{array}\right)\\ \left(\begin{array}{c}-125,-65,24,95\end{array}\right)& \left(\begin{array}{c}-134,-76,-6,65\end{array}\right)\end{array}\right)=\tilde{C}\end{array}$ $\begin{array}{rl}\tilde{D}\tilde{X}+\tilde{Y}\tilde{E}& =\left(\begin{array}{cc}\left(\begin{array}{c}-48,-12,45,78\end{array}\right)& \left(\begin{array}{c}-15,8,60,90\end{array}\right)\\ \left(\begin{array}{c}-55,-27,23,50\end{array}\right)& \left(\begin{array}{c}-80,-47,-8,9\end{array}\right)\end{array}\right)\\ & +\left(\begin{array}{cc}\left(\begin{array}{c}-58,-35,-2,24\end{array}\right)& \left(\begin{array}{c}-64,-30,20,52\end{array}\right)\\ \left(\begin{array}{c}-53,-23,0,36\end{array}\right)& \left(\begin{array}{c}-58,-21,-2,41\end{array}\right)\end{array}\right)\\ & =\left(\begin{array}{cc}\left(\begin{array}{c}-106,-47,43,102\end{array}\right)& \left(\begin{array}{c}-79,-22,80,142\end{array}\right)\\ \left(\begin{array}{c}-108,-50,23,86\end{array}\right)& \left(\begin{array}{c}-138,-68,-10,50\end{array}\right)\end{array}\right)=\tilde{F}\end{array}$

The obtained arbitrary fuzzy solution satisfies the given PCTrFFSME and it is clearly feasible.

The following Example 5.2, was first solved by Kumar et al. [39] and obtained only one fuzzy solution. However, Malkawi et al. [41] considered the same example and obtained two fuzzy solutions. To support the developed method, the same example is considered.

Remark 5.1:

The methods used in solving the example considered by Kumar et al. [39] and Malkawi et al. [41] can only be applied to fuzzy systems with TFNs. Therefore, in the following example, we extend the TFNs to TrFNs by equating the mean values in the TrFNs used, to apply the developed method in this section.

Example 5.2:

Consider the following arbitrary TrFFME: $\left(\begin{array}{cc}(-2,3,3,4)& (-2,2,2,3)\\ (1,2,2,2)& (4,4,4,5)\end{array}\right)\left(\begin{array}{c}(x_{11}^{(1)},x_{11}^{(2)},x_{11}^{(3)},x_{11}^{(4)})\\ (x_{21}^{(1)},x_{21}^{(2)},x_{21}^{(3)},x_{21}^{(4)})\end{array}\right)=\left(\begin{array}{c}(-13,8,8,14)\\ (-14,8,8,14)\end{array}\right)$

Solution: The fuzzy solution to the given TrFFME obtained by Kumar et al. [39] is as follows: $\tilde{X}=\left(\begin{array}{c}(1,2,2,2)\\ (-3,1,1,2)\end{array}\right).$

However, Malkawi et al. [41] were able to obtain two fuzzy solutions for the given TrFFME as follows: $\tilde{X}=\left\{\left(\begin{array}{c}(1,2,2,2)\\ (-3,1,1,2)\end{array}\right),\left(\begin{array}{c}(-23/14,2,2,2)\\ (-15/7,1,1,2)\end{array}\right)\right\}$

The fuzzy solutions to the given TrFFME can be obtained by Algorithm 2 as follows:

Step 1: Convert the arbitrary TrFFME in Equation (1.4) to a reduced non-linear system in Equation (3.5) as follows: $\begin{array}{rl}min(-2x_{11}^{(1)},4x_{11}^{(1)})+min(-2x_{21}^{(1)},3x_{21}^{(1)})& =-13,\end{array}$ $\begin{array}{rl}min(3x_{11}^{(2)},3x_{11}^{(3)})+min(2x_{21}^{(2)},2x_{21}^{(3)})& =8,\end{array}$ $\begin{array}{rl}max(3x_{11}^{(2)},3x_{11}^{(3)})+max(2x_{21}^{(2)},2x_{21}^{(3)})& =8,\end{array}$ $\begin{array}{rl}max(-2x_{11}^{(1)},4x_{11}^{(4)})+max(-2x_{21}^{(1)},3x_{21}^{(4)})& =14,\end{array}$ $\begin{array}{rl}min(x_{11}^{(1)},2x_{11}^{(1)})+min(4x_{21}^{(1)},5x_{21}^{(1)})& =-14,\end{array}$ $\begin{array}{rl}min(2x_{11}^{(2)},2x_{11}^{(3)})+min(4x_{21}^{(2)},4x_{21}^{(3)})& =8,\end{array}$ $\begin{array}{rl}max(2x_{11}^{(2)},2x_{11}^{(3)})+max(4x_{21}^{(2)},4x_{21}^{(3)})& =8,\end{array}$ $\begin{array}{rl}max(x_{11}^{(4)},2x_{11}^{(4)})+max(4x_{21}^{(4)},5x_{21}^{(4)})& =14.\end{array}$

Subject to:

$x_{11}^{(1)}-x_{11}^{(2)}\le 0$,

$x_{11}^{(2)}-x_{11}^{(3)}\le 0$,

$x_{11}^{(3)}-x_{11}^{(4)}\le 0$,

$x_{21}^{(1)}-x_{21}^{(2)}\le 0$,

$x_{21}^{(2)}-x_{21}^{(3)}\le 0$,

$x_{21}^{(3)}-x_{21}^{(4)}\le 0$.

Step 2: Convert the reduced non-linear system to an absolute system using Definition 2.6. $\begin{array}{rl}{x}_{11}^{(1)}+\frac{x_{21}^{(1)}}{2}-3|x_{11}^{(1)}|-\frac{5|x_{21}^{(1)}|}{2}& =-13,\end{array}$ $\begin{array}{rl}\frac{1}{2}(3x_{11}^{(2)}+3x_{11}^{(3)})+\frac{1}{2}(2x_{21}^{(2)}+2x_{21}^{(3)})-\frac{1}{2}|3x_{11}^{(2)}-3x_{11}^{(3)}|-\frac{1}{2}|2x_{21}^{(2)}-2x_{21}^{(3)}|& =8,\end{array}$ $\begin{array}{rl}\frac{1}{2}(3x_{11}^{(2)}+3x_{11}^{(3)})+\frac{1}{2}(2x_{21}^{(2)}+2x_{21}^{(3)})+\frac{1}{2}|3x_{11}^{(2)}-3x_{11}^{(3)}|+\frac{1}{2}|2x_{21}^{(2)}-2x_{21}^{(3)}|& =8,\end{array}$ $\begin{array}{rl}\frac{1}{2}(-2x_{11}^{(1)}+4x_{11}^{(4)})+\frac{1}{2}(-2x_{21}^{(1)}+3x_{21}^{(4)})+\frac{1}{2}|-2x_{11}^{(1)}-4x_{11}^{(4)}|+\frac{1}{2}|-2x_{21}^{(1)}-3x_{21}^{(4)}|& =14,\end{array}$ $\begin{array}{rl}\frac{3x_{11}^{(1)}}{2}+\frac{9x_{21}^{(1)}}{2}-\frac{|x_{11}^{(1)}|}{2}-\frac{|x_{21}^{(1)}|}{2}& =-14,\end{array}$ $\begin{array}{rl}\frac{1}{2}(2x_{11}^{(2)}+2x_{11}^{(3)})+\frac{1}{2}(4x_{21}^{(2)}+4x_{21}^{(3)})-\frac{1}{2}|2x_{11}^{(2)}-2x_{11}^{(3)}|-\frac{1}{2}|4x_{21}^{(2)}-4x_{21}^{(3)}|& =8,\end{array}$ $\begin{array}{rl}\frac{1}{2}(2x_{11}^{(2)}+2x_{11}^{(3)})+\frac{1}{2}(4x_{21}^{(2)}+4x_{21}^{(3)})+\frac{1}{2}|2x_{11}^{(2)}-2x_{11}^{(3)}|+\frac{1}{2}|4x_{21}^{(2)}-4x_{21}^{(3)}|& =8,\end{array}$ $\begin{array}{rl}\frac{3x_{11}^{(4)}}{2}+\frac{9x_{21}^{(4)}}{2}+\frac{|x_{11}^{(4)}|}{2}+\frac{|x_{21}^{(4)}|}{2}& =14.\end{array}$

Subject to:

$x_{11}^{(1)}-x_{11}^{(2)}\le 0$,

$x_{11}^{(2)}-x_{11}^{(3)}\le 0$,

$x_{11}^{(3)}-x_{11}^{(4)}\le 0$,

$x_{21}^{(1)}-x_{21}^{(2)}\le 0$,

$x_{21}^{(2)}-x_{21}^{(3)}\le 0$,

$x_{21}^{(3)}-x_{21}^{(4)}\le 0$.

Step 3: Solve the system of absolute equations and check which solution(s) satisfy the following.

1. $x_{ij}^{(1)}\le x_{ij}^{(2)}\le x_{ij}^{(3)}\le x_{ij}^{(4)}\mathrm{\forall }1\le i\le n,1\le j\le m.$

2. At least one element of $\tilde{X}$ is near-zero TrFN.

Step 4: By solving the system of absolute equations and eliminating the non-fuzzy solutions, the following arbitrary fuzzy solutions are obtained:

Table

Fuzzy Solutions	Non-Fuzzy Solutions
${\tilde{X}}_1=\left(\begin{array}{c}(1,2,2,2)\\ (-3,1,1,2)\end{array}\right),$	${\tilde{X}}_1=\left(\begin{array}{c}(0,2,2,2)\\ (2,1,1,-3)\end{array}\right),$
${\tilde{X}}_2=\left(\begin{array}{c}(-2,2,2,2)\\ (0,1,1,2)\end{array}\right),$	${\tilde{X}}_2=\left(\begin{array}{c}(-1,2,2,2)\\ (2,1,1,-3)\end{array}\right),$
${\tilde{X}}_3=\left(\begin{array}{c}(2,2,2,2)\\ (-3,1,1,2)\end{array}\right),$	${\tilde{X}}_3=\left(\begin{array}{c}(2,2,2,0)\\ (2,1,1,-3)\end{array}\right),$
${\tilde{X}}_4=\left(\begin{array}{c}(0,2,2,2)\\ (-3,1,1,2)\end{array}\right),$	${\tilde{X}}_4=\left(\begin{array}{c}(2,2,2,0)\\ (2,1,1,-3)\end{array}\right),$
${\tilde{X}}_5=\left(\begin{array}{c}(-1,2,2,2)\\ (-3,1,1,2)\end{array}\right),$	${\tilde{X}}_5=\left(\begin{array}{c}(0,2,2,2)\\ (-3,1,1,2)\end{array}\right),$
${\tilde{X}}_6=\left(\begin{array}{c}(-23/14,2,2,2)\\ (-15/7,1,1,2)\end{array}\right).$	${\tilde{X}}_6=\left(\begin{array}{c}(2,2,2,-1)\\ (-3,1,1,2)\end{array}\right).$

Verification of the obtained solution: The first two solutions are verified by Malkawi et al. [41]. Therefore, only the third solution is verified as follows: $\left(\begin{array}{cc}(-2,3,3,4)& (-2,2,2,3)\\ (1,2,2,2)& (4,4,4,5)\end{array}\right)\left(\begin{array}{c}(-23/14,2,2,2)\\ (-15/7,1,1,2)\end{array}\right)=\left(\begin{array}{c}(-13,8,8,14)\\ (-14,8,8,14)\end{array}\right).$

The obtained arbitrary fuzzy solution satisfies the given TrFFME, and it is clearly feasible.

5. Conclusion

In this paper, a new method for solving arbitrary CTrFFSME is proposed based on new arithmetic fuzzy multiplication operations. The CTrFFSME is converted to an equivalent system of non-linear equations, which is converted to a system of absolute equations, where the arbitrary fuzzy solution is obtained by solving that system. In addition, the proposed method is modified and applied to different arbitrary fuzzy equations such as the fully fuzzy matrix equation with triangular and trapezoidal fuzzy numbers. The numerical examples show that the proposed method can find all possible arbitrary fuzzy solutions. Therefore, the proposed method can solve the following fuzzy equations:

1. Couple fully fuzzy Sylvester matrix equations with arbitrary trapezoidal and triangular fuzzy numbers and arbitrary fuzzy solutions.

2. Couple fully fuzzy Lyapunov matrix equations with arbitrary trapezoidal and triangular fuzzy numbers and arbitrary fuzzy solutions.

3. Fully fuzzy Sylvester matrix equations with arbitrary trapezoidal and triangular fuzzy numbers and arbitrary fuzzy solutions.

4. Fully fuzzy Lyapunov matrix equations with arbitrary trapezoidal and triangular fuzzy numbers and arbitrary fuzzy solutions.

5. Fully fuzzy matrix equation with arbitrary triangular and trapezoidal fuzzy numbers and arbitrary fuzzy solutions.

For future work, the proposed method will be extended to a couple of arbitrary generalised trapezoidal fully fuzzy Sylvester matrix equations with arbitrary coefficients.

Human and animal rights

This article does not contain any studies with human participants or animals performed by any of the authors.

Informed consent

Informed consent was obtained from all individual participants included in the study.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Funding

This research received no external funding.

Additional information

Notes on contributors

Ahmed Abdel Aziz Elsayed

Dr. Ahmed Abdel Aziz Elsayed works at the Institute of Applied Technology, Abu Dhabi campus, United Arab Emirates (UAE). He obtained his bachelor's degree in Mathematics and Education in 2007, a Master's in Mathematics in 2018 from the American University in Sharjah, UAE, and a PhD in Mathematics in 2022 from the School of Quantitative Sciences, Universiti Utara Malaysia (UUM). His current research interests include Fuzzy Logic, Numerical Analysis and Linear Algebra.

Nazihah Ahmad

Dr. Nazihah Ahmad works at the School of Quantitative Sciences, Universiti Utara Malaysia (UUM). She obtained her bachelor's degree in Mathematics with honours in 2001, a Master's in Mathematics in 2002, and a PhD in Mathematics in 2009. Her current research interests include fuzzy mathematics, topology, and mathematical modelling in medicine, healthcare, finance, and psychology. She has been awarded several research grants by the Ministry of Education Malaysia and UUM as a principal and co-investigator. The quality of her work has been published in international journals (indexed by SCOPUS or Web of Science), conference proceedings, books and book chapter. She has been teaching mathematics for more than 20 years at the undergraduate and master's level. She has also supervised master and PhD students. Her expertise was acknowledged by her academic peers when she was appointed as an external/internal examiner, article reviewer and external assessor. At national level, she has involved in revising the mathematics syllabus based on the Secondary School Standard Curriculum and in organizing STEM service-learning activities to inspire and strengthen STEM education among young people.

Ghassan Malkawi

Ghassan Malkawi Assistant Professor. Dr. Ghassan Malkawi has been on the mathematics faculty since 2003 at the Higher College of Technology (HCT) in the United Arab Emirates. Where he holds the position of Assistant Professor in the Division of Engineering. Malkawi completed his doctorate from the Universiti Utara Malaysia. His master's degree in applied mathematics from the University of Jordan. Fuzzy linear systems are the primary area of interest for Dr. Malkawi's research.

Appendix A:

Absolute system of equations in Step 3 Example 4.1.

$\begin{array}{rl}& 8x_{11}^{(1)}+6x_{21}^{(1)}+\frac{11y_{11}^{(1)}}{2}-\frac{7y_{12}^{(4)}}{2}-3|x_{11}^{(1)}|-4|x_{21}^{(1)}|-\frac{7|y_{11}^{(1)}|}{2}-\frac{3|y_{12}^{(4)}|}{2}=-119\end{array}$ $\begin{array}{rl}& 8x_{11}^{(2)}+\frac{11x_{21}^{(2)}}{2}+\frac{11y_{11}^{(2)}}{2}-\frac{7y_{12}^{(3)}}{2}-|x_{11}^{(2)}|-\frac{3|x_{21}^{(2)}|}{2}-\frac{3|y_{11}^{(2)}|}{2}-\frac{|y_{12}^{(3)}|}{2}=-51\end{array}$ $\begin{array}{rl}& 8x_{11}^{(3)}+\frac{11x_{21}^{(3)}}{2}+\frac{11y_{11}^{(3)}}{2}-\frac{7y_{12}^{(2)}}{2}+|x_{11}^{(3)}|+\frac{3|x_{21}^{(3)}|}{2}+\frac{3|y_{11}^{(3)}|}{2}+\frac{|y_{12}^{(2)}|}{2}=50\end{array}$ $\begin{array}{rl}& 8x_{11}^{(4)}+6x_{21}^{(4)}+\frac{11y_{11}^{(4)}}{2}-\frac{7y_{12}^{(1)}}{2}+3|x_{11}^{(4)}|+4|x_{21}^{(4)}|+\frac{7|y_{11}^{(4)}|}{2}+\frac{3|y_{12}^{(1)}|}{2}=126\end{array}$ $\begin{array}{rl}& 8x_{12}^{(1)}+6x_{22}^{(1)}+\frac{9y_{11}^{(1)}}{2}+\frac{1}{2}(6y_{12}^{(1)}-4y_{12}^{(4)})-3|x_{12}^{(1)}|-4|x_{22}^{(1)}|-\frac{5|y_{11}^{(1)}|}{2}-\frac{1}{2}|6y_{12}^{(1)}+4y_{12}^{(4)}|=-87\end{array}$ $\begin{array}{rl}& 8x_{12}^{(2)}+\frac{11x_{22}^{(2)}}{2}+4y_{11}^{(2)}+\frac{1}{2}(4y_{12}^{(2)}-2y_{12}^{(3)})-|x_{12}^{(2)}|-\frac{3|x_{22}^{(2)}|}{2}-|y_{11}^{(2)}|-\frac{1}{2}|4y_{12}^{(2)}+2y_{12}^{(3)}|=-25\end{array}$ $\begin{array}{rl}& 8x_{12}^{(3)}+\frac{11x_{22}^{(3)}}{2}+4y_{11}^{(3)}+\frac{1}{2}(-2y_{12}^{(2)}+4y_{12}^{(3)})+|x_{12}^{(3)}|+\frac{3|x_{22}^{(3)}|}{2}+|y_{11}^{(3)}|+\frac{1}{2}|-2y_{12}^{(2)}-4y_{12}^{(3)}|=89\end{array}$ $\begin{array}{rl}& 8x_{12}^{(4)}+6x_{22}^{(4)}+\frac{9y_{11}^{(4)}}{2}+\frac{1}{2}(-4y_{12}^{(1)}+6y_{12}^{(4)})+3|x_{12}^{(4)}|+4|x_{22}^{(4)}|+\frac{5|y_{11}^{(4)}|}{2}+\frac{1}{2}|-4y_{12}^{(1)}-6y_{12}^{(4)}|=160\end{array}$ $\begin{array}{rl}& \frac{1}{2}(-2x_{11}^{(1)}-11x_{11}^{(4)})+\frac{1}{2}(4x_{21}^{(1)}-4x_{21}^{(4)})+\frac{11y_{21}^{(1)}}{2}-\frac{7y_{22}^{(4)}}{2}-\frac{1}{2}|-2x_{11}^{(1)}+11x_{11}^{(4)}|-\frac{1}{2}|4x_{21}^{(1)}+4x_{21}^{(4)}|\\ & -\frac{7|y_{21}^{(1)}|}{2}-\frac{3|y_{22}^{(4)}|}{2}=-125\end{array}$ $\begin{array}{rl}& -7x_{11}^{(3)}+\frac{1}{2}(2x_{21}^{(2)}-3x_{21}^{(3)})+\frac{11y_{21}^{(2)}}{2}-\frac{7y_{22}^{(3)}}{2}-2|x_{11}^{(3)}|-\frac{1}{2}|2x_{21}^{(2)}+3x_{21}^{(3)}|-\frac{3|y_{21}^{(2)}|}{2}-\frac{|y_{22}^{(3)}|}{2}=-65\end{array}$ $\begin{array}{rl}& -7x_{11}^{(2)}+\frac{1}{2}(-3x_{21}^{(2)}+2x_{21}^{(3)})+\frac{11y_{21}^{(3)}}{2}-\frac{7y_{22}^{(2)}}{2}+2|x_{11}^{(2)}|+\frac{1}{2}|-3x_{21}^{(2)}-2x_{21}^{(3)}|\\ & +\frac{3|y_{21}^{(3)}|}{2}+\frac{|y_{22}^{(2)}|}{2}=24\end{array}$ $\begin{array}{rl}& -\frac{13x_{11}^{(1)}}{2}+\frac{1}{2}(-4x_{21}^{(1)}+4x_{21}^{(4)})+\frac{11y_{21}^{(4)}}{2}-\frac{7y_{22}^{(1)}}{2}+\frac{9|x_{11}^{(1)}|}{2}+\frac{1}{2}|-4x_{21}^{(1)}-4x_{21}^{(4)}|\\ & +\frac{7|y_{21}^{(4)}|}{2}+\frac{3|y_{22}^{(1)}|}{2}=95\end{array}$ $\begin{array}{rl}& \frac{1}{2}(-2x_{12}^{(1)}-11x_{12}^{(4)})+\frac{1}{2}(4x_{22}^{(1)}-4x_{22}^{(4)})+\frac{9y_{21}^{(1)}}{2}+\frac{1}{2}(6y_{22}^{(1)}-4y_{22}^{(4)})-\frac{1}{2}|-2x_{12}^{(1)}+11x_{12}^{(4)}|\\ & -\frac{1}{2}|4x_{22}^{(1)}+4x_{22}^{(4)}|-\frac{5|y_{21}^{(1)}|}{2}-\frac{1}{2}|6y_{22}^{(1)}+4y_{22}^{(4)}|=-134\end{array}$ $\begin{array}{rl}& -7x_{12}^{(3)}+\frac{1}{2}(2x_{22}^{(2)}-3x_{22}^{(3)})+4y_{21}^{(2)}+\frac{1}{2}(4y_{22}^{(2)}-2y_{22}^{(3)})-2|x_{12}^{(3)}|-\frac{1}{2}|2x_{22}^{(2)}+3x_{22}^{(3)}|-|y_{21}^{(2)}|\\ & -\frac{1}{2}|4y_{22}^{(2)}+2y_{22}^{(3)}|=-76\end{array}$ $\begin{array}{rl}& -7x_{12}^{(2)}+\frac{1}{2}(-3x_{22}^{(2)}+2x_{22}^{(3)})+4y_{21}^{(3)}+\frac{1}{2}(-2y_{22}^{(2)}+4y_{22}^{(3)})+2|x_{12}^{(2)}|+\frac{1}{2}|-3x_{22}^{(2)}-2x_{22}^{(3)}|+|y_{21}^{(3)}|\\ & +\frac{1}{2}|-2y_{22}^{(2)}-4y_{22}^{(3)}|=-6\end{array}$ $\begin{array}{rl}& -\frac{13x_{12}^{(1)}}{2}+\frac{1}{2}(-4x_{22}^{(1)}+4x_{22}^{(4)})+\frac{9y_{21}^{(4)}}{2}+\frac{1}{2}(-4y_{22}^{(1)}+6y_{22}^{(4)})+\frac{9|x_{12}^{(1)}|}{2}+\frac{1}{2}|-4x_{22}^{(1)}-4x_{22}^{(4)}|\\ & +\frac{5|y_{21}^{(4)}|}{2}+\frac{1}{2}|-4y_{22}^{(1)}-6y_{22}^{(4)}|=65\end{array}$ $\begin{array}{rl}& \frac{17x_{11}^{(1)}}{2}+4x_{21}^{(1)}+4y_{11}^{(1)}-\frac{7y_{12}^{(4)}}{2}-\frac{7|x_{11}^{(1)}|}{2}-2|x_{21}^{(1)}|-3|y_{11}^{(1)}|-\frac{3|y_{12}^{(4)}|}{2}=-106\end{array}$ $\begin{array}{rl}& \frac{17x_{11}^{(2)}}{2}+\frac{9x_{21}^{(2)}}{2}+\frac{9y_{11}^{(2)}}{2}-\frac{7y_{12}^{(3)}}{2}-\frac{3|x_{11}^{(2)}|}{2}-\frac{|x_{21}^{(2)}|}{2}-\frac{|y_{11}^{(2)}|}{2}-\frac{|y_{12}^{(3)}|}{2}=-47\end{array}$ $\begin{array}{rl}& \frac{17x_{11}^{(3)}}{2}+\frac{9x_{21}^{(3)}}{2}+\frac{9y_{11}^{(3)}}{2}-\frac{7y_{12}^{(2)}}{2}+\frac{3|x_{11}^{(3)}|}{2}+\frac{|x_{21}^{(3)}|}{2}+\frac{|y_{11}^{(3)}|}{2}+\frac{|y_{12}^{(2)}|}{2}=43\end{array}$ $\begin{array}{rl}& \frac{17x_{11}^{(4)}}{2}+4x_{21}^{(4)}+4y_{11}^{(4)}-\frac{7y_{12}^{(1)}}{2}+\frac{7|x_{11}^{(4)}|}{2}+2|x_{21}^{(4)}|+3|y_{11}^{(4)}|+\frac{3|y_{12}^{(1)}|}{2}=102\end{array}$ $\begin{array}{rl}& \frac{17x_{12}^{(1)}}{2}+4x_{22}^{(1)}+5y_{11}^{(1)}+\frac{1}{2}(4y_{12}^{(1)}-6y_{12}^{(4)})-\frac{7|x_{12}^{(1)}|}{2}-2|x_{22}^{(1)}|-2|y_{11}^{(1)}|-\frac{1}{2}|4y_{12}^{(1)}+6y_{12}^{(4)}|=-79\end{array}$ $\begin{array}{rl}& \frac{17x_{12}^{(2)}}{2}+\frac{9x_{22}^{(2)}}{2}+\frac{9y_{11}^{(2)}}{2}+\frac{1}{2}(2y_{12}^{(2)}-3y_{12}^{(3)})-\frac{3|x_{12}^{(2)}|}{2}-\frac{|x_{22}^{(2)}|}{2}-\frac{|y_{11}^{(2)}|}{2}-\frac{1}{2}|2y_{12}^{(2)}+3y_{12}^{(3)}|=-22\end{array}$ $\begin{array}{rl}& \frac{17x_{12}^{(3)}}{2}+\frac{9x_{22}^{(3)}}{2}+\frac{9y_{11}^{(3)}}{2}+\frac{1}{2}(-3y_{12}^{(2)}+2y_{12}^{(3)})+\frac{3|x_{12}^{(3)}|}{2}+\frac{|x_{22}^{(3)}|}{2}+\frac{|y_{11}^{(3)}|}{2}+\frac{1}{2}|-3y_{12}^{(2)}-2y_{12}^{(3)}|=80\end{array}$ $\begin{array}{rl}& \frac{17x_{12}^{(4)}}{2}+4x_{22}^{(4)}+5y_{11}^{(4)}+\frac{1}{2}(-6y_{12}^{(1)}+4y_{12}^{(4)})+\frac{7|x_{12}^{(4)}|}{2}+2|x_{22}^{(4)}|+2|y_{11}^{(4)}|+\frac{1}{2}|-6y_{12}^{(1)}-4y_{12}^{(4)}|=142\end{array}$ $\begin{array}{rl}& \frac{1}{2}(-2x_{11}^{(1)}-10x_{11}^{(4)})+\frac{1}{2}(4x_{21}^{(1)}-3x_{21}^{(4)})+4y_{21}^{(1)}-\frac{7y_{22}^{(4)}}{2}-\frac{1}{2}|-2x_{11}^{(1)}+10x_{11}^{(4)}|-\frac{1}{2}|4x_{21}^{(1)}+3x_{21}^{(4)}|\\ & -3|y_{21}^{(1)}|-\frac{3|y_{22}^{(4)}|}{2}=-108\end{array}$ $\begin{array}{rl}& -\frac{11x_{11}^{(3)}}{2}+\frac{1}{2}(3x_{21}^{(2)}-2x_{21}^{(3)})+\frac{9y_{21}^{(2)}}{2}-\frac{7y_{22}^{(3)}}{2}-\frac{3|x_{11}^{(3)}|}{2}-\frac{1}{2}|3x_{21}^{(2)}+2x_{21}^{(3)}|-\frac{|y_{21}^{(2)}|}{2}-\frac{|y_{22}^{(3)}|}{2}=-50\end{array}$ $\begin{array}{rl}& -\frac{11x_{11}^{(2)}}{2}+\frac{1}{2}(-2x_{21}^{(2)}+3x_{21}^{(3)})+\frac{9y_{21}^{(3)}}{2}-\frac{7y_{22}^{(2)}}{2}+\frac{3|x_{11}^{(2)}|}{2}+\frac{1}{2}|-2x_{21}^{(2)}-3x_{21}^{(3)}|\\ & +\frac{|y_{21}^{(3)}|}{2}+\frac{|y_{22}^{(2)}|}{2}=23\end{array}$ $\begin{array}{rl}& -6x_{11}^{(1)}+\frac{1}{2}(-3x_{21}^{(1)}+4x_{21}^{(4)})+4y_{21}^{(4)}-\frac{7y_{22}^{(1)}}{2}+4|x_{11}^{(1)}|+\frac{1}{2}|-3x_{21}^{(1)}-4x_{21}^{(4)}|+3|y_{21}^{(4)}|+\frac{3|y_{22}^{(1)}|}{2}=86\end{array}$ $\begin{array}{rl}& \frac{1}{2}(-2x_{12}^{(1)}-10x_{12}^{(4)})+\frac{1}{2}(4x_{22}^{(1)}-3x_{22}^{(4)})+5y_{21}^{(1)}+\frac{1}{2}(4y_{22}^{(1)}-6y_{22}^{(4)})\\ & -\frac{1}{2}|-2x_{12}^{(1)}+10x_{12}^{(4)}|-\frac{1}{2}|4x_{22}^{(1)}+3x_{22}^{(4)}|-2|y_{21}^{(1)}|-\frac{1}{2}|4y_{22}^{(1)}+6y_{22}^{(4)}|=-138\end{array}$ $\begin{array}{rl}& -\frac{11x_{12}^{(3)}}{2}+\frac{1}{2}(3x_{22}^{(2)}-2x_{22}^{(3)})+\frac{9y_{21}^{(2)}}{2}+\frac{1}{2}(2y_{22}^{(2)}-3y_{22}^{(3)})-\frac{3|x_{12}^{(3)}|}{2}-\frac{1}{2}|3x_{22}^{(2)}+2x_{22}^{(3)}|\\ & -\frac{|y_{21}^{(2)}|}{2}-\frac{1}{2}|2y_{22}^{(2)}+3y_{22}^{(3)}|=-68\end{array}$ $\begin{array}{rl}& -\left(\frac{11x_{12}^{(2)}}{2}\right)+\frac{1}{2}(-2x_{22}^{(2)}+3x_{22}^{(3)})+\frac{9y_{21}^{(3)}}{2}+\frac{1}{2}(-3y_{22}^{(2)}+2y_{22}^{(3)})+\frac{3|x_{12}^{(2)}|}{2}+\frac{1}{2}|-2x_{22}^{(2)}-3x_{22}^{(3)}|\\ & +\frac{|y_{21}^{(3)}|}{2}+\frac{1}{2}|-3y_{22}^{(2)}-2y_{22}^{(3)}|=10\end{array}$ $\begin{array}{rl}& -6x_{12}^{(1)}+\frac{1}{2}(-3x_{22}^{(1)}+4x_{22}^{(4)})+5y_{21}^{(4)}+\frac{1}{2}(-6y_{22}^{(1)}+4y_{22}^{(4)})\\ & +4\left|x_{12}^{(1)}\left|+\frac{1}{2}\right|-3x_{22}^{(1)}-4x_{22}^{(4)}|+2|y_{21}^{(4)}\left|+\frac{1}{2}\right|-6y_{22}^{(1)}-4y_{22}^{(4)}\right|=50\end{array}$