Построение математической модели формирования радиотелефонного канала для передачи аудиоинформации Текст научной статьи по специальности «Математика»

ЧЕГОДАЕВ1 Николай Игоревич ЦЕБРЕНКО2 Константин Николаевич, кандидат технических наук, доцент

ПОСТРОЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ФОРМИРОВАНИЯ РАДИОТЕЛЕФОННОГО КАНАЛА ДЛЯ ПЕРЕДАЧИ АУДИОИНФОРМАЦИИ

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

Ключевые слова: радиоудлинитель, телефонная, линия, дуплексный радиоканал связи, разнесенные частоты, математическая. модель.

This article covers mathematical model of analogue processes under link channel shaping for radio-extension, unit, formulae development under mathematical model equalization; justification of model actuality quoted; there is short description of device operation principles. The subject of design relevance covered.

Keywords: radio extension unit, telephonic line, full-duplex radio path, diversity technique frequencies, mathematical model.

Согласно концепции устойчивого развития сельских территорий России на период до 2020 г., опубликованной Центром экономической безопасности Российской Федерации [1, с. 30], необходимо обеспечить телефонизацию всех сельских населенных пунктов и объектов социальной сферы. Однако прокладка проводных линий связи невозможна в силу действующих экономических или технических ограничений [2]. Задачи телефонизации в таких условиях могут быть решены с помощью систем расширения проводных линий связи, одним из вариантов решения может быть радиоудлинитель телефонной линии. Устройство включает в себя абонентский и базовый блоки, которые находятся на значительном удалении друг от друга. Устройство работает следующим образом. Абонент поднимает трубку стандартного телефонного аппарата, подключенного к абонентскому блоку, и, услышав предварительную готовность, осуществляет выбор линии, далее базовый блок подключает выбранную телефонную линию, и устройство разрешает набор номера в обычном режиме. При входящем вызове со стороны базового блока происходит сканирование линий связи, и при получении звонка происходит пересылка данных о номере линии, с которой пришел вызов на сторону абонента. Информация о выборе линии пересылается с помощью двухтональных многочастоных сигналов набора номера (DTMF) пониженной частоты, чтобы исключить

несанкционированное использование канала устройства. Устройства формируют канал связи между собой посредством стандартных выпускаемых промышленностью радиостанций [2]. Необходимо построить математическую модель устройства для того, чтобы получить представление о процессах, происходящих при передаче аналоговой информации в радиоудлинителе, посредством математического описания. Математическая модель устройства была заявлена на международной заочной конференции «Технические науки: традиции и инновации» [3], где были приведены тезисы, касающиеся ее структуры. Ниже приводится полное математическое описание и обоснование модели, вывод формул и конкретизация параметров.

Цель моделирования радиоудлинителя состоит в теоретическом изучении процесса установления связи между блоками и передачи аудиоинформации при заданных условиях внешней среды, а также внутренних служебных процессов коммутации внутри блоков системы. Необходимо установить математические и аналитические зависимости между параметрами аудиоинформации, приближенно описывающими поведение системы.

Данная модель актуальна в связи с необходимостью показать общую структуру детерминированных сигналов, передача которых является непосредственной задачей ра-

' — Академия маркетинга и социально-информационных технологий ИМСИТ, г. Краснодар, аспирант;

2 — Академия маркетинга и социально-информационных технологий ИМСИТ, г. Краснодар, директор института информационных технологий и инноваций.

Рис. 1. Структура моделируемой системы

диоудлинителя, принцип действия которого заключается в перенесении свойств удаленной телефонной линии к телефонному аппарату абонента посредством дуплексного радиоканала [2]. Модель позволит спрогнозировать прямые и косвенные последствия реализации заданных способов воздействия на объект с учетом изменений отдельных параметров полезного сигнала и ограничивающих факторов передачи.

Для построения модели рассмотрим, как исследуемый объект взаимодействует с внешней средой. Объект имеет свойство непрерывно усложняться, поэтому будем считать его сложной системой различных компонент, взаимосвязанных друг с другом. Таким образом, имеет место системный подход к созданию модели, при котором необходимо выделить систему Б и внешнюю среду Е [5, с. 21]. При этом модель также является системой Б' = Б'(М), рассматриваемой относительно внешней среды Е [5, с.22]. Структура системы приведена на рис. 1.

На основании данных о внешней среде построим модель и выявим ограничения для построения системы [5, с. 24]. Влияние модели внешней среды Е' носит ограничивающий характер в виде рекомендуемых параметров элементов передачи данных, достаточных для функционирования системы в целях решения поставленных в публикации [2] задач. Расчетные данные внешней среды приводятся ниже. Модель адекватна объекту при условии соответствия исчисленных параметров моделируемой системы реальному объекту. Таким образом, в данном случае, кроме математического, имеет место также натурное моделирование, то есть проведение исследования на реальном объекте с последующей обработкой результатов эксперимента на основе теории подобия [5, с. 36], являющейся основанием теории моделирования [5, с.31].

Перейдем от описания объекта к непосредственному составлению модели. Элемент аналоговой коммутации предназначен для формирования аудиоканала связи [2], что указывает на колебательный характер процессов, происходящих в последнем. Сказанное верно и для элемента аутентификации с той разницей, что аутентификация производится с помощью двухтональных многоча-стоных сигналов (БТМР) [2]. Очевидно, что элементы ау-

Рис. 2. Функциональная схема модели аналоговых процессов описываемой системы

тентификации и аналоговой коммутации в общем случае представляют собой источники колебаний, в пределах модели для них могут быть определены сходные функции, поэтому в качестве основы для количественной и качественной описательной оценки происходящих в устройстве процессов, то есть в качестве модели устройства подходит колебательный ЬС-контур [4, с. 50]. Составим функциональную схему математической модели. ункциональная схема модели аналоговых процессов системы приведена на рис. 2.

Уравнение, описывающее колебательный процесс в одном из простейших электрических колебательных контуров [4, с. 50], изображенных на рис. 2, следующее (1):

Ь

d д

-Сд .

(1)

где Ь — значение индуктивности контура; ё.2д/ё.¥ — дифференциал по времени функции, описывающей значения заряда конденсатора; С — значение емкости контура, д — величина заряда, прикладываемого к обкладкам конденсатора. Произведем преобразование выражения (1) к виду, допускающему анализ процессов изменения величины напряжения и^), а также составим уравнения, описывающие поведение сложного сигнала, модулирующую последовательность и уравнение модулированного сигнала.

По определению

Ла

I =--- [4, с.50],

Л

следовательно, из уравнения (1) получим (2):

д(,)С = I ^ . (2)

т

Так как по физическому смыслу величины С, в любой момент времени допустимо равенство и(Ї) = дЩС, то перепишем (2) в виде (3):

/ ч т di

u(t) = L Т- • dt

(3)

/ ч r di

u (0 = A — , dt

T di2

U2( t ) — L2 —— j dt

/ ч r di3

u3(t) = L3 -T ■ dt

(4)

(5)

(6)

u- (t) - L- — + u2 (t) - L2 — = 0.

1W 1 dt 2W 2 dt

(7)

Сгруппируем переменные выражения (7), следующим образом (8):

u ( t) + u2 (t ) - (L — + L2 — ) = 0 .

1W 2 1 dt 2 dt’

(8)

Для дальнейших преобразований в модели в уравнении (8) обозначим

ион (t) = U1 (t) + U2 (t) и Lo

diocu = L di + L dk dt 1 dt 2 dt

ющее [6, с. 136]:

Uam (t ) = Uc (t ) t1 + mS (t ) ],

(10)

где uam (t) — промодулированный сигнал; uc (t) — несущее колебание; S(t) — информационный сигнал; m — коэффициент модуляции, m = const.

Промодулируем сложный основной сигнал, описанный в выражении (9), подставив переменным в выражения (10) выражения (9) и (6), таким образом, получим следующее уравнение для модулированного сигнала (11):

di„„

(t) = L ^3--( ) 3 dt

1 + mL -

dt

(11)

Уравнение (3) справедливо для каждого из электрических колебательных контуров модели, приведенной на рис. 2, однако каждый из колебательных контуров имеет собственную частоту колебаний, поэтому переменные равенств (3), кроме времени ^ различны для каждого из них, поэтому справедливы равенства (4), (5), (6):

Для детектирования сигнала нужно преобразовать уравне ние (11) к уравнению (12) следующим образом [5]:

di

^ / --— гм /

и -

. (t ) = L3 — + L3 — mLo ,W 3 dt 3 dt o

—t

(12)

Кроме того, необходимо учесть соотношение сигнал-шум [6, с. 175] по формуле (13):

С/Ш = P/2rjW,

(13)

Равенства (4), (5), (6) справедливы соответственно для каждого из трех колебательных контуров, изображенных на рис. 2 сверху вниз.

Для математического моделирования сложного сигнала в системе, состоящего из нескольких частот, необходимо сложить сигналы, выраженные равенствами (4) и (5), таким образом, получим равенство (7):

где С/Ш — соотношение сигнал-шум; Р — мощность сигнала, цШ — коэффициент, учитывающий амплитуду шумов и ширину полосы частот.

Для данной модели учитываем, что приемо-передающий элемент, изображенный на рис. 2, генерирует только гауссовы шумы с амплитудой ц/2.

Мощность сигнала является ограничением, накладываемым внешней средой, и имеет квадратичную зависимость от расстояния [7, с. 178]. Запишем выражение для мощности передаваемого сигнала (14):

PD D Л2

р____ г np n

Получим уравнение основного для данной модели сложного сигнала, который для передачи через модель внешней среды необходимо промодулировать. Уравнение (9) выглядит следующим образом:

и (/) = L . (9)

венУ ’ ОСН

Уравнение (9) представляет в модели сложный аналоговый сигнал, который должен быть промодулирован. Выражение (10) для амплитудно-модулированного сигнала следу-

(14)

4п2 S2

где P — суммарная мощность сигнала; РГ — мощность передатчика в момент приема абонентом сигналов, для данной модели примем Рг = const; Dnp, Dn — коэффициенты, определяющие усиление сигнала в направлении к приемнику и от него, для целей данной модели примем равными единице, считаем Dnp Dn = const; Л — длина волны, аналогично вышесказанному, Л = const; S — расстояние между прие-мо-передающими элементами среды E, изображенными на рис. 2, поскольку расстояние между приемо-передающими подсистемами объектав реальности не изменяется, то считаем S = const.

Подставляя выражение для мощности передаваемого сигнала (14) в выражение (13) получаем уравнение (15) для определения соотношения сигнал-шум модели внешней среды E:

PD D X2

C/ш = г . (15)

8п S2nW

БИХ фильтр, или фильтр с бесконечной импульсной характеристикой, в данной модели предназначен для имитации устройств, ответственных за распознавание двухтональных многочастотных сигналов набора номера, и работает по алгоритму Герцеля, подробно описанному в публикации [8].

Опишем математически фильтр с бесконечной импульсной характеристикой, показанный на рис. 4. Изображенный на рисунке БИХ фильтр работает по алгоритму Герцеля. Для вывода уравнений, описывающих математически процессы, происходящие в данном элементе, необходимо, прежде всего, опереться на уравнение БИХ фильтра первого порядка. Функциональная схема БИХ фильтра спектрального отсчета с номером к приведена на рис. 3.

Выведем уравнение для БИХ фильтра, изображенного на рис. 3, используя публикацию [8]. Покажем также передаточную характеристику фильтра. Запишем выражение (16) для дискретного М-точечного преобразования Фурье значений сигнала в(п) при п є (0...Ы — 1):

N-1

пхк N ,

(16)

п=0

где Б(к) — комплексные отсчеты спектра по Фурье; в(п) — значения отсчетов сигнала; WNХк— поворотные коэффициенты, обладающие следующим свойством:

- к х N N

2п

N-1

г; * х лг х X * (п) х ^

п=0

пхк

N

(17)

Умножив выражение (19) на 1 — к X г 1 с учетом ^к + = 2 X ^(^ X к) [12],

получим передаточную характеристику более сложного БИХ фильтра, работающего в соответствии с алгоритмом Герцеля и позволяющего вычислить спектральный отсчет к с меньшей затратой ресурсов, чем фильтр на рис. 4. Значение (20) вещественного коэффициента а для схемы на рис. 4 [8]:

а = 2 х со8(2 х пх —).

(20)

Для схемы, изображенной на рис. 4, значение Б(к) определяется уравнением (21):

= ЄХр(- І-----X П X к) , при к є (0.М — 1).

N

Умножив выражение (16) на кх ж , получим выражение

(17):

5 (к) -

S(к) = ^к х v(N - 1) - v(N - 2).

(21)

В уравнении (21) промежуточные значения у(Ы — 1) и у(Ы — 2) рассчитываются итерационно по формуле (22):

Раскрыв сумму, выделяя спектральные отсчеты с вынесени- V(г) — ^ (г ) + 2 X йХ V (г 1) \>(г 2) .

ем кХ ^ за скобки, как показано в публикации [7], получаем рекуррентное выражение (18) для вычисления спектрального отсчета к на N — 1 шаге вычисления:

(22)

Поворотные коэффициенты уравнения (21) имеют реальную часть

5 (к) = уг ( к ) = г;к х ж X ( s(r ) + л - (к))

(18)

к

где в(г) — отсчет входного сигнала с номером г.

Данное рекуррентное выражение на шаге с номером г = N — 1 позволяет нам получить спектральный отсчет номером к [8]. Выражение (18) представляет собой разностное уравнение фильтра с бесконечной импульсной характеристикой (БИХ фильтра), структурная схема которого приведена на рис. 4. В выражении (18) обозначим через X^) и У^) элементы 8(г) и уг(к), получим передаточную характеристику (19) И^) БИХ фильтра, показанного на рис. 3:

Я = 008(2 х пх —) и мнимую к

ц-к Я = 8Ш( 2 х п х —) части. ж N

С учетом последних исчисляются реальная —

SR = v(N - 1) х W-к г - V(N - 2)

Н (г) =

К

г(г) = ___________

X(г ) 1-К-к X г

-і •

(19)

и мнимая —

S1 — V( N — 1) X WN кI части спектральных отсчетов.

Таким образом, алгоритм Герцеля сводится к итерационному расчету у(г) по формуле (22), где последние итерационные значения помещаются формулу (21) для окончательного

пересчета спектрального отсчета Б(к). Очевидно, что реальные и мнимые части поворотного коэффициента уравнения (21), а также коэффициент а уравнения (22) опираются на характер сигнала. Конкретизируем форму колебательного сигнала для электромагнитного колебательного процесса, описанного в уравнениях (4) — (12).

Электромагнитные колебания в идеальном колебательном контуре, изображенном на рис. 2, подчиняются гармоническому закону. Уравнение (23), описывающее процесс изменения величины и({), выглядит следующим образом:

u(t) = Um c o + ф) ,

T > PzDnp Dn A2 Un Tm > 8n2S2nW

(27)

dt

U cos( ú t + ф ) < 1

тосн V осн i осн ’

и 0 < m < 1,

имеем системы аналогичных уравнений для сигнала с аналогичной системой ограничений.

Система уравнений (28), описывающая характеристики колебательного сигнала модели во времени, выглядит следующим образом:

иосн О) = ) + Щ(? X

(23)

где u(t) — функция изменения величины напряжения; Um — максимальное (амплитудное) значение u(t); ы — угловая частота, ш = 2nf или ш= 2п /Т; ф — сдвиг фаз.

С учетом (23) сложный колебательный процесс, описанный уравнениями (7) и (8), может быть представлен уравнением (24) следующим образом:

и1 (0 + u2(t) = Um c + У1 ) + Um2 cos (a2t+у2). (24)

Выполним преобразование, аналогичное (9), получим (25):

uо-н ) = U тосн C°SKснt + 9 осн )• (25)

По формуле (10) промодулируем сигнал, описанный уравнением (25), получим (26):

Uam (t) = Um3 CO S(®3t + P3 ) [l + MUm0CH COS (fl^ t + + 0CH )]. (26)

Выше было замечено, что влияние модели внешней среды E' носит ограничивающий характер.

Приведем ограничение (27) распространения сигнала, выведенное по формулам (14) и (15), обозначив максимальное значение полезного сигнала, как Um, паразитного — как Un.

L,

diL„ r di r di2 = L — + L—,

1 dt di

dt

dt

/ ч r di3

u (t) = L3 —

am \ ' 3 7,

dt

1 + mLn

di

dt

(28)

Приведем систему ограничений для системы уравнений (28). Ограничения вызваны присутствием материальных объектов, составляющих внешнюю среду Е системы, непосредственно не входящих в модель [4, с. 99]. Соответствующая системе уравнений (28) система ограничений (29) выглядит следующим образом:

di

L,

от < 1 dt ’

di„„.

1 + mLt 0 < m < 1.

dt

> 0,

(29)

Для конкретизации описания формы колебательного сигнала LC-контура модели приведем следующую систему уравнений (30), описывающих гармоническую составляющую сигнала:

иосн (0=и(0+и2(0,

UmocH C0S(®oJ+9ocH ) = U C0S((0lt+(Pl) + Um2 C0S(ffl2í + f2),

Uocu(t) =moc»C0 S(®ocnt +och ) >

Um (t) = Um3 СФ3 + Рз)[ +mmocн Cs(oJ + Poch )]•

(30)

Отношение полезного сигнала к значению шумов должно быть достаточно велико [6, с. 178]. Выражения 1+>0-dt

1 + ти ео8(ю t + ф ) > 0

тосн \ осн тосн’

истинны при выполнении условий di

Ь „„„ —— < 1 ,

Соответствующая ей система ограничений (31) выглядит следующим образом:

Umocn С+ Vocn ) < 1 1 + m U mocn COs(®oJ + 9ocn ) > 0 .

0 < m < 1,

PD D X2U

U ^ г np n n

8n2 S 2r¡W

(31)

Приведем систему уравнений (32) для описания подсистемы, моделирующей обработку сигнала на принимающей стороне:

иначе возникнет перемодуляция сигнала, то есть избыточная модуляция, демодуляция такого сигнала приведет к его искажениям [6, с. 136].

Коэффициент а для уравнения (22) должен изменяться по закону (20), также при моделировании по алгоритму Герце-ля нужно учесть мнимые и реальные части. Таким образом,

\ v(r) = s(r) + 2 x ax v(r -1) - v (r - 2),

(32)

[ ^ (к) = к х у(И - 1) - - 2).

Система ограничивающих условий (33) для данной системы следующая:

к

а = 2 х 008(2 х пх Ш- кЯ = 008(2 х п х N),

<Ш~к I = 8 т(2 х п х к ) , - (33)

Б* = у(N -1) х Ш-кг - у(N - 2),

Б, = у(N -1) х Ш-к1.

При построении математической модели установлены законы получения и распространения детерминированных сигналов, уравнения (28) — (31), при взаимодействии подсистем системы передачи информации радиоудлинителя телефонной линии [2], влияния внешней среды при уточнении способов управления формирования и передачи сигналов при заданных условиях внешней среды. При подготовке данной публикации было составлено математическое описание колебательных процессов передачи аналоговых сигналов через

систему уравнений (28), которое конкретизировано через описание формы сигнала системами гармонических уравнений (30). Также составлена система ограничений (29) и (31) в зависимости от внешней среды; составлена система уравнений (32) с системой условий (33), описывающие процесс детектирования сигнала при аутентификации устройства. Приведенная математическая модель приближенно описывает процессы передачи информации для устройства интеллектуального радиоудлинителя телефонной линии [2]. Модель адекватна реализации полезной модели радиоудлинителя телефонной линии [2] с некоторой степенью приближения, что позволяет оценить параметры передаваемых сложных сигналов по времени и по форме и вычислить их при различных начальных условиях. Модель позволяет вести расчет по параметрам передаваемых сигналов, дает возможность рассчитывать показатели по комплексным отсчетам спектра сигнала при анализе, прогнозировать прямые и косвенные последствия реализации заданных способов воздействия на объект с учетом изменений отдельных параметров полезного сигнала и ограничивающих факторов передачи.

В настоящее время на радиоудлинитель телефонной линии имеется патент [9]

Литература

1. Концепция устойчивого развития сельских территорий Российской Федерации на период до 2Q2Q года./ Центр экономической безопасности [Электронный ресурс] — режим, доступа: http://econsec.ru/files/28.pdf.

2. Чегодаев Н.И. Проблематика научных исследований в области телефонизации удаленных объектов хозяйственной деятельности [Текст]./ Н.И.Чегодаев. Молодой ученый, 2Q11. — № Я — С. 61 — 64.

3. Чегодаев Н.И. Математическая, модель формирования, радиотелефонного канала для. передачи аудиоинформации [Текст]./ Н.И. Чегодаев, К.Н. Цебренко. Технические науки: традиции и инновации: материалы международной заочной научной конференции (г. Челябинск, январь 2Q12 г.). / Под общ. ред. Г. Д. Axметовой. — Челябинск: Два комсомольца,

2Q12. — С. 4Я — SS.

4. Самарский A.A. Математическое моделирование: Идеи. Методы.. Примеры [Текст]./ A.A.Cамарский, АП.Михайлов — М.: Физматлит, 2QQ9 — 32Q с.

5. Советов Б.Я. Моделирование систем. [Текст]./ Б.Я.Советов, С.АЯковлев — М.: Высшая, школа, 2Q1Q. — 343 с.

6. Read. Richard. The Essence of Communications Theory [Text]./ Richard. Read. — Prentice Hall Europe, 2QQ7. — 3Q4 с.

7. Мельник ЮЛ. Основы, радиотехники и радио-технические устройства [Текст]./ ЮЛ. Мельник, Г.В. Стогов —

М.: Радио, 1Я78. — 368с.

8. Aлгоритм Герцеля. Теория, и практика обработки сигналов [Электронный источник] — режим, доступа: http://www.dsplib.ru/content/goertzel/goertzel.html.

Я. Патент, на полезную модель 111732 РФ, МПК 8 HQ4Q7/24. Радиоудлинитель телефонной линии [Текст]./ Н.И. Чегодаев, К.Н. Цебренко; заявитель и патентообладатель — Негосударственное аккредитованное некоммерческое частное образовательное учреждение высшего профессионального образования. «Aкадемия маркетинга и социально-информационных технологий — ИМСИТ», г. Краснодар. — № 2Q111312S1; заявл. 26.Q7.2Q11; опубл. 2Q.12.2Q11.