exponenta event banner

Аналитические выражения и обозначения, используемые в анализе BER

В этом разделе рассматриваются аналитические выражения и обозначения для теоретического анализа, используемые в функциях BER (berawgn, bercoding, berconfint, berfadingberfit, bersync), приложение «Анализ частоты битовых ошибок» и раздел «Методы анализа частоты битовых ошибок».

Общее обозначение

Эта таблица определяет обозначения, используемые в аналитических выражениях в этом разделе.

Описание Примечание
Размер совокупности модуляции

M

Количество битов на символ

k = log2M

Отношение энергии на бит к шуму мощность-спектральная плотность

EbN0

Энергия на символ - шум, мощность - спектральная плотность

EsN0=kEbN0

Частота битовых ошибок (BER)

Свинец

Частота ошибок символов (SER)

Ps

Реальная часть

Re[⋅]

Пол, наибольшее целое число меньше значения, содержащегося в фигурных скобках

В этой таблице описаны термины, используемые для математических выражений в этом разделе.

Функция Математическое выражение
Функция Q

Q (x) =12π∫x∞exp (t2/2) dt

Функция Marcum Q

Q (a, b) =∫b∞texp (t2 + a22) I0 (at) dt

Модифицированная функция Бесселя первого рода порядка

Istart( z) =∑k=0∞ (z/2)

где

Γ (x) =∫0∞e−ttx−1dt

- гамма-функция.

Конфлюэнтная гипергеометрическая функция

F11 (a, c; x) =∑k=0∞ (a) k (c) kxkk!

где символ Почхаммера, (λ) k, определяется как (λ) 0 = 1, (λ) k = λ (λ + 1) (λ + 2) ⋯ (λ + k − 1).

Эта таблица определяет акронимы, используемые в этом разделе.

Акроним Определение
M-PSKМ-арная фазовая манипуляция
DE-M-PSKДифференциально закодированная М-образная фазовая манипуляция
BPSKДвоичная фазовая манипуляция
DE-BPSKДифференциально кодированная двоичная фазовая манипуляция
QPSKЧетвертичная фазовая манипуляция
DE-QPSKДифференциально закодированная квадратурная фазовая манипуляция
OQPSKСмещение квадратурной фазовой манипуляции
DE-OQPSKДифференциально кодированная квадратурная фазовая манипуляция со смещением
M-DPSKМ-арная дифференциальная фазовая манипуляция
M-PAMM-арная импульсная амплитудная модуляция
M-QAMM-арная квадратурная амплитудная модуляция
М-ФСКM-арная частотная манипуляция
MSKМинимальная смена клавиш
M-CPFSKМ-арная непрерывная фазовая частотная манипуляция

Аналитические выражения, используемые в berawgn Приложение для анализа частоты ошибок функций и битов

Эти разделы охватывают основные аналитические выражения, используемые в berawgn функция и приложение Bit Error Rate Analysis.

M-PSK

Из уравнения 8.22 в [2],

Ps=1π∫0 (M 1) δ/Mexp (kEbN0sin2 [λ/M] sin2

Это выражение аналогично, но не строго равно, точному BER (от [4] и уравнению 8.29 от [2]):

Pb = 1k (∑i=1M/2 (wi ') Pi)

где wi '= wi + wM i, wM/2' = wM/2, wi - вес Хэмминга битов, назначенных символу i,

Pi=12π (1− (2i−1) М) exp (−kEbN0sin2 [(2i−1) π/M] sin2θ) dθ−12π (1− (2i+1) М) exp (−kEbN0sin2 [(2i+1) π/M] sin2θ) dθ

Для M-PSK с M = 2, в частности BPSK, применяется это уравнение 5.2-57 из [1]:

Ps = Pb = Q (2EbN0)

Для M-PSK с M = 4, в частности QPSK, применяются следующие уравнения 5.2-59 и 5.2-62 из [1]:

Ps = 2Q (2EbN0) [1 12Q (2EbN0)] Pb = Q (2EbN0)

DE-M-PSK

Для DE-M-PSK с M = 2, в частности DE-BPSK, применяется это уравнение 8.36 из [2]:

Ps = Pb = 2Q (2EbN0) 2Q2 (2EbN0)

Для DE-M-PSK с M = 4, в частности DE-QPSK, применяется это уравнение 8.38 из [2]:

Ps = 4Q (2EbN0) 8Q2 (2EbN0) + 8Q3 (2EbN0) 4Q4 (2EbN0)

Из уравнения 5 в [3],

Pb = 2Q (2EbN0) [1 Q (2EbN0)]

OQPSK

Для OQPSK используйте те же вычисления BER и SER, что и для QPSK в [2].

DE-OQPSK

Для OQPSK используйте те же вычисления BER и SER, что и для DE-QPSK в [3].

M-DPSK

Для M-DPSK применяется это уравнение 8.84 из [2]:

Ps = sin (λ/M) 2π∫−π/2π/2exp ((kEb/N0) (1 cos (λ/M) cos

Это выражение аналогично, но не строго равно точному BER (от [4]):

Pb = 1k (∑i=1M/2 (wi ') Ai)

где wi '= wi + wM i, wM/2' = wM/2, wi - вес Хэмминга битов, назначенных символу i,

Ai=F ((2i+1) πM) −F ((2i−1) πM) F (ψ) = −sinψ4π ∫−π/2π/2exp (−kEb/N0 (1−cosψcost)) 1−cosψcostdt

Для М-ДФСК с М = 2 применяется это уравнение 8,85 из [2]:

Pb = 12exp (EbN0)

M-PAM

Из уравнений 8.3 и 8.7 в [2] и уравнений 5.2-46 в [1],

Ps = 2 (M 1M) Q (6M2 − 1kEbN0)

От [5],

Pb=2Mlog2M×∑k=1log2M ∑i=0 (1 2 k) M 1 {(1) i2k 1M (2k−1−⌊i2k−1M+12 ) Q ((2i + 1) 6log2MM2 − 1EbN0)}

M-QAM

Для квадратного M-QAM k = log2M является четным, поэтому применяется уравнение 8.10 из [2] и уравнения 5.2-78 и 5.2-79 из [1]:

Ps = 4M 1MQ (3M 1kEbN0) 4 (M 1M) 2Q2 (3M − 1kEbN0)

От [5],

Pb=2Mlog2M×∑k=1log2M ∑i=0 (1 2 k) M 1 {(1) i2k 1M (2k−1−⌊i2k−1M+12 ) Q ((2i + 1) 6log2M2 (M − 1) EbN0)}

Для прямоугольной (не квадратной) M-QAM k = log2M нечётно, M = I × J, I = 2k 12, а J = 2k + 12. Так что,

Ps = 4IJ 2I 2JM × Q (6log2 (IJ) (I2 + J2 2) EbN0) 4M (1 + IJ I J) Q2 (6log2 (IJ) (I2 + J2 − 2) EbN0)

От [5],

Pb = 1log2 (IJ) (∑k=1log2IPI (k) +∑l=1log2JPJ (l))

где

PI (k) =2I∑i=0 (1 2 k) I 1 {(1) i2k 1I (2k−1−⌊i2k−1I+12 ) Q ((2i + 1) 6log2 (IJ) I2 + J2 − 2EbN0)}

и

PJ (k) =2J∑j=0 (1−2−l) J−1 {(−1) ⌊j2l−1J  (2l−1−  j2l−1J+12 ) Q ((2j+1) 6log2 (IJ) I2+J2−2EbN0)}

Ортогональный M-FSK с когерентным обнаружением

Из уравнения 8.40 в [2] и уравнения 5.2-21 в [1],

Ps=1−∫−∞∞[Q (q 2kEbN0)] M 112.dexp (q22) dqPb = 2k − 12k − 1Ps

Неоргональные 2-FSK с когерентным обнаружением

Для M = 2 применяются уравнения 5.2-21 в [1] и 8.44 в [2]:

Ps=Pb=Q (Эб (1−Re [ρ]) N0)

start- комплексный коэффициент корреляции, такой, что:

ρ=12Eb∫0Tbs˜1 (t) s˜2 * (t) dt

где s˜1 (t) и s˜2 (t) - сложные низкочастотные сигналы, и

Eb=12∫0Tb's˜1 (t) |2dt=12∫0Tb's˜2 (t) | 2dt

Например, с

s˜1 (t) =  2EbTbej2.df1t, s˜2 (t) = 2EbTbej2.df2t

тогда

ρ=12Eb∫0Tb2EbTbej2πf1t2EbTbe−j2πf2tdt=1Tb∫0Tbej2π (f1 f2) tdt = sin (ΔfTb) ¼ ΔfTbejāΔft

где Δf = f1 − f2.

Из уравнения 8.44 в [2],

    Ре [ρ] = Ре [грех (πΔfTb) πΔfTbejπΔft] =sin (πΔfTb) πΔfTbcos (πΔfTb) =sin (2πΔfTb) 2πΔfTb⇒Pb=Q (Эб (1−sin (2πΔfTb) / (2πΔfTb)) N0)

где h = ΔfTb.

Ортогональный M-FSK с некогерентным обнаружением

Из уравнения 5.4-46 в [1] и уравнения 8.66 в [2],

Ps=∑m=1M−1 (1) m + 1 (M 1 m) 1 m + 1exp [мм + 1kEbN0] Pb = 12MM − 1Ps

Неоргональные 2-FSK с некогерентным обнаружением

Для M = 2 это уравнение 5.4-53 из [1] и это уравнение 8.69 из [2] применяются:

Ps = Pb = Q (a, b) 12exp (a + b2) I0 (ab)

где

a = Eb2N0 (1 1  | 

Предварительно закодированная MSK с когерентным обнаружением

Используйте те же вычисления BER и SER, что и для BPSK.

Дифференциально кодированная MSK с когерентным обнаружением

Используйте те же вычисления BER и SER, что и для DE-BPSK.

MSK с некогерентным обнаружением (оптимальный блок за блоком)

Верхняя граница частоты ошибок из уравнений 10.166 и 10.164 в [6])

Ps=Pb≤12[1−Q (b1, a1) + Q (a1, b1)] + 14 [1 Q (b4, a4) + Q (a4, b4)] + 12e − EbN0

где

a1=EbN0 (1−3−4/π24), b1=EbN0 (1+3−4/π24) a4=EbN0 (1−1−4/π2), b4=EbN0 (1+1−4/π2)

Когерентное обнаружение CPFSK (оптимальный поблочный режим)

Нижняя граница частоты ошибок (из уравнения 5.3-17 в [1]) равна

Ps > KδminQ (EbN0δmin2)

Верхняя граница частоты ошибок:

δmin2>min1≤i≤M−1{2i (1 sinc (2ih))}

где h - индекс модуляции, а Kδmin - количество трактов с минимальным расстоянием.

Pb≅Psk

Аналитические выражения, используемые в berfading Приложение для анализа частоты ошибок функций и битов

В этом разделе рассматриваются основные аналитические выражения, используемые в berfading и приложение Bit Error Rate Analysis.

Примечание

В этой таблице описаны дополнительные обозначения, используемые в аналитических выражениях в этом разделе.

Описание Примечание
Мощность амплитуды замирания rΛ = E [r2], где E[⋅] обозначает статистическое ожидание
Количество разнесенных ветвей

L

Отношение сигнал/шум (SNR) на символ на ветвь

γ¯l = (ΩlEsN0)/L = (ΩlkEbN0)/L

Для идентично распределенных разветвлений разнесения

γ¯ = (ΩkEbN0)/L

Функции генерации момента для каждой ветви разнесения

Для релеевских каналов замирания:

Mγ l (s) = 11

Для каналов замирания Rician:

Mγ l (s) = 1 + K1 + K sγ-1 [Ksγ-1 (1 + K) − sγ-1]

K - отношение энергии в зеркальном компоненте к энергии в диффузном компоненте (линейный масштаб).

Для идентично распределенных разветвлений разнесения Mγ 1 (s) = (s) для всех l.

Эта таблица определяет дополнительные сокращения, используемые в этом разделе.

Акроним Определение
MRCОбъединение максимального отношения
EGCОбъединение с равным коэффициентом усиления

M-PSK с MRC

Из уравнения 9.15 в [2],

Ps=1π∫0 (M 1) π/M∏l=1LMγl (sin2 (λ/M) sin2

Из [4] и [2],

Pb = 1k (∑i=1M/2 (wi ') P _ i)

где wi '= wi + wM i, wM/2' = wM/2, wi - вес Хэмминга битов, назначенных символу i,

P¯i=12π∫0π (1 (2i 1 )/M) ∏l=1LMγl (− 1sin2startsin2 (2i 1) securityM) d −12π∫0π (1 (2i + 1 )/M) ∏l=1LMγl (1sin2startsin2 (2i + 1) δM) d

Для особого случая релеевского замирания с М = 2 (из уравнений C-18 и C-21 и таблицы C-1 в [6]),

Pb=12[1−μ∑i=0L−1 (2ii) (1 мк24) i]

где

μ=γ¯γ¯+1

Если L = 1, то:

Pb = 12 [1 γ fw γ β + 1]

DE-M-PSK с MRC

Для М = 2 (из уравнений 8.37 и 9.8-9.11 в [2]),

Ps=Pb=2π∫0π/2∏l=1LMγl (1sin2start) dθ−2π∫0π/4∏l=1LMγl (− 1sin2start) d

M-PAM с MRC

Из уравнения 9.19 в [2],

Ps = 2 (M 1) Mπ∫0π/2∏l=1LMγl (3/( M2 1) sin2

Из [5] и [2],

Pb=2πMlog2M×∑k=1log2M ∑i=0 (1 2 k) M 1 {(1) i2k 1M (2k−1−⌊i2k−1M+12 ) ∫0π/2∏l=1LMγl ((2i + 1) 23/( M2 − 1) sin2

M-QAM с MRC

Для квадратного M-QAM k = log2M является четным (уравнение 9.21 в [2]),

Ps=4π (1−1M) ∫0π/2∏l=1LMγl (−3/(2 (M−1)) sin2θ) dθ−4π (1−1M) 2∫0π/4∏l=1LMγl (−3/(2 (M−1)) sin2θ) dθ

Из [5] и [2]:

Pb=2πMlog2M×∑k=1log2M ∑i=0 (1 2 k) M 1 {(1) i2k 1M (2k−1−⌊i2k−1M+12 ) ∫0π/2∏l=1LMγl (− (2i + 1) 23/( 2 (M − 1)) sin2

Для прямоугольной (некварной) M-QAM k = log2M нечётно, M = I × J, I = 2k 12, J = 2k + 12, γ ve l = Ωllog2 (IJ) EbN0,

Ps=4IJ−2I−2JMπ 0π/2∏l=1LMγl (−3/(I2+J2−2) sin2θ) dθ−4Mπ (1+IJ−I−J) ∫0π/4∏l=1LMγl (−3/(I2+J2−2) sin2θ) dθ

Из [5] и [2],

Pb=1log2 (IJ) (∑k=1log2IPI (k) + ∑l=1log2JPJ (l)) ПИ (k) =2Iπ i=0 (1−2−k) I−1 {(−1) ⌊i2k−1I  (2k−1−  i2k−1I+12 ) ∫0π/2∏l=1LMγl ((2i+1) 23/(I2+J2−2) sin2θ) } PJ (k) =2Jπ j=0 (1−2−l) J−1 {(−1) ⌊j2l−1J  (2l−1−  j2l−1J+12 ) ∫0π/2∏l=1LMγl (− (2j+1) 23/(I2+J2−2) sin2θ) dθ}

M-DPSK с EGC после обнаружения

Из уравнения 8.165 в [2],

Ps = sin (λ/M) 2π∫−π/2π/21[1−cos (λ/M) cosθ]∏l=1LMγl ([1 cos (λ/M) cos

Из [4] и [2],

Pb = 1k (∑i=1M/2 (wi ') A/i)

где wi '= wi + wM i, wM/2' = wM/2, wi - вес Хэмминга битов, назначенных символу i,

A¯i=F¯ ((2i+1) πM) −F¯ ((2i−1) πM) F¯ (ψ) = −sinψ4π ∫−π/2π/21 (1−cosψcost) ∏l=1LMγl ((1−cosψcost)) dt

Для особого случая релеевского замирания с M = 2 и L = 1 (уравнение 8.173 из [2]),

Pb = 12 (1 + γ

Ортогональное 2-FSK, когерентное обнаружение с MRC

Из уравнения 9.11 в [2],

Ps=Pb=1π∫0π/2∏l=1LMγl (1/2sin2

Для особого случая релеевского замирания (уравнения 14.4-15 и 14.4-21 в [1]),

Ps = Pb = 12L (1 γ 2 + γ p) L∑k=0L−1 (L 1 + kk) 12k (1 + γ 2 + γ b) k

Неоргональное 2-FSK, когерентное обнаружение с MRC

Из уравнений 9.11 и 8.44 в [2],

Ps=Pb=1π∫0π/2∏l=1LMγl ((1 Re [

Для особого случая релеевского замирания с L = 1 (уравнения 20 в [8] и 8.130 в [2]),

Ps = Pb = 12 [1 γ пр (1 Re [

Ортогональный M-FSK, некогерентное обнаружение с EGC

Для релеевского замирания из уравнения 14.4-47 в [1],

Ps​=1−∫0∞1 (1 + γ) L (L 1)! UL 1e U1 + γ (1−e−U∑k=0L−1Ukk!) M − 1dUPb = 12MM − 1Ps

Для замирания Rician из уравнения 41 в [8],

Ps=∑r=1M−1 (1) r + 1e LKγ пр/( 1 + γ пр) (r (1 + γ пр) + 1) L (M 1r) ∑n=0r (L 1) βnrΓ (L + n) Γ (L) [1 + γ пр rr + 1 + пр] nF11 (L + n, L;

где

γ¯r=11+Kγ¯βnr=∑i=n− (L 1) nβi (r 1) (n i)! I [0, (r 1) (L 1)] (i) β00 = β0r = 1βn1 = 1/n! β1r = r

и I [a, b] (i) = 1, если a≤i≤b, и 0 в противном случае.

Неоргональное 2-FSK, некогерентное обнаружение без разнесения

Из уравнения 8.163 в [2],

Ps=Pb=14π −ππ1−ς21+2ςsinθ +ς2Mγ (−14 (1+1−ρ2) (1+2ςsinθ +ς2)) dθ

где

ς=1−1−ρ21+1−ρ2

Аналитические выражения, используемые в bercoding Приложение для анализа частоты ошибок функций и битов

В этом разделе рассматриваются основные аналитические выражения, используемые в bercoding и приложение Bit Error Rate Analysis.

Общее обозначение

В этой таблице описаны дополнительные обозначения, используемые в аналитических выражениях в этом разделе.

ОписаниеПримечание
Отношение энергии на информацию бит к шуму мощность-спектральная плотность

γ b = EbN0

Длина сообщения

K

Длина кода

N

Кодовая скорость

Rc = KN

Блочное кодирование

В этом разделе описывается конкретная нотация для выражений блочного кодирования, где dmin - минимальное расстояние кода.

Мягкое решение

Для BPSK, QPSK, OQPSK, 2-PAM, 4-QAM и предварительно кодированного MSK применяется уравнение 8.1-52 в [1]),

Pb≤12 (2K 1) Q (2γbRcdmin)

Для DE-BPSK, DE-QPSK, DE-OQPSK и DE-MSK,

Pb≤12 (2K 1) [2Q (2γbRcdmin) [1 Q (2γ bRcdmin)]]

Для когерентного обнаружения BFSK применяются уравнения 8.1-50 и 8.1-58 в [1],

Pb≤12 (2K 1) Q (γ bRcdmin)

Для некогерентного обнаружения квадратного закона BFSK применяются уравнения 8.1-65 и 8.1-64 в [1],

Pb≤122K−122dmin−1exp (12γbRcdmin) ∑i=0dmin−1 (12γbRcdmin) i1i!∑r=0dmin−1−i (2dmin − 1r)

Для DPSK,

Pb≤122K−122dmin−1exp (γ bRcdmin) ∑i=0dmin−1 (γ bRcdmin) i1i!∑r=0dmin−1−i (2dmin − 1r)

Тяжелое решение

Для общего линейного блочного кода применяются уравнения 4.3 и 4.4 в [9] и 12.136 в [6],

Pb≤1N∑m=t+1N (m + t) (Нм) pm (1 p) N−mt=⌊12 (дмин − 1) ⌋

Для кода Хэмминга применяются уравнения 4.11 и 4.12 в [9] и 6.72 и 6.73 в [7]

Pb≈1N∑m=2Nm (Нм) pm (1 р) N m = p − p (1 − р) N − 1

Для расширенного кода Голея (24,12) применяются уравнения, 4,17 в [9] и 12,139 в [6]:

Pb≤124∑m=424βm (24 м) пм (1 р) 24 − м

где βm - среднее число ошибок символов канала, которые остаются в исправленном N-кортежном формате, когда канал вызвал m ошибок символов (см. таблицу 4.2 в [9]).

Для кода Рида - Соломона с N = Q 1 = 2q − 1,

Pb≈2q−12q−11N∑m=t+1Nm (Нм) (Ps) m (1 Ps) N − m

Для FSK применяются уравнения 4.25 и 4.27 в [9], 8.1-115 и 8.1-116 в [1], 8.7 и 8.8 в [7] и 12.142 и 12.143 в [6],

Pb≈1q1N∑m=t+1Nm (Нм) (Ps) m (1 Ps) N − m

в противном случае, если log2Q/log2M = q/k = h, где h - целое число (уравнение 1 в [10]),

Ps = 1 (1 − s) h

где s - SER в некодированном канале AWGN.

Например, для BPSK M = 2 и Ps = 1 (1 − s) q, в противном случае Ps задается таблицей 1 и уравнением 2 в [10].

Сверточное кодирование

В этом разделе описывается конкретная нотация для выражений сверточного кодирования, где dfree - свободное расстояние кода, а ad - количество путей расстояния d от пути всех нулей, который впервые сливается с путем всех нулей.

Мягкое решение

Из уравнений 8.2-26, 8.2-24 и 8.2-25 в [1] и 13.28 и 13.27 в [6] применяются:

Pb<∑d=dfree∞adf d) P2 d)

Передаточная функция задается

T (D, N) =∑d=dfree∞adDdNf (d) dT (D, N) dN'N=1=∑d=dfree∞adf (d) Dd

где f (d) - экспонента N как функции d.

Это уравнение дает результаты для BPSK, QPSK, OQPSK, 2-PAM, 4-QAM, предварительно кодированных MSK, DE-BPSK, DE-QPSK, DE-OQPSK, DE-MSK, DPSK и BFSK:

P2 (d) = Pb 'EbN0 = γ bRcd

где Pb - BER в соответствующем некодированном канале AWGN. Например, для BPSK (уравнение 8.2-20 в [1]),

P2 (d) = Q (2γbRcd)

Тяжелое решение

Из уравнений 8.2-33, 8.2-28 и 8.2-29 в [1] и 13.28, 13.24 и 13.25 в [6] применяются,

Pb<∑d=dfree∞adf d) P2 d)

Когда d является нечетным,

P2 (d) =∑k= (d + 1 )/2d (dk) pk (1 − p) d − k

и когда d является четным,

P2 (d) =∑k=d/2+1d (dk) pk (1 p) d k + 12 (dd/2) pd/2 (1 − p) d/2

где p - частота битовых ошибок (BER) в некодированном канале AWGN.

Аналитические выражения, используемые в bersync Приложение для анализа частоты ошибок функций и битов

В этом разделе рассматриваются основные аналитические выражения, используемые в bersync и приложение Bit Error Rate Analysis.

Ошибка синхронизации синхронизации

Для вычисления BER для системы связи с ошибкой синхронизации времени, bersync функция использует эту формулу из [13]:

14πσ ∫− ∞∞ exp (−ξ22σ2) ∫2R (1−2 |) ∞exp (−x22) dxdξ + 122π ∫ 2R∞exp (−x22) дуплекс

где λ - ошибка синхронизации, а R - линейное значение Eb/N0.

Ошибка синхронизации синхронизации

Для вычисления BER для системы связи с ошибкой синхронизации несущей, bersync функция использует эту формулу из [13]:

1πσ 0∞exp (−ϕ22σ2) ∫2Rcosϕ ∞ exp (−y22) dydϕ

где λ - фазовая погрешность R - линейное значение Eb/N0.

См. также

Приложения

Функции

Связанные темы