tale of two formulas REVD
DOCX · 54.3 KB
Open DOCX file
Working note by Phil dated 8.6.13, written while rewriting Section 35 of his spectral theory book. It compares the single pulse train formula P(ω) with the ensemble-average formula <P(ω)>, using the autocorrelation <y_m y_{m+s}>_1 and ensemble <y_m y_{m+s}>. It argues stationarity can apply to finite trains and that the two formulas agree only as N goes to infinity. It ends with an open question about the normalization of r_s. Equations are partly garbled by extraction errors.
AI-written summary; may contain errors. This description is approximate.
Extracted text (machine-read; may contain errors)
Clarification of finite N tale of two formulas PhL 8.6.13
More issues now dealt with in the rewritten (as of Aug 8) Section 35.
(1) In an older FT this was the dichotomy I presented:
(1) one pulse train / autocorrelation, and // see black text below at the start
(2) ensemble.
In the current FT, I describe a completely different dichotomy,
(1) autocorrelation method
(2) double sum method
P(ω) = Ppulse(ω) (!Syntax Error, Irs z-s ) first line of (35.6) "autocorrelation"
P(ω) = Ppulse(ω) (!Syntax Error, Iym z-m )* (!Syntax Error, Iyn z-n ) (35.26)
= Ppulse(ω) (!Syntax Error, I !Syntax Error, I ym yn zm-n ) third line of (35.6) "double sum"
So what happened to the older dichotomy? In Section 35 I have sliced the pie differently. I have separate sections for "one pulse train" and "ensemble of pulse trains". I do show in the infinite case that the single pulse train and ensemble have the same statistics, such as <P(ω)> = P(ω). In fact I have a whole little section on that. For the infinite case, both methods of my first dichotomy always give the same result, so there is no point in computing things both ways (as I once was doing in AMI and Change/Hold). Thus I no longer present that as a dichotomy.
(2) Below I discovered the "gray region" for the first time.
What does stationarity say for my two formulas?
P(ω) = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s one pulse train z ≡ eiωT (35.6)
<P(ω)> = Ppulse(ω) !Syntax Error, I !Syntax Error, I <ymym+s> z-s ensemble of pulse trains (35.5)
It is only a property of <ymym+s>, not of <ymym+s>1 which is always "stationary" by its definition. So in the second equation, stationarity says
<P(ω)> = Ppulse(ω) !Syntax Error, I <ymym+s> z-s ( !Syntax Error, I[1] )
= Ppulse(ω) !Syntax Error, I <ymym+s> z-s
So with stationarity, here are our two formulas:
P(ω) = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s one pulse train
<P(ω)> = Ppulse(ω) !Syntax Error, I <ymym+s> z-s
So this makes pretty clear that for infinite pulse trains, you don't have to do things by "both methods" because they will always be the same. So I need to clean up AMI and Change Hold in that regard.
What happens if the sums are finite??? This is the big question. For a single finite pulse train,
ymym+s = ymym+s θ(m < N)θ(m>-N) θ(m+s < N)θ(m+s>-N)
which is non-vanishing only over the finite gray region shown below
Now consider the type 1 average over a pulse train
<ymym+s>1 = !Syntax Error, I ymym+s
The sum for two positive values of s is shown in the picture as a red line. There are many fewer terms in the sum shown by the upper red line which is for some s close to 2N, so one would expect <ymym+s>1 to be smaller for the upper value of s compared with the lower value of s. But <ymym+s>1 is always independent of m, and we write rs = <ymym+s>1.
Now consider the ensemble average <ymym+s>.
<ymym+s> = !Syntax Error, I ym(i)ym+s(i)
We know that ym(i)ym+s(i) is non vanishing only in the gray area shown above. Imagine now extruding the above picture out of the plane of paper, so we then have a separate (but identical) gray area for each pulse train i. Then the ensemble Σi sum here is a sum perpendicular to the plane of paper and is at a fixed point in the (m,s) plane. For some particular value of m, such as the left blue line, <ymym+s> will be non-vanishing for all s corresponding to the blue line. For a larger value of m, we get the right blue line which has a different range of s. Nothing in this picture precludes the possibility that <ymym+s> be the same for these two values of m and the same value of s, as suggested by the two black dots, each of which represents an ensemble sum perpendicular to the plane of paper.
Conclusion: The assumption of stationarity can be applied to both infinite and finite pulse trains. Nothing blocks the possibility for this assumption when dealing with a finite pulse train.
Now back to the two formulas:
P(ω) = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s one pulse train z ≡ eiωT (35.6)
<P(ω)> = Ppulse(ω) !Syntax Error, I !Syntax Error, I <ymym+s> z-s ensemble of pulse trains (35.5)
For a finite pulse train we can write these as
P(ω) = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s one pulse train z ≡ eiωT (35.6)
<P(ω)> = Ppulse(ω) !Syntax Error, I !Syntax Error, I <ymym+s> z-s ensemble of pulse trains (35.5)
With the stationarity assumption, the lower equation becomes
<P(ω)> = Ppulse(ω) !Syntax Error, I <ymym+s> z-s !Syntax Error, I[1]
= Ppulse(ω) !Syntax Error, I <ymym+s> z-s
so with the stationarity assumption our two formulas become
P(ω) = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s one pulse train z ≡ eiωT (35.6)
<P(ω)> = Ppulse(ω) !Syntax Error, I <ymym+s> z-s ensemble of pulse trains
These formulas are not the same. We have
<ymym+s>1 ≠ <ymym+s> for small N <P(ω)> ≠ P(ω)
<ymym+s>1 ≈ <ymym+s> for large N, <P(ω)> ≈ P(ω)
<ymym+s>1 = <ymym+s> N = ∞ <P(ω)> = P(ω)
Only in the N=∞ case are the two methods the same.
*****************
I thought I was doing well, but now I am confused again. I am losing the main point. There seem to be two formulas, but the second is just the ensemble average of the first. So you can say
P(ω) = Ppulse(ω)!Syntax Error, I rs z-s = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s
= Ppulse(ω) R"(z)
= Ppulse(ω) | Y"(z) |2
<P(ω)> = Ppulse(ω) !Syntax Error, I !Syntax Error, I <ymym+s> z-s
= Ppulse(ω) <R"(z)>
= Ppulse(ω) <| Y"(z) |2>
It just happens that rs = <ymym+s>1 which brings in this funny 1 type average. Stationarity applies only to the second equation which involves Random Variables, and then we have
Ppulse(ω) !Syntax Error, I !Syntax Error, I <ymym+s> z-s = Ppulse(ω) !Syntax Error, I <ymym+s> z-s!Syntax Error, I1
= Ppulse(ω) !Syntax Error, I <ymym+s> z-s
So we then end up with these two equations. The first is exact, the second is exact with the assumption of stationarity
P(ω) = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s
<P(ω)> = Ppulse(ω) !Syntax Error, I<ymym+s> z-s
They are NOT the same equation because they compute different things. It does seem appropriate then to compare the two types of averages as I have done. Now suppose Yn and Ym are independent. Then we have
P(ω) = Ppulse(ω)!Syntax Error, I <ymym+s>1 z-s = Ppulse(ω)R"(z) single pulse train
<P(ω)> = Ppulse(ω) [ <ym2> + !Syntax Error, I<ym>2 z-s ] ensemble of pulse trains
= Ppulse(ω) [ (σ2+μ2) + μ2 !Syntax Error, I z-s ]
Again, for the second line we have assumed both stationarity and independence, whereas in the first line we have assumed nothing at all. Both equations are valid for finite pulse trains as well as ∞ .
What then is the deal with a finite pulse train?
******************************************************
Whole new question! What is the meaning of rs ? I write it for a finite sequence this way
<rs> = !Syntax Error, I <ymym+s>
using my gray area thing. Where does this normalization factor come from? I defined rs for my own personal use in (32.17) and that is why I get R"(z)= | Y"(z) |2 and that in turn is why I got
P(ω) = Ppulse(ω) R"(z). The norm factor must be there to make this all work right.