bersync

Bit error rate (BER) for imperfect synchronization

Syntax

ber = bersync(EbNo,timerr,'timing') ber = bersync(EbNo,phaserr,'carrier')

Description

ber = bersync(EbNo,timerr,'timing') returns the BER of uncoded coherent binary phase shift keying (BPSK) modulation over an additive white Gaussian noise (AWGN) channel with imperfect timing. The normalized timing error is assumed to have a Gaussian distribution. EbNo is the ratio of bit energy to noise power spectral density, in dB. If EbNo is a vector, the output ber is a vector of the same size, whose elements correspond to the different E_b/N₀ levels. timerr is the standard deviation of the timing error, normalized to the symbol interval. timerr must be between 0 and 0.5.

ber = bersync(EbNo,phaserr,'carrier') returns the BER of uncoded BPSK modulation over an AWGN channel with a noisy phase reference. The phase error is assumed to have a Gaussian distribution. phaserr is the standard deviation of the error in the reference carrier phase, in radians.

Examples

collapse all

Calculate Bit Error Rate (BER) for Imperfect Synchronization

Open Live Script

The code below computes the BER of coherent BPSK modulation over an AWGN channel with imperfect timing. The example varies both EbNo and timerr. (When timerr assumes the final value of zero, the bersync command produces the same result as berawgn(EbNo,'psk',2).)

EbNo = [4 8 12];
timerr = [0.2 0.07 0];
ber = zeros(length(timerr),length(EbNo));
for ii = 1:length(timerr)
    ber(ii,:) = bersync(EbNo,timerr(ii),'timerr');
end

Display result using scientific notation.

format short e; ber

ber = 3×3

   5.2073e-02   2.0536e-02   1.1160e-02
   1.8948e-02   7.9757e-04   4.9008e-06
   1.2501e-02   1.9091e-04   9.0060e-09

Switch back to default notation format.

format;

Limitations

The numerical accuracy of this function's output is limited by

Approximations in the analysis leading to the closed-form expressions that the function uses
Approximations related to the numerical implementation of the expressions

You can generally rely on the first couple of significant digits of the function's output.

Limitations Related to Extreme Values of Input Arguments

Inherent limitations in numerical precision force the function to assume perfect synchronization if the value of timerr or phaserr is very small. The table below indicates how the function behaves under these conditions.

Condition	Behavior of Function
`timerr < eps`	`bersync(EbNo,timerr,'timing')` defined as `berawgn(EbNo,'psk',2)`
`phaserr < eps`	`bersync(EbNo,phaserr,'carrier')` defined as `berawgn(EbNo,'psk',2)`

Algorithms

This function uses formulas from [3].

When the last input is 'timing', the function computes

$\frac{1}{4 π σ} \int_{- \infty}^{\infty} \exp (- \frac{ξ^{2}}{2 σ^{2}}) \int_{\sqrt{2 R} (1 - 2 | ξ |)}^{\infty} \exp (- \frac{x^{2}}{2}) d x d ξ + \frac{1}{2 \sqrt{2 π}} \int_{\sqrt{2 R}}^{\infty} \exp (- \frac{x^{2}}{2}) d x$

where σ is the timerr input and R is the value of EbNo converted from dB to a linear scale.

When the last input is 'carrier', the function computes

$\frac{1}{π σ} \int_{0}^{\infty} \exp (- \frac{ϕ^{2}}{2 σ^{2}}) \int_{\sqrt{2 R} \cos ϕ}^{\infty} \exp (- \frac{y^{2}}{2}) d y d ϕ$

where σ is the phaserr input and R is the value of EbNo converted from dB to a linear scale.

Alternatives

As an alternative to the bersync function, invoke the BERTool GUI (bertool) and use the Theoretical tab.

References

[1] Jeruchim, Michel C., Philip Balaban, and K. Sam Shanmugan, Simulation of Communication Systems, Second Edition, New York, Kluwer Academic/Plenum, 2000.

[2] Sklar, Bernard, Digital Communications: Fundamentals and Applications, Second Edition, Upper Saddle River, NJ, Prentice-Hall, 2001.

[3] Stiffler, J. J., Theory of Synchronous Communications, Englewood Cliffs, NJ, Prentice-Hall, 1971.

Documentation