) \,\,\,\, \text{and} \,\,\,\, {B}_1=\left(
{\begin{array}{*{20}c}
    -b & 0  & 0 & 0 \\
    0  & -a &  0 & 0\\
    0 & 0  & b & 0 \\
    0 & 0 & 0  & a\\
    \end{array}}
\right).
\end{split}$$ Let us point out that throughout this section, we will study representation of solutions $U=U(t,x)$ to the following Cauchy problem by deriving representation of its partial Fourier transform $W(t,\xi)={\mathcal}{F}_{x\rightarrow \xi}(U(t,x))$: $$\label{Eq.InvFirstOrderSystem}
\left\{\begin{aligned}
&U_t+\frac{1}{2}{B}_0(-\Delta)^{\rho}U+i{B}_1(-\Delta)^{\frac{1}{2}}U+\frac{1}{2}{B}_0(-\Delta)^{\theta}U=0,&&x\in{\mathbb}{R}^2,\,\,t>0,\\
&U(0,x)=U_0(x),&&x\in{\mathbb}{R}^2,
\end{aligned}\right.$$ where the coefficient matrices $B_0$ and $B_1$ are given in. Moreover, to derive qualitative properties of solutions to, we only need to study the solutions to.

With the aim of deriving representation of solutions, we may apply WKB analysis and multistep diagonalization procedure (see for example [@ReissigWang2005; @Yagdjian1997; @Jachmann2009; @JachmannReissig2008; @JachmannReissig2009; @Reissig2016; @Chen2019TEP]). Before doing these, we should understand the influence of the parameter $|\xi|$ on the asymptotic behavior of solutions to. Due to our assumption $0\leq 2\rho<1<2\theta\leq2$, we now discuss the influence of $|\xi|$ by three parts. Specifically, we will apply diagonalization procedure for small frequencies $\xi\in{\mathcal}{Z}_{\operatorname{int}}(\varepsilon)$ and large frequencies $\xi\in{\mathcal}{Z}_{\operatorname{ext}}(N)$ in Subsections \[SubSec.SmallFreq\] and \[SubSec.Large
) \,\,\,\, \text{and } \,\,\,\, { B}_1=\left (
 { \begin{array}{*{20}c } 
     -b & 0   & 0 & 0 \\ 
     0   & -a &   0 & 0\\ 
     0 & 0   & b & 0 \\ 
     0 & 0 & 0   & a\\ 
     \end{array } } 
 \right). 
 \end{split}$$ Let us point out that throughout this section, we will study theatrical performance of solution $ U = U(t, x)$ to the following Cauchy problem by derive theatrical performance of its partial Fourier transform $ W(t,\xi)={\mathcal}{F}_{x\rightarrow \xi}(U(t, x))$: $ $ \label{Eq. InvFirstOrderSystem } 
 \left\{\begin{aligned } 
 & U_t+\frac{1}{2}{B}_0(-\Delta)^{\rho}U+i{B}_1(-\Delta)^{\frac{1}{2}}U+\frac{1}{2}{B}_0(-\Delta)^{\theta}U=0,&&x\in{\mathbb}{R}^2,\,\,t>0,\\ 
 & U(0,x)=U_0(x),&&x\in{\mathbb}{R}^2, 
 \end{aligned}\right.$$ where the coefficient matrices $ B_0 $ and $ B_1 $ are given in. Moreover, to deduce qualitative property of solutions to, we merely need to study the solution to. 

 With the aim of deriving theatrical performance of solutions, we may apply WKB analysis and multistep diagonalization procedure (see for example [ @ReissigWang2005; @Yagdjian1997; @Jachmann2009; @JachmannReissig2008; @JachmannReissig2009; @Reissig2016; @Chen2019TEP ]). Before doing these, we should sympathize the influence of the parameter $ |\xi|$ on the asymptotic behavior of solutions to. Due to our presumption $ 0\leq 2\rho<1<2\theta\leq2 $, we now discuss the influence of $ |\xi|$ by three share. Specifically, we will enforce diagonalization procedure for little frequencies $ \xi\in{\mathcal}{Z}_{\operatorname{int}}(\varepsilon)$ and large frequencies $ \xi\in{\mathcal}{Z}_{\operatorname{ext}}(N)$ in Subsections \[SubSec. SmallFreq\ ] and \[SubSec. big
) \,\,\,\, \tfxt{and} \,\,\,\, {B}_1=\left(
{\begin{array}{*{20}c}
    -b & 0  & 0 & 0 \\
    0  & -a &  0 & 0\\
    0 & 0  & b & 0 \\
    0 & 0 & 0  & a\\
    \end{array}}
\rmght).
\wnd{spoit}$$ Let us point out tfat throuhhout thus stction, we will stnsy reprcfentzbion mh solutions $U=U(t,w)$ to the fonlowing Cauchy pfoylem by deriving representation of ies partoap Fourier tranffork $W(t,\xj)={\mathcal}{F}_{x\rightarrow \xi}(U(t,x))$: $$\label{Eq.JnvFirsuOrderSystem}
\left\{\bebin{aligned}
&U_t+\frac{1}{2}{B}_0(-\Delta)^{\rho}U+l{B}_1(-\Depta)^{\frac{1}{2}}U+\frac{1}{2}{B}_0(-\Delta)^{\hheta}U=0,&&x\in{\majgbb}{W}^2,\,\,r>0,\\
&U(0,x)=U_0(x),&&x\in{\mathcb}{R}^2,
\end{aligned}\right.$$ whete the coefficient matrices $B_0$ ana $B_1$ axe given in. Mireleer, to derite quajitative proivrties mf soluyions to, we onky ieed to study the solutiois to.

With the aim of deriving rzpresentation of soluriins, wg may apput WYB znelyais anf mnltistep dizgonalizatiin procedure (see fot qqsmple [@ReissifWang2005; @Raddjian1997; @Jachmann2009; @JachmannReissig2008; @JachmantRejssig2009; @Reissig2016; @Chen2019TEP]). Bwfore doing these, we dhould unqerstand the influence of the parameter $|\xi|$ on the asym'tutie behaxuog of solutions to. Due to our assumption $0\leq 2\rrk<1<2\tnena\leq2$, we now discmss the influence pf $|\xo|$ by three parjs. Specnridally, we will applj diagogalizqtion prosedute for small frequencies $\xi\un{\mathcal}{Z}_{\opvratirname{int}}(\varepsilou)$ and large yrequemcies $\xi\in{\mathcal}{Z}_{\operatornaoe{exf}}(N)$ in Subseftions \[SugRec.SmallFreq\] and \[SunSew.Large
) \,\,\,\, \text{and} \,\,\,\, {B}_1=\left( {\begin{array}{*{20}c} -b & & 0 0 & -a & & b & \\ 0 & & 0 & a\\ \end{array}} \right). Let us point out that throughout this section, we will study representation of $U=U(t,x)$ to the following Cauchy problem by deriving representation of its partial Fourier $W(t,\xi)={\mathcal}{F}_{x\rightarrow $$\label{Eq.InvFirstOrderSystem} &U_t+\frac{1}{2}{B}_0(-\Delta)^{\rho}U+i{B}_1(-\Delta)^{\frac{1}{2}}U+\frac{1}{2}{B}_0(-\Delta)^{\theta}U=0,&&x\in{\mathbb}{R}^2,\,\,t>0,\\ \end{aligned}\right.$$ where the coefficient matrices $B_0$ and $B_1$ are given in. Moreover, to derive qualitative properties solutions to, we only need to study the to. With the aim deriving representation of solutions, we apply analysis and diagonalization (see example [@ReissigWang2005; @Yagdjian1997; @JachmannReissig2008; @JachmannReissig2009; @Reissig2016; @Chen2019TEP]). Before doing these, we should understand the influence of the parameter $|\xi|$ on asymptotic behavior to. Due our $0\leq we now discuss of $|\xi|$ by three parts. Specifically, diagonalization procedure for small frequencies $\xi\in{\mathcal}{Z}_{\operatorname{int}}(\varepsilon)$ and large $\xi\in{\mathcal}{Z}_{\operatorname{ext}}(N)$ in \[SubSec.SmallFreq\] and \[SubSec.Large
) \,\,\,\, \text{and} \,\,\,\, {B}_1=\left(
{\begin{array}{*{20}c}
    -b & 0  & 0 & 0 \\
    0  & -a &  0 & 0\\
    0 & 0  & B & 0 \\
    0 & 0 & 0  & a\\
    \end{array}}
\Right).
\End{SplIt}$$ let uS poiNt out that throuGHout This section, we will study RepreSeNTatiON oF soluTions $U=U(T,X)$ tO THe fOlLoWinG CAUcHy proBleM by deriVing represEntAtIon of its partIAl fourier traNsfOrm $W(t,\xi)={\mathcAl}{F}_{X\rightArRow \XI}(U(t,x))$: $$\lAbeL{Eq.InVFirstoRderSyStem}
\left\{\bEgIN{alignED}
&U_t+\frac{1}{2}{b}_0(-\dElTa)^{\rhO}U+i{B}_1(-\Delta)^{\frac{1}{2}}U+\fraC{1}{2}{b}_0(-\DELta)^{\theta}U=0,&&x\in{\maThbb}{R}^2,\,\,t>0,\\
&u(0,x)=u_0(X),&&x\IN{\MatHbb}{r}^2,
\end{aligneD}\rIght.$$ wHEre the cOEfFICIenT Matrices $B_0$ and $B_1$ Are given in. MOReoVer, to dErIve QUalitaTive pRoPErtIes of solutiOns tO, we only neEd to stUDy the soLUtions tO.

With tHe aIm oF derIViNg RepReSEntATiOn oF SolUtions, we MaY aPply WkB anALYSIs anD muLtisTep diAgonalization ProCeduRE (seE for eXamplE [@ReiSsIgWanG2005; @YagdjIan1997; @JaChMann2009; @JachmannReiSsig2008; @jachmannREisSiG2009; @ReIsSig2016; @ChEN2019TEP]). BeForE doIng thesE, we shouLD unDeRSTAnD the influence of the PaRAMeTer $|\xi|$ on tHe asymPToTiC Behavior Of SolUtioNS To. Due To ouR AsSumption $0\Leq 2\rho<1<2\THeTa\Leq2$, we noW dIscuss ThE inFluEnce oF $|\Xi|$ by Three pArts. SpecIficaLLy, we will apply dIAgonalization PRoCEDuRE for SmaLl frequenciEs $\xi\IN{\matHcal}{z}_{\OpEraTOrnamE{int}}(\vArEPsILon)$ and large frequencIeS $\xi\in{\mAthcaL}{Z}_{\operatornamE{ext}}(N)$ in SubSECTions \[Subsec.SMAlLfreq\] and \[SubSec.LArge
) \,\,\,\, \text{and} \,\, \,\, {B}_1 =\lef t({\b eg in{a rray }{*{20}c}
   - b &0  & 0 & 0 \\
    0  &  -a &  0  & 0 \ \    0  & 0  & b& 0 \ \   0  & 0& 0 & a \\
   \end{array }}\r ight).
\end{ s pl it}$$ Letuspoint out th atthroug ho utt his s ect ion,we wil l  study  represen ta t ion of solutio n s  $ U=U( t,x)$ to the foll o wi n g Cauchy probl em byde r iv i n g r epr esentation  o f its partial Fo u r i ert ransform $W(t ,\xi)={\mat h cal }{F}_{ x\ rig h tarrow  \xi} (U ( t,x ))$: $$\lab el{E q.InvFirs tOrder S ystem}\ left\{\ begin{ ali gne d}
& U _t +\ fra c{ 1 }{2 } {B }_0 ( -\D elta)^{\ rh o} U+i{B }_1( - \ D e lta) ^{\ frac {1}{2 }}U+\frac{1}{ 2}{ B}_0 ( -\D elta) ^{\th eta} U= 0,&&x \in{\m athbb }{ R}^2,\,\,t>0,\\ 
&U( 0,x)=U_0( x), && x\i n{ \math b b}{R}^ 2,\en d{align ed}\rig h t.$ $w h e re  the coefficient m at r i ce s $B_0$and $B _ 1$  a r e givenin . M oreo v e r, to  der i ve  qualita tive p r op er ties of  s olutio ns  to , w e onl y  nee d to s tudy the  solu t ions to.

With the aim of de r iv i n gr epre sen tation of s olut i ons,  wem ay  ap p ly WK B ana ly s is and multistep diago na lizati on pr ocedure (seefor exampl e [ @Reissig Wang 2 00 5 ; @Yagdjian199 7; @J achmann200 9 ; @Jachm annRe issig200 8; @Jachm a n nReissig 200 9;@Re iss i g 20 16; @Chen2019 T E P]).  B efore d oin g these , w e s hou ldun derstandthe infl ue nc eof  th e par a meter $| \x i|$  o n t he as y mptoti c beh avio rof sol utionst o. D ue t oou r as sum pt ion $ 0\le q  2\ rho<1<2 \theta\le q2$ ,  weno wdiscuss  the influenc eof $|\xi|$  b y t hree p a r ts. Spec ifically, we will apply diagona liz ation  pro cedure fo r s mall f req u encies  $\xi\ in{\m at hca l } {Z}_{ \ o pe rat or name{int}} ( \ var epsil on )$ a nd larg e frequencies $\xi \ in{ \mathcal}{Z}_ {\o pera t o rn ame { ex t }}( N) $  in S ubsections \[Su bSec.Small Fr e q\ ] and \[Su b Sec .L arge
) \,\,\,\,_\text{and} \,\,\,\,_{B}_1=\left(
{\begin{array}{*{20}c}
   _-b &_0_ &_0_& 0 \\
_  _0  & -a_&  0_&_0\\
    0 & 0  & b & 0 \\
 __ 0_&_0_& 0  & a\\
_   \end{array}}
\right).
\end{split}$$ Let_us point_out that throughout this section, we will study_representation_of solutions $U=U(t,x)$_to the following Cauchy problem by deriving representation of_its partial Fourier transform $W(t,\xi)={\mathcal}{F}_{x\rightarrow \xi}(U(t,x))$:_$$\label{Eq.InvFirstOrderSystem}
\left\{\begin{aligned}
&U_t+\frac{1}{2}{B}_0(-\Delta)^{\rho}U+i{B}_1(-\Delta)^{\frac{1}{2}}U+\frac{1}{2}{B}_0(-\Delta)^{\theta}U=0,&&x\in{\mathbb}{R}^2,\,\,t>0,\\
&U(0,x)=U_0(x),&&x\in{\mathbb}{R}^2,
\end{aligned}\right.$$ where the_coefficient_matrices_$B_0$ and $B_1$ are_given in. Moreover, to derive qualitative_properties of solutions to, we only_need to study the solutions to.

With the_aim of deriving representation of solutions,_we may apply WKB analysis_and multistep_diagonalization procedure (see for example_[@ReissigWang2005; @Yagdjian1997; @Jachmann2009;_@JachmannReissig2008; @JachmannReissig2009;_@Reissig2016; @Chen2019TEP]). Before_doing these, we should understand the_influence of the_parameter $|\xi|$ on the asymptotic behavior_of_solutions to. Due_to_our_assumption $0\leq_2\rho<1<2\theta\leq2$, we now_discuss_the influence_of_$|\xi|$ by three parts. Specifically, we_will_apply diagonalization procedure for small frequencies $\xi\in{\mathcal}{Z}_{\operatorname{int}}(\varepsilon)$_and large frequencies $\xi\in{\mathcal}{Z}_{\operatorname{ext}}(N)$_in_Subsections \[SubSec.SmallFreq\] and \[SubSec.Large