Next Article in Journal
Lower Local Uniform Monotonicity in F-Normed Musielak–Orlicz Spaces
Previous Article in Journal
Satisfiability Threshold of Random Propositional S5 Theories
Previous Article in Special Issue
Uniqueness Results for Some Inverse Electromagnetic Scattering Problems with Phaseless Far-Field Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Covariance and Uncertainty Principle for Dispersive Pulse Propagation

Department of Physics, City University of New York, 695 Park Avenue, New York, NY 10065, USA
Axioms 2024, 13(4), 242; https://doi.org/10.3390/axioms13040242
Submission received: 15 February 2024 / Revised: 27 March 2024 / Accepted: 29 March 2024 / Published: 8 April 2024

Abstract

:
We develop the concept of covariance for waves and show that it plays a fundamental role in understanding the evolution of a propagating pulse. The concept clarifies several issues regarding the spread of a pulse and the motion of the mean. Exact results are obtained for the time dependence of the covariance between position and wavenumber and the covariance between position and group velocity. We also derive relevant uncertainty principles for waves.

1. Introduction

For waves governed by a linear differential equation of the form
n = 0 N x b n n u x n n = 0 N t a n n u t n = 0
one substitutes e i k x i ω t into Equation (1) to obtain [1,2,3,4] the dispersion relation
n = 0 N x b n ( i k ) n n = 0 N t a n ( i ω ) n = 0
Generally, there is more than one solution to the dispersion relation, and the solutions may be complex. Each solution is called a mode. In this paper, we deal with the propagation of a single mode, and we write
ω = ω ( k )
Furthermore, we assume that ω ( k ) is real, which implies that there is no damping. The solution to the wave equation for each mode is [1]
u ( x , t ) = 1 2 π S ( k , 0 ) e i k x i ω ( k ) t d k
where S ( k , 0 ) is the initial spatial spectrum, which is obtained from the initial pulse, u ( x , 0 ) , by way of
S ( k , 0 ) = 1 2 π u ( x , 0 ) e i k x d x
If one defines the time-dependent spectrum, S ( k , t ) , by [5,6,7]
S ( k , t ) = S ( k , 0 ) e i ω ( k ) t
then u ( x , t ) and S ( k , t ) form Fourier transform pairs between x and k for all time
u ( x , t ) = 1 2 π S ( k , t ) e i k x d k
S ( k , t ) = 1 2 π u ( x , t ) e i k x d x
Our aim is to study the dynamics of the modes and, in particular, to study the concept of covariance and uncertainty principles as they apply to wave propagation.
Because of the fact that u ( x , t ) and S ( k , t ) are Fourier transform pairs, one can calculate averages in either the position representation or the k-representation. For example, the moments of the pulse at time t are given by
x n t = x n | u ( x , t ) | 2 d x
but may also be calculated in the k-representation by way of
x n t = S * ( k , t ) X n S ( k , t ) d k
where the position operator, X , in the k-representation, is given by
X = 1 i k
Similarly, the k moments
k n t = k n | S ( k , t ) | 2 d k
may be calculated by way of
k n t = u * ( x , t ) K n u ( x , t ) d x
where K is the wavenumber operator in the position representation
K = 1 i x
Notation: Operators will be denoted by calligraphic letters. Commutators and anticommutators are denoted by the usual notation,
[ A , B ] = AB BA
and
[ A , B ] + = AB + BA
respectively. Indefinite integrals imply integration from to ∞, unless otherwise noted. T i m e -dependent moments, such as, for example, the mean, will be denoted as x t . For the case of time-dependent standard deviation for a quantity x, it will be denoted by σ x 2 ( t ) , and similarly for the covariance.

2. Motion of Mean and Spread

In Appendix A, we show that the mean of a pulse propagates with constant velocity
x t = x 0 + V t
where V is given by
V = W ( k ) k | S ( k , 0 ) | 2 d k
We define, as standard, the group velocity, v ( k ) , by
v ( k ) = W ( k ) k
The velocity, V , is the average of the group velocity at time zero,
V = v ( k ) 0
If we define the group velocity operator by v ( K ) , then V may also calculate the velocity in the position representation by way of
V = u * ( x , 0 ) v ( K ) u ( x , 0 ) d x
Comments:
1. It should be noted that the group velocity, v ( k ) , enters in a straightforward way, namely in the calculation of x t .
2. V not only depends on the group velocity, but also the initial spectrum. Thus, how a medium is disturbed affects the velocity of propagation of the mean and other quantities.
The spread of the pulse is defined in the usual way
σ x 2 ( t ) = x 2 t x t 2
is calculated to be
σ x 2 ( t ) = σ x 2 ( 0 ) + 2 t Cov x v ( 0 ) + t 2 σ v 2 ( 0 )
where
Cov x v ( 0 ) = 1 2 v X + X v 0 v 0 x 0
and σ v is the standard deviation of the group velocity
σ v 2 = v 2 V 2
The properties and meaning of Cov x v ( 0 ) will be discussed in the next section.

3. Covariance for Waves

In standard probability theory, the covariance between two random variables a and b is given by
Cov a b = a b a b
where a b is the first mixed moment between a and b. An important property of covariance is that
σ a 2 σ b 2 Cov a b 2
where σ a and σ b are the standard deviations of a and b .
For the case of operators, we define the covariance by
Cov a b = 1 2 AB + BA A B
and, when the quantities are time-dependent, we shall generally use the notation
Cov a b ( t ) = 1 2 AB + BA t A t B t
Also, we define the first mixed moment by
a b t = 1 2 AB + BA t = 1 2 S * ( k , t ) AB + BA S ( k , t ) d k
To show that the operator definition of covariance, Equation (27), also satisfies Equation (26), we consider the operator
C = A c B
where c is a real number. The standard deviation of C is
σ c 2 = C 2 C 2
= ( A c B ) 2 A c B 2
which simplifies to
σ c 2 = σ a 2 + c 2 σ b 2 c ( AB + AB ) 2 c A B
Using the operator definition of Cov, Equation (27), we have
σ c 2 = σ a 2 + c 2 σ b 2 2 c Cov a b
Now, we take
c = Cov a b σ b 2
and substitute it into Equation (34) to obtain
σ c 2 = σ a 2 σ b 2 Cov a b 2 σ b 2
Since all the standard deviations are positive, we must have
σ a 2 σ b 2 Cov a b 2
This shows that the operator definition of covariance, Equation (27), satisfies the standard condition for covariance, Equation (26).

3.1. Covariance between Position and Group Velocity

For the first mixed moment of position and group velocity, we have
x v t = 1 2 v X + X v t = 1 2 S * ( k , t ) v X + X v S ( k , t ) d k
Writing S ( k , t ) in terms of its amplitude and phase
S ( k , t ) = | S ( k , t ) | e i ψ ( k , t )
we obtain that
x v t = v ( k ) ψ ( k , t ) k | S ( k , t ) | 2 d k
This is true in general. But, for pulse propagation, we have
S ( k , t ) = | S ( k , 0 ) | e i ψ ( k , 0 ) i W ( k ) t
where ψ ( k , 0 ) is the phase of the spectrum at time zero. The phase at time t is, therefore, given by
ψ ( k , t ) = ψ ( k , 0 ) W ( k ) t
Hence,
ψ ( k , t ) k = ψ ( k , 0 ) k v ( k ) t
Substituting Equation (43) into Equation (40), we have
x v t = x v 0 + v 2 t
where
v 2 = v 2 ( k ) | S ( k , 0 ) | 2 d k
The covariance between position and group velocity at time t is, then,
Cov x v ( t ) = x v t V x t
= x v 0 + v 2 t V ( x 0 + V t )
which gives
Cov x v ( t ) = Cov x v ( 0 ) + t σ v 2
This shows that the covariance increases linearly in time. It also shows that, irrespective of the sign of the initial covariance, eventually, the covariance becomes positive and remains so.

3.2. Covariance between Position and Wave Number

Let us first consider the evolution of the first mixed moment
x k t = 1 2 k X + X k t
= 1 2 S * ( k , t ) k X + X k S ( k , t ) d k
Using Equation (39), direct calculation gives
x k t = k ψ ( k , t ) k | S ( k , t ) | 2 d k
This is true in general. But, for Equation (42),
ψ ( k , t ) k = ψ ( k , 0 ) k W ( k ) k t
and, therefore,
x k t = k ψ ( k , 0 ) k W ( k ) k t | S ( k , t ) | 2 d k
which gives
x k t = x k 0 + k v ( k ) 0 t
The covariance between x and k at time t is
Cov x k ( t ) = x k t x t k t
Putting Equation (54) into Equation (55), we obtain
Cov x k ( t ) = Cov x k ( 0 ) + t Cov k v ( 0 )

4. Uncertainty Principle for Waves

For any two operators, A and B , the usual textbook result for the uncertainty product is that [8,9]
σ a 2 σ b 2 1 4 | [ A , B ] | 2
However, there is a stronger version [10,11]
σ a 2 σ b 2 1 4 | [ A , B ] | 2 + Cov a b 2
where
Cov a b = 1 2 AB + BA A B
For completeness, a proof of Equation (58) is given in Appendix B.

4.1. Uncertainty Principle for Position and Wavenumber

For position and wavenumber, we have that
[ X , K ] = i
Therefore, the uncertainty product inequality is
σ x 2 ( t ) σ k 2 ( t ) 1 4 + Cov x k 2 ( t )
Using Equation (56), we may write
σ x 2 ( t ) σ k 2 ( t ) 1 4 + Cov x k 2 ( 0 ) + 2 t Cov k v ( 0 ) Cov x k ( 0 ) + t 2 Cov k v 2 ( 0 )
Since we know the exact standard deviations, given by Equations (22) and (24), we have that
σ x 2 ( t ) σ k 2 ( t ) = σ x 2 ( 0 ) σ k 2 ( 0 ) + 2 t Cov x v ( 0 ) σ k 2 + t 2 σ v 2 σ k 2 ( exact )
Using Equation (37), we have that
σ x 2 ( t ) σ k 2 ( t ) 1 4 + Cov x k 2 ( 0 ) + 2 t Cov k v ( 0 ) Cov x k ( 0 ) + t 2 σ v 2 σ k 2
which shows that the uncertainty principle is an approximate estimation of the uncertainty product, but gives the exact answer for time going to infinity.

4.2. Uncertainty Principle for Position and Group Velocity

For position X and group velocity v ( K ) , the commutator is
[ X , v ( K ) ] = i v ( K )
and, therefore, the uncertainty product for position and group velocity is given by
σ x 2 ( t ) σ v 2 ( t ) 1 4 | v ( K ) | 2 + Cov x v 2 ( t )
Using Equation (48), we have
σ x 2 ( t ) σ v 2 ( t ) 1 4 | v ( K ) | 2 + Cov x v 2 ( 0 ) + 2 t σ v 2 Cov x v ( 0 ) + t 2 σ v 4
The exact uncertainty product is given by
σ x 2 ( t ) σ v 2 ( t ) = σ x 2 ( 0 ) σ v 2 ( 0 ) + 2 t Cov x v ( 0 ) σ v 2 + t 2 σ v 4
In comparing the exact answer, Equation (68) with Equation (67), we see that the time-dependent terms are identical

5. Meaning of Group Velocity and Covariance

What is generally called the group velocity, v ( k ) , is a misnomer as it is not the velocity of the group of waves. The velocity of the group is V. The group velocity indicates what frequencies exist and their individual speeds. To ascertain the contribution of a particular k-frequency, it is weighted by the energy density spectrum | S ( k , 0 ) | 2 . The average of all frequencies is the velocity, V, as given by the integral in Equation (18).
To understand the meaning of covariance for waves, we consider a pulse that, for the moment, is traveling to the right, and has a spatial extension. Let us suppose that the k-frequencies are distributed within the spatial extension of the pulse in a random manner. In such a case, the covariance between the position x and v ( k ) would be zero. On the other hand, let us suppose that the high k-frequencies are concentrated on the right-hand side of the pulse, and the low k-frequencies on the left side of the pulse. In such a case, the covariance between v ( k ) and x would be positive and would indicate the relative relation between high k-frequency velocities and the front of the pulse. Conversely, the covariance would be negative if the high-frequency values are on the left side of the pulse. Hence, the reason that the covariance is fundamental is that it indicates whether the high k- frequencies are on the right or left side of the pulse. If the covariance is positive, the high frequencies are to the right of the low frequencies, but, if it is negative, it means that they are arranged the opposite way. Now, according to Equation (48),
Cov x v ( t ) = Cov x v ( 0 ) + t σ v 2
the evolution of the covariance is such that for large times, it is always positive, irrespective of the sign of the initial covariance
Cov x v ( t ) for t
That is, the pulse evolves so that the high velocity frequencies are at the front of the pulse. The reason is that, if the high frequencies are initially on the left side of the pulse, they will eventually catch up to the front and remain so.

6. Arbitrary Waves

In the above discussion, we have assumed a special type of motion, namely dispersive motion governed by Equation (4). It is of interest to consider what aspects of the previous results apply to an arbitrary pulse where we do not specify a particular equation of motion. Let us suppose the wave u ( x , t ) and its Fourier transform S ( k , t ) are written in terms of their respective amplitudes and phases,
u ( x , t ) = | u ( x , t ) | e i φ ( x , t ) = A ( x , t ) e i φ ( x , t )
S ( k , t ) = | S ( k , t ) | e i ψ ( k , t ) = B ( k , t ) e i ψ ( k , t )
If one calculates the first two moments and standard deviations, one obtains
x t = ψ k B 2 ( k , t ) d k
x 2 t = 1 B B k 2 + ψ k 2 B 2 ( k , t ) d k
σ x | t 2 = 1 B B k 2 + ψ k + x 2 B 2 ( k , t ) d k
and
k t = φ x A 2 ( x , t ) d x
k 2 t = 1 A A x 2 + φ x 2 A 2 ( x , t ) d x
σ k | t 2 = 1 A A x 2 + φ x k 2 A 2 ( x , t ) d x
Now, let us consider the first mixed moment of x and k:
x k = 1 2 X , K +
We calculate x k in the position representation by way of
x k t = 1 2 u * ( x , t ) x K + K x u ( x , t ) d x
resulting in
x k t = x φ x | u ( x , t ) | 2 d x
However, if one calculates it in the k representation using
x k t = 1 2 S * ( k , t ) k X + X k S ( k , t ) d k
then
x k t = k ψ k | S ( k , t ) | 2 d k
is obtained. Equating the two expressions, we have
x φ x | u ( x , t ) | 2 d x = k ψ k | S ( k , t ) | 2 d k
This is an interesting relation because it shows a connection between amplitudes and phase of the pulse and its spectrum.
The covariance, therefore, may be written as
Cov x k = x φ x x k
or
Cov x k = k ψ k x k
This shows that, in the x representation, the wavenumber, k, acts as φ x and that, in the wavenumber representation, x acts as ψ k .
In addition, we speculate that the minimum uncertainty relation takes the form
σ x | t 2 σ k | t 2 1 B B k 2 B 2 ( k , t ) d k 1 A A x 2 A 2 ( x , t ) d x

7. Quasi-Probability Approach

The concept of covariance implies the existence of a joint distribution of the two variables. In the past sections, we have defined it in terms of operators without consideration of a possible joint distribution. In the quantum mechanical case, Cartwright [12] was the first to define covariance in terms of a joint distribution of position and momentum. She used the Wigner distribution [13,14,15] for the joint distribution. In the case of position and wavenumber, the Wigner distribution is defined by
W ( x , k , t ) = 1 2 π u * ( x 1 2 τ , t ) u ( x + 1 2 τ , t ) e i τ k d τ
which, in terms of the spectrum, may be written as
W ( x , k , t ) = 1 2 π S * ( k + 1 2 θ , t ) S ( k 1 2 θ , t ) e i θ x d θ
Because the Wigner and other similar distributions are not manifestly positive, they are often called quasi-distributions. The Wigner distribution satisfies the marginals
W ( x , k , t ) d k = u ( x , t ) 2
W ( x , k , t ) d x = S ( k , t ) 2
Therefore, all functions of position or wavenumber are given correctly by the Wigner distribution. In particular, the moments x n and k n are correctly given by
x n = x n W ( x , k , t ) d k d x
k n = k n W ( x , k , t ) d k d x
Following Cartwright [12], we define the first mixed moment by
x k = x k W ( x , k , t ) d x d k
Straightforward calculation leads to Equation (81) and/or Equation (83).
We now specialize to the case of pulse propagation, where
S ( k , t ) = S ( k , 0 ) e i ω ( k ) t
We have
W ( x , k , t ) = 1 2 π S * ( k + 1 2 θ , 0 ) S ( k 1 2 θ , 0 ) e i θ x e i [ ω ( k + θ / 2 ) ω ( k θ / 2 ) ] t d θ
At t = 0 ,
W ( x , k , 0 ) = 1 2 π S * ( k + 1 2 θ , 0 ) S ( k 1 2 θ , 0 ) e i θ x d θ
and inverting yields
S * ( k + 1 2 θ , 0 ) S ( k 1 2 θ , 0 ) = W ( x , k , 0 ) e i θ x d x
Inserting Equation (98) into Equation (96), we have
W ( x , k , t ) = 1 2 π W ( x , k , 0 ) e i θ ( x x ) e i [ ω ( k + θ / 2 ) ω ( k θ / 2 ) ] t d θ d x
which expresses the Wigner distribution at time t in terms of the Wigner distribution at time zero.
Now, we consider the first mixed moment at time t,
x k t = 1 2 π x k W ( x , k , 0 ) e i θ ( x x ) e i [ ω ( k + θ / 2 ) ω ( k θ / 2 ) ] t d θ d x d x d k
Simplification leads to
x k t = x k W ( x , k , 0 ) d x d k + t k v ( k ) W ( x , k , 0 ) d x d k
and, therefore, we have
x k t = x k 0 + t k v ( k ) 0
which is the answer obtained by the operator method, namely Equation (54).
For the case of position and group velocity, the mixed moment is
x v ( k ) t = x v ( k ) W ( x , k , t ) d x d k
which, for pulse propagation, is given by
x v ( k ) t = 1 2 π x v ( k ) W ( x , k , 0 ) e i θ ( x x ) e i [ ω ( k + θ / 2 ) ω ( k θ / 2 ) ] t d θ d x d x d k
Simplification leads to
x v ( k ) t = x v ( k ) W ( x , k , 0 ) d x d k + t v 2 ( k ) W ( x , k , 0 ) d x d k
and, therefore,
x v ( k ) t = x v ( k ) 0 + t v 2 ( k ) 0
which is the same as Equation (44).

Other Quasi-Distributions

Cartwright [12] pointed out that there are an infinite number of distributions that give the same result as the Winger distribution. All quasi-distributions may be expressed as [10,12,16]
C ( x , k , t ) = 1 4 π 2 u * ( q 1 2 τ , t ) u ( q + 1 2 τ , t ) Φ ( θ , τ ) e i θ x i τ k + i θ q d θ d τ d q
where Φ ( θ , τ ) is the kernel function that characterizes the distribution. If we want the covariance calculated with C ( x , k , t ) to be the same as that with Wigner, that is,
x k C ( x , k , t ) d x d k = x k W ( x , k , t ) d x d k
then the kernel must satisfy
2 θ τ Φ ( θ , τ ) θ , τ = 0
There are an infinite number of distributions that satisfy Equation (109). Among them are the Choi–Williams [17,18], the Margenau–Hill [19], and Rihackek [20] distributions, among others [10].

8. Dual Case

In the previous sections, we considered the physical case where we have a spatial wave at a given time and allow it to evolve. An example is plucking a string and letting go at time zero. The dual case is where, at a fixed position, we create a pulse as a function of time, and examine its evolution as a function of position. Examples are radar, sonar, and a fiber optics pulse. For the fundamental solution, we take e i k x + i ω t instead of e i k x i ω t because that will make the results mathematically and physically similar to the previous case. We sbstitute e i k x + i ω t into Equation (1) to obtain
n = 0 N x b n ( i k ) n n = 0 N t a n ( i ω ) n = 0
In this case, we solve for k as a function of ω
k = k ( ω )
For the case where we generate a signal at position x = 0 , the general solution is then
u ( x , t ) = 1 2 π F ( ω , 0 ) e i k ( ω ) x + i ω t d ω
where F ( ω , 0 ) is the spectrum at x = 0 , given by
F ( ω , 0 ) = 1 2 π u ( 0 , t ) e i ω t d t
Analogous to Equation (6), we define
F ( ω , x ) = F ( ω , 0 ) e i k ( ω ) x
in which case,
u ( x , t ) = 1 2 π F ( ω , x ) e i ω t d ω
F ( ω , x ) = 1 2 π u ( x , t ) e i ω t d t
which shows that u ( x , t ) and F ( ω , x ) form Fourier transform pairs between frequency and time at position x.
Analogous to group velocity, one defines the “group slowness” or the “unit transit time” by
τ ( ω ) = d d ω k ( ω )
All the equations in the previous sections may be written down by analogy. Here, we give just the main ones relevant to the concept of covariance.
The time moments at position x are given by
t n x = t n | u ( x , t ) | 2 d t = F * ( ω , x ) T n F ( ω , x ) d ω
where T is the time operator in the frequency domain
T = i ω
If we write
F ( ω , x ) = | F ( ω , x ) | e i η ( ω , x )
then, in general,
τ t x = 1 2 τ T + T τ x = τ ( ω ) η ω | F ( ω , t ) | 2 d ω
The covariance between time and transit time at position x is given by
Cov t τ ( x ) = t τ x τ x t x
If we write the spectrum in terms of amplitude and phase
F ( ω , x ) = | F ( ω , x ) | e i η ( ω , x ) = | F ( ω , 0 ) | e i η ( ω , 0 ) e i k ( ω ) x
then
η ( ω , x ) = η ( ω , 0 ) k ( ω ) x
and
η ( ω , x ) ω = η ( ω , 0 ) ω τ x
Substituting Equation (125) into Equation (121), the mixed moment is
t τ x = t τ 0 + τ 2 x
from which it follows that
Cov t τ ( x ) = Cov t τ ( 0 ) + x σ τ 2
For the covariance between time and frequency at position x, defined by
Cov t ω ( x ) = t ω x t x ω x
we have
t ω x = 1 2 ω T + T ω x = ω η ( ω , x ) ω | F ( ω , x ) | 2 d ω
and
t ω x = t ω 0 + ω τ x
One obtains
Cov t ω ( x ) = Cov t ω ( 0 ) + x Cov ω τ ( 0 )

Funding

This research received no external funding.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Dynamics of Modes

We first obtain the differential equation for a mode, u ( x , t ) , characterized by ω ( k ) . We differentiate Equation (4) with respect to time to have
i t u ( x , t ) = 1 2 π ω ( k ) S ( k , 0 ) e i k x i W ( k ) t d k
= 1 2 π ω 1 i x S ( k , 0 ) e i k x i ω ( k ) t d k
which we write as
i t u ( x , t ) = ω K u ( x , t )
This equation is of the form of the Schrödinger equation with W ( K ) acting as the Hamiltonian. Since we are taking W to be real, W K is Hermitian, and, further e i W ( K ) t is unitary. The symbolic solution to Equation (A3) is
u ( x , t ) = e i ω ( K ) t u ( x , 0 )
The equation of motion for S ( k , t ) is obtained from Equation (6),
i t S ( k , t ) = ω ( k ) S ( k , t )

Appendix A.1. Dynamical Variable

Starting with
A t = u * ( x , t ) A ( 0 ) u ( x , t ) d x
The usual derivation of the Heisenberg equation of motion leads to [8,9]
A ( t ) = e i ω ( K ) t A ( 0 ) e i ω ( K ) t
and
d A d t = 1 i [ A , ω ( K ) ]

Appendix A.2. Equations of Motion for Position and Wave Number

In Equation (A8), we take A = K , and since K and W ( K ) commute, we have
d K d t = 1 i [ K , ω ( K ) ] = 0
and, therefore,
K ( t ) = K ( 0 )
Also, any function of K ( t ) will not change in time. In particular, all the moments are time-independent:
k n t = k n 0 .
Now, we consider the position operator X . Using Equation (A8), we have
d X d t = 1 i [ X , ω ( K ) ]
The commutator works out to be
[ X , ω ( K ) ] = 1 i d d k ω ( K ) = 1 i v ( K )
Hence,
d X d t = v ( K )
But
v ( K ( t ) ) = v ( K ( 0 ) )
and, therefore,
X ( t ) = X ( 0 ) + v ( K ( 0 ) ) t
Taking the expectation values of both sides, we have
x t = x 0 + V t
where
V = v ( k ) | S ( k , 0 ) | 2 d k
For the second moment, we square Equation (A16) to obtain
X 2 ( t ) = X 2 ( 0 ) + [ X ( 0 ) v ( K ) + v ( K ) X ( 0 ) ] t + v 2 ( K ) t 2
Taking expectation values gives
x 2 t = x 2 0 + t v X + X v 0 + t 2 v 2
where
v 2 = v 2 ( k ) | S ( k , 0 ) | 2 d k
and
v X + X v 0 = S * ( k , 0 ) v X + X v | S ( k , 0 ) d k
For the standard deviation, one obtains
σ x | t 2 = σ x | 0 2 + 2 t Cov x v ( 0 ) + t 2 σ v 2
where
Cov x v ( 0 ) = 1 2 v X + X v 0 v 0 x 0
and where σ v 2 is the dispersion of the group velocity
σ v 2 = ( v ( k ) V ) 2 | S ( k , 0 ) | 2 d k

Appendix B. Proof of the Uncertainty Principle with Covariance

We consider two Hermitian operators A and B and, without loss of generality, we can assume that their means are zero. The uncertainty product in a state ψ ( x ) is
σ a 2 σ b 2 = ψ * ( x ) A 2 ψ ( x ) d x × ψ * ( x ) B 2 ψ ( x ) d x
Since the operators are Hermitian, we have that
σ a 2 σ b 2 = A ψ ( x ) 2 d x × B ψ ( x ) 2 d x
Since the means are zero, the covariance is given by
Cov a b = 1 2 A , B +
Applying the Schwartz inequality to Equation (A27), we have
σ a 2 σ b 2 { A ψ ( x ) } * { B ψ ( x ) } d x 2
= ψ * ( x ) A B ψ ( x ) d x 2
giving that
σ a 2 σ b 2 | A B | 2
The product of two operators may be written as
A B = 1 2 A , B + + 1 2 i A , B / i
where both A , B + and A , B / i are Hermitian. Taking the expectation value of both sides, we have
A B = 1 2 A , B + + 1 2 i A , B / i
= Cov a b + 1 2 i A , B / i
Since A , B + and A , B / i are Hermitian, the expectation values on the right-hand side are real and, therefore,
| A B | 2 = 1 4 A , B 2 + Cov a b 2
Therefore, Equation (A31) gives that
σ a 2 σ b 2 1 4 A , B 2 + Cov a b 2

References

  1. Whitham, G. Linear and Nonlinear Waves; John Wiley and Sons: New York, NY, USA, 1974. [Google Scholar]
  2. Graff, K. Wave Motion in Elastic Solids; Oxford University Press: Oxford, UK, 1975. [Google Scholar]
  3. Jackson, J.D. Classical Electrodynamics; Wiley: Hoboken, NJ, USA, 1992. [Google Scholar]
  4. Morse, P.H.; Ingard, K.U. Theoretical Acoustics; McGraw-Hill: New York, NY, USA, 1968. [Google Scholar]
  5. Cohen, L. Why do wave packets sometimes contract? J. Mod. Opt. 2003, 49, 2365–2382. [Google Scholar] [CrossRef]
  6. Loughlin, P.; Cohen, L. Local properties of dispersive pulses. J. Mod. Opt. 2002, 49, 2645–2655. [Google Scholar] [CrossRef]
  7. Ben-Benjamin, J.S.; Cohen, L. Modes and Noise Propagation in Phase Space. Coherent Opt. Phenom. 2015, 2, 2299. [Google Scholar] [CrossRef]
  8. Bohm, D. Quantum Theory; Prentice-Hall: New York, NY, USA, 1951. [Google Scholar]
  9. Merzbacher, E. Quantum Mechanics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1998. [Google Scholar]
  10. Cohen, L. Time-Frequency Analysis; Prentice-Hall: Upper Saddle River, NJ, USA, 1995. [Google Scholar]
  11. Cohen, L. The Weyl Operator and Its Generalization; Birkhäuser: Basel, Switzerland; Springer: New York, NY, USA, 2013. [Google Scholar]
  12. Cartwright, N. Correlations without joint distributions in quantum mechanics. Found. Phys. 1974, 4, 127–136. [Google Scholar] [CrossRef]
  13. Wigner, E.P. On the quantum correction for thermodynamic equilibrium. Phys. Rev. 1932, 40, 749–759. [Google Scholar] [CrossRef]
  14. Hillery, M.; O’Connell, R.F.; Scully, M.O.; Wigner, E.P. Distribution functions in physics: Fundamentals. Phys. Rep. 1984, 106, 121. [Google Scholar] [CrossRef]
  15. Schleich, W.P. Quantum Optics in Phase Space; Wiley-VCH: Weinheim, Germany, 2001. [Google Scholar]
  16. Lee, H.W. Theory and application of the quantum phase-space distribution functions. Phys. Rep. 1995, 259, 147–211. [Google Scholar] [CrossRef]
  17. Choi, H.; Williams, W. Improved time-frequency representation of multicomponent signals using exponential kernels. IEEE Trans. ASSP 1989, 37, 862–871. [Google Scholar] [CrossRef]
  18. Jeong, J.; Williams, W. Kernel design for reduced interference distributions. IEEE Trans. Signal Process. 1992, 40, 402–412. [Google Scholar] [CrossRef]
  19. Margenau, H.; Hill, R.N. Correlation between Measurements in Quantum Theory. Prog. Theor. Phys. 1961, 26, 722–738. [Google Scholar] [CrossRef]
  20. Rihaczek, W. Signal Energy Distribution in Time and Frequency. IEEE Trans. Inf. Theory 1968, 4, 369–374. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cohen, L. Covariance and Uncertainty Principle for Dispersive Pulse Propagation. Axioms 2024, 13, 242. https://doi.org/10.3390/axioms13040242

AMA Style

Cohen L. Covariance and Uncertainty Principle for Dispersive Pulse Propagation. Axioms. 2024; 13(4):242. https://doi.org/10.3390/axioms13040242

Chicago/Turabian Style

Cohen, Leon. 2024. "Covariance and Uncertainty Principle for Dispersive Pulse Propagation" Axioms 13, no. 4: 242. https://doi.org/10.3390/axioms13040242

APA Style

Cohen, L. (2024). Covariance and Uncertainty Principle for Dispersive Pulse Propagation. Axioms, 13(4), 242. https://doi.org/10.3390/axioms13040242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop