Say that stationary zero-mean Gaussian source x(m)x(m) has
autocorrelation rx(0)=1rx(0)=1, rx(1)=ρrx(1)=ρ, and rx(k)=0rx(k)=0 for k>1k>1.
For a bit rate of R bits per sample, uniformly-quantized PCM
implies a mean-squared reconstruction error of
σ
r
2
|
PCM
=
Δ
2
12
|
Δ
=
2
x
max
/
L
L
=
2
R
=
1
3
x
max
2
σ
x
2
︸
γ
x
σ
x
2
2
-
2
R
=
γ
x
σ
x
2
2
-
2
R
.
σ
r
2
|
PCM
=
Δ
2
12
|
Δ
=
2
x
max
/
L
L
=
2
R
=
1
3
x
max
2
σ
x
2
︸
γ
x
σ
x
2
2
-
2
R
=
γ
x
σ
x
2
2
-
2
R
.
(5)For transform coding, say we choose linear transform
T=t0tt1t=12111-1T=t0tt1t=12111-1
(6)Setting x(m)=x(2m)x(2m-1)tx(m)=x(2m)x(2m-1)t and y(m)=Tx(m)y(m)=Tx(m),
we find that the transformed coefficients have variance
σy02=E{|t0tx(m)|2}=12E|x(2m)+x(2m-1)|2=122rx(0)+2rx(1)=1+ρσy02=E{|t0tx(m)|2}=12E|x(2m)+x(2m-1)|2=122rx(0)+2rx(1)=1+ρ
(7)σy12=E{|t1tx(m)|2}=12E|x(2m)-x(2m-1)|2=122rx(0)-2rx(1)=1-ρσy12=E{|t1tx(m)|2}=12E|x(2m)-x(2m-1)|2=122rx(0)-2rx(1)=1-ρ
(8)and using uniformly-quantized PCM on each coefficient we get
mean-squared reconstruction errors
σq02=(1+ρ)γx2-2R0σq02=(1+ρ)γx2-2R0
(9)σq12=(1-ρ)γx2-2R1.σq12=(1-ρ)γx2-2R1.
(10)We use the same quantizer performance factor γx as before
since linear operations preserve Gaussianity.
For orthogonal matrices T, i.e., T-1=TtT-1=Tt,
we can show that the mean-squared reconstruction error σr2
equals the mean-squared quantization error:
σr2:=1N∑k=0N-1Ex˜(Nm-k)-x(Nm-k)2(hereN=2)=1NE∥x˜(m)-x(m)∥2=1NE∥T-1y˜(m)-x(m)∥2=1NE∥T-1y(m)+q(m)-x(m)∥2=1NE∥T-1Tx(m)+T-1q(m)-x(m)∥2=1NE∥T-1q(m)∥2=1NEqt(m)(T-1)tT-1︸Iq(m)=1NE∥q(m)∥2=1N∑k=0N-1σqk2.σr2:=1N∑k=0N-1Ex˜(Nm-k)-x(Nm-k)2(hereN=2)=1NE∥x˜(m)-x(m)∥2=1NE∥T-1y˜(m)-x(m)∥2=1NE∥T-1y(m)+q(m)-x(m)∥2=1NE∥T-1Tx(m)+T-1q(m)-x(m)∥2=1NE∥T-1q(m)∥2=1NEqt(m)(T-1)tT-1︸Iq(m)=1NE∥q(m)∥2=1N∑k=0N-1σqk2.
(11)Since our 2×22×2 matrix is indeed orthogonal, we have mean-squared
reconstruction error
σr2|TC=12(1+ρ)γx2-2R0+(1-ρ)γx2-2R1σr2|TC=12(1+ρ)γx2-2R0+(1-ρ)γx2-2R1
(12)at bit rate of R0+R1R0+R1 bits per two samples.
Comparing TC to PCM at equal bit rates (i.e. R0+R1=2RR0+R1=2R),
σr2|TCσr2|PCM=12(1+ρ)γx2-2R0+(1-ρ)γx2-2(2R-R0)γx2-2R=(1+ρ)22(R-R0)-1+(1-ρ)22(R0-R)-1.σr2|TCσr2|PCM=12(1+ρ)γx2-2R0+(1-ρ)γx2-2(2R-R0)γx2-2R=(1+ρ)22(R-R0)-1+(1-ρ)22(R0-R)-1.
(13)Figure 2 shows that (i) allocating a higher bit rate to the quantizer
with stronger input signal can reduce the average reconstruction error
relative to PCM, and (ii) the gain over PCM is higher when the input
signal exhibits stronger correlation ρ.
Also note that when R0=R1=RR0=R1=R, there is no gain over PCM—a
verification of the fact that σr2=σq2σr2=σq2 when
T is orthogonal.