Dynamic study of multi-peak solitons and other wave solutions of new coupled KdV and new coupled Zakharov–Kuznetsov systems with their stability

Abstract

In this paper, our aim is to further expand the use of the two-variable $(G^{\prime }/G,1/G)$-expansion approach to a new coupled KdV and Z-K system, which has various significant applications in different fields of applied sciences. The KdV equation, along with shallow-water waves and long internal waves in oceans, basically explains how long, one-dimensional waves propagate in a variety of physical conditions. The study of coastal waves on the basis of the ocean is done using the Zakharov–Kuznetsov (Z-K) equation and this model is utilized to illustrate ion-acoustic wave propagation. By using this method, different forms of analytical solutions of the new coupled KdV (NCKdV) system and the new coupled Z-K (NCZ-K) system, such as solitons, multi-peak solitons, solitary waves, trigonometric, hyperbolic and rational functions and other wave solutions are constructed. The significant features of multi-peak solitons induced by the higher-order effects, including velocity variations, localization or periodicity attenuation and state transitions, are revealed. When the localization disappears then the multi-peak soliton becomes a periodic wave. The constructed solutions are also presented graphically having their applications in engineering, etc. The stability of the solution is examined by utilizing modulation instability. The results obtained show that the proposed technique is universal and efficient. In addition, this technique can also be applied to lots of other new coupled systems arising in other areas of applied sciences.

Keywords:

• New coupled KdV and Z-K systems
• analytical solutions
• two-variable $(G^{\prime }/G1/G)$-expansion technique
• multi-peak solitons
• stability

1. Introduction

Non-linear phenomena play a vital part in applied sciences. Soliton theory heavily depends on the travelling wave solutions of non-linear models in applied sciences [1–4]. The generalized KdV and Zakharov–Kuznetsov (Z-K) equations are essential models for many physical phenomena in partial differential equations (PDEs), including shallow-water waves, non-linear optics, applications in astrophysical and space environments, the interaction between a water wave and a floating ice cover, gravity-capillary waves, shock waves, ion-acoustic waves, hydro-dynamic, stratified waves and many others. The study of explicit analytical solutions to new coupled systems is crucial for illuminating the properties of non-linear phenomena and soliton theory. The travelling wave solutions of these non-linear coupled systems and high-dimensional non-linear equations may enable researchers to thoroughly investigate the aforementioned natural phenomena. Therefore, it is very important to build up non-linear evolution equations and analytical solutions of new coupled systems.

In mathematical physics and the applied sciences, the study of analytical solutions in the explicit form for new coupled systems is crucial for describing the characteristics of non-linear issues and soliton theory. To gain explicit results in the form of solitons, solitary wave and elliptic function solutions, many precise and influential techniques have been developed such as Inverse scattering transform technique [5], ($G^{\prime }/G)$-expansion method [6], the Classical Darboux transformations [7,8], Bäcklund transformation [9–11], expansion scheme, Hirota's bilinear [12,13], Painlevé analysis [14,15], extended tanh scheme [16], auxiliary equation scheme [17], direct algebraic methods [18,19], the mapping and extended mapping methods [20], elliptic function scheme [21], rational expansion method and several others [22,23]. A fantastic progress has been made in the analysis of the solutions, structures, interactions and other properties of solitary waves and solitons, also distinct profound results have been effectively obtained [24,25]. Explicit and numerical solutions were illustrated of KdV and coupled KdV equations by decomposition approaches [26–28]. ‘ The concept of soliton was initiated to show that the solitary waves, which keep their characters unaffected, throughout the proliferation process and after interface, are like particles [29–31]. Because of these characteristics, solitons are believed as standard data bits [32]. Owing to the potential applications of optical solitons in telecommunications and ultrafast signal processing systems, it is the main target of theoretical and experimental analysis [32,33]. As the exact results of the PDEs, the breather plays a central task in semi-conductor quantum wells [34], arrays of micro-mechanical oscillators [35] and junctions of Josephson [36]. Breather specifies that the solutions’ behaviour is periodic in space or time having the property of spatial or temporal localization. Two types of known breathers are the Akhmediev breather of space-periodic and the Kuznetsov-Ma soliton of time-periodic [37–42]. Some other analytical results are the results of rogue waves [?, 42–46]. These rogue wave solutions are confined in together space and time, arise from nowhere and vanish starved of a trace [43,44]. These results may be acquired through the Taylor expansion of the breather results [47–50].

In this paper, the two-variable $(G^{\prime }/G,1/G)$-expansion approach is used to achieve the analytical solutions of the new coupled KdV and Z-K systems. As a result, distinct types of solitons and other wave solutions are obtained. The stability of the solution is examined by utilizing modulation instability (MI). The significant features of multi-peak solitons induced via the higher-order effects, including velocity variations, localization or periodicity attenuation and state transitions, are revealed. When the localization disappears, then the multi-peak soliton becomes a periodic wave. The main advantage of this method is that different forms of analytical solutions of such as solitons, solitary waves, rational solitons, trigonometric function, hyperbolic function and other wave solutions are constructed using this technique. The proposed technique is also universal and efficient. This technique can also be applied to lots of other new coupled systems arising in other areas of applied sciences.

The remaining article's main structure is as follows: Section 2 describes the two-variable $(G^{\prime }/G,1/G)$ technique in detail. In the third Section, the given technique is applied on the NCKdV and NCZ-K systems for producing precise wave results. The analysis of modulation instability is discussed in the fourth section to explain the stability of models. The results discussion and its physical justification are provided in the fifth section. Finally, the sixth section summarizes the entire work.

2. Proposed scheme

In order to arrive at wave results of the aforementioned equations, the $(G^{\prime }/G,1/G)$-expansion technique is described in detail in this section. The following is a full introduction to this strategy from [1,51]:

The aforementioned approach begins with the second-order linear ODE, which is known as follows: (1) $G^{″}(\xi )+\mathrm{\sigma G}(\xi )=\mu ,$(1) by taking (2) $\mathrm{\Psi }=\frac{1}{G},\mathrm{\Phi }=\frac{G^{\prime }}{G},$(2) for the accuracy of the calculations. The variables Φ and Ψ have the derivatives as (3) ${\mathrm{\Phi }}^{\prime }=-{\mathrm{\Phi }}^2+\mu \mathrm{\Psi }-\sigma ,{\mathrm{\Psi }}^{\prime }=-\mathrm{\Phi \Psi }.$(3) There are three types of cases in the general results of Equation (1(1) $G^{″}(\xi )+\mathrm{\sigma G}(\xi )=\mu ,$(1) )

Case-I: $\sigma <0$, the solution is (4) $G(\xi )=A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)+\frac{\mu }{\sigma },$(4) with (5) $\begin{array}{rl}{\mathrm{\Psi }}^2& =\frac{-\sigma }{{\sigma }^2\eta +{\mu }^2}(\sigma -2\mathrm{\Psi \mu }+{\mathrm{\Phi }}^2),\mathrm{where}\\ \eta & =A_1^2-A_2^2.\end{array}$(5) Case-II: $\sigma >0$, the solution is as (6) $G(\xi )=A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)+\frac{\mu }{\sigma },$(6) with (7) ${\mathrm{\Psi }}^2=\frac{\sigma ({\mathrm{\Phi }}^2-2\mu \mathrm{\Psi }+\sigma )}{{\sigma }^2\sigma -{\mu }^2},\mathrm{where}\eta =A_1^2+A_2^2$(7) Case-III: $\sigma =0$, the rational solution is as (8) $G(\xi )=A_1\xi +A_2+\frac{\mu }{2}{\xi }^2,$(8) with (9) ${\mathrm{\Psi }}^2=\frac{{\mathrm{\Phi }}^2-2\mu \mathrm{\Psi }}{A_1^2-2\mu A_2}$(9) In all cases, $A_1$ and $A_2$ are constants.

The following procedures must be taken in order to use this method to obtain precise solutions to the non-linear evolution equation. (10) $P(u,\frac{\mathrm{\partial }u}{\mathrm{\partial }t},\frac{\mathrm{\partial }u}{\mathrm{\partial }x},\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{t}^2},\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2},\dots )=0.$(10) Step 1: Wave transformation $\xi =x-\mathrm{tc}$ is used to convert Equation (10(10) $P(u,\frac{\mathrm{\partial }u}{\mathrm{\partial }t},\frac{\mathrm{\partial }u}{\mathrm{\partial }x},\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{t}^2},\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2},\dots )=0.$(10) ) to ODE (ordinary differential equation); it includes $u(x,t)=u(\xi )$ as (11) $P(u,-c\frac{\mathrm{d}u}{\mathrm{d}\xi },\frac{\mathrm{d}u}{\mathrm{d}\xi },c^2\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},-c\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\dots )=0,$(11) where $u^{(n)}$ is the nth order derivative with respect to ξ.

Step 2: The result could be in the form of Ψ and Φ (12) $u(\xi )=\sum _{j=0}^N{a}_j{\mathrm{\Phi }}^j+\sum _{j=1}^N{b}_j{\mathrm{\Phi }}^{j-1}\mathrm{\Psi },$(12) here G adheres to (1(1) $G^{″}(\xi )+\mathrm{\sigma G}(\xi )=\mu ,$(1) ). constant coefficients $a_j,b_j,\beta ,c$ and α are to be determined. N can be determined by using the balancing principle in (11(11) $P(u,-c\frac{\mathrm{d}u}{\mathrm{d}\xi },\frac{\mathrm{d}u}{\mathrm{d}\xi },c^2\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},-c\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\dots )=0,$(11) ).

Step 3: By considering case 1 as an example, utilizing (12(12) $u(\xi )=\sum _{j=0}^N{a}_j{\mathrm{\Phi }}^j+\sum _{j=1}^N{b}_j{\mathrm{\Phi }}^{j-1}\mathrm{\Psi },$(12) ) into (11(11) $P(u,-c\frac{\mathrm{d}u}{\mathrm{d}\xi },\frac{\mathrm{d}u}{\mathrm{d}\xi },c^2\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},-c\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\dots )=0,$(11) ) while taking into account (3(3) ${\mathrm{\Phi }}^{\prime }=-{\mathrm{\Phi }}^2+\mu \mathrm{\Psi }-\sigma ,{\mathrm{\Psi }}^{\prime }=-\mathrm{\Phi \Psi }.$(3) ) and (5(5) $\begin{array}{rl}{\mathrm{\Psi }}^2& =\frac{-\sigma }{{\sigma }^2\eta +{\mu }^2}(\sigma -2\mathrm{\Psi \mu }+{\mathrm{\Phi }}^2),\mathrm{where}\\ \eta & =A_1^2-A_2^2.\end{array}$(5) ), a polynomial equation in Ψ and Φ is attained and produces a system of algebraic equations.

Step 4: A software package program is used to solve the system. Wave solutions in (11(11) $P(u,-c\frac{\mathrm{d}u}{\mathrm{d}\xi },\frac{\mathrm{d}u}{\mathrm{d}\xi },c^2\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},-c\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\frac{{\mathrm{d}}^{2}u}{\mathrm{d}{\xi }^2},\dots )=0,$(11) ) are constructed as three different types of functions by utilizing resultant values $a_j,b_j,c,\beta ,\alpha ,A_2$ and $A_1$.

Step 5: The solution process is completed by producing results in (10(10) $P(u,\frac{\mathrm{\partial }u}{\mathrm{\partial }t},\frac{\mathrm{\partial }u}{\mathrm{\partial }x},\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{t}^2},\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2},\dots )=0.$(10) ) using $\xi =x-\mathrm{ct}$ (wave transformation) conversely.

3. Applications

3.1. The new coupled KdV system

Qin [40] used a finite-dimensional integrable system to build a new hierarchy of non-linear evolution equations. A novel coupled KdV system [52] is an interesting equation in this hierarchy is as (13) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}u}{\mathrm{\partial }{x}^3}+\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}+\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}v}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}w}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x},\end{array}$(13) where $\lambda ,\alpha ,\beta$ and γ are arbitrary constants. The author in [17] more extended the NCKdV system into a new coupled KP system and deliberated both systems using Hirota's bilinear approach. Consider the travelling wave solution as (14) $\begin{array}{rl}& U(\xi )=u(x,t),V(\xi )=v(x,t),W(\xi )=w(x,t),\\ & \mathrm{where},\xi =\mathrm{\delta x}+\mathrm{\rho t},\end{array}$(14) where δ and ρ are wave number and frequency of results, respectively. Using (14(14) $\begin{array}{rl}& U(\xi )=u(x,t),V(\xi )=v(x,t),W(\xi )=w(x,t),\\ & \mathrm{where},\xi =\mathrm{\delta x}+\mathrm{\rho t},\end{array}$(14) ) in (13(13) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}u}{\mathrm{\partial }{x}^3}+\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}+\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}v}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}w}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x},\end{array}$(13) ) and the ODEs obtained as (15) $\begin{array}{rl}\mathrm{\rho U}& =\beta {\delta }^3{U}^{″}+\mathrm{\alpha \delta }(\mathrm{UV})+\mathrm{\gamma \delta }(\mathrm{VW})+{ϵ}_1,\\ \mathrm{\rho V}& =\beta {\delta }^3{V}^{″}+\mathrm{\lambda \delta }(\mathrm{UW})+{ϵ}_2,\\ \mathrm{\rho W}& =\beta {\delta }^3{W}^{″}+\mathrm{\lambda \delta }(\mathrm{UV})+{ϵ}_3.\end{array}$(15) Using the balancing principle on (15(15) $\begin{array}{rl}\mathrm{\rho U}& =\beta {\delta }^3{U}^{″}+\mathrm{\alpha \delta }(\mathrm{UV})+\mathrm{\gamma \delta }(\mathrm{VW})+{ϵ}_1,\\ \mathrm{\rho V}& =\beta {\delta }^3{V}^{″}+\mathrm{\lambda \delta }(\mathrm{UW})+{ϵ}_2,\\ \mathrm{\rho W}& =\beta {\delta }^3{W}^{″}+\mathrm{\lambda \delta }(\mathrm{UV})+{ϵ}_3.\end{array}$(15) ) and considering the solutions as (16) $\begin{array}{rl}U(\xi )& =a_0+a_1\mathrm{\Phi }(\xi )+a_2\mathrm{\Phi }(\xi {)}^2+b_1\mathrm{\Psi }(\xi )\\ & +b_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ V(\xi )& =c_0+c_1\mathrm{\Phi }(\xi )+c_2\mathrm{\Phi }(\xi {)}^2+d_1\mathrm{\Psi }(\xi )\\ & +d_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ W(\xi )& =e_0+e_1\mathrm{\Phi }(\xi )+e_2\mathrm{\Phi }(\xi {)}^2+f_1\mathrm{\Psi }(\xi )\\ & +f_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ).\end{array}$(16) Utilizing Equations (16(16) $\begin{array}{rl}U(\xi )& =a_0+a_1\mathrm{\Phi }(\xi )+a_2\mathrm{\Phi }(\xi {)}^2+b_1\mathrm{\Psi }(\xi )\\ & +b_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ V(\xi )& =c_0+c_1\mathrm{\Phi }(\xi )+c_2\mathrm{\Phi }(\xi {)}^2+d_1\mathrm{\Psi }(\xi )\\ & +d_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ W(\xi )& =e_0+e_1\mathrm{\Phi }(\xi )+e_2\mathrm{\Phi }(\xi {)}^2+f_1\mathrm{\Psi }(\xi )\\ & +f_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ).\end{array}$(16) ) with (3(3) ${\mathrm{\Phi }}^{\prime }=-{\mathrm{\Phi }}^2+\mu \mathrm{\Psi }-\sigma ,{\mathrm{\Psi }}^{\prime }=-\mathrm{\Phi \Psi }.$(3) ) into (15(15) $\begin{array}{rl}\mathrm{\rho U}& =\beta {\delta }^3{U}^{″}+\mathrm{\alpha \delta }(\mathrm{UV})+\mathrm{\gamma \delta }(\mathrm{VW})+{ϵ}_1,\\ \mathrm{\rho V}& =\beta {\delta }^3{V}^{″}+\mathrm{\lambda \delta }(\mathrm{UW})+{ϵ}_2,\\ \mathrm{\rho W}& =\beta {\delta }^3{W}^{″}+\mathrm{\lambda \delta }(\mathrm{UV})+{ϵ}_3.\end{array}$(15) ), we obtained system of equations in $\rho ,\mu ,{ϵ}_1,{ϵ}_2,{ϵ}_3,a_0,a_1,b_1,$ $c_1,d_1,e_1,f_1,a_2,b_2,c_2,d_2,e_2$ and $f_2$ by equating the coefficients of ${\mathrm{\Psi }}^i{\mathrm{\Phi }}^j$ to zero. After solving the system, the following families of results are produced:

Family 1: (17) $\begin{array}{rl}& a_1=0,a_2=-\frac{6\beta {\delta }^2}{\lambda },b_1=0,b_2=0,\\ & c_0=\frac{a_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-{\alpha }^2)}{2\gamma },c_1=0,\\ & c_2=\frac{3(\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},d_1=0,d_2=0,\\ & e_0=\frac{a_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{2\gamma },\\ & e_1=0,e_2=\frac{3(\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},f_1=0,\\ & f_2=0,\rho =2\delta (a_0\lambda +4\beta {\delta }^2\sigma ),\mu =0,\\ & {ϵ}_1=\frac{\delta (a_0\lambda +2\beta {\delta }^2\sigma )(a_0\lambda +6\beta {\delta }^2\sigma )}{\lambda },\\ & {ϵ}_2=\frac{\left(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha \right)\left(a_0\lambda +2\beta {\delta }^2\sigma \right)\left(a_0\lambda +6\beta {\delta }^2\sigma \right)}{2\mathrm{\gamma \lambda }},\\ & {ϵ}_3=\frac{\left(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha \right)\left(a_0\lambda +2\beta {\delta }^2\sigma \right)\left(a_0\lambda +6\beta {\delta }^2\sigma \right)}{2\mathrm{\gamma \lambda }}.\end{array}$(17) Family 2: (18) $\begin{array}{rl}& a_0=-\frac{6\beta {\delta }^2\sigma }{\lambda },a_1=0,a_2=-\frac{6\beta {\delta }^2}{\lambda },b_1=0,\\ & b_2=0,c_0=\frac{3\beta {\delta }^2\sigma (\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},\\ & c_1=0,c_2=\frac{3\beta {\delta }^2(\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},d_1=0,\\ & d_2=0,e_0=\frac{3\beta {\delta }^2\sigma (\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},\\ & e_1=0,e_2=\frac{3\beta {\delta }^2(\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},f_1=0,\\ & f_2=0,\rho =-4\beta {\delta }^3\sigma ,\\ & \mu =0,{ϵ}_1=0,{ϵ}_2=0,{ϵ}_3=0.\end{array}$(18) Family 3: $\begin{array}{rl}& a_0=\frac{c_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{2\lambda },a_1=0,a_2=-\frac{6\beta {\delta }^2}{\lambda },\\ & b_1=0,b_2=0,c_1=0,\\ & c_2=\frac{3\beta {\delta }^2(\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},d_1=0,d_2=0,\\ & e_0=c_0,e_1=0,\\ & e_2=\frac{3\beta {\delta }^2(\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})}{\mathrm{\gamma \lambda }},f_1=0,f_2=0,\\ & \rho =8\beta {\delta }^3\sigma +c_0\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha c_0\delta ,\mu =0,\\ & {ϵ}_1=\frac{\begin{array}{c}8\beta {\delta }^3\sigma \left(3\beta {\delta }^2\sigma +c_0\left(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha \right)\right)\\ +c_0^2\delta (\alpha \sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+({\alpha }^2+2\mathrm{\gamma \lambda }))\end{array}}{2\lambda },\\ & {ϵ}_2=\frac{\begin{array}{c}12{\beta }^2{\delta }^5{\sigma }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ +\gamma c_0\mathrm{\lambda \delta }(16\beta {\delta }^2\sigma +c_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha ))\end{array}}{2\mathrm{\gamma \lambda }},\\ & {ϵ}_3=\frac{\begin{array}{c}12{\beta }^2{\delta }^5{\sigma }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ +\gamma c_0\mathrm{\lambda \delta }(16\beta {\delta }^2\sigma +c_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha ))\end{array}}{2\mathrm{\gamma \lambda }}.\end{array}$Family 4: $\begin{array}{rl}& a_0=\frac{\begin{array}{c}8\mathrm{\beta \gamma }{\delta }^2\mathrm{\lambda \sigma }-\sqrt{2}\beta {\delta }^2\sigma \sqrt{{\alpha }^2+2\mathrm{\gamma \lambda }-\alpha \sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}}\\ (\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }})\end{array}}{2\gamma {\lambda }^2},\\ & a_1=0,a_2=\frac{6\beta {\delta }^2}{\lambda },b_1=0,b_2=0,\\ & c_0=\frac{\begin{array}{c}\beta {\delta }^2\sigma [\sqrt{2}\sqrt{{\alpha }^2+2\mathrm{\gamma \lambda }-\alpha \sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}}\\ -2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )]\end{array}}{\mathrm{\gamma \lambda }},\\ & c_1=0,c_2=\frac{3\beta {\delta }^2[\alpha -\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}]}{\mathrm{\gamma \lambda }},d_1=0,d_2=0,\\ & e_0=\frac{\begin{array}{c}\beta {\delta }^2\mathrm{\gamma \lambda \sigma }[2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ -\sqrt{2}\sqrt{{\alpha }^2+2\mathrm{\gamma \lambda }-\alpha \sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}}]\end{array}}{{\gamma }^2{\lambda }^2},\\ & e_1=0,e_2=\frac{3\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }},f_1=0,f_2=0,\\ & \rho =\frac{\sqrt{2}\mathrm{\delta \sigma }(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )\sqrt{{\alpha }^2+2\mathrm{\gamma \lambda }-\alpha \sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}}}{\mathrm{\gamma \lambda }},\\ & \mu =0,{ϵ}_1=0,{ϵ}_2=0,{ϵ}_3=0.\end{array}$Case-I: $\sigma <0$, from the first family, the following soliton solution set is obtained in the hyperbolic function form as (19) $\begin{array}{rl}{u}_1(\xi )& =a_0+\frac{\begin{array}{c}6\beta {\delta }^2\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)\\ +A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2\end{array}}{\begin{array}{c}\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)\\ +A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2\end{array}},\\ v_1(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}+\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})\end{array}}{2\gamma },\\ w_1(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}+\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})\end{array}}{2\gamma }.\end{array}$(19) from the second family, the following soliton solution set is also obtained in the hyperbolic function form as (20) $\begin{array}{rl}{u}_2(\xi )& =\frac{6\beta {\delta }^2\sigma (A_1^2-A_2^2)}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ v_2(\xi )& =\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2-A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ w_2(\xi )& =\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2-A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}.\end{array}$(20)

Case-II: $\sigma >0$, from the first and second families, the following soliton solution sets are obtained in the trigonometric function form as (21) $\begin{array}{rl}{u}_3(\xi )& =a_0-\frac{\begin{array}{c}6\beta {\delta }^2\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)\\ -A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2\end{array}}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ v_3(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}-\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2})\end{array}}{2\gamma },\\ w_3(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}-\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2})\end{array}}{2\gamma }.\end{array}$(21) (22) $\begin{array}{rl}{u}_4(\xi )& =-\frac{6\beta {\delta }^2\sigma (A_1^2+A_2^2)}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ v_4(\xi )& =-\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2+A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ w_4(\xi )& =-\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2+A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }{(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right))}^2}.\end{array}$(22) Case-III: $\sigma =0$, from the first and second families, thefollowing soliton solution sets are obtained in the rational function form as (23) $\begin{array}{rl}{u}_5(\xi )& =a_0-\frac{6\beta {\delta }^2{A}_1^2}{\lambda (A_1\xi +A_2){}^2},\\ v_5(\xi )& =\frac{\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )(\frac{a_0}{\beta {\delta }^2}-\frac{6A_1^2}{\lambda (A_1\xi +A_2){}^2})}{2\gamma },\\ w_5(\xi )& =\frac{{\beta }^2{\delta }^4(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )(\frac{a_0}{\beta {\delta }^2}-\frac{6A_1^2}{\lambda (A_1\xi +A_2){}^2})}{2\gamma }.\end{array}$(23) (24) $\begin{array}{rl}{u}_6(\xi )& =\frac{6\beta {\delta }^2(-\frac{A_1^2}{(A_1\xi +A_2){}^2}-\sigma )}{\lambda },\\ v_6(\xi )& =-\frac{\begin{array}{c}3\beta {\delta }^2(A_1^2({\xi }^2\sigma +1)+2A_2{A}_1\mathrm{\xi \sigma }+A_2^2\sigma )\\ (\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\end{array}}{\mathrm{\gamma \lambda }(A_1\xi +A_2){}^2},\\ w_6(\xi )& =-\frac{\begin{array}{c}3\beta {\delta }^2(A_1^2({\xi }^2\sigma +1)+2A_2{A}_1\mathrm{\xi \sigma }+A_2^2\sigma )\\ (\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\end{array}}{\mathrm{\gamma \lambda }(A_1\xi +A_2){}^2}.\end{array}$(24) The more generalized resemble solitons of Equation (13(13) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}u}{\mathrm{\partial }{x}^3}+\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}+\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}v}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}w}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x},\end{array}$(13) ) other families can also be constructed.

3.2. The new coupled Z-K system

The coupled KdV system (13(13) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}u}{\mathrm{\partial }{x}^3}+\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}+\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}v}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}w}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x},\end{array}$(13) ) can extend to the NCZ-K system [52] in the following form (25) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}-\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{wu})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{y}^2})& =0,\end{array}$(25) where $\lambda ,\alpha ,\beta$ and γ are real numbers. The authors in [52] studied this system using the modified extended direct algebraic scheme. By considering the wave transformation as (26) $\begin{array}{rl}U(\xi )& =u(x,t),V(\xi )=v(x,t),W(\xi )=w(x,t),\\ \xi & =\mathrm{\delta x}+\mathrm{\rho t}+\mathrm{\nu y}.\end{array}$(26) where $\delta ,\rho$ and ν are the wave numbers and frequency, respectively. Now Equation (25(25) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}-\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{wu})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{y}^2})& =0,\end{array}$(25) ) is transformed into ODE by using Equation (25(25) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}-\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{wu})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{y}^2})& =0,\end{array}$(25) ) as (27) $\begin{array}{rl}& \mathrm{\rho U}=\mathrm{\alpha \delta UV}+\mathrm{\gamma \delta }(\mathrm{VW})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)U^{″}+{ϵ}_1,\\ & \mathrm{\rho V}=\mathrm{\lambda \delta }(\mathrm{WU})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)V^{″}+{ϵ}_2,\\ & \mathrm{\rho W}=\mathrm{\lambda \delta }(\mathrm{UV})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)W^{″}+{ϵ}_3.\end{array}$(27) On Equation (27(27) $\begin{array}{rl}& \mathrm{\rho U}=\mathrm{\alpha \delta UV}+\mathrm{\gamma \delta }(\mathrm{VW})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)U^{″}+{ϵ}_1,\\ & \mathrm{\rho V}=\mathrm{\lambda \delta }(\mathrm{WU})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)V^{″}+{ϵ}_2,\\ & \mathrm{\rho W}=\mathrm{\lambda \delta }(\mathrm{UV})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)W^{″}+{ϵ}_3.\end{array}$(27) ) using the balancing principle, and presuming solutions as (28) $\begin{array}{rl}U(\xi )& =a_0+a_1\mathrm{\Phi }(\xi )+a_2\mathrm{\Phi }(\xi {)}^2+b_1\mathrm{\Psi }(\xi )\\ & +b_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ V(\xi )& =c_0+c_1\mathrm{\Phi }(\xi )+c_2\mathrm{\Phi }(\xi {)}^2+d_1\mathrm{\Psi }(\xi )\\ & +d_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ W(\xi )& =e_0+e_1\mathrm{\Phi }(\xi )+e_2\mathrm{\Phi }(\xi {)}^2+f_1\mathrm{\Psi }(\xi )\\ & +f_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ).\end{array}$(28) Utilizing Equations (28(28) $\begin{array}{rl}U(\xi )& =a_0+a_1\mathrm{\Phi }(\xi )+a_2\mathrm{\Phi }(\xi {)}^2+b_1\mathrm{\Psi }(\xi )\\ & +b_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ V(\xi )& =c_0+c_1\mathrm{\Phi }(\xi )+c_2\mathrm{\Phi }(\xi {)}^2+d_1\mathrm{\Psi }(\xi )\\ & +d_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ),\\ W(\xi )& =e_0+e_1\mathrm{\Phi }(\xi )+e_2\mathrm{\Phi }(\xi {)}^2+f_1\mathrm{\Psi }(\xi )\\ & +f_2\mathrm{\Psi }(\xi )\mathrm{\Phi }(\xi ).\end{array}$(28) ) with (3(3) ${\mathrm{\Phi }}^{\prime }=-{\mathrm{\Phi }}^2+\mu \mathrm{\Psi }-\sigma ,{\mathrm{\Psi }}^{\prime }=-\mathrm{\Phi \Psi }.$(3) ) into (27(27) $\begin{array}{rl}& \mathrm{\rho U}=\mathrm{\alpha \delta UV}+\mathrm{\gamma \delta }(\mathrm{VW})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)U^{″}+{ϵ}_1,\\ & \mathrm{\rho V}=\mathrm{\lambda \delta }(\mathrm{WU})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)V^{″}+{ϵ}_2,\\ & \mathrm{\rho W}=\mathrm{\lambda \delta }(\mathrm{UV})+\mathrm{\beta \delta }({\delta }^2-{\nu }^2)W^{″}+{ϵ}_3.\end{array}$(27) ), and the system of algebraic equations obtained in $\rho ,\mu ,{ϵ}_1,{ϵ}_2,$ ${ϵ}_3,a_0,a_1,b_1,$ $c_1,d_1,e_1,f_1,a_2,b_2,c_2,d_2,e_2$ and $f_2$ byequating the coefficients of ${\mathrm{\Psi }}^i{\mathrm{\Phi }}^j$ to zero, the following families of solutions by solving these equations are obtained as follows:

Family 1: $\begin{array}{rl}& a_1=0,a_2=\frac{6\beta ({\nu }^2-{\delta }^2)}{\lambda },b_1=0,b_2=0,\\ & c_0=-\frac{a_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{2\gamma },c_1=0,\\ & c_2=\frac{3\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},d_1=0,d_2=0,\\ & e_0=-\frac{a_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{2\gamma },e_1=0,\\ & e_2=\frac{3\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},f_1=0,f_2=0,\\ & \rho =2\delta (a_0\lambda +4\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)),\mu =0,\\ & {ϵ}_1=-\frac{\delta (a_0\lambda +2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2))(a_0\lambda +6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2))}{\lambda },\\ & {ϵ}_2=\frac{\begin{array}{c}\delta (\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )(a_0\lambda +2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2))\\ (a_0\lambda +6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2))\end{array}}{2\mathrm{\gamma \lambda }},\\ & {ϵ}_3=\frac{\begin{array}{c}\delta (\sqrt{({\alpha }^2+4\mathrm{\gamma \lambda })}+\alpha )(a_0\lambda +2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2))\\ (a_0\lambda +6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2))\end{array}}{2\mathrm{\gamma \lambda }}.\end{array}$Family 2: $\begin{array}{rl}& a_0=\frac{2\mathrm{\beta \sigma }({\nu }^2-{\delta }^2)}{\lambda },a_1=0,a_2=\frac{6\beta ({\nu }^2-{\delta }^2)}{\lambda },\\ & b_1=0,b_2=0,\\ & c_0=\frac{\mathrm{\sigma \beta }({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},c_1=0,\\ & c_2=\frac{3\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},d_1=0,d_2=0,\\ & e_0=\frac{\mathrm{\sigma \beta }({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},e_1=0,\\ & e_2=\frac{3\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},f_1=0,f_2=0,\\ & \rho =4\mathrm{\beta \delta \sigma }({\delta }^2-{\nu }^2),\mu =0,{ϵ}_1=0,\\ & {ϵ}_2=0,{ϵ}_3=0.\end{array}$Family 3: $\begin{array}{rl}& a_1=0,a_2=\frac{6\beta ({\delta }^2-{\nu }^2)}{\lambda },b_1=0,b_2=0,\\ & c_0=\frac{a_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{2\gamma },c_1=0,\\ & c_2=-\frac{3\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},d_1=0,\\ & d_2=0,e_0=\frac{a_0(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{2\gamma },e_1=0,\\ & e_2=\frac{3\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},\\ & f_1=0,f_2=0,\rho =8\mathrm{\beta \delta \sigma }({\delta }^2-{\nu }^2)-2a_0\mathrm{\delta \lambda },\mu =0,\\ & {ϵ}_1=\frac{\delta (2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)-a_0\lambda )(6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)-a_0\lambda )}{\lambda },\\ & {ϵ}_2=-\frac{\begin{array}{c}\delta (\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )(2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)-a_0\lambda )\\ (6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)-a_0\lambda )\end{array}}{2\mathrm{\gamma \lambda }},\\ & {ϵ}_3=\frac{\begin{array}{c}\delta (\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )(2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)-a_0\lambda )\\ (6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)-a_0\lambda )\end{array}}{2\mathrm{\gamma \lambda }}.\end{array}$Family 4: $\begin{array}{rl}& a_0=\frac{2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)}{\lambda },a_1\to 0,a_2\to \frac{6\beta ({\delta }^2-{\nu }^2)}{\lambda },\\ & b_1=0,b_2=0,\\ & c_0=-\frac{\mathrm{\sigma \beta }({\nu }^2-{\delta }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},c_1=0,\\ & c_2=-\frac{3\beta ({\nu }^2-{\delta }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},\\ & d_1=0,d_2=0,\\ & e_0=\frac{\mathrm{\sigma \beta }({\nu }^2-{\delta }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},e_1=0,\\ & e_2=\frac{3\beta ({\nu }^2-{\delta }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )}{\mathrm{\gamma \lambda }},f_1=0,f_2=0,\\ & \rho =4\mathrm{\beta \delta \sigma }({\delta }^2-{\nu }^2),\mu =0,{ϵ}_1=0,\\ & {ϵ}_2=0,{ϵ}_3=0.\end{array}$Case-1: $\sigma <0$, from the first family, the following solution set is obtained in the hyperbolic function form as

(29) $\begin{array}{rl}{u}_1(\xi )& =a_0+\frac{6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)(A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ v_1(\xi )& =\frac{\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )(-\frac{a_0}{\beta {\delta }^2-\beta {\nu }^2}-\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})}{2\gamma },\\ w_1(\xi )& =\frac{\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )(-\frac{a_0}{\beta {\delta }^2-\beta {\nu }^2}-\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})}{2\gamma }.\end{array}$(29) From the second family, the following solution set is obtained in the hyperbolic function form as

(30) $\begin{array}{rl}{u}_2(\xi )& =\frac{2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)(2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ v_2(\xi )& =-\frac{\sigma (2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ w_2(\xi )& =-\frac{\sigma (2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}.\end{array}$(30) Case-2: $\sigma >0$, from the first and second families, the following solution sets are obtained in the trigonometric function form as

(31) $\begin{array}{r}\begin{array}{rl}{u}_3(\xi )& =a_0-\frac{6\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)(A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ v_3(\xi )& =\frac{\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )(\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}-\frac{a_0}{\beta {\delta }^2-\beta {\nu }^2})}{2\gamma },\\ w_3(\xi )& =\frac{\beta ({\delta }^2-{\nu }^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}+\alpha )(\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}-\frac{a_0}{\beta {\delta }^2-\beta {\nu }^2})}{2\gamma }.\end{array}\end{array}$(31) (32) $\begin{array}{r}\begin{array}{rl}{u}_4(\xi )& =\frac{2\beta ({\nu }^2-{\delta }^2)(\frac{3\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}+\sigma )}{\lambda },\\ v_4(\xi )& =\frac{\sigma (-2A_2{A}_1\sin \left(2\xi \sqrt{\sigma }\right)+A_1^2(\cos \left(2\xi \sqrt{\sigma }\right)+2)-A_2^2(\cos \left(2\xi \sqrt{\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ w_4(\xi )& =\frac{\sigma (-2A_2{A}_1\sin \left(2\xi \sqrt{\sigma }\right)+A_1^2(\cos \left(2\xi \sqrt{\sigma }\right)+2)-A_2^2(\cos \left(2\xi \sqrt{\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}.\end{array}\end{array}$(32) Case-III: $\sigma =0$, from the first and second families, the following solution sets are obtained in the rational function form as

(33) $\begin{array}{r}\begin{array}{rl}{u}_5(\xi )& =a_0+\frac{6A_1^2\beta ({\nu }^2-{\delta }^2)}{\lambda (A_1\xi +A_2){}^2},\\ v_5(\xi )& =\frac{\frac{6A_1^2(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\lambda (A_1\xi +A_2){}^2}-\frac{a_0(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\beta ({\delta }^2-{\nu }^2)}}{2\gamma },\\ w_5(\xi )& =\frac{\frac{6A_1^2(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\lambda (A_1\xi +A_2){}^2}-\frac{a_0(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\beta ({\delta }^2-{\nu }^2)}}{2\gamma }.\end{array}\end{array}$(33) (34) $\begin{array}{r}\begin{array}{rl}{u}_6(\xi )& =\frac{2\beta ({\nu }^2-{\delta }^2)(\frac{3A_1^2}{(A_1\xi +A_2){}^2}+\sigma )}{\lambda },\\ v_6(\xi )& =\frac{(A_1^2({\xi }^2\sigma +3)+2A_2{A}_1\mathrm{\xi \sigma }+A_2^2\sigma )(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\xi +A_2){}^2},\\ w_6(\xi )& =\frac{(A_1^2({\xi }^2\sigma +3)+2A_2{A}_1\mathrm{\xi \sigma }+A_2^2\sigma )(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\xi +A_2){}^2}.\end{array}\end{array}$(34) The more generalized resemble solitons of Equation (25(25) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}-\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{wu})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{y}^2})& =0,\end{array}$(25) ) other families can also be constructed.

4. Modulation stability analysis

Various non-linear PDEs exhibit scattering and non-linear consequences as a result of uncertainty in the steady-state modulation. Utilizing linear stability analysis, the modulation instability of NCKdV and NCZ-K systems is investigated [36,38].

4.1. The new coupled KdV system

The following is the form of steady-state solution (SSS) of the NCKdV system (35) $\begin{array}{rl}& \sqrt{P}+{\mathrm{\Phi }}_1(x,t)e^{\mathrm{P\delta ϵt}}=u(x,t),\\ & \sqrt{P}+{\mathrm{\Phi }}_2(x,t)e^{\mathrm{P\delta ϵt}}=v(x,t),\\ & \sqrt{P}+{\mathrm{\Phi }}_3(x,t)e^{\mathrm{P\delta ϵt}}=w(x,t).\end{array}$(35) Here, normalized optical power P. Evolution of the perturbation $\mathrm{\Phi }(x,t)$ is verified by using the linear stability (LS) analysis. We can linearize by putting Equations (35(35) $\begin{array}{rl}& \sqrt{P}+{\mathrm{\Phi }}_1(x,t)e^{\mathrm{P\delta ϵt}}=u(x,t),\\ & \sqrt{P}+{\mathrm{\Phi }}_2(x,t)e^{\mathrm{P\delta ϵt}}=v(x,t),\\ & \sqrt{P}+{\mathrm{\Phi }}_3(x,t)e^{\mathrm{P\delta ϵt}}=w(x,t).\end{array}$(35) ) into (13(13) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}u}{\mathrm{\partial }{x}^3}+\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}+\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}v}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}w}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x},\end{array}$(13) ), we have (36) $\begin{array}{rl}& \frac{\mathrm{\partial }{\mathrm{\Phi }}_1}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_1-\sqrt{P}\frac{\mathrm{\partial }}{\mathrm{\partial }x}(\alpha {\mathrm{\Phi }}_1+(\alpha +\gamma ){\mathrm{\Phi }}_2+\gamma {\mathrm{\Phi }}_3)\\ & +\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_1}{\mathrm{\partial }{x}^3}=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Phi }}_2}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_2-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Phi }}_1+{\mathrm{\Phi }}_3)-\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_2}{\mathrm{\partial }{x}^3}=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Phi }}_3}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_3-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Phi }}_1+{\mathrm{\Phi }}_2)-\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_3}{\mathrm{\partial }{x}^3}=0.\end{array}$(36) It is supposed that the solution of Equation (36(36) $\begin{array}{rl}& \frac{\mathrm{\partial }{\mathrm{\Phi }}_1}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_1-\sqrt{P}\frac{\mathrm{\partial }}{\mathrm{\partial }x}(\alpha {\mathrm{\Phi }}_1+(\alpha +\gamma ){\mathrm{\Phi }}_2+\gamma {\mathrm{\Phi }}_3)\\ & +\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_1}{\mathrm{\partial }{x}^3}=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Phi }}_2}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_2-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Phi }}_1+{\mathrm{\Phi }}_3)-\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_2}{\mathrm{\partial }{x}^3}=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Phi }}_3}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_3-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Phi }}_1+{\mathrm{\Phi }}_2)-\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_3}{\mathrm{\partial }{x}^3}=0.\end{array}$(36) ) has as (37) $\begin{array}{rl}{\mathrm{\Phi }}_1(x,t)& ={\rho }_1{e}^{(\mathrm{kx}-\mathrm{\omega t})},{\mathrm{\Phi }}_2(x,t)={\rho }_2{e}^{(\mathrm{kx}-\mathrm{\omega t})},\\ {\mathrm{\Phi }}_3(x,t)& ={\rho }_3{e}^{(\mathrm{kx}-\mathrm{\omega t})},\end{array}$(37) where k and ω are the wave number and frequency of perturbation, respectively. Utilizing Equations (37(37) $\begin{array}{rl}{\mathrm{\Phi }}_1(x,t)& ={\rho }_1{e}^{(\mathrm{kx}-\mathrm{\omega t})},{\mathrm{\Phi }}_2(x,t)={\rho }_2{e}^{(\mathrm{kx}-\mathrm{\omega t})},\\ {\mathrm{\Phi }}_3(x,t)& ={\rho }_3{e}^{(\mathrm{kx}-\mathrm{\omega t})},\end{array}$(37) ) into (36(36) $\begin{array}{rl}& \frac{\mathrm{\partial }{\mathrm{\Phi }}_1}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_1-\sqrt{P}\frac{\mathrm{\partial }}{\mathrm{\partial }x}(\alpha {\mathrm{\Phi }}_1+(\alpha +\gamma ){\mathrm{\Phi }}_2+\gamma {\mathrm{\Phi }}_3)\\ & +\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_1}{\mathrm{\partial }{x}^3}=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Phi }}_2}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_2-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Phi }}_1+{\mathrm{\Phi }}_3)-\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_2}{\mathrm{\partial }{x}^3}=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Phi }}_3}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Phi }}_3-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Phi }}_1+{\mathrm{\Phi }}_2)-\beta \frac{{\mathrm{\partial }}^3{\mathrm{\Phi }}_3}{\mathrm{\partial }{x}^3}=0.\end{array}$(36) ), the dispersion relation (DR) is acquired as (38) $\begin{array}{rl}\omega & =\frac{1}{2}(2\mathrm{\delta Pϵ}-2\beta k^3-\mathrm{\alpha k}\sqrt{P}-\mathrm{\lambda k}\sqrt{P}\\ & \pm \sqrt{k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)}).\end{array}$(38) The dispersion relations in (38(38) $\begin{array}{rl}\omega & =\frac{1}{2}(2\mathrm{\delta Pϵ}-2\beta k^3-\mathrm{\alpha k}\sqrt{P}-\mathrm{\lambda k}\sqrt{P}\\ & \pm \sqrt{k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)}).\end{array}$(38) ) indicate that steady-state stability varies upon the wave number, self-phase modulation and stimulating Raman scattering. For wave numbers k, the velocity dispersion ω is real and the steady state is stable along the small perturbation if $k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)>0$. when $k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)<0$, it becomes unstable, i.e. ω is imaginary as the perturbation builds exponentially. It can be easily observed MI(modulation instability) when $k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)<0$. These conditions indicate the increased rate of MI achieved spectrum $g(k)$ can be stated as (39) $g(k)=\sqrt{k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)}.$(39)

4.2. The new coupled Z-K system

The following form is the SSS of the NCZ-K system (40) $\begin{array}{rl}& \sqrt{P}+{\mathrm{\Psi }}_1(x,y,t)e^{\mathrm{P\delta ϵt}}=u(x,y,t),\\ & \sqrt{P}+{\mathrm{\Psi }}_2(x,y,t)e^{\mathrm{P\delta ϵt}}=v(x,y,t),\\ & \sqrt{P}+{\mathrm{\Psi }}_3(x,y,t)e^{\mathrm{P\delta ϵt}}=w(x,y,t).\end{array}$(40) Here, the normalized optical power is P. Evolution of the perturbation $\mathrm{\Psi }(x,y,t)$ is verified by using the LS analysis. We are using Equations (25(25) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}-\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{wu})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{y}^2})& =0,\end{array}$(25) ) and (40(40) $\begin{array}{rl}& \sqrt{P}+{\mathrm{\Psi }}_1(x,y,t)e^{\mathrm{P\delta ϵt}}=u(x,y,t),\\ & \sqrt{P}+{\mathrm{\Psi }}_2(x,y,t)e^{\mathrm{P\delta ϵt}}=v(x,y,t),\\ & \sqrt{P}+{\mathrm{\Psi }}_3(x,y,t)e^{\mathrm{P\delta ϵt}}=w(x,y,t).\end{array}$(40) ) and linearizing (41) $\begin{array}{rl}& \frac{\mathrm{\partial }{\mathrm{\Psi }}_1}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_1-\sqrt{P}\frac{\mathrm{\partial }}{\mathrm{\partial }x}(\alpha {\mathrm{\Psi }}_1+(\alpha +\gamma ){\mathrm{\Psi }}_2+\gamma {\mathrm{\Psi }}_3)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_1}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_1}{\mathrm{\partial }{y}^2})=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Psi }}_2}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_2-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Psi }}_1+{\mathrm{\Psi }}_3)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_2}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_2}{\mathrm{\partial }{y}^2})=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Psi }}_3}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_3-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_3}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_3}{\mathrm{\partial }{y}^2})=0.\end{array}$(41) It is supposed that the solution of Equation (41(41) $\begin{array}{rl}& \frac{\mathrm{\partial }{\mathrm{\Psi }}_1}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_1-\sqrt{P}\frac{\mathrm{\partial }}{\mathrm{\partial }x}(\alpha {\mathrm{\Psi }}_1+(\alpha +\gamma ){\mathrm{\Psi }}_2+\gamma {\mathrm{\Psi }}_3)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_1}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_1}{\mathrm{\partial }{y}^2})=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Psi }}_2}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_2-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Psi }}_1+{\mathrm{\Psi }}_3)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_2}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_2}{\mathrm{\partial }{y}^2})=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Psi }}_3}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_3-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_3}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_3}{\mathrm{\partial }{y}^2})=0.\end{array}$(41) ) has as (42) $\begin{array}{rl}{\mathrm{\Psi }}_1(x,y,t)& ={\rho }_1{e}^{({\sigma }_{1}x+{\sigma }_{2}y-\mathrm{\nu t})},\\ {\mathrm{\Psi }}_2(x,y,t)& ={\rho }_2{e}^{({\sigma }_{1}x+{\sigma }_{2}y-\mathrm{\nu t})},\\ {\mathrm{\Psi }}_3(x,y,t)& ={\rho }_3{e}^{({\sigma }_{1}x+{\sigma }_{2}y-\mathrm{\nu t})},\end{array}$(42) where ν and ${\sigma }_1,{\sigma }_2$ are the frequency of perturbation and wave numbers, respectively. By using Equations (42(42) $\begin{array}{rl}{\mathrm{\Psi }}_1(x,y,t)& ={\rho }_1{e}^{({\sigma }_{1}x+{\sigma }_{2}y-\mathrm{\nu t})},\\ {\mathrm{\Psi }}_2(x,y,t)& ={\rho }_2{e}^{({\sigma }_{1}x+{\sigma }_{2}y-\mathrm{\nu t})},\\ {\mathrm{\Psi }}_3(x,y,t)& ={\rho }_3{e}^{({\sigma }_{1}x+{\sigma }_{2}y-\mathrm{\nu t})},\end{array}$(42) ) and (41(41) $\begin{array}{rl}& \frac{\mathrm{\partial }{\mathrm{\Psi }}_1}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_1-\sqrt{P}\frac{\mathrm{\partial }}{\mathrm{\partial }x}(\alpha {\mathrm{\Psi }}_1+(\alpha +\gamma ){\mathrm{\Psi }}_2+\gamma {\mathrm{\Psi }}_3)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_1}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_1}{\mathrm{\partial }{y}^2})=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Psi }}_2}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_2-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Psi }}_1+{\mathrm{\Psi }}_3)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_2}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_2}{\mathrm{\partial }{y}^2})=0,\\ & \frac{\mathrm{\partial }{\mathrm{\Psi }}_3}{\mathrm{\partial }t}+\mathrm{P\delta ϵ}{\mathrm{\Psi }}_3-\sqrt{P}\lambda \frac{\mathrm{\partial }}{\mathrm{\partial }x}({\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2)\\ & -\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_3}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^2{\mathrm{\Psi }}_3}{\mathrm{\partial }{y}^2})=0.\end{array}$(41) ), the following dispersion relation is obtained (43) $\begin{array}{rl}\nu & =\frac{1}{2}(-2\beta {\sigma }_1^3-2\beta {\sigma }_2^2{\sigma }_1-\alpha \sqrt{P}{\sigma }_1+2\mathrm{\delta Pϵ}-\lambda \sqrt{P}{\sigma }_1\\ & \pm \sqrt{{\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2}).\end{array}$(43) The dispersion relations in (43(43) $\begin{array}{rl}\nu & =\frac{1}{2}(-2\beta {\sigma }_1^3-2\beta {\sigma }_2^2{\sigma }_1-\alpha \sqrt{P}{\sigma }_1+2\mathrm{\delta Pϵ}-\lambda \sqrt{P}{\sigma }_1\\ & \pm \sqrt{{\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2}).\end{array}$(43) )) indicate that the steady-state stability varies upon the wave number, self-phase modulation and stimulating Raman scattering. For all wave numbers ${\sigma }_1,{\sigma }_2$, the velocity dispersion ν is real and the steady state is stable along the small perturbation if ${\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2>0$. It turns out to be unstable if ${\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2<0$, i.e. ν is imaginary as the perturbation builds exponentially. It can easily be observed modulational instability (MI) if ${\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2<0$. These conditions indicate that the increased rate of MI gain spectrum $g(k)$ can be conveyed as (44) $g(k)=\sqrt{{\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2}.$(44)

5. Discussion and physical interpretation of results

The acquired results in this article are dissimilar of the gained results of several academics because Equation (1(1) $G^{″}(\xi )+\mathrm{\sigma G}(\xi )=\mu ,$(1) ) is dissimilar from the existing methods. By assigning precise values of parameters, different families of solutions have been obtained for ODE (1(1) $G^{″}(\xi )+\mathrm{\sigma G}(\xi )=\mu ,$(1) ). The new coupled KdV system was resolved by Hirota's bilinear approach [17]. In [52], modified extended direct algebraic approach has been utilized to get precise wave results of (13(13) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}u}{\mathrm{\partial }{x}^3}+\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}+\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}v}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uw})}{\mathrm{\partial }x},\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}& =\beta \frac{{\mathrm{\partial }}^{3}w}{\mathrm{\partial }{x}^3}+\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x},\end{array}$(13) ) and (25(25) $\begin{array}{rl}\frac{\mathrm{\partial }u}{\mathrm{\partial }t}-\alpha \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\gamma \frac{\mathrm{\partial }(\mathrm{vw})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}u}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }v}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{wu})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}v}{\mathrm{\partial }{y}^2})& =0,\\ \frac{\mathrm{\partial }w}{\mathrm{\partial }t}-\lambda \frac{\mathrm{\partial }(\mathrm{uv})}{\mathrm{\partial }x}-\beta \frac{\mathrm{\partial }}{\mathrm{\partial }x}(\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{x}^2}+\frac{{\mathrm{\partial }}^{2}w}{\mathrm{\partial }{y}^2})& =0,\end{array}$(25) ) in travelling wave solutions in trignometric and elliptic functions forms. In [53], the researchers found Soliton solutions for a KdV equation and a generalized Hirota–Satsuma CKdV equation. By a natural decomposition method, a new solution of coupled KdV equation has been obtained in [54]. Sine-Gordon and Modified Kudryashov methods were used to evaluate dual-mode Hirota-Satsuma coupled KdV equations in [55]. The main advantage of this method is that different forms of analytical solutions of such as solitons, solitary waves, rational solitons, trigonometric function, hyperbolic function and other wave solutions are constructed using this technique. If $A_1=0$ or $A_2=0$ in hyperbolic function solutions, then we can obtain an analytical one-soliton solution; otherwise, we can get two soliton solutions. So, in this work, several innovative results have been accomplished, which never have been stated earlier.

In Figure 1, by specifying suitable values of parameters, the wave solution (19(19) $\begin{array}{rl}{u}_1(\xi )& =a_0+\frac{\begin{array}{c}6\beta {\delta }^2\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)\\ +A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2\end{array}}{\begin{array}{c}\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)\\ +A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2\end{array}},\\ v_1(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}+\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})\end{array}}{2\gamma },\\ w_1(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}+\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})\end{array}}{2\gamma }.\end{array}$(19) ) in disparate constitution is represented as follows: Figure 1(A,B,E) describes the bright soliton and their 2D contour plot figures (Figure 1(B,D ,F)) correspondingly. In Figure 2, by specifying  proper values of parameters, the wave solution (20(20) $\begin{array}{rl}{u}_2(\xi )& =\frac{6\beta {\delta }^2\sigma (A_1^2-A_2^2)}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ v_2(\xi )& =\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2-A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ w_2(\xi )& =\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2-A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}.\end{array}$(20) ) is shown in disparate forms as follows: Figure 2(A,C,E) are solions in multi-peak forms of dissimilar amplitude and their 2D contour plots in figures (Figure 2(B,D ,F)) correspondingly. In Figure 3, by specifying suitable values of parameters, the solution (21(21) $\begin{array}{rl}{u}_3(\xi )& =a_0-\frac{\begin{array}{c}6\beta {\delta }^2\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)\\ -A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2\end{array}}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ v_3(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}-\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2})\end{array}}{2\gamma },\\ w_3(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}-\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2})\end{array}}{2\gamma }.\end{array}$(21) ) is shown in different forms as follows: Figure 3(A,C,E) are solitons in periodic form of dissimilar amplitude and their 2D contourplots in Figure 3(B,D ,F), correspondingly.

Figure 1. Appropriate parameter values to result (19(19) $\begin{array}{rl}{u}_1(\xi )& =a_0+\frac{\begin{array}{c}6\beta {\delta }^2\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)\\ +A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2\end{array}}{\begin{array}{c}\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)\\ +A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2\end{array}},\\ v_1(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}+\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})\end{array}}{2\gamma },\\ w_1(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}+\frac{6\sigma (A_2\sinh \left(\xi \sqrt{-\sigma }\right)+A_1\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2})\end{array}}{2\gamma }.\end{array}$(19) ) are illustrated as follows: Figures (A,B,E) describe the bright soliton and their 2D contour plot figures (B,D,F), respectively, at $\delta =1,\gamma =0.5,\beta =1.5,A_1=1.5,A_2=0.1,\sigma =-1,\lambda =1,\alpha =0,a_0=1$.

Figure 2. Appropriate parameter values to result (20(20) $\begin{array}{rl}{u}_2(\xi )& =\frac{6\beta {\delta }^2\sigma (A_1^2-A_2^2)}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ v_2(\xi )& =\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2-A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ w_2(\xi )& =\frac{3\mathrm{\sigma \beta }{\delta }^2(A_1^2-A_2^2)(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}.\end{array}$(20) ) are depicted as follows: Figures (A,C,E) are multi-peak solitons with various amplitudes and their 2D contourplots in figures (B,D,F) respectively, at $\delta =1,\gamma =0.75,\beta =0.5,A_1=0.5,A_2=0.1,\sigma =-0.5,\lambda =1$.

Figure 3. Appropriate parameter values to the result (21(21) $\begin{array}{rl}{u}_3(\xi )& =a_0-\frac{\begin{array}{c}6\beta {\delta }^2\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)\\ -A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2\end{array}}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ v_3(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}-\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2})\end{array}}{2\gamma },\\ w_3(\xi )& =\frac{\begin{array}{c}\beta {\delta }^2(\sqrt{{\alpha }^2+4\mathrm{\gamma \lambda }}-\alpha )\\ (\frac{a_0}{\beta {\delta }^2}-\frac{6\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{\lambda (A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2})\end{array}}{2\gamma }.\end{array}$(21) ) are depicted as follows: Figures (A,C,E) are periodic solitons with various amplitudes and their 2D contourplots in figures (B,D,F), respectively, at $\delta =1,\gamma =2,\beta =3,A_1=1.5,A_2=0.1,\sigma =1,\lambda =1,\alpha =0,a_0=1$.

In Figure 4, by specifying suitable values of parameters, the wave solution (30(30) $\begin{array}{rl}{u}_2(\xi )& =\frac{2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)(2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ v_2(\xi )& =-\frac{\sigma (2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ w_2(\xi )& =-\frac{\sigma (2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}.\end{array}$(30) ) is illustrated in disparate shapes as follows: Figure 4(A,C,E) are solitons in multi-peak forms and their 2D contourplot in Figure 4(B,D ,F) correspondingly. In Figure 5, by specifying suitable parameter values, the solution (32(32) $\begin{array}{r}\begin{array}{rl}{u}_4(\xi )& =\frac{2\beta ({\nu }^2-{\delta }^2)(\frac{3\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}+\sigma )}{\lambda },\\ v_4(\xi )& =\frac{\sigma (-2A_2{A}_1\sin \left(2\xi \sqrt{\sigma }\right)+A_1^2(\cos \left(2\xi \sqrt{\sigma }\right)+2)-A_2^2(\cos \left(2\xi \sqrt{\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ w_4(\xi )& =\frac{\sigma (-2A_2{A}_1\sin \left(2\xi \sqrt{\sigma }\right)+A_1^2(\cos \left(2\xi \sqrt{\sigma }\right)+2)-A_2^2(\cos \left(2\xi \sqrt{\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}.\end{array}\end{array}$(32) ) is illustrated in various shapes as follows: Periodic solitons with varying amplitudes are shown in Figure 5(A,C, E), and their 2D contourplots are shown in Figure 5(B,D, F), respectively. Figure 6(A) depicts the DR between the ω and k of (38(38) $\begin{array}{rl}\omega & =\frac{1}{2}(2\mathrm{\delta Pϵ}-2\beta k^3-\mathrm{\alpha k}\sqrt{P}-\mathrm{\lambda k}\sqrt{P}\\ & \pm \sqrt{k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)}).\end{array}$(38) ) and the DR between ν and wave numbers (${\omega }_1,{\omega }_2$) of (43(43) $\begin{array}{rl}\nu & =\frac{1}{2}(-2\beta {\sigma }_1^3-2\beta {\sigma }_2^2{\sigma }_1-\alpha \sqrt{P}{\sigma }_1+2\mathrm{\delta Pϵ}-\lambda \sqrt{P}{\sigma }_1\\ & \pm \sqrt{{\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2}).\end{array}$(43) ) is demonstrated in Figure 6(B), respectively.

Figure 4. Appropriate parameter values to the result (30(30) $\begin{array}{rl}{u}_2(\xi )& =\frac{2\mathrm{\beta \sigma }({\delta }^2-{\nu }^2)(2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))}{\lambda (A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ v_2(\xi )& =-\frac{\sigma (2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2},\\ w_2(\xi )& =-\frac{\sigma (2A_2{A}_1\sinh \left(2\xi \sqrt{-\sigma }\right)+A_1^2(\cosh \left(2\xi \sqrt{-\sigma }\right)+2)+A_2^2(\cosh \left(2\xi \sqrt{-\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sinh \left(\xi \sqrt{-\sigma }\right)+A_2\cosh \left(\xi \sqrt{-\sigma }\right)){}^2}.\end{array}$(30) ) are followed as follows: Figures (A,C,E) are the solitons of multi-peak and their 2D contourplot in figures (B,D,F), respectively, at $\delta =0.3,\gamma =0.8,\beta =1,A_1=0.5,A_2=0.1,\sigma =-0.1,\lambda =1,\alpha =0,a_0=1,y=0.1,\nu =1$.

Figure 5. Appropriate parameter values to result (32(32) $\begin{array}{r}\begin{array}{rl}{u}_4(\xi )& =\frac{2\beta ({\nu }^2-{\delta }^2)(\frac{3\sigma (A_1\cos \left(\xi \sqrt{\sigma }\right)-A_2\sin \left(\xi \sqrt{\sigma }\right)){}^2}{(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}+\sigma )}{\lambda },\\ v_4(\xi )& =\frac{\sigma (-2A_2{A}_1\sin \left(2\xi \sqrt{\sigma }\right)+A_1^2(\cos \left(2\xi \sqrt{\sigma }\right)+2)-A_2^2(\cos \left(2\xi \sqrt{\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2},\\ w_4(\xi )& =\frac{\sigma (-2A_2{A}_1\sin \left(2\xi \sqrt{\sigma }\right)+A_1^2(\cos \left(2\xi \sqrt{\sigma }\right)+2)-A_2^2(\cos \left(2\xi \sqrt{\sigma }\right)-2))(\sqrt{{\beta }^2{({\delta }^2-{\nu }^2)}^2({\alpha }^2+4\mathrm{\gamma \lambda })}+\mathrm{\alpha \beta }({\delta }^2-{\nu }^2))}{\mathrm{\gamma \lambda }(A_1\sin \left(\xi \sqrt{\sigma }\right)+A_2\cos \left(\xi \sqrt{\sigma }\right)){}^2}.\end{array}\end{array}$(32) ) are illustrated as follows: Figures (A,C,E) are periodic solitons with dissimilar amplitude and their 2D Contourplots in figures (B,D,F) respectively, at $\delta =0.3,\gamma =0.8,\beta =1,A_1=0.5,A_2=0.1,\sigma =-0.1,\lambda =1,a_0=1,y=0.1,\nu =1$.

Figure 6. The DR between frequency(ω) and wave number (k) of (38(38) $\begin{array}{rl}\omega & =\frac{1}{2}(2\mathrm{\delta Pϵ}-2\beta k^3-\mathrm{\alpha k}\sqrt{P}-\mathrm{\lambda k}\sqrt{P}\\ & \pm \sqrt{k^{2}P({\alpha }^2+2\mathrm{\alpha \lambda }+8\mathrm{\gamma \lambda }+{\lambda }^2)}).\end{array}$(38) ) is shown in (A) and DR between frequency(ν) and wave numbers (${\omega }_1,{\omega }_2$) of (43(43) $\begin{array}{rl}\nu & =\frac{1}{2}(-2\beta {\sigma }_1^3-2\beta {\sigma }_2^2{\sigma }_1-\alpha \sqrt{P}{\sigma }_1+2\mathrm{\delta Pϵ}-\lambda \sqrt{P}{\sigma }_1\\ & \pm \sqrt{{\alpha }^{2}P{\sigma }_1^2+2\mathrm{\alpha \lambda P}{\sigma }_1^2+8\mathrm{\gamma \lambda P}{\sigma }_1^2+{\lambda }^{2}P{\sigma }_1^2}).\end{array}$(43) ) is shown in (B).

6. Conclusion

We effectively applied the proposed technique to the NCKdV and the NCZ-K systems in this work. The KdV equation, along with shallow-water waves and long internal waves in oceans, basically explains how long, one-dimensional waves propagate in a variety of physical conditions. The study of coastal waves on the basis of the ocean is done using the Zakharov–Kuznetsov (Z-K) equation and this model is utilized to illustrate ion-acoustic wave propagation. As a result, different forms of analytical solutions of the new coupled KdV (NCKdV) system and new coupled Z-K (NCZ-K) system, such as solitons, multi-peak solitons, solitary waves, trigonometric, hyperbolic and rational function solutions and other wave solutions, are constructed in the explicit form using the proposed scheme. By giving the parameters, the appropriate values, novel structures of constructed solutions of these models are represented. Graphical representations of the physical structures of a few obtained results are efficient for conveying the intricate physical existences of both models. The stability of the solution is examined by utilizing modulation instability. The results obtained show that the proposed technique is universal and efficient, and can be applied on lots of other new coupled systems.

Disclosure statement

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