Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Math / Spectral Theory Book / support

line code verifications

DOCX · 92.0 KB
Open DOCX file

Working notes dated 5.12.13 by Phil, in the support files for his spectral theory book. He compares his power spectral density formulas for unipolar and bipolar NRZ, RZ, Manchester and AMI codes with those of Smith and Xiong, converting from ω to f and setting p = 1/2. He finds and fixes a factor-of-4 error in the Manchester result. NRZI and hold/change codes remain to be checked.

AI-written summary; may contain errors.

Extracted text (machine-read; may contain errors)
Verification of Line Codes PhL 5.12.13 1. The NRZ waveforms For my unipolar NRZ I get x = ω/ω1 <P(ω)> = (V2/ω1) sinc2(πx) [ p(1-p) + p2 δ(x) ] // unipolar NRZ (36.3) and for my bipolar NRZ I get <P(ω)> = (V2/ω1) sinc2(πx) [ 4p(1-p) + (1-2p)2δ(x) ] // bipolar NRZ (36.6) If I set p = 1/2 in these results I get <P(ω)> = (V2/ω1) sinc2(πx) [1/4 + 1/4 δ(x) ] // unipolar NRZ (36.3) <P(ω)> = (V2/ω1) sinc2(πx) // bipolar NRZ (36.6) Smith defines his NRZ in terms of voltages V and -V, so his is bipolar NRZ. My result in f space is <P(f)>/2π = (V2/ω1) sinc2(πx) = V2T1/(2π) * sinc2( πx) <P(f)> = V2T1 sinc2( πx) = V2T1 sinc2(π fT1) x = ω/ω1 = T1/T = fT1 This is in agreement with Smith which I now quote and this agrees also with Xiong which I quote Going back to my unipolar, I get for p = 1/2 <P(ω)> = (V2/ω1) sinc2(πx) [1/4 + 1/4 δ(x) ] // unipolar NRZ (36.3) <P(f)>/2π = V2T1/(2π) * sinc2( πx) [ 1/4 + 1/4 δ(x) ] <P(f)> = V2T1 sinc2( π fT1) [ 1/4 + 1/4 δ(fT1) ] / (2π) = V2T1 sinc2( π fT1) [ 1/4 + 1/4 (1/T1)δ(f) ] = (1/4) V2T1 sinc2( π fT1) + 1/4 V2T1(1/T1)δ(f) = (1/4) V2T1 sinc2( π fT1) + (1/4) V2δ(f) This agrees exactly with Xiong, so I am happy with my NRZ verifications. 2. The RZ waveforms My RZ is unipolar and the pulse fills the left half of the gap as in my Fig 36.8. My result is <P(ω)> = (V/2)2 (1/ω1) sinc2( x) [ p(1-p) + p2δ(x) + p2!Syntax Error, I δ(x - m) ] (36.9) and with p = 1/2 this becomes [ recall x = fT1] <P(ω)> = (V/2)2 (1/ω1) sinc2( x) [ 1/4 + 1/4 δ(x) +1/4!Syntax Error, I δ(x - m) ] (36.9) Now write δ(x) = δ(fT1) = (1/T1)δ(f) δ(x-m) = δ(fT1-m) = (1/T1)δ(f - m/T1) so my result can be written <P(f)>/2π = (V/2)2 (T1/2π) sinc2( fT1) [ 1/4 + 1/4 (1/T1)δ(f) +1/4 (1/T1)!Syntax Error, Iδ(f - m/T1) ] Now I have to use, for m odd, sinc2( m) = sin2(πm/2)/(πm/2)2 = (2/πm)2 and then <P(f)> = (V/2)2 T1 { 1/4 sinc2( fT1) + 1/4 (1/T1)δ(f) + (2/πm)2 1/4 (1/T1)!Syntax Error, Iδ(f - m/T1) } Xiong's unipolar RZ seems to be the same as my waveform at least, But here is his result but what is Rb ? Earlier he says it is 1/T. Well, his form does not show explicitly that the even lines are missing. Let's start over with my OTHER form where I have [ p(1-p) + p2!Syntax Error, I δ(x - m) ] = [1/4 + 1/4!Syntax Error, I δ(x - m) ] = [ 1/4 + 1/4 (1/T1) !Syntax Error, I δ(f - m/T1) ] Then my expression says <P(f)>/2π = (V/2)2 (T1/2π) sinc2( fT1) [ 1/4 + 1/4 (1/T1) !Syntax Error, I δ(f - m/T1) ] <P(f)> = (V2/16) T1 sinc2( fT1) [ 1 + (1/T1) !Syntax Error, I δ(f - m/T1) ] and this DOES agree with Xiong's result. Fine. Verification complete. 2. Manchester Line Code My result here is <P(ω)> = (V/2)2 (1/ω1) sin2( x) sinc2( x) [4p(1-p) + (2p-1)2!Syntax Error, I δ(x - m) ] which with p = 1/2 becomes <P(ω)> = (V/2)2 (1/ω1) sin2( x) sinc2( x) <P(f)>/2π = (V/2)2 T1/2π sin2( x) sinc2( x) <P(f)> = (V/2)2 T1 sin2( x) sinc2( x) Here is the Xiong result for Manchester, where x = fT1 His result seems to be 4 times larger than mine! Yes, I found a bug in my work starting near 36.11, am fixing it now. Fix is A/2 = V. I did a complete repair just now. My result got 4x larger and now matches Xiong, thank you Mr. X! So I have now verified Manchester with Xiong for p = 1/2. Question: How does Xiong handle negative frequency in his results? He says absolutely nothing about it! I just have to assume that his spectrum is for f in (-∞.∞) which is the way my stuff works and we are in agreement with the above examples./ 3. AMI Code My pulse is height V, and my result is this <P(ω) > =  Ppulse(ω) a = (1-2p) . (37.11) I go on to say that with p = 1/2 the result is <P(ω) > =  Ppulse(ω) sin2(πx) p = 1/2 (37.12) <P(ω) > =  V2 (1/ω1) sinc2(πx) sin2(πx) p = 1/2 x = ω/ω1 (37.13) Then my result says that <P(f)>/2π = V2 T1/2π sinc2(πx) sin2(πx) x = fT1 <P(f)> = V2 T1 sinc2(πfT1) sin2(πfT1) Xiong approaches this subject on page 35. He calls AMI the name AMI-NRZ pseudoternary since it does in fact have 3 levels. Eventually he quotes this result and we have verification for this difficult case as well!! Hurray! Still to go: NRZI and hold/change linecodes. On 5.13.13 I integrated all these Xiong verifications into spectral doc, and I comment more on the Bennet Davey verification and its factor of 2 and I maintain their credit for the whole discussion.