Учебное пособие 2026
.pdfВ квадратичных преобразователях КВ вычисляются
|
|
|
|
|
, |
(8.13) |
||
z0,i |
= |
|
u012 ,i + u022 |
,i |
||||
z |
= |
|
|
|
. |
(8.14) |
||
u2 |
+ u2 |
|
||||||
1,i |
|
|
11,i |
12,i |
|
|
|
Рис. 8.9
Решение о переданном символе принимается после поступления элемента сигнала с ОФТ по импульсам тактовой
синхронизации демодулятора: если z0 > z1 , то сдвиг фазы 0, а
иначе он равен π .
Зависимости z0 и z1 от номера текущего периода i представлены на рис. 8.10. Функция z0 показана сплошной линией, а z1 - пунктиром, а прямоугольными импульсами отмечено изменение начальной фазы ϕ . В начале наблюдается переходной процесс заполнения регистров сдвига.
100
Рис. 8.10
Как видно, если сдвиг фаз между колебаниями принимаемого и предшествующего элементов равен π , то максимума по завершении приема элемента достигает величина z1 ,
при этом z0 = 0. Если фазы этих элементов одинаковы, то максимум имеет место в канале z0 , а z1 = 0. Таким образом,
решение о принятом информационном символе принимается по наибольшему значению z0 или z1 .
Неравенство z0 > z1 при подстановке (8.12) – (8.14) можно записать в виде: при
y0,i y0,i−N + y1,i y1,i−N > 0, |
(8.15) |
принимается решение о сдвиге фазы 0, а иначе – π . Соответствующая (8.15) структурная схема демодулятора показана на рис. 8.11.
На рис. 8.12 сплошной линией показана зависимость нормированных значений zi /2N (8.15) от номера текущего
периода i/N без помех при N = 256, пунктиром показано изменение фазы сигнала. В моменты формирования решения результат аналогичен показанному на рис. 8.10.
101
Рис. 8.11
Рис. 8.12
8.4. Некогерентный демодулятор четырехпозиционного сигнала с ОФМ
У сигнала с четырехпозиционной ОФМ сдвиг фаз между соседними символами может принимать значения 0 , π / 2, π и 3π / 2 . Структурная схема демодулятора показана на рис. 8.13.
102
Рис. 8.13
По показанным на рис. 8.13 величинам y0,i , y1,i и за-
держанным на |
N |
|
периодов (на длительность элемента) их |
||||||||||||||||
значениям y0,(i−N) , y1,(i−N) |
в четырех умножителях УМН вы- |
||||||||||||||||||
числяются произведения |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
a = y0,i y0,(i−N), |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
b = y |
|
y |
|
, |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
1,i |
|
1,(i−N) |
|
|
|
|
|
|
|
(8.16) |
||
|
|
|
|
|
|
c |
= y |
|
y |
0,(i−N) |
, |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
1,i |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
d |
= y |
|
y |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,i |
1,(i−N) |
|
|
|
|
|
|
|
|
||
далее в сумматорах СУМ и вычитателях ВЫЧ определя |
ются |
||||||||||||||||||
значения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z0 |
= a+b+c−d = y0,i y0,(i−N) |
+ y1,i y1,(i−N) + y1,i y0,(i−N) − y0,i y1,(i−N), |
(8.17) |
||||||||||||||||
z |
= a+b−c+d |
= y |
0,i |
y |
0,(i |
−N) |
+ y |
y |
|
|
− y |
y |
0,(i−N) |
+ y |
0,i |
y |
, |
||
1 |
|
|
|
1,i |
1,(i−N) |
|
1,i |
|
|
1,(i−N) |
|
|
по которым в решающих устройствах РУ0 и РУ1 (компарато-
103
рах) принимаются решения о значениях принимаемых четырехпозиционных информационных символов sИ со значения-
ми 00, 01, 10, 11.
При отсутствии шума по окончании информационного символа y0,i , y1,i принимают значения
y0,i = 2NScos(ϕ), |
|
(8.18) |
|
y1,i = 2NSsin(ϕ), |
|
||
|
|
||
а для предыдущего символа соответственно |
|
||
y0,i−N = 2NScos(ϕ + a π |
/2), |
(8.19) |
|
y1,i−N = 2NSsin(ϕ + a π /2). |
|||
|
В результате из (8.17)
a = y0,i y0,(i−N) = (2NS)2 cos(ϕ)cos(ϕ + a π /2),b = y1,i y1,(i−N) = (2NS)2 sin(ϕ)sin(ϕ + a π /2),c = y1,i y0,(i−N) = (2NS)2 sin(ϕ)cos(ϕ + a π /2),d = y0,i y1,(i−N) = (2NS)2 cos(ϕ)sin(ϕ + a π /2),
тогда получим
|
|
|
|
|
|
|
|
2 |
|
|
a +b = y0,i y0,(i−N) |
+ y1,i y1,(i−N) |
= (2NS) |
cos(a π /2), |
|||||||
|
||||||||||
|
|
y |
|
− y |
|
y |
= (2NS)2 sin(a π /2). |
|||
c−d = y |
0,(i−N) |
0,i |
||||||||
|
1,i |
|
|
1,(i−N) |
|
|
|
(8.20)
(8.21)
идалее
z0 = a+b+c−d = (2NS)2 cos(a π /2)+(2NS)2 sin(a π /2),(8.22)z1 = a+b−c+ d = (2NS)2 cos(a π /2)−(2NS)2 sin(a π /2).
104
По окончании приема символа величины z0 и z1 не
зависят от начальной фазы сигнала и равны ±(2NS)2 при всех значениях информационных символов.
На рис. 8.14 представлены зависимости от номера i/N принятого периода сигнала нормированных откликов z0 и z1
на входе решающих устройств при отсутствии шума и
N = 256.
Рис. 8.14
Пунктиром показаны изменения фаз элементов сигна-
ла, деленных на 2π . При целочисленных i/N принимаются решения
105
|
00→сдвигфаз 0 |
при |
z0 |
> 0 и z1 > 0, |
|
||||
|
01→сдвигфаз π /2 |
при |
|
z0 |
< 0 и z1 |
> 0, |
|
||
|
|
(8.23) |
|||||||
sИ = |
10→сдвигфаз π |
при |
z |
|
< 0 и z < 0, |
||||
|
0 |
|
|||||||
|
|
|
|
|
1 |
|
|
||
11→сдвигфаз 3π /2 |
при |
|
|
z |
0 |
> 0 и z |
< 0. |
|
|
|
|
|
|
|
|
1 |
|
|
9. ЦИФРОВЫЕ ИЗМЕРИТЕЛИ
9.1. Интегратор
Цифровой интегратор выполняет операцию численного интегрирования сигнала (функции) s(t) на заданном времен-
ном инетервале TИНТ от t −TИНТ до текущего значения t,
t |
|
F(t) = ∫s(τ)dτ , |
(9.1) |
t−TИНТ
Эти устройства применяется в различных устройствах управления и обработки сигналов.
Интервал интегрирования TИНТ разбивается на N −1 отрезков с шагом ∆,
N = |
TИЗМ |
, |
(9.2) |
|
∆ |
|
|
N = 2n , n ‒ целое число.
На интервале TИНТ в моменты ti = ∆ i с помощью АЦП формируются дискретные отсчеты si = s(ti ), i - номер
отсчета, по которой численно определяется интеграл (9.1). Применяются [118, 119] разнообразные методы вычисления определенного интеграла:
‒ прямоугольников
106
|
|
t |
|
N−1 |
|
|
|
|
|
|
F(t) = ∫s(τ)dτ = ∆∑si−k ; |
(9.3) |
|||||
|
|
t−TИНТ |
|
k=0 |
|
|
|
|
‒ трапеций |
|
|
|
|
|
|
||
F(t)= |
∫s(τ)dτ = ∆(si +si−N+1)+∆∑si−k = ∆∑si−k − ∆(si +si−N+1); (9.4) |
|||||||
|
t |
|
N−2 |
|
N−1 |
|
|
|
|
t−TИНТ |
2 |
k=1 |
|
k=0 |
2 |
|
|
‒ парабол |
|
|
|
|
|
|
|
|
|
t |
|
|
N−2 |
|
|
N−2 |
|
F(t) = |
|
|
|
|
2∆ ∑si−2k , (9.5) |
|||
∫s(τ)dτ = ∆ (si + si−N+1)+ 4∆ ∑si−2k+1 + |
||||||||
|
|
|
|
2 |
|
2 |
|
|
t−TИНТ |
3 |
3 |
k=1 |
3 k=1 |
||||
|
|
|
|
|
|
|
и другие.
Метод парабол (Симпсона) обеспечивает наименьшую
погрешность при заданном числе интервалов N −1, которое в простых вычислениях обычно выбирается от 30 до 50. Общим свойством рассмотренных методов является снижение по-
грешности с ростом N , особенно при N >>1000, однако при этом возникают вычислительные трудности и необходимо использовать быстрые цифровые алгоритмы.
Наиболее прост интегратор на основе по формулы прямоугольников (9.3), структурная схема интегратора показана на рис. 9.1.
Рис. 9.1
107
Рассмотрим интегрирование гармонического сигнала с единичной амплитудой и частотой F вида
s(t) = cos(2πFt), |
|
(9.6) |
|||
аналитическое значение интеграла равно |
|
||||
F0(t) = ∫t |
s(τ)dτ = |
∫t |
cos(2πFτ)dτ = |
(9.7) |
|
t−TИНТ |
t−TИНТ |
||||
|
= sin(2πFt)−sin(2πF(t −TИНТ )).
На рис. 9.2 показаны результаты имитационного моделирования работы интегратора методом прямоугольников при
F =1000Гц, N = 2048 (n =11), частоте дискретизации АЦП fКВ =1МГц (∆ =1мкс) длительности интервала интегрирования TИНТ = 2,048мс. На рис. 9.2а сплошной линией показан результат интегрирования F(t) (9.3), а точечной – его точное значение F0(t) (9.7). В начале в течение TИНТ возникает переходной процесс заполнения регистров сдвига, а затем F(t) практически совпадает с F0(t) (рис. 9.2б), график погрешности F(t)− F0(t) показан на рис. 9.2в (менее 0,5%).
Рис. 9.2
Цифровой интегратор (рис. 9.1) для обеспечения малой погрешности требует, чтобы частота дискретизации была много больше частоты сигнала.
108
Значительно лучшие результаты обеспечивает интегратор на основе метода Симпсона. Структурная схема интегратора показана на рис. 9.3. Входной сигнал x(t), поступает на
вход АЦП, который в моменты времени ti , определяемые
ГТИ 3 ( i – порядковый номер), с интервалом времени τ формирует отсчеты xi входного сигнала. В текущий момент ti ,
k > N, обрабатывается (интегрируется) выборка отсчетов объемом
N = 2m+2 +1 (9.8)
(m – целое число). Выбор N (9.8) необходим для организации быстрой вычислительной процедуры в каналах накопления отсчетов КНО0 и КНО1.
Значение интеграла на интервале времени от ti −(N −1)τ до ti вычисляется по формуле парабол (9.5) в со-
ответствии с выражением
ti |
|
|
τ |
2m+1 |
|
2m+1−1 |
|
|
|
|
|
|||
Fi = ∫x(t)dt ≈ |
|
xi + 4∑xi−2 j+1 + 2 ∑xi−2 j + xi−N+1 |
= |
|
||||||||||
ti −(N−1)τ |
|
|
3 |
j=1 |
|
j=1 |
|
|
|
|
(9.9) |
|||
|
|
τ |
|
2m+1 |
|
2m+1 |
|
|
|
|
||||
|
|
|
|
|
|
|
||||||||
|
= |
|
3 |
xi + 4 |
∑xi−2 j+1 + 2∑xi−2 j − xi−N+1 |
|
|
|
||||||
|
|
|
|
j=1 |
|
j=1 |
|
|
|
|
|
|||
с абсолютной погрешностью R, равной |
|
|
|
|
|
|||||||||
|
|
|
|
R ≤ |
|
(T1)5 A |
, |
|
|
(9.10) |
||||
где |
|
|
|
180[2(N −1)]4 |
|
|
|
|||||||
|
|
|
|
|
|
|
dx4 (t) dt4 |
|
. |
|
|
(9.11) |
||
|
|
|
|
A= |
max |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
t [ti −T1,ti ] |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
Структурная схема интегратора показана на рис. 9.3.
109