Electromagnetic Propagation in Multi-Mode Random Media
Electromagnetic Propagation in Multi-Mode Random Media
Table of Contents
Electromagnetic Propagation in
Multi-Mode Random Media
Author(s): Harrison E. Rowe
Series Editor(s): Kai Chang
ISBN: 0471200700 (Electronic) 0471110035 (Print)
Copyright © 1999 John Wiley & Sons, Inc.
Chapter 1: Introduction
1-3 Harrison E. Rowe
Summary PDF Full Text (Size: 42K)
Published Online: 30 Oct 2001
DOI 10.1002/0471200700.ch1
[Wiley InterScience Home Page] [Personal Home Page] [Journal Finder] [Book Finder] [Search Wiley InterScience] [Reference Works] [Help]
[Contact Us] [Logout]
Copyright © 1999-2001 by John Wiley & Sons, Inc. All rights reserved.
Electromagnetic Propagation
in Multi-Mode Random Media
Electromagnetic
Propagation
in Multi-Mode
Random Media
HARRISON E. ROWE
A WILEY-INTERSCIENCE PUBLICATION
JOHN WILEY & SONS, INC.
NEW YORK / CHICHESTER / WEINHEIM / BRISBANE / SINGAPORE / TORONTO
Designations used by companies to distinguish their products are often claimed as trademarks. In all
instances where John Wiley & Sons, Inc., is aware of a claim, the product names appear in initial
capital or ALL CAPITAL LETTERS. Readers, however, should contact the appropriate companies for
more complete information regarding trademarks and registration.
Copyright © 1999 by John Wiley & Sons, Inc. All rights reserved.
No part of this publication may be reproduced, stored in a retrieval system or transmitted in any form
or by any means, electronic or mechanical, including uploading, downloading, printing, decompiling,
recording or otherwise, except as permitted under Sections 107 or 108 of the 1976 United States
Copyright Act, without the prior written permission of the Publisher. Requests to the Publisher
for permission should be addressed to the Permissions Department, John Wiley & Sons, Inc.,
605 Third Avenue, New York, NY 10158-0012, (212) 850-6011, fax (212) 850-6008, E-Mail:
[email protected].
This publication is designed to provide accurate and authoritative information in regard to the subject
matter covered. It is sold with the understanding that the publisher is not engaged in rendering
professional services. If professional advice or other expert assistance is required, the services of a
competent professional person should be sought.
ISBN 0-471-20070-0
For more information about Wiley products, visit our web site at www.Wiley.com.
To Alicia,
Amy, Elizabeth, Edward, and Alison,
and to the memory of
Stephen O. Rice
Contents
1 Introduction 1
References 3
4 Examples—White Coupling 49
4.1 Introduction 49
vii
viii CONTENTS
10 Conclusion 161
References 162
x CONTENTS
References 214
Index 231
Electromagnetic Propagation
in Multi-Mode Random Media
Electromagnetic Propagation in Multi-Mode Random Media. Harrison E. Rowe
Copyright © 1999 John Wiley & Sons, Inc.
Print ISBN 0-471-11003-5; Electronic ISBN 0-471-20070-0
Shiban K. Koul
MICROWAVE DEVICES, CIRCUITS AND THEIR INTERACTION l Char/es A. lee and
G. Conrad Da/man
ADVANCES IN MICROSTRIP AND PRINTED ANTENNAS l Kai-Fong Lee and Wei Chen (eds.)
OPTICAL FILTER DESIGN AND ANALYSIS: A SIGNAL PROCESSING APPROACH l
C. K. Madsen and 1. H. Zhao
OPTOELECTRONIC PACKAGING l A. R. Mickelson, N. R. Basavanhaly, and Y. C. Lee feds.)
ANTENNAS FOR RADAR AND COMMUNICATIONS: A POLARIMETRIC APPROACH l
Harold Mott
INTEGRATED ACTIVE ANTENNAS AND SPATIAL POWER COMBINING l julio A. Navarro and
Kai Chang
FREQUENCY CONTROL OF SEMICONDUCTOR LASERS l Motoichi Ohtsu (ecf.)
SOLAR CELLS AND THEIR APPLICATIONS l Larry 0. Par&in (ed.)
ANALYSIS OF MULTICONDUCTOR TRANSMISSION LINES l Clayton R. Paul
INTRODUCTION TO ELECTROMAGNETIC COMPATIBILITY l Clayton R. Paul
INTRODUCTION TO HIGH-SPEED ELECTRONICS AND OPTOELECTRONICS .
Leonard M. Riaziat
NEW FRONTIERS IN MEDICAL DEVICE TECHNOLOGY l Arye Rosen and Hare/ Rosen (eds.)
ELECTROMAGNETIC PROPAGATION IN MULTI-MODE RANDOM MEDIA l Harrison E. Rowe
NONLINEAR OPTICS l E. C. Sauter
InP-BASED MATERIALS AND DEVICES: PHYSICS AND TECHNOLOGY l Osamu Wada
and Hideki Hasegawa (eds.)
FREQUENCY SELECTIVESURFACE AND GRID ARRAY l T. K. Wu (ed.)
ACTIVE AND QUASI-OPTICAL ARRAYS FOR SOLID-STATE POWER COMBINING l
Robert A. York and Zoya B. PopovZ (eds.)
OPTICAL SIGNAL PROCESSING, COMPUTING AND NEURAL NETWORKS l Francis T. S. Yu
and Suganda lutamulia
Electromagnetic Propagation in Multi-Mode Random Media. Harrison E. Rowe
Copyright © 1999 John Wiley & Sons, Inc.
Print ISBN 0-471-11003-5; Electronic ISBN 0-471-20070-0
232 INDEX
nondegenerate,83 propagation,seeRandompropagation
two modes,37-38 parameters
examples,60-63 separationbetweenguides,2
Powersand cross-powers,31 straightness,2,5
multi-mode, 33-34,50-52 Randompropagationparameters,88-9 1
two modes,32,93-95 correlationlength,89
Propagationparameters,6. Seealso
Randompropagationparameters Schelkunoff,5
multi-mode, 15 Square-wavecoupling:
attenuationandphase,18 coupling function, 130,215
constant,17,22, 100 coupling spectra:
example,68 binary independent,133,218
degenerate,16,206 binary Markov, 136137,2 18
two modes,6-7 five-level, band-passMarkov, 141,
attenuationandphase,6-7, 11-12, 224
22,24-26,38,40 six-level, low-passMarkov, 140,219
constant,9-10,22,30,49,100 examples:
degenerate,8,203 four modes:
six-level, low-passMarkov
Randomcoupling coefficients,seealso sections,141-142
Square-wavecoupling two modes:
almost-whitespectra,1OO-101 binary independentsections,
correlation length, 13, 101 132-134
Gaussianor Poisson,2 1 binary Markov sections,135, 138
generalspectra: inputs, 130
band-pass,107
low-pass,100 Transferfunction(s),23-24
white spectrum,2 1,3 1,49 covariance:
Randomoptical thickness,2. SeeaZso multi-mode,43
Multi-layer coatings two modes,42
Randomparameters: Transmissionstatistics,1,2, 6. Seealso
coupling,seeRandomcoupling Averagetransferfunctions;
coefficients Coupledpowerequations;Impulse
guide width, 2 response,statistics;Power
layer thickness,seeMulti-layer coatings fluctuations;Square-wave
optical thickness,seeMulti-layer coupling;Transferfunctions,
coatings covariance;Multi-layer coatings
Electromagnetic Propagation in Multi-Mode Random Media. Harrison E. Rowe
Copyright © 1999 John Wiley & Sons, Inc.
Print ISBN 0-471-11003-5; Electronic ISBN 0-471-20070-0
CHAPTER ONE
Introduction
REFERENCES
CHAPTER TWO
2.1. INTRODUCTION
In the lossless case powers of different modes add, and Pz repre-
sents the total power in the guide. More generally, Pz is simply the
sum of the mode powers computed individually. Substituting Equa-
tion (2.3) into (2.4), differentiating, and substituting Equation (2.1)
for the resulting derivatives, we obtain
dPz
= −2α0 zP0 z − 2α1 zP1 z (2.5)
dz
This result states that each mode contributes to the decrease in
Pz in proportion to the product of its attenuation constant and
the power it carries. For the lossless case, α0 = α1 = 0, Equation
(2.5) yields conservation of power.
Equations (2.1)–(2.2) approximate the response of a multi-mode
guide in which only two forward-traveling modes are significant, and
other modes may be neglected. We denote the signal (desired) mode
by I0 and the spurious (undesired) mode by I1 . For the ideal guide
(without imperfections) 0 and 1 are constant, independent of z,
and cz = 0. Since both modes travel in the forward +z direction
α0 , α1 , β0 , and β1 are all positive. We assume that the signal mode
has lower loss; α0 ≤ α1 . Random imperfections are modeled by re-
garding cz 0 z and 1 z as stationary random processes with
appropriate statistics. As a particular example, a two-mode guide
with random two-dimensional straightness deviations has constant
propagation parameters 0 and 1 and a coupling coefficient in-
versely proportional to the radius of curvature Rz of the guide
8 COUPLED LINE EQUATIONS
axis,
C
cz = (2.6)
Rz
the constant C being a function of the two modes, determined by
electromagnetic theory. In the special case of “single-mode” fiber,
I0 and I1 represent the two polarizations of the dominant mode and
0 = 1 .
The coupled line equations can describe a directional coupler
where backward modes are neglected. For this case, random im-
perfections are modeled by taking the propagation parameters to be
stationary random processes. As a simple example, cz = constant
and the propagation parameters have equal mean values, 1 z =
2 z .
Alternatively, Equations (2.1) and (2.2) can be applied to a trans-
mission medium with a reflected wave. Denoting the forward wave
by I0 and the backward wave by I1 , 0 = −1 = α + jβ with positive
α and β; c represents the reflection coefficient. However, we will not
use these relations in our later treatment of multi-layer devices, but
will rather employ direct methods.
Direct solution of the coupled line equations for arbitrary cz
and/or 0 z 1 z is not possible. The matrix methods we employ
to determine the transmission statistics for random guides require
the solutions described in the remainder of this chapter.
(2.8)
2.3. EXACT SOLUTIONS 9
In this case, only the integrals of z and cz matter, their de-
tailed functional behavior being unimportant. These results are a
consequence of the fact that the two modes have identical attenu-
ation and phase parameters; it doesn’t matter how the coupling is
distributed.
Equation (2.8) yields directly the wave matrix for delta-function
coupling. Set
0 +1
Tz = e− 2 z /K+ − K−
√ √ √ √
−K− e/2z + K+ e−/2z e/2z − e−/2z
· √ √ √ √
e/2z − e−/2z K+ e/2z − K− e−/2z
(2.14)
10 COUPLED LINE EQUATIONS
√
1±
K± = −j K+ K− = −1 (2.15)
2c0 /
1 c /
= j 0√ (2.16)
K+ − K−
√
= 1 − 2c0 /2 (2.17)
= 0 − 1 (2.18)
Let us approach a discrete coupler, Equation (2.9), by setting
c
c0 = (2.19)
z
with fixed c, and taking the limit as z → 0 (and c0 → ∞) in Equa-
tions (2.12)–(2.18). Then [6],
1 1
lim = (2.22)
z→0 K+ − K− 2
0 +1
lim e− 2 z
= 1 (2.23)
z→0
where
kz
ck = czdz (2.26)
k−1z
The kth section is illustrated in Figure 2.1, which shows cz and
cδ z; this approximation places all of the coupling at the end of the
section.
The corresponding matrix description becomes
−γ
I0 kz cos ck j sin ck e 0k 0 I0 k − 1z
= · ·
I1 kz j sin ck cos ck 0 e−γ1k I1 k − 1z
(2.27)
where ck is given by Equation (2.26) and the γk ’s give the total
attenuation and phase shift in the length z:
kz kz kz
γ0k = 0 zdz = α0 zdz + j β0 zdz
k−1z k−1z k−1z
(2.28)
kz kz kz
γ1k = 1 zdz = α1 zdz + j β1 zdz
k−1z k−1z k−1z
k z
area = c(z) dz = ck
(k – 1) z
c(z)
c (z)
ck (z – k z)
(k – 1) z k z
1 2π
z z (2.30)
α0 − α1 β0 − β1
2.5. PERTURBATION THEORY 13
where
z z
γ0 z = 0 xdx γ1 z = 1 xdx (2.32)
0 0
where
z
γz = xdx z = 0 z − 1 z (2.34)
0
14 COUPLED LINE EQUATIONS
In this case, two alternative forms for m00 and m11 are as follows:
z x
m00 = 1 − cxe+x dx cye−y dy
0 0
z z−ζ
=1− e+ζ dζ cxcx + ζdx
0 0
1 z z
=1− cxcye+x−y dxdy (2.39)
2 0 0
z x
m11 = 1 − cxe−x dx cye+y dy
0 0
z z−ζ
=1− e−ζ dζ cxcx + ζdx
0 0
1 z z
=1− cxcye−x−y dxdy (2.40)
2 0 0
2.6. MULTI-MODE COUPLED LINE EQUATIONS 15
Each of these three expressions for m00 and for m11 has its uses. The
first yields a physical interpretation of these perturbation results. The
expected value of the middle expressions introduces the covariance
of the coupling coefficient for guides with random coupling, and
constant propagation parameters.
Equations (2.35)–(2.40) yield useful approximations for small cou-
pling. They are the first terms of infinite series given in Appendix A.
These series converge rapidly when
z
cxdx 1 (2.41)
0
The above relations for two modes are readily extended to many for-
ward modes [2–6] using matrix notation as follows. Equations (2.1)
and (2.2), the coupled line equations, become
where
I T z = I0 z I1 z I2 z ··· (2.43)
0 z 0 0 ···
0 1 z 0 ···
z = (2.44)
0 0 2 z ···
0 C01 C02 ···
C01 0 C12 ···
C= C = CT = C∗ (2.45)
C02 C12 0 ···
Define
Pz = Pi z = I T z · I ∗ z (2.47)
i
Then,
dPz
= I T z · I ∗ z + I T z · I ∗ z (2.48)
dz
Substituting Equation (2.42) into (2.48), the multi-mode generaliza-
tion of Equation (2.5) becomes
dPz
= − I T z · z + ∗ z · I ∗ z = − 2 αi zPi z (2.49)
dz i
For the lossless case all αi = 0 and Equation (2.49) yields conser-
vation of power.
Exact solutions corresponding to those of Section 2.3 for two
modes are readily found for the multi-mode case. In the degenerate
case, all modes have identical propagation constants,
where the wave matrix for the multi-mode discrete coupler is given
by
S = ejcC (2.56)
and
kz
k = zdz (2.65)
k−1z
where
z
z = xdx (2.68)
0
This result holds for small coupling, i.e., from Equations (A.22) and
(A.24)
z
cxdxC 1 (2.71)
0
and
z
0 xdx 0
z =
0
z (2.73)
0 1 xdx
0
REFERENCES
CHAPTER THREE
3.1. INTRODUCTION
= 0 − 1 = α + jβ
(3.8)
α = α0 − α1 β = β0 − β1
+C0 f >0
C = C0 sgn f = C0 a positive constant. (3.18)
−C0 f <0
gα τ of Equations (3.16) and (3.17) is therefore the normalized im-
pulse response of the idealized dispersionless guide with frequency-
independent attenuation constants and coupling.
While Equation (3.23) is strictly true only for an idealized guide,
it will provide a useful approximation for a physical guide in many
cases. This approximation requires the frequency variation of α,
C, and of the group delay of the two modes to be small over the
frequency range of interest.
The general properties of G0 α β and gα τ are described in
Appendix B. These properties are reflected in the statistical results
to follow.
are satisfied:
1 2π
z z (3.25)
α β
By Equations (3.1) and (3.2), cz is white with zero mean; by Equa-
tion (3.3), the different sections remain independent for any z.
Therefore Equation (3.24) may be identified with Equation (D.8).
Define the average complex mode amplitudes as
I i z = Ii z G i z = Gi z (3.26)
by Equation (D.14), the expected value of Equation 3.24 yields
I 0 kz e−0 z cos ck e−1 z jsin ck I 0 k − 1z
= ·
I 1 kz e−0 z jsin ck e−1 z cos ck I 1 k − 1z
(3.27)
For sufficiently small z,1
1 1
cos ck ≈ 1 − ck2 = 1 − S0 z
2 2 (3.28)
sin ck ≈ ck = 0
Substituting Equation (3.28) into Equation (3.27) and taking the
limit as z → 0,
S
I 0 z = − 0 + 0 I 0 z
2
(3.29)
S0
I 1 z = − 1 + I 1 z
2
1
From Equation (2.26), since cz has zero mean,
kz
ck = cz dz = 0
k−1z
where I0 0 and I1 0 are, respectively, the sinusoidal signal and
spurious mode inputs. By Equation (3.6), the expected normalized
responses are
S0
G 0 z = G0 z = e− 2 z G0 0
S0
(3.31)
− z
G 1 z = G1 z = e 2 G 0
1
1
ejck C ≈ I + jck C − ck2 C2
2 (3.39)
−z
e ≈ I − z
where I is the unit matrix. The first factor on the right-hand side of
Equation (3.37) becomes
1 1
ejck C · e−z ≈ I − z − ck2 C2 = I − z − S0 C2 z (3.40)
2 2
The substitution
0 1
C= C2 = I (3.42)
1 0
3.3. AVERAGE TRANSFER FUNCTIONS 29
1 1 2 2
I 1 z = − S0 C02 C12 · I 0 z − 1 + S0 C01 + C12 I 1 z
2 2
1
− S0 C01 C02 · I 2 z (3.45)
2
1 1
I 2 z = − S0 C01 C12 · I 0 z − S0 C01 C02 · I 1 z
2 2
1 2 2
− 2 + S0 C02 + C12 I 2 z (3.46)
2
the spurious mode averages I 1 z = I1 z and I 2 z = I2 z
depart from zero for a coherent input signal.
We use coherent inputs to the coupled line equations to study
transfer function statistics, here the expected transfer functions.
For a coherent sinusoidal input with unit power to mode i we set
Ii 0 = 1, a deterministic quantity representing a complex wave ejωt .
Sources such as light-emitting diodes (LEDs) produce incoherent or
partially coherent inputs to multi-mode optical fibers. We model an
incoherent sinusoidal input of unit power with constant amplitude
and random phase as Ii 0 = ejθ , with θ a random variable uni-
formly distributed from 0 < θ < 2π, representing the complex wave
ejωt ejθ ; here the expected input is zero, I0 0 = 0 The present
results, Equations (3.29), (3.41), or (3.44)–(3.46), show that all ex-
pected mode amplitudes are zero for incoherent inputs; for coherent
or partially coherent inputs the expected mode amplitudes depart
from zero.
The present results simplify in the nondegenerate case, where no
two modes have the same propagation constants. We show in Section
4.6.1 that for large βij = βi − βj the coupling between expected
mode amplitudes approaches zero, as in the two-mode case.
The expected complex transfer functions give little information
about the transmission behavior of guides with random coupling. Ex-
act results for higher order statistics—average powers, power fluctua-
tions, transfer function vs. frequency statistics—are calculated below
using Kronecker products.
For two modes with constant propagation parameters, take the Kro-
necker product of Equation (3.24) with its complex conjugate, to
yield by Equation (C.7)
I0 kz2 I0 k − 1z 2
I0 kzI1∗ kz I k − 1z I1∗ k − 1z
= M ⊗ M · 0
∗
∗ ∗
I1 kzI0 kz I1 k − 1z I0 k − 1z
I1 kz2 I1 k − 1z 2
(3.48)
3.4. COUPLED POWER EQUATIONS 31
where C and are given by Equations (3.35) and (3.36) and I is the
N × N unit matrix. Equation (3.61) reduces to the two-mode case
by the substitution of Equation (3.42).
Even the two-mode case involves a disagreeable amount of hand
calculation, requiring the expected values of 16 terms of a 4 × 4
matrix. It is evident that greater numbers of modes require the use
of a computer algebra system capable of programmable symbolic
and numeric algebraic, matrix, and calculus operations. The N-mode
coupled power equations of Equation (3.61) comprise a set of N 2
differential equations; N of these are real, N 2 − N occur in complex
conjugate pairs.
The three-mode case of Equation (3.43) leads to a 9 × 9 matrix
with 81 terms, yielding 9 differential equations relating the average
powers and cross powers, each with 9 terms on the right-hand side.
These have been obtained with MAPLE; the first two of these 9
equations are given below:
2 2
P 0 z = − S0 C01 + C02 + 0 + ∗0 P 0 z − 5S0 C02 C12 P 01 z
2
− 5S0 C01 C12 P 02 z − 5S0 C02 C12 P 10 z + S0 C01 P 1 z
+ S0 C01 C02 P 12 z − 5S0 C01 C12 P 20 z + S0 C01 C02 P 21 z
2
+ S0 C02 P 2 z
P 01 z = −5S0 C02 C12 P 0 z
2 2 2
− S0 C01 + 5C12 + 5C02 + 0 + ∗1 P 01 z
2
− 5S0 C01 C02 P 02 z + S0 C01 P 10 z − 5S0 C02 C12 P 1 z
+ S0 C01 C12 P 12 z + S0 C01 C02 P 20 z − 5S0 C01 C12 P 21 z
+ S0 C02 C12 P 2 z
··· ··· ···················································
(3.62)
3.5. POWER FLUCTUATIONS 35
It is unlikely that anyone will wish to deal with even three modes
by hand. Subsequent examples given in Chapters 4, 6, 7, and 8
demonstrate that MAPLE is capable of treating significant multi-
mode problems.
The coupled power equations of this section apply equally well
to coherent, incoherent, or partially coherent inputs; their solutions
yield the mode powers as a function of distance z along the guide in
terms of the input mode powers, whatever the state of coherence at
the input. For coherent or partially coherent inputs, the modes will
be correlated and will contain coherent components given by the
results of Section 3.3; but the total mode powers and cross-powers
depend only on their input values.
Strict decoupling of the average powers and cross-powers occurs
only in the two-mode case. However, if all propagation constants
are distinct (i.e., no degenerate modes), the powers and the cross-
powers are approximately decoupled. In this approximation, the av-
erage mode powers are described by N differential equations [2],
rather than the N 2 differential equations of Equation (3.61). Fur-
ther discussion is given in Section 4.6.2.
P ijk% z = Pijk% z = Ii zIj∗ zIk zI%∗ z (3.64)
2
Pi z 2 = Pi2 z − Pi z = P iiii z − P 2i z (3.68)
P i z was found in Section 3.4; we determine P iiii z in this section,
and through Equation (3.68), the power fluctuation. From Equation
(3.34),
Define
− I ⊗ I ⊗ I ⊗ ∗ − I ⊗ I ⊗ ⊗ I
∗
− I ⊗ ⊗ I ⊗ I − ⊗ I ⊗ I ⊗ I · 4 z (3.72)
QT z = P 0000 z P 0011 z P 0101 z P 1010 z P 1111 z
(3.73)
38 GUIDES WITH WHITE RANDOM COUPLING
−4α0 0 0 0 0
0 −2α+ 0 0 0
N2 = 0 0 −2α+ + jβ 0 0
0 0 0 −2α+ − jβ 0
0 0 0 0 −4α1
(3.77)
where
α+ = α0 + α1 β = β0 − β1 (3.78)
3.6. TRANSFER FUNCTION STATISTICS 39
G0 β + σ
· (3.82)
G01 β + σ k−1z
where the + signs are taken when β and β + σ have the same
sign (i.e., both frequencies are positive, as in the usual case, or both
are negative), the − signs when they have opposite signs, in accord
with Equations (3.15) and (3.18).
The analysis follows that of Section 3.4:
R0 z G0 β + σG∗0 β
R01 z G0 β + σG∗01 β
= (3.83)
∗
R10 z G01 β + σG0 β
R1 z σ
G01 β + σG∗01 β z
The initial conditions are obtained from Equation (3.14) and (3.79):
The four equations again decouple into two sets of two equations,
in the two-mode case.
For the cross-moments, Equations (3.84) and (3.86) yield
Signal and spurious modes remain uncorrelated for all z, for zero
spurious mode input.
The transfer function covariances are defined in Equation (3.83)
as
The solutions yield exact expressions for the normalized signal and
signal-spurious mode transfer function covariances, for white cou-
pling, as follows:
with C and given by Equations (3.35) and (3.36). Define the trans-
fer function co- and cross-variances as
Ri zxy = Ii zx Ii∗ zy
(3.95)
Rij zxy = Ii zx Ij∗ zy = R∗ji zyx i = j
zxy = S0 C ⊗ C − 21 I ⊗ C2 − 21 C2 ⊗ I
− I ⊗ y∗ − x ⊗ I · zxy (3.97)
Consider a causal filter with real impulse response ɡt and trans-
fer function f :
∞
ɡt = f ej2πft df (3.98)
−∞
The expected value and covariance of the transfer function are in-
dependent of the frequency f . The spectral density of the transfer
function is −t, where t is the Fourier transform of ν:
∞
t = νej2πtν dν t = ∗ t (3.104)
−∞
Let the impulse response of the overall filter be t, with envelope
t = t + jt (3.106)
t is the Hilbert transform of t:
1 1 ∞ τ
t = + t = dτ (3.107)
πt π −∞ t − τ
3.7. IMPULSE RESPONSE STATISTICS 45
3.8. DISCUSSION
REFERENCES
CHAPTER FOUR
Examples—White Coupling
4.1. INTRODUCTION
49
50 EXAMPLES—WHITE COUPLING
1 0 < θ0 < 2π
wθ0 = 2π (4.4)
0 otherwise.
Then,
I0 0 = 0 I0 0 2 = 1 (4.5)
Ii 0 = 1 0 ≤ i ≤ N − 1 (4.10)
Then,
α0 = α1 = 0 (4.21)
Assume input only to the signal mode; i.e., from Equation (4.1) with
N = 2,
I1 0 = 0 (4.22)
The coupled power equations for the lossless case are, from Equa-
tions (3.53) and (4.21),
The initial conditions, from Equations (4.7), are the same for both
coherent and incoherent signal-mode inputs; for N = 2,
1 1 −2S0 z
P 0 z = + e
2 2
(4.26)
1 1
P 1 z = − e−2S0 z
2 2
The results are plotted in Figure 4.1.
The cross-powers are given by Equation (3.54); in the lossless case
they become
Therefore,
For signal only input the two modes are uncorrelated at the input,
and remain so for all z; this is true for both coherent and incoherent
signal mode inputs.
For a coherent signal input with zero spurious mode input, Equa-
tions (4.1) and (4.2), the signal mode contains a coherent component
given by Equation (4.23), plus a random component. The spurious
54 EXAMPLES—WHITE COUPLING
1.0
0.8
P0 (z)
0.6
0.4
P1 (z)
0.2
0
0 0.5 1 1.5 2 2.5 3
S0 z
FIGURE 4.1. Signal and spurious mode powers vs. distance, for coherent
or incoherent signal input.
1.0
0.8
P0 (z)
0.6
0.4
P0 (z) incoherent
0.2
P0 (z) coherent
0
0 0.5 1 1.5 2 2.5 3
S0 z
FIGURE 4.2. Coherent and incoherent signal mode powers vs. distance, for
coherent signal input.
2
β β −S0 1− β
S z
2
(4.35)
· 1− −1+j e 0
S0 S0
2
β β +S0 1− β
S z
2
+ 1− +1−j e 0
S0 S0
The following special cases are of interest. For equal propagation
constants,
P 01 z = P 10 z = 1 β = 0 (4.36)
this result is readily verified by elementary means. Next,
P 01 z = P ∗10 z = 1 + 1 − j · S0 ze−S0 z β = S0 (4.37)
4.2. COUPLED POWER EQUATIONS—TWO-MODE CASE 57
=0
S0
1
0.5
0.8
0.6
P01 (z)
5
P10 (z)
50 –
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3
S0 z
FIGURE 4.3. Cross-powers vs. distance, for coherent signal and spurious
mode inputs.
58 EXAMPLES—WHITE COUPLING
ers, the two cross-power second moments, and the spurious power
second moment.
z is given by the solution of Equation (3.75),
and N2 is given by Equation (3.77) and (3.78) for the lossless case as
0 0 0 0 0
0 0 0 0 0
N2 = 0 0 −j2β 0 0 (4.46)
0 0 0 +j2β 0
0 0 0 0 0
The signal and spurious mode power variances are given by Equation
(3.68):
where P 0 z and P 1 z are given in Section 4.2 for various inputs.
60 EXAMPLES—WHITE COUPLING
=0
S0
0.35
0.30
2
0.25
10 –
0.20
P0, 1 (z)
0.10
0.15
0.0
0
0 0.5 1 1.5 2 2.5 3
S0 z
FIGURE 4.4. Standard deviations of signal and spurious mode powers vs.
distance, with differential propagation constant as a parameter, for coher-
ent or incoherent signal input.
100 –
0.6
0.5
10
0.4 2
P0, 1 (z)
0.3
=1
S0
0.2
0.1
=0
S0
0
0 0.5 1 1.5 2 2.5 3
S0 z
FIGURE 4.5. Standard deviations of signal and spurious mode powers vs.
distance, with differential propagation constant as a parameter, for coher-
ent signal and spurious mode inputs.
4.3.4. Discussion
The power fluctuations in the signal and spurious modes are equal
in each of the three cases of Figures 4.4–4.6. This is required by
conservation of power, Equation (2.5) in the lossless case. For large
z, the mode powers have standard deviations comparable to their
average values.
4.3. POWER FLUCTUATIONS—TWO-MODE GUIDE 63
0.6
20 –
2
1
0.5
=0
S0
0.4
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3
S0 z
FIGURE 4.6. Standard deviations of signal and spurious mode powers vs.
distance, with differential propagation constant as a parameter, for inco-
herent signal and spurious mode inputs.
where
S z · e−S0 z 1 − τ I 2S z τ1 − τ
0 ≤ τ < 1
0
P0ac τ = τ 1 0
0 otherwise.
(4.56)
1.3
1.0
S0 z = 10
3
P0ac ( )
0.5
0.5
0
0 0.5 1.0
As an example, for S0 z = 628, P0ac τ has a peak value P0ac 05 ≈ 10
and width ≈ 0032 at the 368% points, where P0ac 05 ± 0016 ≈
368; the coherent component has become insignificant.
As S0 increases, more terms of the summation in Equation (B.1)
become important. The partial impulse response corresponding to
Equation (B.2) is given in Equation (E.33), which must be integrated
over all possible sets of 2n transitions; for the great majority of these
the signal will have spent about half of the time traveling in the signal
and in the spurious modes. Consequently, the output will be a narrow
pulse containing half the input power, with delay corresponding to
the average velocity of signal and spurious modes. This behavior is
described by the results of Equations (4.58) and (4.59).
We apply the results of Sections 3.3 and 3.4, for the average transfer
functions and for the coupled power equations, to four lossless cou-
pled modes. The large number of parameters necessary—6 coupling
coefficients, 4 propagation constants, and the spectral density of the
geometric parameter—do not allow a compact treatment of the gen-
eral case, as was done for the two-mode case. Consequently, we con-
sider as a specific example TE modes in a symmetric slab waveguide,
with particular numerical parameters for dielectric constants, thick-
ness, and frequency/wavelength given in Appendix F. The average
transfer functions are given by Equation (3.41) as
1
z = − − S0 C2 · z z = Iz (4.60)
2
The dielectric slab is illustrated in Figures F.1 and F.3, with index
1.01, thickness 123 · 10−5 m, and wavelength 10−6 m, given in Equa-
tions (F.14) and (F.17). The propagation and coupling matrices are
68 EXAMPLES—WHITE COUPLING
too small to appear in Figure 4.8. The average signal transfer func-
tion appears to decay exponentially, similar to the two-mode case.
For incoherent signal input, all four modes have zero coherent
component.
70 EXAMPLES—WHITE COUPLING
0.8
0.6
I0 (z)
0.4
0.2
0
0 50 100 150 200
z - meters
FIGURE 4.8. Average signal transfer function magnitude, for coherent
single-mode input.
The average mode powers for the four modes are shown in Fig-
ure 4.9.
4.5. COUPLED POWER EQUATIONS—FOUR-MODE CASE 71
1.0
0.8
P0 (z)
0.6
0.4
P3 (z)
0.2
P2 (z)
P1 (z)
0
0 20 40 60 80 100
z - meters
FIGURE 4.9. Mode powers vs. distance, for coherent or incoherent signal
input.
For large z, the total power divides equally between all four
modes. If the input signal is coherent, P 0 z contains a coherent
component I0 z, whose magnitude is shown in Figure 4.8, with
power given by I0 z 2 . Modes 1 and 3 have zero coherent com-
ponent; the coherent component of mode 2 is too small to appear
on the scale of Figure 4.9.
72 EXAMPLES—WHITE COUPLING
The magnitudes of the expected values of the four modes are shown
in Figure 4.10.
Each mode appears to decay exponentially, with its own decay
constant. Further study shows that the coupling between all four
expected mode amplitudes is weak, as in the preceding section.
The initial conditions for the coupled power equations, Equation
(4.61), for four-mode coherent input are
T 0 = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1!
16 ones
(4.74)
The average powers, P 0 z P 1 z P 2 z and P 3 z, depart only
infinitesimally from their initial values of 1, by less than ±4 × 10−6
over the range 0 < z < 30, and imperceptibly for greater z.
Each contains a coherent component, given by the square mag-
nitude of the corresponding average mode amplitude, shown on
Figure 4.10.
T 0 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1
(4.76)
0.8
0.6
0.4
I0 (z)
I1 (z)
0.2
I3 (z)
I2 (z)
0
0 50 100 150 200
z - meters
FIGURE 4.10. Average mode amplitudes, for four-mode coherent input.
all four modes have zero coherent component. Equations (4.76) and
(4.61) yield
For large z the average mode powers are equal, with zero coherent
components. Since the initial conditions satisfy this asymptotic be-
havior, the coupled power equations yield no further change in the
average powers.
74 EXAMPLES—WHITE COUPLING
1 " #
Gz = − S0 e+z · C2 · e−z · Gz (4.80)
2
$
N−1 $
N−1
C2 ij = Ci−1k Ckj−1 = Ci−1k Cj−1k 1 ≤ i j ≤ N (4.81)
k=0 k=0
$
N−1 $
N−1
C2 ii = Ci−1k Cki−1 = 2
Ci−1k 1 ≤ i ≤ N (4.82)
k=0 k=0
" # $
N−1
e+z · C2 · e−z ii
= C2 ii = 2
Ci−1k 1 ≤ i ≤ N (4.84)
k=0
3
The matrix elements are numbered from 1 1 to N N, while the nonzero coupling
coefficients in Equation 2.45, with unequal subscripts, are numbered from 0 1 to
N − 2 N − 1.
76 EXAMPLES—WHITE COUPLING
"i−1j−1 z
$
N−1
=e Ci−1k Ckj−1 1 ≤ i j ≤ N (4.85)
k=0
where
In this case, the exponential factor ejβ02 z = ejβ0 −β2 z has approxi-
mately 106 periods along the horizontal axis of Figure 4.8. Therefore
the asymptotic approximation of Equation (4.89) is a good approxi-
mation to the exact solution shown in Figure 4.8, obtained by solving
4.6. NONDEGENERATE CASE—APPROXIMATE RESULTS 77
Rij z = R∗ji z = Rij z Rij z = Gi zG∗j z (4.94)
T
&
z = R0 R01 · · · R0N−1 R10 R1 · · · R1N−1 ···
'
··· RN−10 RN−11 · · · RN−1 (4.96)
N as follows:
T
& T T T
'
z = 1 z 2 z ··· N z (4.106)
T
&
i z = Ri−10 z Ri−11 z ··· Rii−2 z
'
Ri−1 z Ri−1i z ··· Ri−1N−1 z (4.107)
$
N
1 1
i z = S0 M1 ik − M2 ik − M3 ik · k z 1 ≤ i ≤ N
k=1
2 2
(4.109)
Appendix H gives the asymptotic form for the M matrices and their
submatrices.
First, substitute Equations (H.5)–(H.7) into the ith row of Equa-
tion (4.109):
$
N
2
Ri−1 z∞ = S0 Ci−1k−1 e2αi−1k−1 z Rk−1 z∞ − Ri−1 z∞
k=1
1 ≤ i ≤ N (4.110)
Using Equations (4.79), (4.93), and (3.50), and relabeling the sum-
mation index in Equation (4.110), we have the final desired result:
$
N−1
2
$
N−1
2
P i z∞ = − 2αi + S0 Cik P i z∞ + S0 Cik P k z∞
k=0 k=0
0 ≤ i ≤ N − 1 (4.111)
80 EXAMPLES—WHITE COUPLING
S0 $
N "
2 2
#
Ri−1j−1 z∞ = − Ci−1k−1 + Cj−1k−1 Ri−1j−1 z∞ (4.114)
2 k=1
4.6. NONDEGENERATE CASE—APPROXIMATE RESULTS 81
For the lossless two-mode case, Equation (3.42) yields the following
analytic result for the β/S0 = ∞ curve of Figure (4.3):
where
" ∗ ∗ #
4 z = ez ⊗ e z ⊗ ez ⊗ e z · z (4.118)
(4.121)
where
with solutions
1 1 −6z
P 0 z = P 1 z = − e (4.127)
3 3
Equations (4.125) and (4.127) yield the β/S0 = ∞ curves for Fig-
ures 4.4 and for Figures 4.5 and 4.6, respectively.
84 EXAMPLES—WHITE COUPLING
4.6.4. Discussion
The asymptotic approximation provides major simplification in the
solutions for transfer function statistics of multi-mode transmission
systems with white coupling. In the two-mode case, the differential
equations for power fluctuations are reduced in number from 16 to
3, and provide a simple closed-form solution.
The physical example of Appendix F, a dielectric slab waveguide
supporting four modes, suggests that this approximation will apply
to many practical guides. Departures from these approximate results
can occur only with strict degeneracy between one or more pairs of
modes, i.e. with mode pairs having identical propagation constants.
Under these special conditions additional non-zero terms in the var-
ious matrices are readily included, with corresponding modifications
to the corresponding differential equations.
4.7. DISCUSSION
REFERENCES
CHAPTER FIVE
5.1. INTRODUCTION
I0 (z)
I1 (z)
β0 z = β + b0 z
(5.6)
β1 z = β + b1 z
5.2. STATISTICAL MODEL 89
cz = c (5.10)
1
λb = (5.12)
νb
90 DIRECTIONAL COUPLER WITH WHITE PROPAGATION PARAMETERS
1. Independent sections.
1
z λb = (5.13)
νb
3. Small coupling.
c z 1 (5.16)
2B0 νb β2 (5.17)
νb β (5.18)
1
c νb (5.19)
z
Equations (5.13) and (5.15) yield
1
2B0 νb (5.20)
z
Equations (5.18) and (5.20) yield Equation (5.17).
5.3. AVERAGE TRANSFER FUNCTIONS 91
z and
We have defined for convenience the auxiliary quantities G 0
G1 z as follows:
I0 z z
G 0
= e−jβz (5.23)
I1 z z
G 1
From Equations (D.8) and (D.14), for small z the difference equa-
tions of Equation (5.25) are approximated by the differential
equations
B0
G 0 z −2 jc G 0 z
= · (5.26)
B0
G 1 z jc − G 1 z
2
Assume the initial conditions of Equation (5.4) and Figure 5.1:
I 0 0 =
G 0 0 = 1 I 1 0 =
G 1 0 = 0 (5.27)
B0 > 0, I0 π/2c = 0, but this of course does not imply com-
plete power transfer, since the angles I0 π/2c and I1 π/2c
are random variables.
For zero coupling, c = 0, Equation (5.29) becomes
I 0 z 1
= e− jβ+B0 /2 z (5.30)
I 1 z 0
For this case the solution to Equations (5.1), (5.2), (5.4), and (5.6) is
−j z b xdx
I0 z e 0 0
= e−jβz (5.31)
I1 z 0
with
I0 z z
G 0
Iz =
Gz = (5.34)
I1 z z
G 1
1 1 −2cz
P 0 z = 1 − P 1 z = + e 1 + 2cz B0 = 4c (5.39)
2 2
P 01 z = P ∗10 z = −jcze−2cz B0 = 4c (5.40)
5.5. DISCUSSION
B0
=
c
1.0
10
0.8
0.6
P0 (z) = 1 – P1 (z)
1
0.4
0.2
0
0 1 2 3 4 5
2 · cz
FIGURE 5.2. Average power vs. distance for excited mode of symmetric
coupler.
good approximation for the almost white case. The analysis uses
Kronecker product methods, similar to those of Chapter 3 for white
coupling and constant propagation parameters. For two modes with
white coupling, the average powers were decoupled from the cross-
powers. Here the cross-powers are essential in the analysis.
These results are appropriate for investigating the necessary tol-
erances in two-branch couplers or switches, for which B0 /c 1, and
REFERENCES 97
REFERENCES
CHAPTER SIX
6.1. INTRODUCTION
99
100 GUIDES WITH GENERAL COUPLING SPECTRA
S( )
6.2.2. N Modes
Equations (2.66) for constant propagation parameters, and Equa-
tions (6.3) and (6.4) yield
1 1 1 2π
z z 0 ≤ i j ≤ N − 1
νc αij νc βij
(6.20)
C = L · · L−1 (6.21)
where
λ0
0
λ1
= (6.22)
0 λN−1
Then,
ck C = L · ck · L −1 (6.23)
6.2. ALMOST-WHITE COUPLING SPECTRA 105
ejck λ0
0
ejck λ1
−1
ejck C =L· ·L (6.24)
0 ejck λN−1
ejck C
1 − 21 λ20 S0z
0
1 − 21 λ21 S0z
−1
=L· ·L
0 1 2
1 − 2 λN−1 S0z
λ20 0
1 λ21
−1
= I − S0zL · ·L
2
0 λ2N−1
1
= I − S0zC2 S0zλ2i max 1 (6.26)
2
106 GUIDES WITH GENERAL COUPLING SPECTRA
S( )
2 c
where
z z−ζ
H0 = − e+ζ dζ cxcx + ζdx
0 z−z
z
H01 = cxe+x dx
z−z
z (6.31)
H10 = cxe−x dx
z−z
z z−ζ
H1 = − e−ζ dζ cxcx + ζdx
0 z−z
For N modes,
where
z
G =I +j cxex · C · e−x dx
z−z
z x
− cxex · C · e−x dx cyey · C · e−y dy (6.37)
z−z z−z
1 ≤ i j ≤ N i = j (6.38)
N z
Gii = 1 − Ci−1k−1 Ck−1i−1 ei−1k−1 ζ dζ
k=1 0
z−ζ
· cxcx + ζdx 1 ≤ i ≤ N (6.39)
z−z
where C denotes the matrix norm, Equation (A.25). The restric-
tion corresponding to Equation (6.33) becomes
∞
SνdνzC 1 (6.41)
−∞
6.3. GENERAL COUPLING SPECTRA—LOSSLESS CASE 109
αij = 0 0 ≤ i j ≤ N − 1 (6.42)
where we use the notation of Equation (3.26) for the average com-
plex mode amplitudes, I i z = Ii z. By Equation (6.4), Rc ζ is
narrow compared to 1 − ζ/z, and we have approximately
−jβ0 z 1 β β
I 0 z = e 1− S + jS z I 0 z − z (6.45)
2 2π 2π
−jβ1 z 1 β β
I 1 z = e 1− S − jS z I 1 z − z (6.46)
2 2π 2π
Sν is the Hilbert transform of the coupling spectrum Sν of Equa-
tion (6.2),
1 1 1 ∞ Sµ
Sν = Sν = dµ (6.47)
π ν π −∞ ν − µ
where
z ζ
B11 = B44 = 1 − z Rc ζ 1 − e+jβζ dζ
−z z
z (6.51)
ζ
B14 = B41 = z Rc ζ 1 − e+jβζ dζ
−z z
The average powers P 0 , P 1 are uncoupled from the cross-powers
P 01 P 10 ; we confine our current interest to the average powers.
Equation (6.4) permits neglect of the factor 1 − ζ /z in Equa-
tion (6.51), and the approximation of Appendix I yields the coupled
power equations for two modes:
β β
P 0 z −S 2π S
2π P 0 z
= · (6.52)
P 1 z β β P 1 z
S −S
2π 2π
6.3.2. N Modes
We retain the assumption of Equation (6.4) requiring guide sections
to be approximately independent:
z λc (6.53)
6.3. GENERAL COUPLING SPECTRA—LOSSLESS CASE 111
N z ζ
Gii = 1 − Ci−1k−1 Ck−1i−1 z Rc ζ 1 − ejβi−1k−1 ζ dζ
k=1 0 z
1 ≤ i ≤ N (6.57)
Gij∞ = 0 i = j (6.59)
Changing notation,
with elements
N
1 βi−1k−1 βi−1k−1
Hi = 2
Ci−1k−1 S
+ jS (6.63)
k=1
2 2π 2π
where
M = G ⊗ G∗ (6.70)
M ij = Gij G∗ (6.72)
N
i z = M ij · j z − z (6.73)
j=1
1
Gij represents the ij element of the matrix G, previously denoted by the slightly
different notation of Equation (6.38).
114 GUIDES WITH GENERAL COUPLING SPECTRA
N
Ri−1 z∞ = M ii ii ∞ Ri−1 z − z∞ + M ij ij ∞ Rj−1 z − z∞
j=1
j=i
(6.81)
Substituting Equations (6.76) and (6.80),
N
2
βi−1k−1
Ri−1 z∞ = 1 − Ci−1k−1 S z Ri−1 z − z∞
k=1
2π
N
2 βi−1j−1
+ Ci−1j−1 S z Rj−1 z − z∞ (6.82)
j=1
2π
where
∞
S0 ν = Rc ζe− α ζ e−j2πνζ dζ (6.87)
−∞
∞
S1 ν = Rc ζe+ α ζ e+j2πνζ dζ (6.88)
−∞
and we assume that the signal mode has lower loss, i.e., α = α0 −
α1 < 0. Note that S0 ν = S1 ν, in contrast to the lossless case.
Rc ζ must decay more rapidly than e+ α ζ grows for large ζ , in
order for Equation (6.88) to hold.
For the lossy coupled power equations, Equation (6.51) generalize
as follows. For the diagonal elements,
z
−2α0 z − α ζ ζ +jβζ
B11 = e 1 − z Rc ζe 1− e dζ
−z z
z (6.89)
−2α1 z + α ζ ζ +jβζ
B44 = e 1 − z Rc ζe 1− e dζ
−z z
Since Rc ζ is narrow by Equation (6.4), we have the usual approx-
imations, neglecting the factor 1 − ζ /z in the integrands of
Equation (6.89):
−2α0 z β
B11 = e 1 − S0 z
2π
(6.90)
−2α1 z β
B44 = e 1 − S1 z
2π
where S0 and S1 are given by Equations (6.87) and (6.88). Again,
the product Rc ζe+ α ζ must decay for large ζ in order for this
result for B44 to remain valid.
6.4. GENERAL COUPLING SPECTRA—LOSSY CASE 117
z 1 − e−2αz1− z
ζ
1 − e−2αz1− z 1 − e+2αz1− z
ζ ζ
≈
2αz −2αz
ζ
≈ 1− α z 1 (6.92)
z
−2α1 zβ
B14 = e S1 z
2π
α z 1 (6.93)
−2α0 z β
B41 = e S0 z
2π
6.4.2. N Modes
The lossy two-mode results of Section 6.4.1 are readily extended to
N modes in the nondegenerate case. The lossless result of Equation
(6.66) for the average transfer functions becomes
where
∞
Sik ν = Rc ζeαik ζ ej2πνζ dζ (6.97)
−∞
becomes
2
M ij ij ∞ = Ci−1j−1 · e2αi−1j−1 z z
z 1 − e−2αi−1j−1 z1− 2π
ζ
−αi−1j−1 ζ jβi−1j−1 ζ
· Rc ζe e dζ
−z 2αi−1j−1 z
i = j (6.98)
2 2αi−1j−1 z
βi−1j−1
M ij ij ∞ = Ci−1j−1 ·e Sj−1i−1 z
2π
αi−1j−1 z 1 (6.99)
N
2
M ii ii ∞ = 1 − z Ci−1k−1
k=1
z
αi−1k−1 ζ ζ
· Rc ζe 1− ejβi−1k−1 ζ dζ (6.100)
−z z
N
2
βi−1k−1
M ii ii ∞ = 1 − z Ci−1k−1 Si−1k−1 (6.101)
k=1
2π
These results do not require small loss; for the low-loss case
Si−1k−1 is again approximated by S of Equation (6.2).
120 GUIDES WITH GENERAL COUPLING SPECTRA
2
N−1
βik
P i z∞ = 1 − Cik Sik z e−2αi z P i z − z∞
k=0
2π
N−1 z
2
+ Cik z Rc ζe−αik ζ ejβik ζ
k=0 −z
1 − e−2αik z1− 2π
ζ
−2αk z
· dζ e P k z − z∞ (6.103)
2αik z
αik z 1 0 ≤ i k ≤ N − 1 (6.105)
6.5. DISCUSSION
Rc ζ = 1 − ζ ζ ≤ 1 (6.106)
α = 1 (6.107)
REFERENCES
CHAPTER SEVEN
7.1. INTRODUCTION
2/ζc
Sν = Rc 0 (7.2)
1/ζc + 2πν2
2
4πν
Sν = Rc 0 (7.3)
1/ζc + 2πν2
2
123
124 FOUR-MODE GUIDE WITH EXPONENTIAL COUPLING COVARIANCE
z = 0005 m. (7.9)
We take
Rc 0 = 5 × 10−10 m., Rc 0 = 25 × 10−19 m2 (7.13)
625 × 10−15
Sν = (7.15)
15625 × 10+8 + 2πν2
πν × 10−18
Sν = (7.16)
15625 × 10+8 + 2πν2
1 3 2 βik
I i z∞ = Ii 0 exp − z C S 0 ≤ i ≤ 3 (7.17)
2 k=0 ik 2π
Finally, the angle of Equation (6.66) for real (zero phase) inputs is
1 3 2 βik
I i z∞ = −z βi + C S (7.19)
2 k=0 ik 2π
7.4. DISCUSSION
The results of the present chapter illustrate the coupled power equa-
tions of Chapter 6 for the particular dielectric slab waveguide of
Appendix F with straightness deviation spectrum of Equation (7.15).
Other results of prior chapters for white imperfections may be simi-
larly extended to general spectra.
Such results are only approximate; they require the imperfections
to be small. For the present problem this requirement is given in
Equation (7.12). However, such limits are only heuristic; no rigorous
bounds on errors are available.
Exact transmission statistics for non-white coupling or propaga-
tion parameters are available only for random square-wave imper-
fections. Chapter 8 derives results for random square-wave coupling.
Comparison with these exact results permits validation of the ap-
proximate transmission statistics of Chapter 6.
Electromagnetic Propagation in Multi-Mode Random Media. Harrison E. Rowe
Copyright © 1999 John Wiley & Sons, Inc.
Print ISBN 0-471-11003-5; Electronic ISBN 0-471-20070-0
CHAPTER EIGHT
Random Square-Wave
Coupling
8.1. INTRODUCTION
c(z)
c1
c2
c4
z
2 3 4
c3
1 c /β
= k√ (8.6)
K+ − K−
√
= 1 + 2ck /β2 (8.7)
β = β0 − β1 (8.8)
0 23225 0 44966
23225
9
0 68737 0 −2
C = 10 m (8.12)
0 68737 0 133083
44966 0 133083 0
Equations (J.21) and (J.22) give the coupling statistics for indepen-
dent, symmetric, binary coupling. The coupling spectrum is given by
Equation (J.24), and plotted in Figure 8.2.
The small-coupling condition of Equations (6.32) and (6.33)
becomes
a 1 (8.15)
Equation (6.52) yields the following approximation for the signal
power with input given by Equation (8.2):
1 1 −2S β z
P 0app z = + e 2π (8.16)
2 2
For n sections, total guide length z = n, Equations (8.3), (D.14),
and (D.17) yield the following exact solutions for the inputs of Equa-
tion (8.2):
I 0 n I0 n 1
= =Tck n · coherent input.
I 1 n I1 n 0
(8.17)
I 0 n I0 n 0
= = incoherent input. (8.18)
I 1 n I1 n 0
P 0 n 1
P 01 n 0
= Tck ⊗ T∗ ck n ·
P 10 n 0
P 1 n 0
coherent or incoherent input, (8.19)
8.2. TWO MODES—BINARY INDEPENDENT SECTIONS 133
1.0
0.8
0.6
1
S( )
a2
0.4
0.2
0
0 0.5 1 1.5 2
FIGURE 8.2. Coupling spectrum: binary independent square-wave cou-
pling.
1.0
0.8
approximate
0.6 exact
P0
0.4
0.2
0
0 1 2 3 4
z/
FIGURE 8.3. Signal power vs. normalized distance; two modes; binary in-
dependent square-wave coupling. a = 05, β = 02π.
8.3. TWO MODES—BINARY MARKOV SECTIONS 135
1 1 −3242·a2 z
P 0app z = + e (8.23)
2 2
The exact coupled power equations for unit signal input, Equation
(8.2), are given from Equations (D.58) and (D.59):
P 0 n
n−1
P 01 n 1 02A1 ⊗ A1∗ 08A1 ⊗ A1∗
= I I ·
P 10 n 2 08A2 ⊗ A2∗ 02A2 ⊗ A2∗
P 1 n
1
A1 ⊗ A1∗ 0
· ∗
·
(8.24)
A2 ⊗ A2 0
0
136 RANDOM SQUARE-WAVE COUPLING
2
1
S( )
a2
0
0 0.5 1 1.5 2
FIGURE 8.4. Coupling spectrum: binary Markov square-wave coupling.
p = 02
1.6
1.4
1.2
1.0
1
S( )
a2
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2
FIGURE 8.5. Coupling spectrum: binary Markov square-wave coupling.
p = 08
exact
0.8
approximate
0.6
P0
0.4
0.2
0
0 2 4 6 8 10 12 14
z/
FIGURE 8.6. Signal power vs. normalized distance; two modes; binary
Markov square-wave coupling. a = 02, β = π.
We choose
d = 10−6 m, (8.30)
= 2 × 10−5 m. (8.31)
Then the spectral width of Sν in Figure 8.7 is 104 m−1 , roughly com-
parable to that of the low-pass example in Chapter 4, in Equation
(4.66). However, the low-pass spectral density, S0 = 5 × 10−16 m3 ,
is much greater than that of Equation (4.67). The root mean square
(rms) straightness deviation is
∞
c 2 z = Sνdν = 2236 × 10−6 m, (8.32)
−∞
also much larger than that of Equation (4.68), but still smaller than
the guide width of 123×10−5 m from Equation (F.17). The left-hand
140 RANDOM SQUARE-WAVE COUPLING
25
20
15
1
S( )
d 2
10
0
0 0.2 0.4 0.6 0.8 1
FIGURE 8.7. Coupling spectrum: six-level Markov square-wave coupling.
10
6
1
S( )
d 2
0
0 0.2 0.4 0.6 0.8 1
FIGURE 8.8. Coupling spectrum: five-level Markov square-wave coupling.
Equation (8.10) as
Ai = Tai = ej−+dai C (8.35)
where the propagation and coupling matrices and C are given by
Equations (8.11) and (8.12). These approximate and exact results
are plotted in Figure 8.9.
The small-coupling condition of Equation (6.41) is significantly
violated by Equation (8.33); the difference between the two results
in Figure 8.9 is not surprising.
0.8
approximate
0.6
I0 (z)
0.4
exact
0.2
0
0 2 ×10–4 4×10–4 6×10–4
z - meters
FIGURE 8.9. Coherent signal component vs. distance; four modes; six-level
Markov square-wave coupling.
8.5. DISCUSSION 143
8.5. DISCUSSION
Exact results for the coupled line equations with random coupling or
propagation parameters are available for the limiting case of white
imperfection spectra. General coupling spectra, treated in Chapter
6, require the use of both Kronecker product methods, to find the
overall transmission statistics of cascaded sections, and perturbation
theory, to determine the performance of a typical section. These
results are therefore only approximate; heuristic conditions for small
coupling are given in Equations (6.33) and (6.41) for two modes and
N modes, respectively.
The present chapter has extended Kronecker product methods to
Markov matrices, permitting the treatment of correlated sections.
The choice of square-wave coupling permits exact treatment of the
individual sections, and therefore yields exact transmission statistics
for square-wave coupling with a variety of spectra. These results of-
fer partial validation of the intuitive small-coupling conditions of
Section 6.3.
The present results suggest the following physical interpretation
of the small-coupling requirement: perturbation theory must apply
in a correlation length of the coupling coefficient cz.
The final example, of Section 8.4, represents an approach to a
more realistic model of a random coupling function. The example
chosen requires numerical computing with quite large matrices—
128 × 128 in these calculations for the expected transfer functions.
These results are within the capabilities of MAPLE on a Pentium
machine running at 233 MHz. Significant reduction in computing
time is achieved by taking advantage of the special properties of the
matrices in this particular example.
Higher order statistics require correspondingly increased comput-
ing resources. Determining the mode powers would require matrix
products 512 × 512 in size.
Finally, our choice of square-wave coupling for the individual sec-
tions permits exact results. However, the Markov–Kronecker product
methods may be applied to other coupling functions for the individ-
ual sections, greatly enlarging the scope of these methods. The in-
dividual sections would no longer permit exact treatment, as in the
present case, but would have to be solved by perturbation theory
or other numerical methods.
Electromagnetic Propagation in Multi-Mode Random Media. Harrison E. Rowe
Copyright © 1999 John Wiley & Sons, Inc.
Print ISBN 0-471-11003-5; Electronic ISBN 0-471-20070-0
CHAPTER NINE
9.1. INTRODUCTION
Ei
d1 d2 dk d
Et
ni n1 n2 ... nk ... n nt
Er
0 z1 z2 zk – 1 zk z – 1 z
R + T = 1 (9.3)
λ is the free-space wavelength; the n’s and d’s are shown in Fig-
ure 9.1.
The circuit-theoretic model of a dominant-mode transmission line
with random parameters is closely related to the present problem
[7]. In both, input and output are related by matrix products, per-
mitting Kronecker product methods to be applied. However, unlike
the forward-mode case of preceding chapters, the outputs are known
while the inputs are not. Therefore direct analysis yields the statistics
for only the reciprocal transfer function 1/T . We desire statistics for
T and for R, corresponding to mixed boundary conditions at input
and output, i.e., a unit forward wave at the input and zero backward
wave at the output.
A heuristic solution for the average powers in two randomly cou-
pled modes traveling in opposite directions has been based on physi-
cal arguments [8, 9]. A general approach, based on the present meth-
ods, provides detailed results for various transfer function statistics
of multi-layer coatings. This is possible because useful devices of
the class studied here must exhibit behavior close to that of their
theoretical design. We combine perturbation theory with Kronecker
product methods to obtain accurate, general results for the trans-
mission and reflection statistics of multi-layer optical coatings whose
layers exhibit random optical thickness.
For convenience set the output electric field equal to unity in Fig-
ure 9.1 and Equations (9.1) and (9.4):
Et = 1 (9.6)
dk = dk0 + δk (9.9)
where dk0 is the design thickness of the kth layer and δk is its ran-
dom thickness error. For the ideal design, substitute dk = dk0 into
Equation (9.5); denote the corresponding results of Equations (9.7)
and (9.2) as R0 T0 R0 and T 0 for the ideal device.
Next, let the different layers have random thickness variations
δk of Equation (9.9). Equations (9.8) and (9.5) become
1
1
Ei + ei 1
ni
1 1 1
= · Mk · · (9.10)
Er + er 2 1 k=1 nt −nt 0
1 −
ni
where
j sin φk
cos φk nk 2π
Mk = φk = n d + δk (9.11)
λ k k0
jnk sin φk cos φk
Er + er 1
R= T = (9.12)
Ei + ei Ei + ei
9.2. MATRIX ANALYSIS 149
2 2
Er + er 1 nt
R= T = (9.13)
Ei + ei Ei + ei ni
2
1 e ei
T = 1− i + (9.14)
Ei Ei Ei
ei = Ei + ei − Ei (9.15)
2
1 3 E +e
T = 3− Ei + ei + i 2 i (9.16)
Ei Ei Ei
2
E + er 3 E +e
R= r 3− E + ei + i 2 i (9.17)
Ei Ei i Ei
6 8
T = − Re Ei Ei∗ + e∗i
Ei 2 Ei 4
2 Ei + ei 2 nt
+ 6
Re Ei2 Ei∗ + e∗i 2
+ (9.18)
Ei Ei 4 ni
150 MULTI-LAYER COATINGS WITH RANDOM OPTICAL THICKNESS
6 8
R = Er + er 2 · 2
− Re Ei Ei∗ + e∗i
Ei Ei 4
2 Ei + ei 2
+ Re Ei2 Ei∗ + e∗i 2
+ (9.19)
Ei 6 Ei 4
ni = nt = 1 (9.29)
k 1 2 3 4 5 6 7 8 9 10 11 12 13
nk 4 1.35 4 4 1.35 4 1.35 4 1.35 4 4 1.6 4
dk0 (nm)2 31.2 92.6 31.2 31.2 92.6 31.2 92.6 31.2 92.6 31.2 31.2 78.1 31.2
π π π π π π π π π π π π π
φk
2 2 2 2 2 2 2 2 2 2 2 2 2
(9.30)
1
δ < δmax
p δ = 2δmax (9.31)
0 δ > δmax
2
These values are approximate. The exact values, given by dk = 125/nk correspond-
ing to φk = π/2 for this design, are used throughout the following calculations.
154 MULTI-LAYER COATINGS WITH RANDOM OPTICAL THICKNESS
1.0
0.8
0.6
T0
0.4
0.2
0
480 490 500 510 520
– nm
9.4.2. Transmittance
Random thickness errors cause the average transmittance T to
depart from the design value T 0 of Figure 9.2. Figure 9.3 shows
this deviation T 0 − T , determined from the expected value of
the second-order series expansion of Equation (9.18) and the re-
sults of Section 9.3, with random thickness errors having parameters
of Equations (9.31) and (9.34).
The in-band transmittance has been reduced by a small amount
near the band center, and increased slightly at the band edges. Other
results, not given here, show that the ratio T /T 0 is very close to 1
over the entire band 475–525 nm. Therefore the tolerances have not
reduced the out-of-band discrimination significantly. Calculation of
the in-band signal distortion requires in addition the two-frequency
transmission statistics, discussed in Section 9.4.3.
T0 – T
0.008
0.006
0.004
0.002
0
480 490 500 510 520
– nm
–0.002
T = T − T (9.36)
0.010
0.008
0.006
0.004
0.002
0
490 495 500 505 510
– nm
FIGURE 9.4. Transmission transfer function variance.
with T1 T2∗ given by the expected value of Equation (9.20). The
two-frequency correlation coefficient is given by
T1 · T2∗
C12 = (9.39)
T1 2 · T2 2
158 MULTI-LAYER COATINGS WITH RANDOM OPTICAL THICKNESS
9.5. DISCUSSION
REFERENCES
CHAPTER TEN
Conclusion
The main thrust of the present work is the use of Kronecker prod-
uct methods in the calculation of transmission statistics of multi-
mode systems with random coupling. Such systems exhibit linearity
between complex mode amplitudes; Kronecker products yield linear
relations between corresponding moments of the mode amplitudes.
These methods require that the system must be regarded as a
cascade of sections, described by a matrix product. This descrip-
tion is natural for systems such as multi-layer coatings. For systems
with continuous coupling, a matrix description must be artificially
introduced.
For forward-traveling modes with white coupling, i.e., resulting
from imperfections with a very short correlation length, the sys-
tem may be divided into sections of infinitesimal length; exact linear
differential equations for the mode powers or for higher moments
are obtained. For more general spectra, the sections must be long
enough for the coupling in different sections to be approximately
statistically independent, and the coupling must be small enough for
perturbation theory to be applied to the individual sections; the re-
sulting coupled power or other differential equations for the higher
moments are only approximate.
Multi-layer coatings offer an important two-mode system with
modes traveling in opposite directions. A combination of Kronecker
product methods with perturbation theory yields approximate trans-
fer function statistics for systems with performance close to that of
the ideal design.
161
162 CONCLUSION
REFERENCES
APPENDIX A
in order to obtain p00 and p10 . Let G0n z and G1n z be the
nth approximation to the solution of Equation (2.33). Let the initial
approximation be given simply by the initial conditions of Equation
163
164 SERIES SOLUTION FOR THE COUPLED LINE EQUATIONS
G01 z = 1
z
G11 z = j cxe−γx dx
0
z x
G02 z = 1 − cxe+γx dx cye−γy dy
0 0
z (A.4)
G12 = j cxe−γx dx
0
∞
p00 = 1 + p00n z
n=1
(A.8)
∞
p10 = p10n z
n=1
Note that:
The first and second summations of Equations (A.6) and (A.8) con-
tain only even and odd terms, respectively.
We obtain the following bounds on the terms of Equations (A.6)
and (A.8):
z
≤ 0
cx dx n
n even.
p00n z n! (A.10)
= 0 n odd.
= 0
n even.
p10n z z
cx dx n − (A.11)
≤ 0
e
z
0 αxdx
n odd.
n!
Let Equation (A.10) hold for some even n. Then from the middle
relation of Equation (A.7),
z x
x
cy dy n
0
p10n+1 z ≤ cx e− 0 αydy
dx
0 n!
e−
z
0 αydy z x n x
≤ cy dy d cy dy
n! 0 0 0
z n+1
0
cx dx z
= e− 0 αxdx
(A.12)
n + 1!
166 SERIES SOLUTION FOR THE COUPLED LINE EQUATIONS
where
Define
z
γij z = i z − j zdz (A.20)
0
∞
Gz = G0 + pn z (A.22)
n=1
where
z
pn z = j cxKx · pn−1 xdx p0 z = G0 (A.23)
0
REFERENCES
APPENDIX B
General Transmission
Properties of Two-Mode Guide
Set
s = x1 − x2 + x3 − x4 + · · · − x2n (B.4)
s = Lτ (B.6)
in
Equation
(B.5), substitute into Equation (B.1), and interchange
and ; then
1
∞
G0 α β = 1 + L fn Lτ eαLτ ejβLτ dτ (B.7)
0 n=1
Finally, define
∞
f Lτ = fn Lτ (B.8)
n=1
to yield
1
G0 α β = 1 + f LτeαLτ ejβLτ (B.9)
0
= A (B.12)
1 1
Aα β = Aα β
π β
1 ∞ Aα σ
= dσ (B.13)
π −∞ β − σ
+ A 1 B
A= A = lim Aα σdσ (B.14)
B→∞ 2B −B
1 1
G0 α β = G0 0 β (B.16)
2π α − jβ
1 α
Aα β = A0 β (B.17)
π α + β2
2
1 α
α β = 0 β (B.18)
π α + β2
2
1 ∞ 1
dβ = 1 (B.19)
πα −∞ 1 + β 2
α
1 1
lim 2
= δβ (B.20)
α→0 πα
1 + β α
2πn
βn = (B.21)
L
Then,
∞
2πn sin π βL − n
G0 α β = ejπ 2π
βL
2π
G0 α −1n
n=−∞ L βL
π 2π − n
(B.22)
REFERENCES
APPENDIX C
Kronecker Products
A ⊗ B = B ⊗ A (C.4)
A ⊗ B ⊗ C = A ⊗ B ⊗ C = A ⊗ B ⊗ C (C.5)
A + B ⊗ C + D = A ⊗ C + A ⊗ D + B ⊗ C + B ⊗ D (C.6)
A ⊗ B · C ⊗ D = A · C ⊗ B · D (C.7)
A1 ⊗ A2 ⊗ A3 ⊗ · · · · B1 ⊗ B2 ⊗ B3 ⊗ · · ·
· C1 ⊗ C2 ⊗ C3 ⊗ · · · · · · ·
= A1 · B1 · C1 · · · · ⊗ A2 · B2 · C2 · · · ·
⊗ A3 · B3 · C3 · · · · ⊗ · · · (C.8)
A · ai = αi ai B · bj = βj bj (C.11)
Define γk and ck by
A ⊗ B · ck = γk ck (C.12)
γk = αi βj ck = ai ⊗ bj (C.13)
a · bT = a ⊗ bT = bT ⊗ a (C.14)
REFERENCES
APPENDIX D
Expected Values of
Matrix Products
is
a11 a12 ···
A =
a21 a22 ··· (D.6)
Let Mk denote the wave matrix describing the kth section:
where I0 is the vector input, and the Mk are N × N random
matrices.
Assume different Mk are independent. Any element of Mk
is statistically independent of any element of Mr, k = r; however,
different elements of any particular matrix are not independent, e.g.,
see Equation (3.24). Denote the elements of Mk by mij k; then
mij k is independent of mpq r if k = r for any i j p q The th
element of Equation (D.9) is
I n = mk1 nmk1 k2 n − 1 · · · mkn−1 kn Ikn 0 (D.10)
k1 ···kn
D.1. INDEPENDENT MATRICES 181
n
= Mn + 1 − k · I0 (D.12)
k=1
follows from Equation (D.12). Finally, if the Mk have the same
statistics,
n
In ⊗ I ∗ n = Mn + 1 − k ⊗ M∗ n + 1 − k · I0 ⊗ I ∗ 0
k=1
(D.15)
average output:
n
In ⊗ I ∗ n = Mn + 1− k ⊗ M∗ n + 1− k · I0 ⊗ I ∗ 0
k=1
(D.16)
If the Mk have the same statistics,
In ⊗ I ∗ n = M ⊗ M∗ n · I0 ⊗ I ∗ 0
(D.17)
Ik ⊗ I ∗ k = M ⊗ M∗ · Ik − 1 ⊗ I ∗ k − 1
We write out Equation (D.17) explicitly for the two-mode case:
I0 n2 I0 02
I0 nI1∗ n I 0I1∗ 0
= M ⊗ M∗ n · 0 (D.18)
∗ ∗
I1 nI0 n I1 0I0 0
I1 n2 I1 02
where
m11 2 m11 m∗12 m12 m∗11 m12 2
m11 m∗21 m11 m∗22 m12 m∗21 m12 m∗22
M ⊗ M =
∗
∗ ∗ ∗ ∗
(D.19)
m21 m11 m21 m12 m22 m11 m22 m12
m21 2 m21 m∗22 m22 m∗21 m22 2
The input vector on the right-hand side of Equation (D.18) will nor-
mally be simplified. For example, if the signal mode I0 has determin-
istic unit input and the spurious mode I1 has zero input, the input
vector will be
1
0
I0 ⊗ I ∗ 0 = I0 ⊗ I ∗ 0 =
(D.20)
0
0
Higher order statistics are obtained in a similar manner. The ma-
trices in the analysis increase in size as the order and/or the number
of modes increases.
D.2. MARKOV MATRICES 183
a1 a2 · · · aL (D.21)
A1 A2 · · · AL (D.22)
Then,
pk = p (D.28)
From Equations (D.27) and (D.28), the solution for the unit eigen-
value of
P − I · p = 0 (D.29)
a1 a1 0 ··· 0
a2 0 a2 ··· 0
a=
aD =
(D.30)
aL 0 0 ··· aL
p1 p1 0 ··· 0
p2 0 p2 ··· 0
p=
pD =
(D.31)
pL 0 0 ··· pL
1T = 1 1 ··· 1 (D.32)
L elements
D.2. MARKOV MATRICES 185
Therefore,
MnM1 = ai PMn = ai M1 = aj aj (D.34)
ij
(D.45)
Now
(D.47)
Similarly,
A · pj = Aj pj (D.48)
where I0 is the vector input, and the other matrix quantities are
given in Equations (D.39)–(D.43).
Second- (and higher) order statistics are again found by use of
Kronecker products. For second-order, make the following substitu-
tions in the above analysis:
Mk −→ Mk ⊗ M∗ k
Ai −→ Ai ⊗ Ai∗
(D.51)
N −→ N 2
A −→ A⊗
Equations (D.39)–(D.43) become
1 0 ··· 0
0 1 ··· 0
I N 2 ×N 2 =
(D.52)
0 0 0 1
A1 ⊗ A1∗ 0 0 0
0 A2 ⊗ A2∗ ··· 0
A⊗
LN 2 ×LN 2
=
(D.54)
0 0 ··· AL ⊗ AL∗
p1 I
p2 I
pLN 2 ×N 2 = pL×1 ⊗ I N 2 ×N 2 =
(D.55)
pL I
D.2. MARKOV MATRICES 189
Note that A⊗
LN 2 ×LN 2
of Equation (D.54) is not equal to the Kronecker
product of ALN×LN of Equation (D.41) with itself. Then the expected
value of Equation (D.15) yields
n−1
1 1 − pA1 pA1 A1
In = I I · · · I0
2 pA2 1 − pA2 A2
(D.58)
In ⊗ I ∗ n
n−1
1 1 − pA1 ⊗ A1∗ pA1 ⊗ A1∗
= I I ·
2 pA2 ⊗ A2∗ 1 − pA2 ⊗ A2∗
A1 ⊗ A1∗
· · I0 ⊗ I ∗ 0 (D.59)
A2 ⊗ A2∗
REFERENCES
APPENDIX E
ɡt = ɡ t
∗
ɡt = 0 t < 0 (E.2)
∞
f = ɡte−j2πft dt (E.3)
0
f = (E.7)
191
192 TIME- AND FREQUENCY-DOMAIN STATISTICS
= (E.10)
2 1 F F
− = lim f f df df
F→∞ 2F2 −F −F
−2 + 2
1 F F
= lim f − f df df − 2
F→∞ 2F2 −F −F
2
= ∞ − = 0 (E.11)
Therefore,
2 ∗
ν = + ac ν ac ν = ac f + νac f (E.13)
where
t is given by Equation (3.107) in terms of the convolution
operator :
1 1 ∞ τ
t = t = dτ (E.18)
πt π −∞ t − τ
1 1
+ t = t + j
t ⇐⇒ + f = 1 + sgn f f (E.19)
2 2
1 1
− t = t − j
t ⇐⇒ − f = 1 − sgn f f (E.20)
2 2
where sgn is the signum function:
1 f > 0
sgn f = (E.21)
−1 f < 0
t ⇐⇒ 4
2
+ f − f = 4+ f + ∗ f
f0 +B−f +f /2
=4 a + f ∗ ada
f0 −B−f −f /2
(E.22)
194 TIME- AND FREQUENCY-DOMAIN STATISTICS
( f )
–( f ) +( f )
f
– f0 – B – f0 – f0 + B f0 – B f0 f0 + B
Since
2
sin 2πBt
lim 2B = δt (E.26)
B→∞ 2πBt
t
2
≈ 8B · t B large, (E.27)
t = 4
2
+
t− t
∞ ∞
= τ1 τ2 τ1 τ2
−∞ −∞
t
2
= t
2
∞
=
−∞
2
τ t − τ dτ
2
= 2
t t
2
(E.30)
1
C0 dxdx · j sgn f ⇐⇒ −C0 dxdx (E.32)
πt
196 TIME- AND FREQUENCY-DOMAIN STATISTICS
REFERENCES
APPENDIX F
2π
k0 = propagation constant in free space. (F.2)
λ0
2π
ki = propagation constant in medium of index ni . (F.3)
λi
i λ
ni = = 0 n1 > n2 (F.4)
0 λi
The propagation constants of TE modes traveling in the z direction,
indexed by the subscript ν, are given by
βν = n21 k20 − κ2ν = k21 − κ2ν (F.5)
197
198 SYMMETRIC SLAB WAVEGUIDE—LOSSLESS TE MODES
n2
n1
d
z
n2
n1 > n2
where
γ= n21 − n22 k20 − κ2 = k21 − k22 − κ2 (F.6)
κd γ
tan = even modes, (F.7)
2 κ
κd κ
tan =− odd modes. (F.8)
2 γ
tan d
2
5
0
0 2 4 6 8 10 12 d
–5
–10
max d = 3.5
n2
n1 c(z)
d z
n2
n1 > n2
For this case, κmax d = 3 5π, there are four TE modes, two even
and two odd. More generally, for N − 1π < κmax d < Nπ there are
N TE modes, N + 1/2 even and N/2 odd, where x represents
the greatest integer contained in x.
Finally, consider the slab guide with random straightness deviation
cz shown in Figure F.3. The coupling coefficients between two
forward-traveling TE modes are given as follows:
2κµ κν εµν
Cµν = − 1/2 (F.12)
1 1
βµ βν d+ d+
γµ γν
1 one mode even and the other odd.
εµν = (F.13)
0 both modes even or both odd.
Cµν is to be substituted into Equations (3.35) and (3.13); cz in
Equation (3.9) represents the straightness deviation of the random
guide, shown in Figure F.3.
F.2. EXAMPLE
2π
κmax = 1 012 − 1 m−1 (F.15)
10−6
Take
κmax d = 3 5π (F.16)
3 5 × 10−6
d= √ = 1 2343 × 10−5 m. (F.17)
2
2 1 01 − 1
γ 3 5π2
= −1 (F.18)
κ κd2
Equation (F.18) with Equations (F.7) and (F.8) yield the following
solutions for the eigenvalues, illustrated in Figure F.21 :
Substituting these values into Equations (F.5), (F.6), and (F.12), the
mode propagation constants and coupling coefficients to be substi-
tuted into Equations (3.35) and (3.36) are
6 3424 0 0 0
6
0 6 3316 0 0 −1
= j × 10 m (F.20)
0 0 6 3142 0
0 0 0 6 2921
1
We number the eigenvalues 0 1 2 3 · · · to correspond with the notation used
throughout the text.
202 SYMMETRIC SLAB WAVEGUIDE—LOSSLESS TE MODES
0 2 3225 0 4 4966
2 3225
9
0 6 8737 0 −2
C = 10 m (F.21)
0 6 8737 0 13 3083
4 4966 0 13 3083 0
REFERENCE
APPENDIX G
0 = 1 = jβ (G.1)
Cz = 0 (G.6)
and variance
z z
C 2 z = cxcydx dy
0 0
z z
= S0 dy δx − ydx = S0 z (G.7)
0 0
1
P 0 z = P0 z = I0 z 2 = cos2 Cz =1 + e−2S0 z
2
(G.9)
2 2 1 −2S0 z
P 1 z = P1 z = I1 z = sin Cz = 1 − e
2
the mode powers (but not the complex mode amplitudes) are deter-
ministic in the degenerate case, and consequently have zero variance,
in agreement with the results of Section 4.3.2 for the β = 0 case.
For two-mode incoherent input,
z = I (G.18)
= α + jβ (G.19)
with I the unit matrix of Equation (2.52). Then from Equation (2.53)
ments, as in the prior results for white coupling. For simplicity, the
following treatment is restricted to the large z case. Then, Equation
(G.32) yields
APPENDIX H
M1 = V ⊗ V ∗ (H.1)
∗ ∗
V = ez · C · e−z V ∗ = e z · C · e− z (H.2)
1 ≤ i j ≤ N (H.3)
M1 ij = Vij V ∗ = e i−1j−1 z
Ci−1j−1 V ∗ 1 ≤ i j ≤ N (H.4)
N−1
M2 ii kk = C2 kk = 2
Ck−1 1 ≤ i k ≤ N (H.6)
=0
N−1
M3 ii kk = C2 ii = 2
Ci−1 1 ≤ i k ≤ N (H.7)
=0
APPENDIX I
Differential Equations
Corresponding to
Matrix Equations
dy
= Dy (I.3)
dz
The solution to Equation (I.3), evaluated at the points z = kz, is
ln M
D= (I.5)
z
211
212 DIFFERENTIAL EQUATIONS CORRESPONDING TO MATRIX EQUATIONS
M ≈ 1 + Ez E z 1 (I.6)
D ≈ E (I.7)
d
yz = D · yz (I.10)
dz
Its solution is
eDz = M (I.12)
Diagonalizing
λM1 0 ··· 0
0 λM2 ··· 0
M = LM · M · −1
LM M =
(I.13)
0 0 ··· λMN
I.2. MATRIX CASE 213
λD1 0 ··· 0
0 λD2 ··· 0
D = LD · D · LD−1 D =
(I.14)
0 0 ··· λDN
Then,
eλD1 z 0 ··· 0
0 eλD2 z ··· 0 −1
Dz
e = LD ·
·L
D (I.15)
0 0 ··· eλDN z
ln λMi
LD = LM λDi = i = 1 · · · N (I.16)
z
M ≈ I + Ez E z 1 (I.19)
D ≈ E (I.20)
REFERENCE
APPENDIX J
Random Square-Wave
Coupling Statistics
J.1. INTRODUCTION
The coupling function of Equation (8.1) and Figure 8.1 may be writ-
ten as follows:
z 1
cz = ck rect − −k (J.1)
2
where
z 1 1 k < z < k + 1
rect − −k = (J.2)
2 0 otherwise.
ck = 0 (J.5)
Then Sν, the spectral density of the coupling cz defined in Equa-
tions (3.2) or (6.2), is given by
2
sin πν
∞
Sν = φc κe−j2πνκ (J.6)
πν κ=−∞
∞
∞
ck = dak ak = dn hk − n = hndk−n (J.9)
n=−∞ n=−∞
dk = ±1 dk = 0 (J.10)
∞
ϕh κ = hk + κhk (J.12)
k=−∞
J.2. BINARY SECTIONS 217
Let the ck of Equations (8.1) and (J.1) take on the values ai ,
1 ≤ i ≤ L, with unconditional probabilities
Pck = ai cn = aj k = n (J.19)
218 RANDOM SQUARE-WAVE COUPLING STATISTICS
J.2.1. Independent
For independent ck ,
J.2.2. Markov
Let p be the probability that adjacent ck are different:
κ 0 1 2 3 4 ≥5
(J.28)
ϕLP κ 5 4 3 2 1 0
2
2 sin πν
4
Sν = d 5+2 5 − κ cos2πνκ (J.29)
πν κ=1
We omit the closed form for the summation in this equation, subse-
quently using MAPLE to evaluate it directly.
The ck for this model are correlated, pseudo-Gaussian random
variables, with the following probability distribution:
ak −5 −3 −1 1 3 5
ck = dak 1 5 10 10 5 1 (J.30)
Pak
32 32 32 32 32 32
dk−2 dk−1 dk dk+1 dk+2 . The state can take on the (vector) val-
ues Si 1 ≤ i ≤ 32 given in the following table:
i Si aiLP
1 −1 −1 −1 −1 −1 −5
2 −1 −1 −1 −1 1 −3
3 −1 −1 −1 1 −1 −3
4 −1 −1 −1 1 1 −1
5 −1 −1 1 −1 −1 −3
6 −1 −1 1 −1 1 −1
7 −1 −1 1 1 −1 −1
8 −1 −1 1 1 1 1
9 −1 1 −1 −1 −1 −3
10 −1 1 −1 −1 1 −1
11 −1 1 −1 1 −1 −1
12 −1 1 −1 1 1 1
13 −1 1 1 −1 −1 −1
14 −1 1 1 −1 1 1
15 −1 1 1 1 −1 1
16 −1 1 1 1 1 3 (J.32)
17 1 −1 −1 −1 −1 −3
18 1 −1 −1 −1 1 −1
19 1 −1 −1 1 −1 −1
20 1 −1 −1 1 1 1
21 1 −1 1 −1 −1 −1
22 1 −1 1 −1 1 1
23 1 −1 1 1 −1 1
24 1 −1 1 1 1 3
25 1 1 −1 −1 −1 −1
26 1 1 −1 −1 1 1
27 1 1 −1 1 −1 1
28 1 1 −1 1 1 3
29 1 1 1 −1 −1 1
30 1 1 1 −1 1 3
31 1 1 1 1 −1 3
32 1 1 1 1 1 5
The left column i indexes the states; the middle column lists the state
vectors Si ; the final column aiLP , calculated from Si and Equation
(J.31), determines the coupling constant ck of the kth section as
ck = daiLP STATEk = Si (J.33)
for the ith state at time k.
J.3. MULTI-LEVEL MARKOV SECTIONS 221
aD = d
a 0
D11
(J.35)
0 a D22
(J.36)
222 RANDOM SQUARE-WAVE COUPLING STATISTICS
−3
−1
0
−1
1
−1
1
1
a
3
D22 =
−1
1
1
3
1
3
3
0 5
(J.37)
Then since dk+2 can take on only the values ±1, STATEK can have
only the two possible values
and
dk−2 dk−1 dk dk+1 dk+2
κ 0 1 2 3 4 ≥5
(J.45)
ϕBP κ 4 2 −1 −2 −1 0
2
sin πν
2
Sν = d 4 + 22 cos2πν − cos4πν
πν
− 2 cos6πν − cos8πν (J.46)
Equation (J.32) remains valid with the exception of the last column,
relabeled aiBP , calculated from the Si of Equations (J.32) and (J.47):
i 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
aiBP 0 2 2 4 0 2 2 4 −2 0 0 2 −2 0 0 2
i 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
aiBP −2 0 0 2 −2 0 0 2 −4 −2 −2 0 −4 −2 −2 0
(J.48)
REFERENCES
APPENDIX K
The electric and magnetic fields for the linearly polarized wave in
Figure 9.1 consist of forward and backward waves:
Ez = E + z + E − z Hz = H + z + H − z (K.1)
where
E + z E − z
H + z = H − z = − (K.2)
Z Z
and Z is the characteristic impedance of the dielectric medium:
µ0
Z= (K.3)
ε
The forward and backward waves have the following z dependence:
E + z = E + e−jβz E − z = E − e+jβz (K.4)
where the propagation constant is
2πn
β= (K.5)
λ
λ is the free-space wavelength,
1
λ= √ (K.6)
f µ0 ε0
227
228 MATRIX FOR A MULTI-LAYER STRUCTURE
2π
φk = βk dk = n d (K.16)
λ k k
The electric and magnetic fields are continuous across the bound-
aries at 0 z1 · · · z Therefore,
Ez0 − Ez +
= Mk · (K.17)
Z0 Hz0 − k=1 Z0 Hz +