Eigenfunction Matching for a Semi-Infinite Floating Elastic Plate

Introduction

We show here a solution for a semi-infinite Floating Elastic Plate on Finite Depth. The problem was solved by Fox and Squire 1994 but the solution method here is slightly different. The simpler theory for a Dock describes many of the ideas here in more detail.

Equations

We consider the problem of small-amplitude waves which are incident on a semi-infinite floating elastic plate occupying water surface for [math]\displaystyle{ x\gt 0 }[/math]. The submergence of the plate is considered negligible. We assume that the problem is invariant in the [math]\displaystyle{ y }[/math] direction, although we allow the waves to be incident from an angle. We also assume that the plate edges are free to move at each boundary, although other boundary conditions could easily be considered using the methods of solution presented here. We begin with the Frequency Domain Problem for a semi-infinite Floating Elastic Plates in the non-dimensional form of Tayler 1986 (Dispersion Relation for a Floating Elastic Plate). We also assume that the waves are normally incident (incidence at an angle will be discussed later).

[math]\displaystyle{ \begin{matrix} \Delta \phi = 0, \;\;\; -h \lt z \leq 0, \end{matrix} }[/math] [math]\displaystyle{ \partial_z \phi = 0, \;\;\; z = - h, }[/math] [math]\displaystyle{ \partial_z\phi=\alpha\phi, \,\, z=0,\,x\lt 0, }[/math] [math]\displaystyle{ \beta \partial_x^4\partial_z \phi - \left( \gamma\alpha - 1 \right) \partial_z \phi + \alpha\phi = 0, \;\; z = 0, \;\;\; x \geq 0, }[/math]

where [math]\displaystyle{ \alpha = \omega^2 }[/math], [math]\displaystyle{ \beta }[/math] and [math]\displaystyle{ \gamma }[/math] are the stiffness and mass constant for the plate respectively. The free edge conditions at the edge of the plate imply

[math]\displaystyle{ \partial_x^3 \partial_z\phi = 0, \;\; z = 0, \;\;\; x = 0, }[/math] [math]\displaystyle{ \partial_x^2 \partial_z\phi = 0, \;\; z = 0, \;\;\; x = 0, }[/math]

Method of solution

Eigenfunction expansion

We will solve the system of equations using an Eigenfunction Matching Method. The method was developed by Fox and Squire 1994. The separation of variables for the left hand region where there is open water, and the incident wave, is described in Eigenfunction Matching for a Semi-Infinite Dock and we consider here only the separation of variables in the plate covered region.

Separation of variables under the Plate

The potential velocity can be written in terms of an infinite series of separated eigenfunctions under each elastic plate, of the form

[math]\displaystyle{ \psi = \frac{\cos(\kappa (z+h))}{\cos\kappa h}.\, }[/math]

If we apply the boundary conditions given we obtain the Dispersion Relation for a Floating Elastic Plate

[math]\displaystyle{ \kappa \tan{(\kappa h)}= -\frac{\alpha}{\beta \kappa^{4} + 1 - \alpha\gamma} }[/math]

Solving for [math]\displaystyle{ \kappa }[/math] gives a pure imaginary root with positive imaginary part, two complex roots (two complex conjugate paired roots with positive imaginary part in most physical situations), an infinite number of positive real roots which approach [math]\displaystyle{ {n\pi}/{h} }[/math] as [math]\displaystyle{ n }[/math] approaches infinity, and also the negative of all these roots (Dispersion Relation for a Floating Elastic Plate) . We denote the two complex roots with positive imaginary part by [math]\displaystyle{ \kappa_(-2) }[/math] and [math]\displaystyle{ \kappa_(-1) }[/math], the purely imaginary root with positive imaginary part by [math]\displaystyle{ \kappa_{0} }[/math] and the real roots with positive imaginary part by [math]\displaystyle{ \kappa_(n) }[/math] for [math]\displaystyle{ n }[/math] a positive integer. The imaginary root with positive imaginary part corresponds to a reflected travelling mode propagating along the [math]\displaystyle{ x }[/math] axis. The complex roots with positive imaginary parts correspond to damped reflected travelling modes and the real roots correspond to reflected evanescent modes.

Expressions for the potential velocity

We now expand the potential in the two regions using the separation of variables solution. The potential [math]\displaystyle{ \phi }[/math] can now be expressed as the following sum of eigenfunctions: be expanded as

[math]\displaystyle{ \phi(x,z)= e^{-k_{0}x}\phi_{0}(z) + \sum_{m=0}^{\infty}a_{m}e^{k_{m}x}\phi_{m}(z), \;\;x\lt 0 }[/math]

and

[math]\displaystyle{ \phi(x,z)=\sum_{m=-2}^{\infty}b_{m} e^{-\kappa_{m}x}\psi_{m}(z), \;\;x\gt 0 }[/math]

where [math]\displaystyle{ a_{m} }[/math] and [math]\displaystyle{ b_{m} }[/math] are the coefficients to be determined.

We define

[math]\displaystyle{ \int\nolimits_{-h}^{0}\phi_{m}(z)\phi_{n}(z) d z=A_{m}\delta_{mn} }[/math]

and

[math]\displaystyle{ \int\nolimits_{-h}^{0}\phi_{n}(z)\psi_{m}(z) d z=B_{mn} }[/math]

Expressions for [math]\displaystyle{ A_{mn} }[/math] and [math]\displaystyle{ B_{mn} }[/math] can be found in Eigenfunction Matching for a Semi-Infinite Dock.

An infinite dimensional system of equations

The potential and its derivative must be continuous across the transition from open water to the plate covered region. Therefore, the potentials and their derivatives at [math]\displaystyle{ x=0 }[/math] have to be equal. We also truncate the sum at [math]\displaystyle{ N }[/math] being careful that we have two extra modes on the plate covered region to satisfy the edge conditions. We obtain

[math]\displaystyle{ \phi_{0}\left( z\right) + \sum_{m=0}^{N} a_{m} \phi_{m}\left( z\right) =\sum_{m=-2}^{N}b_{m}\psi_{m}(z) }[/math]

and

[math]\displaystyle{ -k_{0}\phi_{0}\left( z\right) +\sum _{m=0}^{N} a_{m}k_{m}\phi_{m}\left( z\right) =-\sum_{m=-2}^{N}b_{m}\kappa_{m}\psi _{m}(z) }[/math]

for each [math]\displaystyle{ n }[/math]. We solve these equations by multiplying both equations by [math]\displaystyle{ \phi_{l}(z)\, }[/math] and integrating from [math]\displaystyle{ -h }[/math] to [math]\displaystyle{ 0 }[/math] to obtain:

[math]\displaystyle{ A_{0}\delta_{0l}+a_{l}A_{l} =\sum_{m=0}^{N}b_{m}B_{ml}+b_{-1}B_{-1l}+b_{-2}B_{-2l} }[/math]

and

[math]\displaystyle{ -k_{0}A_{0}\delta_{0l}+a_{l}k_{l}A_l =-\sum_{m=0}^{N}b_{m}\kappa_{m}B_{ml}+b_{-1}\kappa_{-1}B_{-1l}+b_{-2}\kappa_{-2}B_{-2l} }[/math]

We separate the terms for [math]\displaystyle{ m=-1 }[/math] and [math]\displaystyle{ m=-2 }[/math] because [math]\displaystyle{ l \in [0,N] }[/math]. If we multiply the first equation by [math]\displaystyle{ k_l }[/math] and subtract the second equation we obtain

[math]\displaystyle{ 2k_{0}A_{0}\delta_{0l} =\sum_{m=0}^{N}b_{m}(k_l + \kappa_{m})B_{ml}+b_{-1}(k_l+\kappa_{-1})B_{-1l}+b_{-2}(k_l+\kappa_{-2})B_{-2l} }[/math]

Finally, we need to apply the conditions at the edge of the plate to give us two further equations,

[math]\displaystyle{ \partial_x^2\partial_z\phi = - \sum_{m=-2}^{N}b_{m} \kappa_m^3 \tan\kappa_m h = 0 }[/math]

and

[math]\displaystyle{ \partial_x^3\partial_z\phi = \sum_{m=-2}^{N}b_{m} \kappa_m^4 \tan\kappa_m h = 0 }[/math]

These equations will permit us to calculate [math]\displaystyle{ b_{-1} }[/math] and [math]\displaystyle{ b_{-2} }[/math], knowing [math]\displaystyle{ b_m }[/math], [math]\displaystyle{ m \in [0,N] }[/math].

Numerical Solution

To solve the system of equations previously defined we set the upper limit of [math]\displaystyle{ l }[/math] to be [math]\displaystyle{ N }[/math], as we already said before. In terms of matrix, we obtain

[math]\displaystyle{ \begin{bmatrix} \begin{bmatrix} A_0&0 \quad \cdots&0\\ 0&&\\ \vdots&A_l\delta_{0l}&\vdots\\ &&0\\ 0&\cdots \quad 0 &A_N \end{bmatrix} & \begin{bmatrix} -B_{-20}&-B_{-10}\\ &\\ \vdots&\vdots\\ &\\ -B_{-2N}&-B_{-1N} \end{bmatrix} & \begin{bmatrix} -B_{00}&\cdots&-B_{0N}\\ &&\\ \vdots&-B_{ml}&\vdots\\ &&\\ -B_{N0}&\cdots&-B_{NN} \end{bmatrix} \\ \begin{bmatrix} 0&\cdots&0\\ 0&\cdots&0 \end{bmatrix} & \begin{bmatrix} \kappa_{-2}^3\tan\kappa_{-2}h&\kappa_{-1}^3\tan\kappa_{-1}h\\ \kappa_{-2}^4\tan\kappa_{-2}h&\kappa_{-1}^4\tan\kappa_{-1}h \end{bmatrix} & \begin{bmatrix} \kappa_0^3\tan\kappa_0h&\cdots&\kappa_N^3\tan\kappa_Nh\\ \kappa_0^4\tan\kappa_0h&\cdots&\kappa_N^4\tan\kappa_Nh \end{bmatrix} \\ \begin{bmatrix} 0&\cdots&0\\ \vdots&\ddots&\vdots\\ 0&\cdots&0 \end{bmatrix} & \begin{bmatrix} (k_0+\kappa_{-2})B_{-20}&(k_0+\kappa_{-1})B_{-10}\\ \vdots&\vdots\\ (k_N+\kappa_{-2})B_{-2N}&(k_N+\kappa_{-1})B_{-1N} \end{bmatrix} & \begin{bmatrix} (k_0 + \kappa_0) \, B_{00}&\cdots&(k_N + \kappa_{0}) \, B_{0N}\\ &&\\ \vdots&(k_l + \kappa_{m}) \, B_{ml}&\vdots\\ &&\\ (k_0 + \kappa_N) \, B_{N0}&\cdots&(k_N + \kappa_{N}) \, B_{NN}\\ \end{bmatrix} \end{bmatrix} \begin{bmatrix} a_{0} \\ a_{1} \\ \vdots \\ a_M \\ \\ b_{-2}\\ b_{-1}\\ \\ b_{0}\\ b_1 \\ \vdots \\ \\ b_M \end{bmatrix} = \begin{bmatrix} - A_{0} \\ 0 \\ \vdots \\ \\ 0 \\ \\ 0 \\ 0 \\ \\ 2k_{0}A_{0} \\ 0 \\ \vdots \\ \\ 0 \end{bmatrix} }[/math]

We then simply need to solve the linear system of equations.

Equations

We consider the problem of small-amplitude waves which are incident on a semi-infinite floating elastic plate occupying water surface for [math]\displaystyle{ x\gt 0 }[/math]. The submergence of the plate is considered negligible. We assume that the problem is invariant in the [math]\displaystyle{ y }[/math] direction, although we allow the waves to be incident from an angle. We also assume that the plate edges are free to move at each boundary, although other boundary conditions could easily be considered using the methods of solution presented here. We begin with the Frequency Domain Problem for a semi-infinite Floating Elastic Plates in the non-dimensional form of Tayler 1986 (Dispersion Relation for a Floating Elastic Plate)

[math]\displaystyle{ \begin{matrix} \left(\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial z^2} - k_y^2\right) \phi = 0, \;\;\;\; \mbox{ for } -h \lt z \leq 0, \end{matrix} }[/math] [math]\displaystyle{ \begin{matrix} \frac{\partial \phi}{\partial z} = 0, \;\;\;\; \mbox{ at } z = - h, \end{matrix} }[/math] [math]\displaystyle{ \begin{matrix} \left( \beta \left(\frac{\partial^2}{\partial x^2} - k^2_y\right)^2 - \gamma\alpha + 1\right)\frac{\partial \phi}{\partial z} - \alpha\phi = 0, \;\;\;\; \mbox{ at } z = 0, \;\;\; x \geq 0, \end{matrix} }[/math] [math]\displaystyle{ \begin{matrix} \frac{\partial \phi}{\partial z} - \alpha\phi = 0, \;\;\;\; \mbox{ at } z = 0, \;\;\; x \leq 0, \end{matrix} }[/math]

where [math]\displaystyle{ \alpha = \omega^2 }[/math], [math]\displaystyle{ \beta }[/math] and [math]\displaystyle{ \gamma }[/math] and the stiffness and mass constant for the plate. The conditions at the edge of the plate os

[math]\displaystyle{ \begin{matrix} \left(\frac{\partial^3}{\partial x^3} - (2 - \nu)k^2_y\frac{\partial}{\partial x}\right) \frac{\partial\phi}{\partial z}= 0, \;\;\;\; \mbox{ at } z = 0, \;\;\; \mbox{ for } x = 0, \end{matrix} }[/math] [math]\displaystyle{ \begin{matrix} \left(\frac{\partial^2}{\partial x^2} - \nu k^2_y\right)\frac{\partial\phi}{\partial z} = 0, \;\;\;\;\mbox{ at } z = 0, \;\;\; \mbox{ for } x = 0. \end{matrix} }[/math]

Method of solution

Eigenfunction expansion

We will solve the system of equations using an Eigenfunction Matching Method. The method was developed by Fox and Squire 1994. The separation of variables for the left hand region where there is open water is described in Eigenfunction Matching for a Semi-Infinite Dock and we consider here only the separation of variables in the plate covered region.

Separation of variables under the Plate

The potential velocity can be written in terms of an infinite series of separated eigenfunctions under each elastic plate, of the form [math]\displaystyle{ \phi = e^{\kappa_\mu x} \cos(k_\mu(z+h)). }[/math] If we apply the boundary conditions given we obtain the Dispersion Relation for a Floating Elastic Plate

[math]\displaystyle{ \begin{matrix} k_\mu\tan{(k_\mu h)}= & -\frac{\alpha}{\beta_\mu k_\mu^{4} + 1 - \alpha\gamma_\mu} \end{matrix} }[/math]

Solving for [math]\displaystyle{ k_\mu }[/math] gives a pure imaginary root with positive imaginary part, two complex roots (two complex conjugate paired roots with positive imaginary part in all physical situations), an infinite number of positive real roots which approach [math]\displaystyle{ {n\pi}/{h} }[/math] as [math]\displaystyle{ n }[/math] approaches infinity, and also the negative of all these roots (Dispersion Relation for a Floating Elastic Plate) . We denote the two complex roots with positive imaginary part by [math]\displaystyle{ k_\mu(-2) }[/math] and [math]\displaystyle{ k_\mu(-1) }[/math], the purely imaginary root with positive imaginary part by [math]\displaystyle{ k_\mu(0) }[/math] and the real roots with positive imaginary part by [math]\displaystyle{ k_\mu(n) }[/math] for [math]\displaystyle{ n }[/math] a positive integer. The imaginary root with positive imaginary part corresponds to a reflected travelling mode propagating along the [math]\displaystyle{ x }[/math] axis. The complex roots with positive imaginary parts correspond to damped reflected travelling modes and the real roots correspond to reflected evanescent modes. In a similar manner, the negative of these correspond to the transmitted travelling, damped and evanescent modes respectively. The coefficient [math]\displaystyle{ \kappa_\mu }[/math] is

[math]\displaystyle{ \kappa_\mu(n) = \sqrt{k_\mu(n)^2 + k_y^2}, }[/math]

where the root with positive real part is chosen or if the real part is negative with negative imaginary part.

Expressions for the potential velocity

We now expand the potential in the two regions using the separation of variables solution. We always include the two complex and one imaginary root under the plate and truncate the expansion at [math]\displaystyle{ M }[/math] real roots on both sides. The potential [math]\displaystyle{ \phi }[/math] can now be expressed as the following sum of eigenfunctions:

[math]\displaystyle{ \phi \approx \left\{ \begin{matrix} { Ie^{-\kappa_{0}x}\frac{\cos(k_0(z+h))}{\cos(k_0h)} }+ { \qquad \qquad \sum_{n=0}^{M}R_n e^{\kappa_n x} \frac{\cos(k_n(z+h))}{\cos(k_n h)} },& \mbox{ for } x \lt 0,\\ { \sum_{n=-2}^{M}T_n e^{-\kappa^\mu_n x} \frac{\cos(k^\mu_n(z+h))}{\cos(k^\mu_n h)} } \mbox{ for } x \gt 0, \end{matrix} \right. }[/math]

Note that the coefficients are normalised by the potential at the free surface rather than at the bottom surface.

Solving via eigenfunction matching

To solve for the coefficients, we require as many equations as we have unknowns. We derive the equations from the free edge conditions and from imposing conditions of continuity of the potential and its derivative in the [math]\displaystyle{ x }[/math]-direction at the plate-water boundary. We impose the latter condition by taking inner products with respect to the orthogonal functions [math]\displaystyle{ \cos \frac{m\pi}{h}(z+h) }[/math], where [math]\displaystyle{ m }[/math] is a natural number. The main reason we use these functions is because this is the method used in the solution for multiple plates. However, if would be better in this simple problem to use the eigenfunctions for the water.

Taking inner products leads to the following equations

[math]\displaystyle{ \begin{matrix} { \int_{-h}^0 \phi(0^-,z)\cos \frac{m\pi}{h}(z+h) \, dz } &=& { \int_{-h}^0 \phi(0^+,z)\cos \frac{m\pi}{h}(z+h) \, dz }\\ { \int_{-h}^0 \frac{\partial\phi}{\partial x}(0^-,z) \cos \frac{m\pi}{h}(z+h) \, dz } &=& { \int_{-h}^0 \frac{\partial\phi}{\partial x}(0^+,z) \cos \frac{m\pi}{h}(z+h) \, dz } \end{matrix} }[/math]

where [math]\displaystyle{ m\in[0,M] }[/math] and [math]\displaystyle{ \phi_\mu }[/math] denotes the potential under the [math]\displaystyle{ \mu }[/math]th plate, i.e. the expression for [math]\displaystyle{ \phi }[/math] valid for [math]\displaystyle{ l_\mu \lt x\lt r_\mu }[/math]. The remaining equations to be solved are given by the two edge conditions satisfied at both edges of each plate

[math]\displaystyle{ \begin{matrix} { \left(\frac{\partial^3}{\partial x^3} - (2 - \nu)k_y^2\frac{\partial}{\partial x}\right)\frac{\partial\phi_\mu}{\partial z} } &=&0, & \mbox{ for } z = 0 \mbox{ and } x = l_\mu,r_\mu,\\ { \left(\frac{\partial^2}{\partial x^2} - \nu k_y^2\right)\frac{\partial\phi_\mu}{\partial z} } &=&0, & \mbox{ for } z = 0 \mbox{ and } x = l_\mu,r_\mu. \end{matrix} }[/math]

We will show the explicit form of the linear system of equations which arise when we solve these equations. Let [math]\displaystyle{ {\mathbf T} }[/math] be a column vector given by [math]\displaystyle{ \left[T_{-2}, . . ., T_M\right]^{{\mathbf T}} }[/math] and [math]\displaystyle{ {\mathbf R} }[/math] be a column vector given by [math]\displaystyle{ \left[R_0 . . . R_M)\right]^{{\mathbf T}} }[/math].

The equations which arise from matching at the boundary between the water and plate are

[math]\displaystyle{ \begin{matrix} I{\mathbf C} + {\mathbf M}^{+} {\mathbf R} ={\mathbf M}^{-} {\mathbf T}\\ -\kappa_1(0)I\mathbf{C} + {\mathbf N}^{+} {\mathbf R} = {\mathbf N}^{-} {\mathbf T}. \end{matrix} }[/math]

where

[math]\displaystyle{ \begin{matrix} { {\mathbf M}^{+}(m,n) = \int_{-h}^0 \frac{\cos \left(k_n (z+h)\right)}{\cos \left(k_n h\right)} \cos \left(\frac{m\pi}{h}(z+h)\right) \, dz}, \\ { {\mathbf M}^{-}(m,n) = \int_{-h}^0 \frac{\cos \left(k^{\mu}_n (z+h)\right)}{\cos \left(k^{\mu}_n h\right)} \cos \left(\frac{m\pi}{h}(z+h)\right) \, dz } \end{matrix} }[/math]

and [math]\displaystyle{ {\mathbf N}^{+} }[/math] and [math]\displaystyle{ {\mathbf N}^{-} }[/math] are given by

[math]\displaystyle{ \begin{matrix} {\mathbf N}^{+}(m,n)= -\kappa_n{\mathbf M}^{+}(m,n),\\ {\mathbf N}^{-}(m,n)= \kappa^f_\mu{\mathbf M}^{-}(m,n). \end{matrix} }[/math]

[math]\displaystyle{ \mathbf{C} }[/math] is a [math]\displaystyle{ (M+1) }[/math] vector which is given by

[math]\displaystyle{ {\mathbf C}(m)=\int_{-h}^0 \frac{\cos (k_1(0)(z+h))}{\cos (k_1(0)h)} \cos \left(\frac{m\pi}{h}(z+h)\right)\, dz. }[/math]

The integrals in the above equation are each solved analytically. Now, for all but the first and [math]\displaystyle{ \Lambda }[/math]th plate, the edge equation becomes

[math]\displaystyle{ \begin{matrix} {\mathbf E}^{+}_{T_\mu} {\mathbf T}_\mu + {\mathbf E}^{+}_{R_\mu} {\mathbf R}_\mu = 0,\\ {\mathbf E}^{-}_{T_\mu} {\mathbf T}_\mu + {\mathbf E}^{-}_{R_\mu} {\mathbf R}_\mu = 0. \end{matrix} }[/math]

The first and last plates only require two equations, because each has only one plate edge. The equation for the first plate must be modified to include the effect of the incident wave. This gives us

[math]\displaystyle{ \begin{matrix} I \left( \begin{matrix} {\mathbf E}^{+}_{T_1}(1,0)\\ {\mathbf E}^{+}_{T_1}(2,0) \end{matrix} \right) + {\mathbf E}^{+}_{R_1} {\mathbf R}_1 = 0,\\ \end{matrix} }[/math]

and for the [math]\displaystyle{ \Lambda }[/math]th plate we have no reflection so

[math]\displaystyle{ \begin{matrix} {\mathbf E}^{-}_{T_\mu} {\mathbf T}_\mu = 0.\\ \end{matrix} }[/math]

[math]\displaystyle{ {\mathbf E}^{+}_{T_\mu} }[/math], [math]\displaystyle{ {\mathbf E}^{+}_{R_\mu} }[/math], [math]\displaystyle{ {\mathbf E}^{-}_{T_\mu} }[/math] and [math]\displaystyle{ {\mathbf E}^{-}_{R_\mu} }[/math] are 2 by M+3 matrices given by

[math]\displaystyle{ \begin{matrix} {\mathbf E}^{-}_{T_\mu}(1,n) = (\kappa_\mu(n)^2 - (2 - \nu)k_y^2)(k_{\mu}(n)\kappa_\mu(n)\tan{(k_{\mu}(n)h)}),\\ {\mathbf E}^{+}_{T_\mu}(1,n) = (\kappa_\mu(n)^2 - (2 - \nu)k_y^2)(k_{\mu}(n)\kappa_\mu(n)e^{-\kappa_\mu(n)(r_\mu - l_\mu)}\tan{(k_{\mu}(n)h)}),\\ {\mathbf E}^{-}_{R_\mu}(1,n) = (\kappa_\mu(n)^2 - (2 - \nu)k_y^2)(-k_{\mu}(n)\kappa_\mu(n)e^{\kappa_\mu(n)(l_\mu - r_\mu)}\tan{(k_{\mu}(n)h)}),\\ {\mathbf E}^{+}_{R_\mu}(1,n) = (\kappa_\mu(n)^2 - (2 - \nu)k_y^2)(-k_{\mu}(n)\kappa_\mu(n)\tan{(k_{\mu}(n)h)}),\\ {\mathbf E}^{-}_{T_\mu}(2,n) = (\kappa_\mu(n)^2 - \nu k_y^2)(-k_{\mu}(n)\tan{(k_{\mu}(n)h)}),\\ {\mathbf E}^{+}_{T_\mu}(2,n) = (\kappa_\mu(n)^2 - \nu k_y^2)(-k_{\mu}(n)e^{-\kappa_\mu(n)(r_\mu - l_\mu)}\tan{(k_{\mu}(n)h)}),\\ {\mathbf E}^{-}_{R_\mu}(2,n) = (\kappa_\mu(n)^2 - \nu k_y^2)(-k_{\mu}(n)e^{\kappa_\mu(n)(l_\mu - r_\mu)}\tan{(k_{\mu}(n)h)}),\\ {\mathbf E}^{+}_{R_\mu}(2,n) = (\kappa_\mu(n)^2 - \nu k_y^2)(-k_{\mu}(n)\tan{(k_{\mu}(n)h)}).\\ \end{matrix} }[/math]

Now, the matching matrix is a [math]\displaystyle{ (2M+6)\times(\Lambda-1) }[/math] by [math]\displaystyle{ (2M+1)\times(\Lambda -1) }[/math] matrix given by

[math]\displaystyle{ {\mathbf M} = \left( \begin{matrix} {\mathbf M}^{+}_{R_1} & -{\mathbf M}^{-}_{T_2} & -{\mathbf M}^{-}_{R_2} & 0 & 0 & & 0 & 0 & 0 \\ {\mathbf N}^{+}_{R_1} & -{\mathbf N}^{-}_{T_2} & -{\mathbf N}^{-}_{R_2} & 0 & 0 & & 0 & 0 & 0 \\ 0 & {\mathbf M}^{+}_{T_2} & {\mathbf M}^{+}_{R_2} & -{\mathbf M}^{-}_{T_3} & -{\mathbf M}^{-}_{R_3} & & 0 & 0 & 0 \\ 0 & {\mathbf N}^{+}_{T_2} & {\mathbf N}^{+}_{R_2} & -{\mathbf N}^{-}_{T_3} & -{\mathbf N}^{-}_{R_3} & & 0 & 0 & 0 \\ & & \vdots & & & \ddots & \\ 0 & 0 & 0 & 0 & 0 & & {\mathbf M}^{+}_{T_{\Lambda - 1}} & {\mathbf M}^{+}_{R_{\Lambda - 1}} & -{\mathbf M}^{-}_{ T_{\Lambda}} \\ 0 & 0 & 0 & 0 & 0 & & {\mathbf N}^{+}_{T_{\Lambda - 1}} & {\mathbf N}^{+}_{R_{\Lambda - 1}} & -{\mathbf N}^{-}_{T_{\Lambda }} \\ \end{matrix} \right), }[/math]

the edge matrix is a [math]\displaystyle{ (2M+6)\times(\Lambda-1) }[/math] by [math]\displaystyle{ 4(\Lambda-1) }[/math] matrix given by

[math]\displaystyle{ {\mathbf E} = \left( \begin{matrix} {\mathbf E}^{+}_{R_1} & 0 & 0 & 0 & 0 & & 0 & 0 & 0 \\ 0 & {\mathbf E}^{+}_{T_2} & {\mathbf E}^{+}_{R_2} & 0 & 0 & & 0 & 0 & 0 \\ 0 & {\mathbf E}^{-}_{T_2} & {\mathbf E}^{-}_{R_2} & 0 & 0 & & 0 & 0 & 0 \\ 0 & 0 & 0 & {\mathbf E}^{+}_{T_3} & {\mathbf E}^{+}_{R_3} & & 0 & 0 & 0 \\ 0 & 0 & 0 & {\mathbf E}^{-}_{T_3} & {\mathbf E}^{-}_{R_3} & & 0 & 0 & 0 \\ & & \vdots & & & \ddots & \\ 0 & 0 & 0 & 0 & 0 & & {\mathbf E}^{+}_{T_{\Lambda-1}} & {\mathbf E}^{+}_{R_{\Lambda-1}} & 0 \\ 0 & 0 & 0 & 0 & 0 & & {\mathbf E}^{-}_{T_{\Lambda-1}} & {\mathbf E}^{-}_{R_{\Lambda-1}} & 0 \\ 0 & 0 & 0 & 0 & 0 & & 0 & 0 & {\mathbf E}^{-}_{ T_\Lambda} \end{matrix} \right), }[/math]

and finally the complete system to be solved is given by

[math]\displaystyle{ \left( \begin{matrix} {\mathbf M}\\ {\mathbf E}\\ \end{matrix} \right) \times \left( \begin{matrix} {\mathbf R}_1\\ {\mathbf T}_2\\ {\mathbf R}_2\\ {\mathbf T}_3\\ {\mathbf R}_3\\ \vdots\\ {\mathbf T}_{\Lambda-1}\\ {\mathbf R}_{\Lambda-1}\\ {\mathbf T}_{\Lambda} \end{matrix} \right) = \left( \begin{matrix} -I{\mathbf C}\\ \kappa_{1}(0)I{\mathbf C}\\ 0\\ \vdots\\ -IE^{+}_{T_1}(1,0)\\ -IE^{+}_{T_1}(2,0)\\ 0\\ \vdots \end{matrix} \right). }[/math]

The final system of equations has size [math]\displaystyle{ (2M+6)\times (\Lambda - 1) }[/math] by [math]\displaystyle{ (2M+6)\times (\Lambda - 1) }[/math].