7th International Conference on Recent Advances in Pure and Applied Mathematics Goddess of Bodrum Isis Hotel, Bodrum/Muğla, TURKEY September 25-28, 2020
Proceeding Book of ICRAPAM (2020)
Editör Ekrem SAVAS
Associate Editors Mahpeyker OZTURK, Emel Aşıcı Veli CAPALI
Date of Publication: 20.12.2020
ISBN Number:978-625-409-857-4
1 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Goddess of Bodrum Isis Hotel Bodrum/Muğla, TURKEY ICRAPAM CONFERENCE PROCEEDING
PREFACE
International Conference on Recent Advances in Pure and Applied Mathematics (ICRAPAM 2020) was held in Bodrum Muğla, Turkey, at the Goddess of Bodrum Isis Hotel Bodrum/Muğla Turkey, from September 25-28, 2020. It was the 7th edition of such conferences. The chairman of the Organizing Committee of ICRAPAM 2020 was Associate Prof. Mahpeyker Ozturk, and the Scientic Committee consisted of mathematicians from 20 countries. 200 participants from 130 countries attended the conference and 120 papers have been presented, including 6 plenary lectures. The conference was devoted to almost all fields of mathematics and variety of its applications. This issue of the proceeding contains 12 papers presented at the conference and selected by the usual editorial procedure of scientific committee. We would like to express our gratitude to the authors of articles published in this issue and to the referees for their kind assistance and help in evaluation of contributions. I would like to thank to the following my colleagues and students who helped us at every stage of International Conference on Recent Advances in Pure and Applied Mathematics (ICRAPAM 2020).
Editor: Ekrem SAVAS
Usak University, Usak – Turkey
1 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Goddess of Bodrum Isis Hotel Bodrum/Muğla, TURKEY ICRAPAM CONFERENCE PROCEEDING
AIM OF THE CONFERENCE
International Conference on Recent Advances in Pure and Applied Mathematics (ICRAPAM 2020) is aimed to bring researchers and professionals to discuss recent developments in both theoretical and applied mathematics and to create a professional knowledge exchange platform between mathematicians.
SCOPUS
Prospective authors are invited to submit their one-page abstracts on the related, but not limited, following topics of interest:
Numerical Analysis, Ordinary and Partial Differential Equations, Scientific computing, Boundary Value Problems, Approximation Theory, Sequence Spaces and Summability, Real Analysis, Functional Analysis, Fixed Point Theory, Optimization, Geometry, Computational Geometry, Differential Geometry, Applied Algebra, Combinatorics, Complex Analysis, Flow Dynamics, Control, Mathematical modelling in scientific disciplines, Computing Theory, Numerical and Semi-Numerical Algorithms, Game Theory, Operations Research, Optimization Techniques, Fuzzy sequence spaces, Symbolic Computation, Fractals and Bifurcations, Analysis and design tools, Cryptography, Number Theory and Mathematics Education, Finance Mathematics, Fractional Dynamics, Fuzzy systems and fuzzy control, Dynamical systems and chaos, Biomathematics & modeling. Soft Computing, Cryptology & Security Analysis, Image Processing, etc.
PROCEEDING BOOK
The full texts contained in this proceeding book contain all oral presentations in ICRAPAM 2020 Conference.
2 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Goddess of Bodrum Isis Hotel Bodrum/Muğla, TURKEY ICRAPAM CONFERENCE PROCEEDING
COMMITTEE
Honorary Chair
Prof. Dr. Ekrem SAVAŞ Rector of Uşak University, Turkey
Prof. Dr. Vatan KARAKAYA Rector of Ahi Evran University, Turkey
Prof. Dr. Hamdullah ŞEVLİ Rector of Yüzüncü Yıl University, Turkey
Prof. Dr. Billy E. RHOADES Indiana University Bloomington, USA
Prof. Dr. Gradimir V. MILOVANOVİĆ Serbian Academt OF Science and Arts,SERBIA
Scientific Committee
Prof. Dr. Abdelmejid Bayad Prof. Dr. Pratulananda Das (France) (India) Prof. Dr. Abedallah M. Rababah Prof. Reza Saadati (Jordan) (Iran)
Prof. Dr. Alberto Manuel Tavares Prof. Roman Dwilewicz (Portugal) (Usa)
Prof. Dr. Amiran Gogatishvili Prof. Dr. S.sadiq Basha (Czech Republic) (India)
Prof. Dr. Ants Aasma Tallinn Prof. Dr. Taras Banakh (Estonia) (Ukraınıan)
Prof. Dr. Hüseyin Çakallı Prof.dr. Wasfi Shatanawi (Turkey) (Jordan)
Prof. Dr. Hemen Dutta Prof. Dr. Wutiphol Sintunavarat (India) (Thailand) Prof. Dr. Emine Misirli Prof. Dr. Varga Kalantarov (Turkey) (Turkey) Prof. Dr. Erdal Karapinar Prof. Dr. Vasil Angelov (Taiwan) (Bulgaria) Prof. Dr. Fahreddin Abdullayev Prof. Dr. Vasile Berinde (Turkey) (Romania)
Prof. Dr. Fairouz Tchier Prof. Dr. Yilmaz Şimşek (Saudi Arabia) (turkey)
3 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Goddess of Bodrum Isis Hotel Bodrum/Muğla, TURKEY ICRAPAM CONFERENCE PROCEEDING
Prof. Dr. Ibrahim Çanak Prof. Dr. Zoran D. Mitrovic (Turkey) (Vietnam) Prof. Dr. Oktay Duman Prof. Dr. Eman Samir Bhaya (Turkey) (Iraq) Prof. Dr. Ishak Altun Prof. Dr. Ayhan Şerbetçi (Turkey) (Turkey) Prof. Dr. Jeff Connor Prof. Dr. Ismail Ekincioğlu (Usa) (Turkey) Prof. Dr. Lubomira Softova Prof. Dr. Erhan Deniz (Italy) (Turkey) Prof. Dr. Manuel De La Sen Prof. Dr. Rifat Çolak (Spain) (Turkey) Prof. Dr. Martin Bohner Prof. Dr. Ammar Isam Edress (Usa) (Iraq) Prof. Dr. Marc Lassonde Prof. Dr. Sadek Bouroubi (Guadeloupe) (Algeria) Prof.dr. Mehmet Dik Prof. Dr. Hamlet Guliyev (Usa) (Azerbaijan) Prof.dr. Metin Başarir Prof. Dr. Akbar Aliyev (Turkey) (Azerbaijan) Prof.dr. Mikail Et Assoc.prof. Yusuf Zeren (Turkey) (Turkey) Prof. Dr. M. Mursaleen Assist. Prof.dr. Arzu Akgül (India) (Turkey) Prof.dr. Naseer Shahzad Assoc. Prof. Dr. Şükran Konca (Saudi Arabia) (Turkey) Prof.dr. Necip Şimşek Assoc.prof. Dr. Ulaş Yamanci (Turkey) (Turkey) Prof. Dr. Poom Kumam Assist. Prof. Dr. Vuqar Mehrabov (Thailand) (Azerbaijan)
4 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Goddess of Bodrum Isis Hotel Bodrum/Muğla, TURKEY ICRAPAM CONFERENCE PROCEEDING
Organizing Committee
Prof. Dr. Debasis Giri Assoc. Prof. Dr. Emel Aşici (India) (Turkey) Prof. Dr. Maria Zelster Assoc. Prof. Dr. Mahpeyker Öztürk (Estonia) (Turkey) Prof. Dr. Mehmet Gurdal Assoc. Prof. Dr. Moosa Gabeleh (Turkey) (Iran) Prof. Dr. Ram Mohapatra Assoc. Prof. Dr. Murat Kirişçi (Usa) (Turkey) Prof. Dr. Richard F. Patterson Assoc.prof. Dr. Narin Petrot (Usa) (Thailand) Prof. Dr. Werner Varnhorn Assoc prof. Dr. Olivier Olela Otafudu (Germany) (Witwatersrand South Africa) Prof. Dr. Hari Mohan Srivastava Assoc. Prof. Dr. Rahmet Savaş (Canada) (Turkey) Prof. Dr. Naim Braha Assist. Prof. Dr. Veli Çapali (Chuang) (Turkey) Prof. Dr. Mujahid Abbas Assist. Prof. Dr. Banu Güntürk (Pakistan) (Turkey) Prof. Dr. Ali Akhmedov Assist. Prof. Dr. Serkan Araci (Azerbaijan) (Turkey) Prof. Dr. Ziyatkhan Aliyev Dr. Abdurrahman Büyükkaya (Azerbaijan) (Turkey) Assoc. Prof. Dr. Azhar Hussain Dr. Melek Erişçi Büyükkaya (Pakistan) (Turkey) Assoc. Prof. Bancha Panyanak Dr. Rabia Savaş (Thailand) (Turkey) Assoc. Prof. Dr. Edixon M. Rojas Dr. Sefa Anil Sezer (Colombia) (Turkey)
5 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Goddess of Bodrum Isis Hotel Bodrum/Muğla, TURKEY ICRAPAM CONFERENCE PROCEEDING
INVITED SPEAKERS
6 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Goddess of Bodrum Isis Hotel Bodrum/Muğla, TURKEY ICRAPAM CONFERENCE PROCEEDING
7 /
CONTENT
Preface 2
Committee 3
Invited Speakers 6
Content 8
Some investigations on nullnorms on bounded lattices Emel Aşıcı 9
On solving the symmetric regularized long-wave equation by (1 / Gꞌ)- expansion method 12 Gizem Aydin, Esin Ilhan, Haci Mehmet Baskonus, Hasan Bulut The Bounds for the Length Between Dirichlet and the Periodic Eigenvalues of Hill’s Equation with Symmetric Single Well Potential 19 Elif Başkaya On Solving the (2+1)-Dimensional Zoomeron Equation by (1/Gꞌ) -Expansion Method 24 Hasan Bulut,Gizem Aydin
Some Comparisons of BLUPs under General Linear Random Effects Models Melek Eriş Büyükkaya, Nesrin Güler 29
On Generalized Statistical Convergence via Ideal in Cone Metric Spaces 36 Isıl Açık Demirci, Mehmet Gürdal
I-Statistical Limit Superior and I-Statistical Limit Inferior of Triple Sequences Mualla Birgül Huban, Mehmet Gürdal and Ekrem Savaş 42
I-Statistically Localized Sequences of Weighted g via Modulus Functions in 2-Normed Space Mualla Birgül Huban, Mehmet Gürdal and Ekrem Savaş 50
Green’s Function for Hill’s Equation with Symmetric Single Well Potential 57 Ayşe Kabataş N-dimensional bound state solutions of the hyperbolic potential function in approximate analytic form 63 Aysel Özfidan Almost Lacunary P- Bounded Variation and Matrix Transformation 69 Rabia Savaş
A new sequence spaces defined by bounded variation in n- normed spaces 73 Ekrem Savaş
Hankel matrices for blind image deconvolution 78 Belhaj Skander, Alsulami Abdulrahman
Asymptotically Lacunary Statistical Equivalent Functions on Time Scales 84 Bayram Sözbir, Selma Altundağ
8 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Bodrum Isis Hotel in Mugla, TURKEY ICRAPAM CONFERANCE PROCEEDING
Some investigations on nullnorms on bounded lattices
Emel A¸sıcı1 1Department of Software Engineering, Faculty of Technology, Karadeniz Technical University, 61830 Trabzon, Turkey
Keywords Abstract: In this paper, we study an order induced by nullnorms on bounded lattices. Bounded lattice, We study on the set of incomparable elements with respect to the F-partial order which Nullnorm, L is denoted by KF . Also, we give some examples for clarity. Finally, we obtain some conclusions for t-norms and t-conorms on bounded lattices.
1. Introduction and motivation
Triangular norms (briefly t-norms) and triangular conorms (briefly t-conorms) were introduced by Schweizer and Sklar [19] in the study of probabilistic metric spaces as a special kind of associative functions defined on the unit interval [0,1]. Although the t-norms and t-conorms were strictly defined on the unit interval [0,1], they were mostly studied on bounded lattices. Nullnorms and t-operators were introduced in [15] and [8], respectively, which are also generalizations of t-norms and t-conorms. It has been pointed out that nullnorms and t-operators are equivalent [16]. Nullnorms are generalization of triangular norms and triangular conorms with a zero element in the interior of the unit interval and have to satisfy an additional condition. Such nullnorms are interesting not only from a theoretical point of view, but also for their applications, since they have been proved to be useful in several fields like expert systems, neural networks, fuzzy quantifiers. Also, nullnorms that has to be used as aggregators or in fuzzy logic maintain as many logical properties as possible. A natural order for semigroups was defined in [17]. Similarly, a partial order defined by means of t-norms on a bounded lattice was introduced in [13]: x T y :⇔ T(`,y) = x for some ` ∈ L, where L is a bounded lattice, x,y of a bounded lattice L and T is a t-norm on L. This partial order T is called a T-partial order on L. Then, an order F induced by a nullnorm on a bounded lattice L was defined in [6]. Also the set of incomparable elements with respect to the F-partial order for any nullnorm on [0,1] was defined and properties of the introduced set was investigated. In this paper, we investigate some properties of an order induced by nullnorms on bounded lattices. The paper is organized as follows. We shortly recall some basic notions in Section 2. In Section 3, we study on the set of incomparable elements L with respect to the F-partial order which is denoted by KF . Then, we give example and counter example. Finally, by using Proposition 3.2, we obtain some conclusions for t-norms and t-conorms on bounded lattices.
2. Preliminaries
A lattice [2, 3, 7, 9–11] is a partially ordered set (L,≤) in which each two element subset {x,y} has an infimum, denoted as x ∧ y and a supremum, denoted as x ∨ y. A bounded lattice (L,≤,0,1) is a lattice that has the bottom and top elements written as 0 and 1, respectively. Given a bounded lattice (L,≤,0,1) and a,b ∈ L, if a and b are incomparable, in this case, we use the notation a k b. We denote the set of elements which are incomparable with a by Ia. So Ia = {x ∈ L | x k a}. Given a bounded lattice (L,≤,0,1) and a,b ∈ L, a ≤ b, a subinterval [a,b] of L is defined as [a,b] = {x ∈ L | a ≤ x ≤ b} Similarly, [a,b) = {x ∈ L | a ≤ x < b}, (a,b] = {x ∈ L | a < x ≤ b} and (a,b) = {x ∈ L | a < x < b}. Definition 2.1. [4, 18] Let (L,≤,0,1) be a bounded lattice. A t-norm T is a binary operation on L which is commutative, associative, increasing with respect to both variables and it satisfies T(x,1) = x for all x ∈ L.
Example 2.2. [14] The following are the four basic t-norms TM, TP, TL and TD on the real unit interval [0,1] given by, respectively,
∗ Corresponding author: [email protected] 9 TM(x,y) = min(x,y), TP(x,y) = xy, TL(x,y) = max(x + y − 1,0), ( 0 if (x,y) ∈ [0,1)2 , TD(x,y) = min(x,y) otherwise. Extremal t-norms T∧ and TW on L are defined as follows, respectively: x y = 1 , T∧(x,y) = x ∧ y TW (x,y) = y x = 1, 0 otherwise.
Definition 2.3. [18] Let (L,≤,0,1) be a bounded lattice. A triangular conorm S is a binary operation on L which is commutative, associative, increasing with respect to both variables and it satisfies S(x,0) = x for all x ∈ L.
Example 2.4. [3, 14] The following are the four basic t-conorms SM, SP, SL and SD on the real unit interval [0,1] given by, respectively, SM(x,y) = max(x,y), SP(x,y) = x + y − xy, SL(x,y) = min(x + y,1), ( 1 if (x,y) ∈ (0,1]2 , SD(x,y) = max(x,y) otherwise. Similarly, extremal t-conorms S∨ and TW on L are defined on bounded lattice. Example 2.5. [14] The t-norm T nM on [0,1] is defined as follows: ( 0 x + y ≤ 1 , T nM(x,y) = min(x,y) otherwise. is called nilpotent minimum t-norm. The t-norm T ∗ on [0,1] is defined as follows: ( 0 (x,y) ∈ (0,k)2 , T ∗(x,y) = 0 < k < 1 min(x,y) otherwise. Definition 2.6. [5, 8, 12] A nullnorm is a binary operator F : L2 → L which is commutative, associative, non-decreasing in each variable and there exists some element a ∈ L such that F(x,0) = x for all x ≤ a, F(x,1) = x for all x ≥ a. Clearly, F is a t-norm if a = 0 and a t-conorm a = 1. It is easy to show that F(x,a) = a for all x ∈ L. Therefore, a ∈ L is the zero (absorbing) element for F. It is clear that F(1,0) = a. In this study, for the sake of brevity, the set [0,a) × (a,1] ∪ (a,1] × [0,a) for a ∈ L\{0,1} is denoted by Da, i.e., Da = [0,a) × (a,1] ∪ (a,1] × [0,a) for a ∈ L\{0,1}. Proposition 2.7. [12] Let (L,≤,0,1) be a bounded lattice, a ∈ L \{0,1} and F be a nullnorm with zero element a on L. Then ∗ 2 (i) S = F |[0,a]2 : [0,a] → [0,a] is a t-conorm on [0,a]. ∗ 2 (ii) T = F |[a,1]2 : [a,1] → [a,1] is a t-norm on [a,1]. Definition 2.8. [13] Let L be a bounded lattice and T be a t-norm on L. The order defined as follows is called a T− partial order (triangular order) for t-norm T:
x T y :⇔ T(`,y) = x for some ` ∈ L. Definition 2.9. [6] Let (L,≤,0,1) be a bounded lattice and F be a nullnorm with zero element a on L. Define the following relation, for x,y ∈ L, as if x,y ∈ [0,a] and there exist k ∈ [0,a] such that F(x,k) = y or, x F y :⇔ if x,y ∈ [a,1] and there exist ` ∈ [a,1] such that F(y,`) = x or, (1) if (x,y) ∈ L∗ and x ≤ y, ∗ where Ia = {x ∈ L | x k a} and L = [0,a] × [a,1] ∪ [0,a] × Ia ∪ [a,1] × Ia ∪ [a,1] × [0,a] ∪ Ia × [0,a] ∪ Ia × [a,1] ∪ Ia × Ia.
Note: The partial order F in (1) is called F-partial order on L.
Lemma 2.10. [6] Let (L,≤,0,1) be a bounded lattice. For all nullnorms F and all x ∈ L it holds that 0 F x, x F x and x F 1.
Proposition 2.11. [6] Let (L,≤,0,1) be a bounded lattice and F be a nullnorm on L. If x F y for any x,y ∈ L, then x ≤ y. 10 L 3. About the set KF on any bounded lattice
L In this section we study on the set of incomparable elements with respect to the F-partial order which is denoted by KF . We give Proposition 3.2 and by using Proposition 3.2, we obtain some conclusions for t-norms and t-conorms on bounded lattices.
L Definition 3.1. [1] Let F be a nullnorm on (L,≤,0,1) with zero element a and let KF be defined by L KF = {x ∈ L\{0,1} | for some y ∈ L\{0,1}, [x < y implies x F y]or
[y < x implies y F x]orx k y}. Proposition 3.2. Let (L,≤,0,1) be a bounded lattice and F be a nullnorm on L with zero element a. If there exist two L elements in L which are incomparable with each other, then KF 6= /0. Example 3.3. Consider the lattice (L = {0, p,m,a,n,1},≤,0,1) whose lattice diagram is displayed in Figure 1. Consider
Figure 1. The order ≤ on L the nullnorm on L with zero element a defined as follows:
x ∨ y ,(x,y) ∈ [0,a]2 y ,x = 1 and y ≥ a F(x,y) = x ,y = 1 and x ≥ a a ,otherwise
It is clear that KF = {m,a,n}. L Remark 3.4. The converse of Proposition 3.2 is not correct. That is, although the set KF 6= /0, it need not be the case that all elements in L are incomparable. To illustrate this claim, we shall give the following example.
2 1 Example 3.5. Consider the nullnorm F1 : [0,1] → [0,1] with zero element defined by 2 2
2 max(x,y) ,(x,y) ∈ [ , 1 ] 0 2 1 1 2 F1 (x,y) = 2 ,(x,y) ∈ [ 2 ,1) ∪ D 1 2 2 min(x,y) ,otherwise.
1 Then, KF1 = ( 2 ,1), but since L = [0,1] is a chain all elements are comparable with each other in L. Let us show that 2 1 1 K 1 . Let x ∈ KF1 . We will show that x ∈ ( ,1). Suppose that x ∈/ ( ,1). Since x ∈ KF1 , there exists an element F1 =( 2 ,1) 2 2 2 2 2 y ∈ (0,1) such that x < y and x F1 y or y < x and y F1 x. 2 2 1 1 Let x < y and x F1 y. Since x ∈/ ( 2 ,1), we have x ≤ 2 . Since x < y, we have that max(x,y) = y. By the definition of F1 , 2 2 we obtain that y = max(x,y) = F1 (x,y) 2
Then, it is obtained that x F1 y, a contradiction. 2 If y < x, then similarly it can be obtained y x, a contradiction. So, we have x ∈ ( ,1 ). Thus, it is obtained that F1 21 2 11 1 KF1 ⊆ ( 2 ,1). 2 1 1 Let x ∈ ( 2 ,1) and x < y < 1. Then x F1 y. Suppose that x F1 y. Then, there exists an element k ∈ [ 2 ,1] such that 2 2 1 2 x = F1 (y,k). Since x 6= y, it is not possible k = 1. Since (y,k) ∈ [ ,1) , we have that 2 2 1 x = F1 (y,k) = , 2 2
1 a contradiction. So, we have that x F1 y, x ∈ KF1 . So it is obtained that ( 2 ,1) ⊆ KF1 . Therefore, we have that 2 2 2 1 KF1 = ( 2 ,1). 2 Corollary 3.6. Let (L,≤,0,1) be a bounded lattice and T be a t-norm on L. If there exist two elements in L which are L incomparable with each other, then KT 6= /0. Corollary 3.7. Let (L,≤,0,1) be a bounded lattice and S be a t-conorm on L. If there exist two elements in L which are L incomparable with each other, then KS 6= /0. References
[1] A¸sıcı,E. 2019. An extension of the ordering based on nullnorms. Kybernetika, 55(2) (2019), 217-232. [2] A¸sıcı,E. 2019. The equivalence of uninorms induced by the U-partial order. Hacettepe Journal of Mathematics and Statistics. 48 (2019), 439-450. [3] A¸sıcı,E., Mesiar, R. 2020. New constructions of triangular norms and triangular conorms on an arbitrary bounded lattice. International Journal of General Systems. 49(2) (2020), 143-160. [4] A¸sıcı,E., Mesiar, R. 2021. On the construction of uninorms on bounded lattices. Fuzzy Sets and Systems. doi: 10.1016/j.fss.2020.02.007, in press. [5] A¸sıcı,E., Mesiar, R. 2019. Direct product of nullnorms on bounded lattices. Journal of Intelligent and Fuzzy Systems. 36(6) (2019), 5745-5756. [6] A¸sıcı,E. 2017. An order induced by nullnorms and its properties. Fuzzy Sets and Systems. 325 (2017), 35-46. [7] Birkhoff, G. 1967. Lattice Theory. 3 rd edition, Providence, 1967. [8] Calvo, T. De Baets, B., Fodor, J. 2001. The functional equations of Frank and Alsina for uninorms and nullnorms. Fuzzy Sets and Systems. 120 (2001) 385-394. [9] Çaylı, G.D. 2020. Nullnorms on bounded lattices derived from t-norms and t-conorms. Information Sciences. 512 (2020), 1134-1154. [10] Çaylı, G.D. 2019. Alternative approaches for generating uninorms on bounded lattices. Information Sciences. 488 (2019), 111-139. [11] Çaylı, G.D. 2019. New methods to construct uninorms on bounded lattices. International Journal of Approximate Reasoning. 115 (2019), 254-264. [12] Karaçal, F., Ince, M.A., Mesiar, R. 2015. Nullnorms on bounded lattices. Information Sciences. 325(2015), 227-236. [13] Karaçal, F., Kesicioglu,˘ M.N. 2011. A T-partial order obtained from t-norms. Kybernetika. 47(2011) 300-314. [14] Klement E.P, Mesiar, R., Pap, E. 2000. Triangular Norms. Kluwer Academic Publishers, Dordrecht 2000. [15] Mas, M., Mayor, G., Torrens, J. 1999. t-operators. International Journal of Uncertainty, Fuzziness and Knowledge- Based Systems. 7 (1999) 31-50. [16] Mas, M., Mayor, G., Torrens, J. 2002. The distributivity condition for uninorms and t-operators. Fuzzy Sets and Systems. 128 (2002) 209-225. [17] Mitsch, H. 1986. A natural partial order for semigroups, Proceedings of the American Mathematical Society, 97 (1986) 384-388. [18] Saminger, S. 2006. On ordinal sums of triangular norms on bounded lattices. Fuzzy Sets and Systems. 157 (2006), 1403-1416. [19] Schweizer, B., Sklar, A. Statistical metric spaces. Pacific Journal of Mathematics. 10 (1960), 313-334.
12 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Bodrum Isis Hotel in Mugla, TURKEY ICRAPAM CONFERENCE PROCEEDING
1 On solving the symmetric regularized long-wave equation by (G 0 ) expansion method
Gizem Aydin1, Esin Ilhan2, Haci Mehmet Baskonus∗3, Hasan Bulut1 1Department of Mathematics, Firat University, Elazig, Turkey 2Kirsehir Ahi Evran University, Kirsehir, Turkey 3Faculty of Education, Harran University, Sanliurfa, Turkey
1 Keywords Abstract: In this article, the ( G0 ) expansion method is considered to find new wave The symmetric regularized solutions for nonlinear partial differential equations named symmetric regularized long- long-wave equation, wave equation. Many new traveling wave solutions to the given equation are found. 1 The ( G0 ) expansion method, In addition, 2D and 3D graphs are plotted by selecting the values of the appropriate Analytical solutions, parameters. Traveling wave solutions.
1. Introduction
It is well known that nonlinear partial differential equations (NLPDEs) define a wide variety of phenomena not only in physics but also in biology, applied science, chemistry, engineering, and many other fields. The investigation of moving wave solutions for NLPDEs plays an important role in the investigation of nonlinear physical phenomena. In recent years, many scientists have applied various methods to obtain analytical and numerical solutions for the nonlinear partial G 1 differential equations in mathematical physics, such as the ( G0 ) expansion method [1-2], the ( G0 ) expansion method [3-4], the tanh-coth method [5], the sine-Gordon expansion method [6-8], the Laplace perturbation method [9], the extended simple equation method [10], the integral approximation [11], the hyperbolic tangent function expansion method [12], the generalized expansion method [13], the Bernoulli sub-equation method [14-17], the Hirota’s bilinear transformation method [18], the extended sinh-Gordon method[19-20], the modified simple equation method[21], a direct algebraic method[22], an inverse spectral transform technique[23], the multiplier method [24], the generalized Kudryashov method [25], the sechp −tanhp functions method [26], the homotopy method [27-28], and many other numerical and analytical methods 1 [29-40]. The main purpose of this work is to apply the ( G0 ) expansion method to the symmetric regularized long-wave (SRLW) equation, u2 u − u + ( ) − u = 0, (1) tt xx 2 xx xxtt which was first introduced by Seyler and Fenstermacher in 1983 [41] from a weakly nonlinear analysis of the cold-electron plasma equations suitable for a strongly magnetized non-relativistic electron beam such that the fluid motion is restricted to one direction [41,42]. The SRLW equation represents a weekly nonlinear ion-acoustic and space-charge waves, and the real-valued u(x,t) refers to the dimensionless fluid velocity with a decay condition limx→∞ u = 0 . The SRLW equation is explicitly symmetry in the and derivatives and is quite similar to the regularized-long-wave equation that represents shallow-water waves and plasma drift waves [41, 42]. The radial basis functions collocation method has been used to study the single solitary wave solution, the interaction of two positive solitary waves, and the clash of two solitary waves [43, 44].The existence and uniqueness of solutions for the SRLW have been established, e.g., in [45,46]. On the other hand, many authors tried to find the approximation and exact solution to the SRLW equation, such as the radial basis functions collocation method[47]. The modified cubic B-Splines has been implemented, then very accurate solutions in different settings of parameters have been provided [48]. The four-level linear implicit finite difference has been used G to obtain some solitary wave solutions [49]. The ( G0 ) -expansion method has been utilized to find exact solutions to the SRLW equation [50]. The Lie symmetry approach along with the simplest equation and exp-function methods are used to obtain solutions of the SRLW equation [51]. The existence and nonlinear stability of periodic traveling wave solutions of the cnoidal type of the SRLW equation have been investigated in [52].We organized the paper as follows: In Section 2, 1 the basic concepts of the ( G0 ) -expansion method are addressed. In Section 3, the application of the used method to the equation under consideration is presented. The Conclusion is given in Section 4.
∗ Corresponding author: [email protected] 12 /
1 2. The( G0 ) expansion method
1 In this portion, the general details of the ( G0 ) - expansion method are provided. Let us consider a three-variable general form of the nonlinear partial differential equation (PDE) as
P(u,ut , uux, uuy, uuxx,···) = 0, (2) transform equation (2) with the aid of u(x,y,t) = U(ξ),ξ = κ1x+κ2y−κ3t in which κ1, κ2 and κ3 are non-zero constants. After transformation, we have a nonlinear ordinary differential equation (ODE) as follows
N(U, U0, U00, U000,....) = 0. (3)
The solution of (3) shall be taken in the following form
m 1 i U(ξ) = a0 + ∑ ai( 0 ) . (4) i=1 G
Here G = G(ξ) is a function that verifies the following second-order ordinary differential equation
G00 + λG0 + µ = 0, (5) here a0,a1,a2,··· ,λ, µ are constants and m is a term of balance. The balance term is a fixed number obtained between the linear term and the non-linear term in any nonlinear ODE that has the highest order. This number is written in place of (4) and then the derivatives required for the solution are obtained. Here we gain a system of algebraic equations by equalizing 1 i the coefficients of ( G0 ) ,(i = 1,2,3,.....m) terms to zero. This algebraic equation scheme is solved automatically or with the aid of a computer program depending on the degree of difficulty. These solutions (4) given by the solution function and moving wave transformation is written in place (2) given by the solution of the PDE is found [41].
The general solution of the ODE (5) could be reported as below
µξ G(ξ) = − −C eλξ +C , (6) λ 1 2 here C1 and C2 are fixed. Depending on the variable ξ the derivative of the (6) solution function is taken and necessary arrangements are made, we have 1 1 = (7) G0 µ λξ − λ +C1e When we transform the algebraic expression provided in (7) to a trigonometric function may write as
1 λ 0 = . (8) G −µλC1(Cosh(ξλ) − Sinh(ξλ))
1 3. Implementation of ( G0 ) expansion method
1 In this sub-section, we apply the ( G0 ) expansion method to the SRLW equation. We start by inserting u(x,t) = U(ξ),ξ = x − ct, (9) into Eq. (9), we get
U2 (c2 − 1)U00 − c( )00 − c2U0000 = 0. (10) 2 Twice integrate the Eq. (10) and setting to zero of constants of integrations, we have c (c2 − 1)U − U2 − c2U00 = 0. (11) 2 Balancing among U2 and U00 in Equation (11), we get .In this case, the Eq. (4) takes the form 1 1 U(ξ) = a + a ( ) + a ( )2, (12) 0 1 G0 2 G0
13 /
1 i using Eq. (12) together with its derivative of Eq. (11). Then, summing the ( G0 ) coefficients with the likely powers and setting every summation to zero, we get a system of algebraic equations as
1 ca2 ( )0 : −a + c2a − 0 = 0, G0 0 0 2 1 ( )1 : −a + c2a − c2λ 2a − ca a = 0, G0 1 1 1 0 1 1 ca2 ( )2 : −3c2λ µa − 1 − a + c2a − 4c2λ 2a − ca a = 0, G0 1 2 2 2 2 0 2 1 ( )3 : −2c2µ2a − 10c2λ µa − ca a = 0, G0 1 2 1 2 1 ca2 ( )4 : −6c2µ2a − 2 = 0, (13) G0 2 2 we solve the above system of equations by the computer package program, we achieve the following families of solution. 12λ µ 12µ2 1 Case-1 When a0 = 0, a1 = −√ , a2 = −√ , c = √ , λ 2−1 λ 2−1 λ 2−1 we can gain the following soliton solution
12µ2 u(x,t) = √ 2 2 µ 1 1 λ − 1 −C1cosh(λ(x − √ )) +C1sinh(λ(x − √ )) λ λ 2−1t λ 2−1t (14) 12λ µ + . √ 2 2 µ 1 1 λ − 1 −C1cosh(λ(x − √ )) +C1sinh(λ(x − √ )) λ λ 2−1t λ 2−1t
2λ 2 12λ µ 12µ2 1 Case-2 When a0 = √ , a1 = √ , a2 = √ , c = −√ , 1+λ 2 1+λ 2 1+λ 2 1+λ 2 we obtain the following soliton solution
2λ 2 12µ2 u(x,t) =√ + 1 + 2 √ 2 λ 2 1 1 µ 1 + λ C1cosh(λ(x + √ )) −C1sinh(λ(x + √ )) − 1+λ 2t 1+λ 2t λ (15) 12λ µ + . √ 2 2 1 1 µ 1 + λ C1cosh(λ(x + √ )) −C1sinh(λ(x + √ )) − 1+λ 2t 1+λ 2t λ √ √ 2 2 2 −1+c Case-3 When a0 = 0, a1 = −12 −1 + c , a2 = 12cµ , λ = c , we have 12cµ2 u(x,t) = √ √ 2 cµ c2−1(x−ct) c2−1(x−ct) √ +C1sinh( ) −C1cosh( ) c2−1 c c √ (16) 12 c2 − 1µ + . √ √ 2 cµ c2−1(x−ct) c2−1(x−ct) √ −C1sinh( ) +C1cosh( ) c2−1 c c √ √ 2 2 2 2 c −1 Case-4 When a0 = −( c ), a1 = −12 +c − 1µ, a2 = −12cµ , λ = c , we get 2 12cµ2 u(x,t) = − + 2c − √ √ 2 c c2−1(x−ct) c2−1(x−ct) cµ C1cosh( ) −C1cosh( ) − √ c c c2−1 √ (17) 12 c2 − 1µ − . √ √ 2 c2−1(x−ct) c2−1(x−ct) cµ C1cosh( ) −C1cosh( ) − √ c c c2−1
14 /
, Figure 1. The topological type of solution to the SRLW equation plotted when λ = −0.5, µ = 0.3,C1 = 0.4 and t = 1 for 2D figure.
, Figure 2. The non-topological type of solution to the SRLW equation plotted when c = 3, µ = 4,C1 = 0.2 and t = 1 for 2D figure.
, Figure 3. The singular type of solution to the SRLW equation plotted when c = 3, µ = 4,C1 = 0.2 and t = 1 for 2D figure.
15 /
, Figure 4. The non-topological type of solution to the SRLW equation plotted when c = −0.9, µ = 0.3,C1 = 0.2and t = 1 for 2D figure.
4. Conclusion
1 In this study, we have successfully obtained some new exact traveling wave solutions by using the ( G0 ) -expansion method to the SRLW equation, which has a very important place in applied sciences. All solutions obtained in this article are 1 verified by the PDE. Chand and Malik [50] have used the ( G0 ) -expansion method and kink-antikink shaped, bell shaped soliton and periodic wave solutions to the SRLW equation have been constructed. In Ref. [51], the Lie symmetry approach, the exp-function method and the simplest equation method have been utilized to construct some solutions to the studied equation. Comparing our results with the results reported in Refs. [51,52], it is concluded that our results are new and different. In the solutions presented, constants are given specific values to obtain graphics that represent the wave. Figures (1-4) symbolize the traveling wave behaviors of the governing model. As these constants change, the position of speed, amplitude, and wavelength of the wave may change.
References
G [1] A. Yokus, M. Tuz, An application of a new version of ( G0 )-expansion method,In AIP Conference Proceedings 2017 Jan 27 (Vol. 1798, No. 1, p. 020165). AIP Publishing LLC. [2] A.Yokus, D. Kaya, Traveling wave solutions of some nonlinear partial differential equations by using extended (G’/G)-expansion method, Istanbul Ticaret Universitesi Fen Bilimleri Dergisi 14, no. 28 (2015): 85. [3] A. Yokus, Comparison of Caputo and conformable derivatives for time-fractional Korteweg-de Vries equation via the finite difference method, Int. J. Mod. Phys. B, 2018. 1 [4] A. Yokus and H. Durur, Complex hyperbolic traveling wave solutions of KuramotoSivashinsky equation using ( G0 ) expansion method for nonlinear dynamic theory, Balikesir Universitesi Fen Bilim. Enstitusu Derg., vol. 21, no. 2, pp. 590-599,2019. [5] A.Irshad,andS.T.Mohyud-Din,Solitary wave solutions for Zoomeron equation, Walailak J. Sci. Technol., vol. 10, no. 2, pp. 201-208, 2012. [6] K. K. Ali, R. Yilmazer, and H. Bulut, Analytical Solutions to the Coupled Boussinesq-Burgers Equations via Sine-Gordon Expansion Method, in 4th International Conference on Computational Mathematics and Engineering Sciences(CMES-2019), p. 233. [7] H. F. Ismael, H. Bulut, and H. M. Baskonus, Optical soliton solutions to the Fokas-Lenells equation via sine-Gordon G0 expansion method and (m+ ( G )))-expansion method, Pramana - J. Phys., vol. 94, no. 1, Dec.2020. [8] K. K. Ali, R. Yilmazer, A. Yokus, and H. Bulut,Analytical solutions for the (3+ 1)-dimensional nonlinear extended quantum Zakharov-Kuznetsov equation in plasma physics, Phys. A Stat. Mech. its Appl., p. 124327,2020. [9] M. Yavuz, N. Ozdemir, and H. M. Baskonus, Solutions of partial differential equations using the fractional operator involving Mittag-Leffler kernel, Eur. Phys. J.Plus, 2018. [10] E. M. E. Zayed, R. M. A. Shohib, and A.G. Al-Nowehy, On solving the (3+ 1)-dimensional NLEQZK equation and the (3+ 1)-dimensional NLmZK equation using the extended simplest equation method, Comput. Math. with Appl.,2019. [11] H. F. Ismael, Linear Stability of Thin Liquid Films flows down on an Inclined Plane using Integral Approximation, Raf. J. Comp. Math’s, vol. 11, no. 2, pp. 13-24,2014.
16 /
[12] S. Liu, Z. Fu, S. Liu, and Q. Zhao, Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations, Phys. Lett. Sect. A Gen. At. Solid State Phys.,2001. [13] R. Sabry, W. M. Moslem, F. Haas, S. Ali, and P. K. Shukla, Nonlinear structures:Explosive, soliton, and shock in a quantum electron-positron-ion magnetoplasma,Phys. Plasmas,2008. [14] T. A. Sulaiman, A. Yokus, N. Gulluoglu, and H. M. Baskonus, Regarding the numerical solutions of the Sharma- Tasso-Olver equation, in ITM Web of Conferences,2018, vol. 22, p. 1036. [15] H. H. Abdulkareem, H. F. Ismael, E. S. Panakhov, and H. Bulut, Some Novel Solutions of the Coupled Whitham- Broer-Kaup Equations, in International Conference on Computational Mathematics and Engineering Sciences, 2019, pp. 200-208. [16] K. K. Ali, R. Yilmazer, H. M. Baskonus, and H. Bulut,Modulation instability analysis andanalytical solutions to the system of equations for the ionsoundand Langmuir waves, Phys. Scr., 2020. [17] H. F. Ismael and H. Bulut, On the Solitary Wave Solutions to the (2+ 1)-Dimensional Davey-Stewartson Equations, in International Conference on Computational Mathematics and Engineering Sciences, 2019, pp.156-165. [18] R. Hirota and J. Satsuma, Soliton solutions of a coupled Korteweg-de Vriesequation,Phys. Lett. A, vol. 85, no. 8-9, pp. 407-408, 1981. [19] W. Gao, H. F. Ismael, S. A. Mohammed, H. M. Baskonus, and H. Bulut, Complex and Real Optical Soliton Properties of the Paraxial Non-linear Schrodinger Equation in Kerr Media With M-Fractional, Front. Phys., vol. 7, Nov.2019. [20] H. Dutta, H. Gunerhan, K. K. Ali, and R. Yilmazer, Exact soliton solutions to the cubic-quartic nonlinear Schrodinger equation with conformable derivative," Front. Phys., vol. 8, p. 62,2020. [21] K. Khan and M. Ali Akbar, "Traveling wave solutions of the (2 + 1)-dimensional Zoomeron equation and the Burgers equations via the MSE method and the Exp- function method, Ain Shams Eng. J.,2014. [22] A. R. Seadawy, Stability analysis for Zakharov-Kuznetsov equation of weakly nonlinear ion-acoustic waves in a plasma, Comput. Math. with Appl., vol. 67, no. 1, pp. 172-180, 2014. [23] F. Calogero and A. Degasperis, Nonlinear evolution equations solvable by the inverse spectral transform.-I, Nuovo Cim. B Ser. 11,1976. [24] T. Motsepa, C. Khalique, and M. Gandarias, Symmetry Analysis and Conservation Laws of the Zoomeron Equation, Symmetry (Basel).,2017. [25] E. Aksoy, A. C. Cevikel, and A. Bekir, Soliton solutions of (2+1)-dimensional time- fractional Zoomeron equation, Optik (Stuttg).,2016. [26] M. Alquran and K. Al-Khaled, Mathematical methods for a reliable treatment of the (2+1)-dimensional Zoomeron equation, Math. Sci.,2013. [27] M. Yavuz, Dynamical Behaviors of Separated Homotopy Method Defined by Conformable Operator, Konuralp J. Math., vol. 7, no. 1, pp.1-6. [28] M. A. Yousif, B. A. Mahmood, K. K. Ali, and H. F. Ismael, Numerical simulation using the homotopy perturbation method for a thin liquid film over an unsteady stretching sheet, Int. J. Pure Appl. Math., vol. 107, no. 2, 2016. [29] H.F.Ismael and K.K.Ali, MHD cass on flow over a nunsteady stretching sheet, Adv. Appl. Fluid Mech., vol. 20, no. 4, pp. 533-541, Oct. 2017. [30] H. F. Ismael, Carreau-Casson fluids flow and heat transfer over stretching plate with internal heat source/sink and radiation, Int. J. Adv. Appl. Sci. J., vol. 6, no. 2, pp. 81- 86, 2017. [31] A. Yokus, "Numerical Solutions of Time Fractional Korteweg-de Vries Equation and Its Stability Analysis," Commun. Fac. Sci. Univ. Ankara Ser. A1 Math. Stat., vol. 68, no. 1, pp. 353-361, 2019. [32] K. K. Ali and A. Varol, "Weissenberg and Williamson MHD Flow over a stretching surface with thermal radiation and chemical reaction,"2019. [33] A. Zeeshan, H. F. Ismael, M. A. Yousif, T. Mahmood, and S. U. Rahman, "Simultaneous effects of slip and wall stretching/shrinking on radiative flow of magneto nanofluid through porous medium," J. Magn., vol. 23, no. 4, pp. 491-498, 2018. [34] A. Yokus, B. Kuzu, and U. Demiroglu, "Investigation of solitary wave solutions for the (3+ 1)-dimensional Zakharov- Kuznetsov equation," Int. J. Mod. Phys. B, vol. 33, no. 29, p. 1950350, 2019. [35] K. K. Ali, H. F. Ismael, B. A. Mahmood, and M. A. Yousif, "MHD Casson fluid with heat transfer in a liquid film over unsteady stretching plate," Int. J. Adv. Appl. Sci., vol. 4, no. 1, pp. 55-58, 2017.
17 /
[36] W. Gao, H. F. Ismael, A. M. Husien, H. Bulut, and H. M. Baskonus, "Optical Soliton Solutions of the Cubic-Quartic Nonlinear Schrodinger and Resonant Nonlinear Schrodinger Equation with the Parabolic Law," Appl. Sci., vol. 10, no. 1, p. 219, Dec. 2019. [37] M. A. Yousif, M. Hatami, and H. F. Ismael, "Heat Transfer Analysis of MHD Three Dimensional Casson Fluid Flow Over a Porous Stretching Sheet by DTM-Pade," Int. J. Appl. Comput. Math., vol. 3, no. S1, pp. 813-828, 2017. [38] H. F. Ismael and N. M. Arifin, "Flow and heat transfer in a maxwell liquid sheet overa stretching surface with thermal radiation and viscous dissipation," JP J. Heat Mass Transf., vol. 15, no. 4, 2018. [39] A. Yokus and S. Gulbahar, "Numerical Solutions with Linearization Techniques of the Fractional Harry Dym Equation," Appl. Math. Nonlinear Sci.,2019. [40] W. Gao and H. F. I. and H. B. and H. M. Baskonus, "Instability modulation for the (2+1)-dimension paraxial wave equation and its new optical soliton solutions in Kerr media," Phys. Scr.,2019. [41] Seyler, C. E., and D. L. Fenstermacher. "A symmetric regularized?long?wave equation." The Physics of fluids 27, no. 1 (1984): 4-7. [42] Dong, Xuanchun. "Numerical solutions of the symmetric regularized-long-wave equation by trigonometric integrator pseudospectral discretization." arXiv preprint arXiv:1109.0764 (2011). [43] Albert, John. "On the decay of solutions of the generalized Benjamin-Bona-Mahony equation." Journal of Mathemati- cal Analysis and Applications 141, no. 2 (1989): 527-537. [44] Amick, Charles J., Jerry L. Bona, and Maria E. Schonbek. "Decay of solutions of some nonlinear wave equations." Journal of Differential Equations 81, no. 1 (1989): 1-49. [45] Bo-Ling, Guo. "The spectral method for symmetric regularized wave equations." Journal of Computational Mathe- matics 5, no. 4 (1987): 297-306. [46] Lin, Chen. "Stability and instability of solitary waves for generalized symmetric regularized-long-wave equations." Physica D: Nonlinear Phenomena 118, no. 1-2 (1998): 53-68. [47] Kaplan, Ayse Gull, and Yilmaz Dereli. "Numerical solutions of the symmetric regularized long wave equation using radial basis functions." Comput. Model. Eng. Sci 84 (2012): 423-438. [48] Mittal, R. C., and A. Tripathi. "Numerical solutions of symmetric regularized long wave equations using collocation of cubic B-splines finite element." Int. J.Computational Methods in Engineering Science and Mechanics 16, no. 2 (2015): 142-150. [49] Yimnet, S., B. Wongsaijai, Thaned Rojsiraphisal, and K. Poochinapan. "Numerical implementation for solving the symmetric regularized long wave equation." Applied Mathematics and Computation 273 (2016): 809-825. [50] F. Chand and A. K. Malik, "Exact Traveling Wave Solutions of Some Nonlinear Equations Using G’/G -expansion Method methods The SRLW equation," International Journal of Nonlinear Science, vol. 14, no. 4, pp. 416-424,2012. [51] Mhlanga, Isaiah Elvis, and Chaudry Masood Khalique. "Exact solutions of the symmetric regularized long wave equation and the Klein-Gordon-Zakharov equations." In Abstract and Applied Analysis, vol. 2014. Hindawi, 2014. [52] Brango, Carlos Banquet. "The symmetric regularized-long-wave equation: Well-posedness and nonlinear stability." Physica D: Nonlinear Phenomena 241, no. 2 (2012): 125-133.
18 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Bodrum Isis Hotel in Mugla, TURKEY ICRAPAM CONFERANCE PROCEEDING
The Bounds for the Length Between Dirichlet and the Periodic Eigenvalues of Hill’ s Equation with Symmetric Single Well Potential
Elif Ba¸skaya∗1 1Karadeniz Technical University, Science Faculty, Department of Mathematics, Turkey
Keywords Abstract: A symmetric single well potential on [0,a] is defined as symmetric with Hill’ s equation, respect to the midpoint a/2 and nonincreasing on [0,a/2]. In a symmetric hydrogen bond Symmetric single well that occurs in many biological structures such as DNA and water, the proton free energy potential, landscape is a symmetric single well potential. Besides, especially in recent years, since Dirichlet eigenvalues, quantum mechanic has gained importance, there are a lot of studies on eigenvalues of The periodic eigenvalues Hill’ s equation and Schrödinger’ s operator with symmetric single well potential such as anharmonic oscillator. The eigenvalues of these equations represent excitation energy and eigenfunctions are named as wavefunction in physics. In this study, the upper and lower bounds for the length between Dirichlet and the periodic eigenvalues are determined for Hill’ s equation when the potential is symmetric single well.
1. Introduction
In this study, we deal with
y00 (x) + [λ − q(x)]y(x) = 0, x ∈ [0,a] (1) where λ is a real parameter; q(x) is a real- valued, continuous and periodic function of periodic a. This equation is known as the standart form of Hill’ s equation. Let y1 and y2 be the linearly independent solutions of Hill’ s equation and initial conditions be as following for these solutions:
0 0 y1 (0,λ) = y2 (0,λ) = 1, y1 (0,λ) = y2 (0,λ) = 0.
D(λ) that is an analytic function of λ , is defined by
0 D(λ) = y1 (a,λ) + y2 (a,λ) is called the discriminant of Hill’ s equation. The periodic problem of (1) is defined with boundary conditions y(0) = y(a) 0 0 and y (0) = y (a). The eigenvalues of this problem are zeros of D(λ) − 2 and shown by {λn}. The semi-periodic problem of (1) is described with boundary conditions y(0) = −y(a) and y0 (0) = −y0 (a). The eigenvalues of this problem are the zeros of D(λ) + 2 and shown by {µn}. The periodic and semi-periodic eigenvalues {λn} and {µn} occur in the order [1]:
−∞ < λ0 < µ0 ≤ µ1 < λ1 ≤ λ2 < µ2 ≤ µ3 < ...
Besides, the eigenvalues of Hill’ s equation with y(0) = y(a) = 0 are denoted by Λn and the eigenvalues of Hill’ s equations 0 0 with y (0) = y (a) = 0 are denoted by νn. These conditions are named as Dirichlet and Neumann boundary conditions, respectively. It is also proven in [1] that for n = 0,1,2,...
µ2n ≤ Λ2n ≤ µ2n+1, λ2n+1 ≤ Λ2n+1 ≤ λ2n+2, (2) and µ2n ≤ ν2n+1 ≤ µ2n+1, λ2n+1 ≤ ν2n+2 ≤ λ2n+2. (3)
All solutions of Hill’ s equation are bounded in (−∞,+∞) for λ ∈ (λ2n, µ2n)∪(µ2n+1,λ2n+1), while all nontrivial solutions
∗ Corresponding author: [email protected] 19 / of Hill’ s equation are unbounded in (−∞,+∞) for λ ∈ (−∞,λ0) ∪ (µ2n, µ2n+1) ∪ (λ2n+1,λ2n+2).
A symmetric single well potential on [0,a] is defined as symmetric for the midpoint a/2 and nonincreasing on [0,a/2]. For example, in a symmetric hydrogen bond that occurs in many biological structures such as DNA and water, the proton free energy landscape is a symmetric single well potential. Besides, especially in recent years, since quantum mechanic has gained importance, there are a lot of studies on eigenvalues of Hill’ s equation and Schrödinger’ s operator with symmetric single well potential such as anharmonic oscillator. The eigenvalues of these equations represent excitation energy and eigenfunctions are named as wavefunction in physics.
Instability intervals for Hill’ s equation with symmetric single well potential have been investigated by many authors. Some important studies of them are [2], [4], [5], [6], [7], [8]. Especially, we refer to [4] includes the length for Λ0 − λ0. In this study, the upper and lower bounds for the length between Dirichlet and the periodic eigenvalues are determined for Hill’ s equation when the potential is symmetric single well. We note that a symmetric single well potential on [0,a] means a continuous function q(x) on [0,a] which is symmetric about x = a/2 and non-increasing on [0,a/2], that is, q(x) = q(a − x) mathematically.
2. Main Results Firstly, we emphasize that q0 (x) exists since a monotone function on an interval I is differentiable almost everywhere on I [3]. Our analysis is based on the following theorem of [8]: The periodic eigenvalues of (1) satisfy, as n → ∞
a/2 Z 1/2 2(n + 1)π a 0 4(n + 1)π λ n+ = − q (t)sin t dt 2 1 a 8(n + 1)2π2 a 0 a/2 a/2 2 Z Z (4) a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 64(n + 1)3π3 0 0 + o n−3 and
a/2 Z 1/2 2(n + 1)π a 0 4(n + 1)π λ n+ = + q (t)sin t dt 2 2 a 8(n + 1)2π2 a 0 a/2 a/2 2 Z Z (5) a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 64(n + 1)3π3 0 0 + o n−3. The semi-periodic eigenvalues of (1) satisfy, as n → ∞
a/2 Z 1/2 (2n + 1)π a 0 2(2n + 1)π µ n = − q (t)sin t dt 2 a 2(2n + 1)2π2 a 0 a/2 a/2 2 Z Z (6) a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 8(2n + 1)3π3 0 0 + o n−3 and
a/2 Z 1/2 (2n + 1)π a 0 2(2n + 1)π µ n+ = + q (t)sin t dt 2 1 a 2(2n + 1)2π2 a 0 a/2 a/2 2 Z Z (7) a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 8(2n + 1)3π3 0 0 + o n−3. The purpose of this study is to prove the following theorem: 20 /
Theorem 2.1. Let q(x) be a symmetric single well potential on [0,a]. Then, the bounds for the length between Dirichlet and the periodic eigenvalues as n → ∞
a/2 2 Z (4n + 1)π 1 0 2(2n + 1)π Λ2n − λ2n ≥ − q (t)sin t dt a2 (2n + 1)π a 0 a/2 Z 1 0 4nπ − q (t)sin t dt 2nπ a 0 a/2 a/2 " # Z Z a a 2 0 0 − − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 4(2n + 1)2π2 16n2π2 0 0 + o n−2 and
a/2 2 Z (4n + 1)π 1 0 2(2n + 1)π Λ2n − λ2n ≤ + q (t)sin t dt a2 (2n + 1)π a 0 a/2 Z 1 0 4nπ − q (t)sin t dt 2nπ a 0 a/2 a/2 " # Z Z a a 2 0 0 − − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 4(2n + 1)2π2 16n2π2 0 0 + o n−2. Proof. Firstly, if we take the squares of the eigenvalues in (5)-(7), we find that
a/2 2 2 Z 4(n + 1) π 1 0 4(n + 1)π λ2n+2 = + q (t)sin t dt a2 2(n + 1)π a 0 a/2 a/2 Z Z (8) a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 16(n + 1)2π2 0 0 + o n−2,
a/2 2 2 Z (2n + 1) π 1 0 2(2n + 1)π µ2n = − q (t)sin t dt a2 (2n + 1)π a 0 a/2 a/2 Z Z (9) a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 4(2n + 1)2π2 0 0 + o n−2 and
a/2 2 2 Z (2n + 1) π 1 0 2(2n + 1)π µ2n+1 = + q (t)sin t dt a2 (2n + 1)π a 0 a/2 a/2 Z Z (10) a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 4(2n + 1)2π2 0 0 + o n−2.
By (2), we obtain the bounds for Λ2n − λ2n that
µ2n − λ2n ≤ Λ2n − λ2n ≤ µ2n+1 − λ2n. (11) 21 /
By rewriting λ2n+2 from (8) for 2n, we get that
a/2 2 2 Z 4n π 1 0 4nπ λ2n = + q (t)sin t dt a2 2nπ a 0 a/2 a/2 Z Z a 2 0 0 − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 16n2π2 0 0 + o n−2. By using the last equation and (9)-(10), we calculate that
a/2 2 Z (4n + 1)π 1 0 2(2n + 1)π µ2n+1 − λ2n = + q (t)sin t dt a2 (2n + 1)π a 0 a/2 Z 1 0 4nπ − q (t)sin t dt 2nπ a 0 a/2 a/2 " # Z Z a a 2 0 0 − − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 4(2n + 1)2π2 16n2π2 0 0 + o n−2 and
a/2 2 Z (4n + 1)π 1 0 2(2n + 1)π µ2n − λ2n = − q (t)sin t dt a2 (2n + 1)π a 0 a/2 Z 1 0 4nπ − q (t)sin t dt 2nπ a 0 a/2 a/2 " # Z Z a a 2 0 0 − − aq (a) + 2a q(t)q (t)dt − 4 tq(t)q (t)dt 4(2n + 1)2π2 16n2π2 0 0 + o n−2. By using last equations and (11), we prove the theorem.
Example 2.2. Consider an eigenvalue problem is named as anharmonic oscillator
y00 (x) + [λ − q(x)]y(x) = 0, x ∈ [0,π]
1 π 4 1 π 2 where q(x) = 4 x − 2 + 2 x − 2 and extended by periodicity. Since q(x) has mean value zero in the reference theorem, q(x) can taken as follows:
1 π 4 1 π 2 π2 π4 q(x) = x − + x − − − . 4 2 2 2 24 320 In this case, by calculating integral terms in the theorem, we get as n → ∞
1 3 2 3 3 Λ2n − λ2n ≥(4n + 1) − 4π + 16π n + 4π + 16π n + π − 2π 16(2n + 1)4π 1 − 2π3 + 8πn2 − 3π 128n4π " # 1 1 π4 − − 112π4 + 1920π2 + 8960 4(2n + 1)2 16n2 6451200 + o n−2 and
22 /
1 3 2 3 3 Λ2n − λ2n ≤(4n + 1) + 4π + 16π n + 4π + 16π n + π − 2π 16(2n + 1)4π 1 − 2π3 + 8πn2 − 3π 128n4π " # 1 1 π4 − − 112π4 + 1920π2 + 8960 4(2n + 1)2 16n2 6451200 + o n−2.
3. Discussion and Conclusion
We obtain the following conclusion from the main theorem: Conclusion: If q(x) is a constant, the length between Dirichlet and the periodic eigenvalues as n → ∞
2 " # (4n + 1)π 1 1 2 2 Λ2n − λ2n = − − a q (a) a2 4(2n + 1)2π2 16n2π2 + o n−2. This conclusion holds since q0 (x) = 0
Acknowledgment
The author is grateful to Prof. Dr. Haskız CO¸SKUNfor supervisor.
References
[1] Eastham, M. S. P. 1973. The spectral theory of periodic differential equations. Scottish Academic Press, Edinburgh and London, chapter 2 and 3. [2] Ashbaugh, M. and Benguria, R. 1989. Optimal lower bound for the gap between the first two eigenvalues of one- dimensional Schrödinger operators with symmetric single-well potentials and related results. Proc. Amer. Math. Soc., 105(1989), 419-424. [3] Haaser, N. B. and Sullivian, J. A. 1991. Analysis. Van Nostrand Reinhold Co., New York. [4] Huang, M.- J. 1997. The first instability interval for Hill equations with symmetric single-well potentials. Proc. Amer. Math. Soc., 125(1997), no.3, 775-778. [5] Horvath, M. 2002. On the first two eigenvalues of Sturm-Liouville Operators. Proc. Amer. Math. Soc., 131(2002), no.4, 1215-1224. [6] Huang, M.- J. and Tsai, T.- M. 2009. The eigenvalue gap for one-dimensional Schrödinger operators with symmetric potentials. Proc. Roy. Soc. Edinburg, 139(2009), 359-366. [7] Chen, D.- Y. and Huang, M.- J. 2011. Comparison theorems for the eigenvalue gap of Schrödinger operators on the real line. Ann. Henri Poincare, 13(2011), 85-101. [8] Co¸skun,H., Ba¸skaya,E. and Kabata¸s,A. 2019. Instability intervals for Hill’s equation with symmetric single well potential. Ukr. Mat. Zh., 71(2019), no.6, 858-864.
23 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Bodrum Isis Hotel in Mugla, TURKEY ICRAPAM CONFERENCE PROCEEDING
1 On solving the (2+1)-dimensional Zoomeron equation via (G 0 )- expansion method
Hasan Bulut∗1, Gizem Aydin2 1Department of Mathematics, Firat University, Elazig, Turkey 2Department of Mathematics, Firat University, Elazig, Turkey
1 Keywords Abstract: In this article, ( G0 )-expansion method is applied on (2+1)-dimensional Zoomeron equation, Zoomeron equation to seek some new traveling wave solutions. Many operations have 1 ( G0 )-expansion method. been performed with the package program code for the accuracy of these solutions. The gained solutions in this paper are verified by plugging them back to the studied partial differential equation. The exact solutions of the suggested equation are obtained via using 1 the ( G0 )-expansion method. We are confident that the solutions will be useful for scientists to work in this area.
1. Introduction
Nonlinear partial differential equations (PDEs) have an important place in physics and mathematics. Such equations are mathematical models for processes in the aspects of physics, genetics, chemicals and science. Computer programs like Matlab, Maple, and Mathematica are used to gain the general solutions of PDEs. The general solutions of these equations provide researchers with information about the character of physical events. To understand the nature of physical phenomena, researchers have studied many methods, some of these methods can be listed as follows: an inverse spectral 1 transform technique[1], the Laplace perturbation method [2],( G0 )-expansion method [3], direct algebraic method [4], the hyperbolic tangent function expansion method [5], the generalized expansion method [6], sine-Gordon expansion method (SGEM) [7], the Hirota‘s bilinear transformation method [8], the Bernoulli sub-equation function method [9][10], sechp − tanhp functions methods [11], the modified simple equation method [12], generalized Kudryashov method G0 [13],( G )-expansion method [14], Tanh-Coth method [15], a multiplier method [16], and more [17-23].The nonlinear Zoomeron equation in (2+1)-dimension is given by [1] can be written as: u u xy − xy + 2 u2 = 0, (1) u xy u xx xt here u(x,y,t) represent the amplitude of the relative wave mode.
1 2. ( G0 )- Expansion Method
1 In this portion, the general details of the ( G0 )-expansion method are provided. Let us consider a three-variable general form of the nonlinear PDE as O(u,ut ,uux,uuy,uxx,....) = 0, (2) and transform equation (2) with u(x,y,t) = u(ξ) = u,ξ = κ1x + κ2y − κ3t,in which κ1,κ2 and κ3 are non-zero constants. After conversion, we’ll have a nonlinear ordinary differential equation (ODE)
S(u,u0,uu00,u000,.....) = 0. (3)
The solution of (3) shall be taken in the following form
m 1 i u(ξ) = a0 + ∑ ai( 0 ) . (4) i=1 G
Here G = G(ξ) is a function that verifies the following second-order ordinary differential equation
G00 + λG0 + µ = 0, (5)
∗ Corresponding author: hbulut@firat.edu.tr 24 / here a0,a1,a2,...,λ, µ are constants and m is a term of balance. The term balancing is a fixed number obtained between the linear term and the non-linear term in any nonlinear ODE that have highest order. This number is written in place of (4) 00 1 and then the derivatives required for the solution are obtained. These derivatives are taken as G = −λG − µ.( G0 ) is a polynomial and homogeneous equation. Here we gain a system of algebraic equations by equalizing the coefficients of 1 i ( Gi ) ,i = 1,2,3,.... terms to zero. This algebraic equation scheme is solved automatically or with the aid of a computer program depending on the degree of difficulty. These solutions (4) given by the solution function and moving wave transformation is written in place (2 ) given by the solution of the PDE is found [24]. The general solution of the ODE (5) could be reported as below. µξ G = G(ξ) = − − c e−λξ + c , (6) λ 1 2 here c1 and c2 are fixed. Depending on the variable the derivative of the (6) solution function is taken and necessary arrangements are made then 1 1 = , (7) G0 µ λξ − λ + c1e is obtained. When we transform the algebraic expression provided in (7) to a trigonometric function with c1 = A, then it may write as 1 λ = . (8) G0 −µλA(Cosh(ξλ) − Sinh(ξλ))
1 3. Implementation ( G0 )- Expansion Method
1 Now, we use the ( G0 ) -expansion method on (2 + 1) dimensional Zoomeron equation. Consider the (2 + 1) dimensional Zoomeron equation given by Eq. (1). If we apply
u(x,y,t) = U(ξ),ξ = κ1x + κ2y − κ3t, (9) using a wave transformation 00 00 00 00 U U 00 κ κ κ2 − κ3κ + 2κ κ u2 = 0, (10) 1 2 3 u 1 2 u 1 3 we acquire the nonlinear ODE. Now integrating the above equation two times with belong to ξ , we have
2 2 00 3 κ1κ2(κ3 − κ1 )U − κ1κ2U + εU = 0, (11) we get the equation. Here ε is the scalar of integration. Balancing between the U00 and U3 in Equation (11), one have m = 1. In this case, the equation (4) becomes 1 u(ξ) = a + a ( ), (12) 0 1 G0
1 i the solution (12) together with its derivative inserting into Eq. (11) and the necessary arrangements are made, Gi ) ,i = 1,2,3....,m coefficients of the term polynomial terms equal to zero 1 ( )0 : εa − κ κ a3, G0 0 1 3 0 1 ( )1 : εa − κ3λ 2κ a + κ λ 2κ κ2a − 6κ κ a2a , G0 1 1 2 1 1 2 3 1 1 3 0 1 1 ( )2 : −3κ3λ µκ a + 3κ λ µκ κ2a − 6κ κ a a2, G0 1 2 1 1 2 3 1 1 3 0 1 1 ( )3 : −2κ3µ2κ a + 2κ µ2κ κ2a − 2κ κ a3, G0 1 2 1 1 2 3 1 1 3 1 (13) algebraic equation system is obtained. The system of equations (13) given by the computer package program is solved. √ √ 1 √ √ 1 √ 4 4 3 2 i ε λκ2 i 2 εµκ2 2ε+κ1 λ κ2 Case-1 When a0 = − , a1 = − , κ3 = − √ √ , one could gain the following √ 1 1 √ 1 1 κ1λ κ2 4 3 2 4 4 3 2 4 2κ1 (2ε+κ λ κ2) λκ1 (2ε+κ1 λ κ2) complex soliton solution
√ √ 1 √ √ 1 4 4 i ε λκ2 i 2 εµκ2 u(x,y,t) = − 1 − √ 1 , √ 1 1 µ 4 3 2 4 4 3 2 4 2κ1 (2ε + κ λ κ2) λκ1 (2ε + κ1 λ κ2) − λ +C1cosh[λξ] −C1sinh[λξ]
25 / where ξ = κ1x + κ2y − κ3t. √ √ 1 i ε λκ 4 u(x,t) = − 2 √ 1 1 4 3 2 4 2κ1 (2ε + κ λ κ2) √ √ 1 i 2 εµκ 4 − √ 2 √ √ 1 1 3 2 3 2 4 3 2 µ t 2ε+κ1 λ κ2 t 2ε+κ1 λ κ2 λκ (2ε + κ λ κ2) 4 − +C1cosh[λ(κ1x + κ2y − √ √ )] −C1sinh[λ(κ1x + κ2y − √ √ )] 1 1 λ κ1λ κ2 κ1λ κ2 providing ε > 0, κ1 > 0, λ > 0, κ2 > 0. √ √ 1 √ √ 1 √ 4 4 3 2 ε λκ2 2 εµκ2 2ε+κ1 λ κ2 Case-2 When a0 = − , a1 = − , κ3 = − √ √ , we can gain the following real √ 1 1 √ 1 1 κ1λ κ2 4 3 2 4 4 3 2 4 2κ1 (2ε+κ λ κ2) λκ1 (2ε+κ1 λ κ2) soliton solution √ √ 1 √ √ 1 4 4 ε λκ2 2 εµκ2 u(x,y,t) = − 1 − √ 1 , √ 1 1 µ 4 3 2 4 4 3 2 4 2κ1 (2ε + κ λ κ2) λκ1 (2ε + κ1 λ κ2) − λ +C1cosh[λξ] −C1sinh[λξ] where ξ = κ1x + κ2y − κ3t. √ √ 1 ε λκ 4 u(x,y,t) = − 2 √ 1 1 4 3 2 4 2κ1 (2ε + κ λ κ2) √ √ 1 2 εµκ 4 − √ 2 √ . √ 1 1 3 2 3 2 4 3 2 µ t 2ε+κ1 λ κ2 t 2ε+κ1 λ κ2 λκ (2ε + κ λ κ2) 4 − +C1cosh[λ(κ1x + κ2y − √ √ )] −C1sinh[λ(κ1x + κ2y − √ √ )] 1 1 λ κ1λ κ2 κ1λ κ2 √ √ 2 2 2 2 2 2 λ κ2(−κ1 +κ3 ) λ κ2(−κ1 +κ3 ) 1 2 2 2 Case-3 When a0 = √ , a1 = √ , ε = κ1λ κ2(−κ + κ ), we have the following real soliton 2 κ3 λ κ3 2 1 3 solution q q 2 2 2 2 2 2 λ κ2(−κ1 + κ3 ) λ κ2(−κ1 + κ3 ) u(x,y,t) = − √ − √ µ , 2 κ3 λ κ3(− λ +C1cosh[λξ] −C1sinh[λξ]) where ξ = κx + νy − ωt. q 2 2 2 2µ λ κ2(−κ1 + κ3 ) −1 + µ−C λcosh[λ(κ x+κ y−κ t)]−C λsinh[λκ x+κ y−κ t] u(x,y,t) = 1 √ 1 2 3 1 1 2 3 , 2 κ3
2 2 2 providing λ κ2(−κ + κ ) > 0,κ3 > 0. 1 3√ √ √ ε 2 εµ 2 √ √ √ √ √ ε Case-4 When a0 = − ,a1 = − κ λ κ , κ2 = 2 2 2 , 2 κ1 κ3 1 3 κ1λ (−κ1 −κ3 ) we have the following real soliton solution √ √ √ ε 2 εµ u(x,y,t) = −√ √ √ − √ √ µ , 2 κ1 κ3 κ1λ κ3(− λ +C1cosh[λξ] −C1sinh[λξ]) where ξ = κ1x + κ2y − κ3t. √ 2µ ε −1 + 2ε 2ε µ−C1λcosh λ xκ1−tκ3+ +C1sinh λ xκ1−tκ3+ κ λ2(−κ2−κ2) κ λ2(−κ2−κ2) u(x,y,t) = 1√ √1 √3 1 1 3 , 2 κ1 κ3 providing κ1 > 0,κ3 > 0,ε > 0.
4. Conclusions
1 We have implemented the( G0 ) -expansion method successfully to the Zoomeron equation. The gained solutions are novel compared to other results presented in article research [25-27]. The findings of these studies would be useful in describing the physical definitions of certain non-linear models in research. This approach can also be extended to many other nonlinear equations. 26 /
References
[1] F. Calogero and A. Degasperis, Nonlinear evolution equations solvable by the inverse spectral transform.-I, Nuovo Cim. B Ser. 11, 1976. [2] M. Yavuz, N. Ozdemir, and H. M. Baskonus, Solutions of partial differential equations using the fractional operator involving Mittag-Leffler kernel, Eur. Phys. J. Plus, 2018. [3] K. K. Ali, H. Dutta, R. Yilmazer, and S. Noeiaghdam, On the new wave behaviors of the Gilson-Pickering equation, Front. Phys., vol. 8, p. 54, 2020. [4] A. R. Seadawy, Stability analysis for Zakharov-Kuznetsov equation of weakly nonlinear ion-acoustic waves in a plasma, Comput. Math. with Appl., 2014. [5] S. Liu, Z. Fu, S. Liu, and Q. Zhao, Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations, Phys. Lett. Sect. A Gen. At. Solid State Phys., 2001. [6] R. Sabry, W. M. Moslem, F. Haas, S. Ali, and P. K. Shukla, Nonlinear structures: Explosive, soliton, and shock in a quantum electron-positron-ion magnetoplasma, Phys. Plasmas, 2008. [7] G. Yel, H. M. Baskonus, and H. Bulut, Novel archetypes of new coupled Konno-Oono equation by using sine-Gordon expansion method, Opt. Quantum Electron., 2017. [8] R. Hirota and J. Satsuma, Soliton solutions of a coupled Korteweg-de Vries equation, Phys. Lett. A, vol. 85, no. 8-9, pp. 407-408, 1981. [9] H. H. Abdulkareem, H. F. Ismael, E. S. Panakhov, and H. Bulut, Some Novel Solutions of the Coupled Whitham- Broer-Kaup Equations, in International Conference on Computational Mathematics and Engineering Sciences, 2019, pp. 200-208. [10] H. F. Ismael and H. Bulut, On the Solitary Wave Solutions to the (2+ 1)-Dimensional Davey-Stewartson Equations, in International Conference on Computational Mathematics and Engineering Sciences, 2019, pp. 156-165. [11] M. Alquran and K. Al-Khaled, Mathematical methods for a reliable treatment of the (2+1)-dimensional Zoomeron equation, Math. Sci., 2013. [12] K. Khan and M. Ali Akbar, Traveling wave solutions of the (2 + 1)-dimensional Zoomeron equation and the Burgers equations via the MSE method and the Exp-function method, Ain Shams Eng. J., 2014. [13] E. Aksoy, A. C. Cevikel, and A. Bekir, Soliton solutions of (2+1)-dimensional time-fractional Zoomeron equation, Optik (Stuttg)., 2016. [14] R. Abazari, The solitary wave solutions of Zoomeron equation, Appl. Math. Sci, vol. 5, no. 59, pp. 2943-2949, 2011. [15] A. IRSHAD and S. T. MOHYUD-DIN, Solitary wave solutions for Zoomeron equation, Walailak J. Sci. Technol., vol. 10, no. 2, pp. 201-208, 2012. [16] T. Motsepa, C. Khalique, and M. Gandarias, Symmetry Analysis and Conservation Laws of the Zoomeron Equation, Symmetry (Basel)., 2017. [17] H. F. Ismael, H. Bulut, and H. M. Baskonus, Optical soliton solutions to the Fokas-Lenells equation via sine-Gordon G0 expansion method and (m+ ( G )) -expansion method, Pramana - J. Phys., vol. 94, no. 1, Dec. 2020. [18] H. Dutta, H. Gunerhan, K. K. Ali, and R. Yilmazer, Exact soliton solutions to the cubic-quartic nonlinear Schrodinger equation with conformable derivative, Front. Phys., vol. 8, p. 62, 2020. [19] H. F. Ismael, Carreau-Casson fluids flow and heat transfer over stretching plate with internal heat source/sink and radiation, Int. J. Adv. Appl. Sci. J., vol. 6, no. 2, pp. 81-86, 2017. [20] W. G. and H. F. I. and H. B. and H. M. Baskonus, Instability modulation for the (2+1)-dimension paraxial wave equation and its new optical soliton solutions in Kerr media, Phys. Scr., 2019. [21] W. Gao, H. F. Ismael, A. M. Husien, H. Bulut, and H. M. Baskonus, Optical Soliton Solutions of the Cubic-Quartic Nonlinear Schro¨dinger and Resonant Nonlinear Schrodinger Equation with the Parabolic Law, Appl. Sci., vol. 10, no. 1, p. 219, Dec. 2019. [22] C. K. Kuo and B. Ghanbari, Resonant multi-soliton solutions to new (3+1)-dimensional Jimbo-Miwa equations by applying the linear superposition principle, Nonlinear Dyn., 2019. [23] W. Gao, H. F. Ismael, S. A. Mohammed, H. M. Baskonus, and H. Bulut, Complex and real optical soliton properties of the paraxial nonlinear Schrodinger equation in Kerr media with M-fractional, Front. Phys., vol. 7, p. 197, 2019. [24] A. Yokus, An expansion method for finding traveling wave solutions to nonlinear pdes," 2015.
27 /
G0 [25] K. Khan, M. Ali Akbar, M. Abdus Salam, and M. Hamidul Islam, A note on enhanced ( G )-expansion method in nonlinear physics, Ain Shams Eng. J., 2014. [26] O. Guner and A. Bekir, Bright and dark soliton solutions for some nonlinear fractional differential equations, Chinese Phys. B, 2016. [27] A. Bekir, F. Tascan, and O¨. U¨ nsal, Exact solutions of the Zoomeron and Klein-Gordon-Zakharov equations, J. Assoc. Arab Univ. Basic Appl. Sci., 2015.
28 7Th INTERNATIONAL CONFERENCE ON RECENT ADVANCES IN PURE AND APPLIED MATHEMATICS (ICRAPAM 2020) September 25-28 2020 at Bodrum Isis Hotel in Mugla, TURKEY ICRAPAM CONFERANCE PROCEEDING
Some Comparisons of BLUPs under General Linear Random Effects Models
Melek ERI˙S¸ BÜYÜKKAYA∗1, Nesrin GÜLER2 1Karadeniz Technical University, Faculty of Science, Department of Statistics and Computer Sciences, TURKEY 2Sakarya University, Faculty of Political Sciences, Department of Econometrics, TURKEY
Keywords Abstract: In the study, covariance matrices’ comparison of predictors is considered in BLUP, the context of general linear random-effects models. Assume a general linear random- Covariance matrix, effects model with its sub-sample models having correlated random-effects. Our aim is to Inertia, give some comparisons for covariance matrices of the Best Linear Unbiased Predictors Linear random-effects model, (BLUPs) of unknown vectors under sub-sample models by using some well-known rank Rank and inertia formulas for matrices and their Moore-Penrose generalized inverses.
1. Introduction
A general linear random-effects model, defined by
M : y = Xβ + ε, β = Aα + γ, (1)
nx nxp px where y ∈ R 1 is a vector of observable response variables, X ∈ R is a known matrix of arbitrary rank, β ∈ R 1 is pxk kx an unknown random vector, A ∈ R is a known matrix of arbitrary rank, α ∈ R 1 is a vector of fixed but unknown px nx parameters, γ ∈ R 1 is a vector of unobservable random variables, and ε ∈ R 1 is an unobservable vector of random errors. We can divide y, X and ε in M as y1 X1 ε1 X1 ε1 y2 X2 ε2 X2 ε2 : = + = (A + ) + (2) M . . β . . α γ . . . . . . yn Xn εn Xn εn and then we obtain sub-sample models of general linear random-effects model of M as follows
M 1 : y1 = X1β + ε1 = X1 (Aα + γ) + ε1,
M 2 : y2 = X2β + ε2 = X2 (Aα + γ) + ε2, (3) . . M n : yn = Xnβ + εn = Xn (Aα + γ) + εn.
Models M i, i = 1,...,n are transformed version of model M . They are obtained from pre-multiplying the model M by the matrices T1 = In1 0 ··· 0 , T2 = 0 In2 ··· 0 ,..., Tn = 0 0 ··· Inn , (4) respectively. For the models M i, i = 1,...,n, we will assume the following general assumptions on the expectation vector and covariance matrix of random vectors γ γ Σ00 Σ01 ··· Σ0n ε1 ε1 Σ10 Σ11 ··· Σ1n E = 0, cov = := Σ, (5) . . . . .. . . . . . . . εn εn Σn0 Σn1 ··· Σnn
∗ Corresponding author: [email protected] 29 /
(n+p)×(n+p) where Σ ∈ R is a positive semi-definite matrix and all the entries of Σ are known. Under the assumptions in (5), we can obtain the following formulas of the expectations and dispersion of yi
0 E(yi) = Xbiα and cov(yi,yi) = D(yi) = XeiΣXei := Vi, (6) where we use the following expressions for convenience of the representation of above formulas
Xbi = XiA and Xei = [Xi Ti], i = 1,...,n. (7)
We assume that the sub-sample models Mi are consistent, i.e., yi ∈ C [Xbi Vi] holds with probability 1; see, [17, p. 282]. It is known that if the general linear random-effects model M is consistent, then the sub-models are consistent; see, [24].
In order to establish some general results on simultaneous predictions of all unknown parameters under the models, we can consider the following general linear function of fixed effects and random-effects
φ i = Fα + Gγ + HiTiε, i = 1,...,n, (8)
sxk sxp sxni with E(φ i) = Fα for given matrices F ∈ R , G ∈ R and Hi ∈ R . From the assumptions in (5),
0 0 D(φ i) = JiΣJi and cov(φ i,yi) = JiΣXei := Ci (9) is obtained, where Ji = [GHiTi], i = 1,...,n. The vector φ i in (8) is said to be predictable under M i if there exists an sxn Li ∈ R i such that D(Liyi − φ i) = min and E (Liyi − φ i) = 0, i = 1,...,n, (10) holds in the Löwner partial ordering, the linear statistic Li is defined to be the best linear unbiased predictor (BLUP) of φ i under M i and is denoted by
Liyi = BLUPM i (φ i) = BLUPM i (Fα + Gγ + HiTiε). (11)
This is well-known definition of the BLUP of φ i which is originated from [3]. If G = 0 and Hi = 0, (11) reduces the best linear unbiased estimator (BLUE) of Fα under Mi.
Consideration of the models M and Mi separately or simultaneously for making prediction on unknown vectors is meaningful because they have the same random coefficient vectors. Predictors of unknown vectors have different algebraic expressions and different properties under these models. Then it is natural to consider the comparison problem among co- variance matrices of predictors because of having different performances under the models. Covariance matrices of BLUPs are usually used as a comparison criteria to determine optimal predictors among other types of unbiased predictors because of their minimum covariance requirement in the Löwner partial ordering. In this study, we give variety of inequalities and equalities for comparison of covariance matrices of BLUPs under two sub-sample models of M by using rank and inertia formulas. We also present results for comparisons of covariance matrices of BLUEs under considered models. As further reference for comparison of covariance matrix of predictors/estimators, we may mention [5, 6, 22, 23, 25]. For more details on inertias and ranks of symmetric matrices and relations between inertias and Löwner partial ordering of symmetric matrices; see, e.g., [16, 18, 19, 26]. For more related work on BLUP under linear mixed models; see, e.g., [2, 4, 7–14, 20, 21].
mxn 0 + Let R stand for the collection of all m × n real matrices, A , r (A), C (A), and A denote the transpose, the rank, mxn the column space, and the Moore-Penrose generalized inverse of a matrix A ∈ R , respectively, Im denote the identity + ⊥ + + matrix of order m. Furthermore, let PA = AA , EA = A = Im −AA , FA = In −A A stand for the orthogonal projectors. mxm The number of positive and negative eigenvalues of symmetric matrix A ∈ R counted with multiplicities are denoted by i+ (A) and i− (A), called the positive and negative inertia of A, respectively, and also for brief i± (A) denotes the both numbers. It is easy to see that r (A) = i+ (A)+i− (A). A 0, A ≺ 0, A 0 and A 0 mean that A is a symmetric positive definite, negative definite, positive semi-definite, and negative semi-definite matrix, respectively.
2. Preliminaries
In the study of the equality and inequality for covariance matrices of the BLUPs of unknown parameters under general linear random-effects model, it is important to simplify all kinds of matrix expressions involving Moore–Penrose generalized inverses of matrices. One of the powerful tool for simplifying complicated matrix expressions is rank formulas for partitioned matrices. The following rank equalities for partitioned matrices are well known; see [15].
m×n m×k l×n l×k Lemma 2.1. Let A ∈ R , B ∈ R , C ∈ R , and D ∈ R . Then,
r[A B] = r(A) + r(EAB) = r(B) + r(EBA), (12) 30 /
A r = r(A) + r(CF ) = r(C) + r(AF ), (13) C A C AB r = r(B) + r(C) + r(E AF ), (14) C 0 B C AB r = r(A) + r(D − CA+B) if (B) ⊆ (A) and (C0) ⊆ (A0). (15) CD C C C C The following lemmas were given by [18], which will be used to establish rank and inertias of matrices and to characterize equalities and inequalities of BLUPs’ covariances matrices. m×n m×m Lemma 2.2. Let A, B ∈ R , or, A, B ∈ R be symmetric matrices. Then, the following results hold. A = B ⇔ r(A − B) = 0, (16)
A B ⇔ i+(A − B) = m and A ≺ B ⇔ i−(A − B) = m, (17) A < B ⇔ i−(A − B) = 0 and A 4 B ⇔ i+(A − B) = 0. (18) m×m n×n m×n Lemma 2.3. Let A ∈ R and B ∈ R be symmetric matrices and Q ∈ R . Then, + i±(A ) = i±(A), i±(−A) = i∓(A), (19) AQ A −Q −AQ i = i = i , (20) ± Q0 B ± −Q0 B ∓ Q0 −B A 0 0 Q 0 Q i = i (A) + i (B), i = i = r(Q). (21) ± 0 B ± ± + Q0 0 − Q0 0 m×m n×n m×n Lemma 2.4. Let A ∈ R and D ∈ R be symmetric matrices and B ∈ R . Then, AB AB i = r(B) + i (E AE ) and r = 2r(B) + r(E AE ). (22) ± B0 0 ± B B B0 0 B B In particular, AA0 B AA0 B i = r[AB], i = r(B), (23) + B0 0 − B0 0 AA0 B AA0 BC r = r[AB] + r(B), r = r[ABC] + r(B), (24) B0 0 B0 0 0 AB i = i (A) + i (D − B0A+B) if (B) ⊆ (A). (25) ± B0 D ± ± C C For the following lemma, see [20].
Lemma 2.5. Let Mi be as given in (3), i = 1,...,n. The parameter vector φ i in (8) is predictable by yi under Mi if and only if 0 0 C (F ) ⊆ C (Xbi). (26)
In this case, a linear statistic Liyi is the BLUP for φ i if and only if the equation ⊥ ⊥ Li[Xbi ViXbi ] = [FCiXbi ] (27) is satisfied. Then, the general expression of Li and the corresponding BLUPMi (φ i) can be written as ⊥ ⊥ + ⊥ ⊥ BLUPMi (φ i) = Liyi = [FCiXbi ][Xbi ViXbi ] + Ui[Xbi ViXbi ] yi, (28) s×n where Ui ∈ R i is arbitrary. Furthermore, the following results hold. ⊥ ⊥ ⊥ a) r[Xbi ViXbi ] = r[Xbi Vi], C [Xbi ViXbi ] = C [Xbi Vi], and C (Xbi) ∩ C (ViXbi ) = {0},
b) Li is unique if and only if r[Xbi Vi] = ni,
c) BLUPMi (φ i) is unique if and only if yi ∈ C [Xbi Vi] holds with probability 1,
d) BLUPMi (φ i) satisfies + + 0 D[BLUPMi (φ i)] = (ZiWi )Vi(ZiWi ) , (29) + 0 cov(BLUPMi (φ i),φ i) = ZiWi Ci, (30) + + 0 D[φ i − BLUPMi (φ i)] = (ZiWi Xei − Ji)Σ(ZiWi Xei − Ji) , (31) ⊥ ⊥ where Wi = [Xbi ViXbi ], Zi = [FCiXbi ] and Ji = [GHiTi].
It is noteworthy that the requirement in (26) also means Fα is an estimable parametric function under Mi, which is a well-known definition for estimability of a parametric function, see, e.g., [1]. 31 /
3. Main Results
In this section, we give links between the covariance matrix of BLUP of φ i under two sub-sample models of general linear random-effects model M .
Theorem 3.1. Assume that φ k and φ l are predictable under sub-sample models M k and M l, respectively, k,l = 1,...,n, k 6= l, and let the BLUP equation be as given in (28). Also, denote
0 Vk 0 Ck 0 Xbk 0 0 −Vl Cl −Xbl 0 0 0 N = Ck −Cl JkΣJk − JlΣJl −FF . 0 0 0 −Xbl −F 0 0 0 0 Xbk 0 F 0 0 Then, i+ cov φ k − BLUP M k (φ k) − cov φ l − BLUP M l (φ l) = i+ (N) − r Xbk Vk − r Xbl , (32)