Б.К. Чостковский
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ФОРМИРОВАНИЯ ОБОБЩЕННЫХ ПАРАМЕТРОВ
КАЧЕСТВА НЕРЕГУЛЯРНЫХ КАБЕЛЕЙ СВЯЗИ В СТОХАСТИЧЕСКОЙ ПОСТАНОВКЕ
Определены обобщенные параметры качества кабелей связи и построены детерминированная и
стохастическая модели, отражающие взаимосвязь обобщенных параметров кабелей и вероятностных характеристик параметров, формируемых в ходе управляемых технологических процессов.
Выбор обобщенных параметров и постановка задачи. Проводные каналы связи создаются с помощью направляющих систем, в качестве которых используются двухпроводные цепи, коаксиальные и симметричные [1]. Промышленность изготавливает кабели магистральной, зоновой или сельской связи. Для соединения антенных систем с радиочастотной аппаратурой выпускаются как радиочастотные, так и коаксиальные и симметричные кабели. Связь с движущимися объектами осуществляется с помощью коаксиальных радиочастотных излучающих кабелей [2]. Для построения локальных вычислительных систем используются симметричные радиочастотные кабели с парной скруткой (LAN-кабели).
Реальные кабели связи изготавливаются на технологическом оборудовании, подверженном влиянию случайных возмущающих воздействий. Поэтому, несмотря на наличие встроенных локальных систем автоматической стабилизации (САС), большинство технологических режимных параметров и контролируемых “на проход” параметров кабельного изделия имеют непостоянные по длине кабеля значения и являются случайными функциями координаты длины кабеля х. Линии с переменными по длине параметрами принято называть нерегулярными. Как правило, нерегулярность параметров кабеля по длине невелика и составляет единицы, а чаще десятые и сотые доли процента от номинальных значений. Следовательно, кабели являются слабонерегулярными линиями. Тем не менее, слабое непостоянство размеров и свойств кабеля по длине решительным образом влияет на характер распространения электромагнитных волн и параметры кабеля, среди которых необходимо выделить обобщенный параметр, количественно характеризующий уровень нерегулярности кабеля и, соответственно, его качество.
Количество организуемых по кабелю каналов связи и длина усилительного участка (или тракта) определяются первичными и вторичными параметрами передачи. В связи с тем, что передача сигнала в кабелях связи осуществляется по многим двухпроводным цепям, расположенным в непосредственной близости друг от друга, наблюдается взаимное влияние между цепями, приводящее к появлению переходных шумов и разговоров, которые характеризуются с помощью первичных и вторичных параметров влияния [1, 3].
Как следует из общего решения уравнений распространения электромагнитных волн, в линиях передачи существуют падающие и отраженные волны [1]. Важнейшим параметром, характеризующим качество линии, является отношение амплитуды отраженной волны к амплитуде падающей волны Г5х в любом заданном сечении линии. Это отношение принято называть
коэффициентом отражения [1, 2]. Обоснованность его выбора в качестве обобщенного параметра обусловливается тем, что остальные параметры передачи и влияния как коаксиальных, так и симметричных кабелей, определяются также решением обобщенных телеграфных уравнений и их величины коррелированы, так как отражают степень нерегулярности кабеля [1-3].Так, коэффициент затухания выражается через коэффициент стоячей волны, определяемый величиной Г_дх, и показано, что при нерегулярности периодического характера вызванное ею приращение собственного затухания кабеля Dal прямо пропорционально квадрату модуля входного коэффициента отражения [4].
Следует отметить, что в технических условиях на кабель оговаривается неравномерность волнового сопротивления, выраженная величиной коэффициента стоячей волны напряжения КСВн, более удобной при измерениях. Величина КСВн определяется величиной модуля входного коэффициента отражения Г:
КСВн = . (1)
1 -Г
Причем в технических условиях оговаривается как максимально допустимое, так и среднее КСВн в заданной полосе частот.
Достижение минимального уровня нерегулярности кабеля возможно только при оптимальном управлении его параметрами, формируемыми в ходе технологического процесса. При изготовлении низкочастотных кабелей (до 10 кГц) хороший результат дает оптимизация по типовым критериям оптимизации - максимуму быстродействия или минимуму дисперсии ошибки. Лучшие результаты дает использование регуляторов, робастных к скорости технологического процесса и межконтурных терминальных регуляторов [5].
Однако в высокочастотных кабелях проявляется эффект интерференции гармоник передаваемого по кабелю высокочастотного сигнала с соответствующими гармониками остаточной нерегулярности кабеля, обусловленной динамическими ошибками комплекса систем автоматической стабилизации, которым оснащено технологическое оборудование [4].
Проведенные исследования показали, что использование робастных регуляторов в САС первичных параметров передачи коаксиального кабеля, благодаря интерференции гармоник сигнала с остаточной нерегулярностью, приведет к получению кабеля с существенно худшей характеристикой ГВх (/), чем без автоматической стабилизации, то есть в режиме ручного управления [6].
Поэтому оптимизация всего комплекса локальных САС должна основываться на использовании математической модели, отражающей взаимосвязь обобщенного параметра, характеризующего качество нерегулярного кабеля, с управляемыми параметрами кабеля и технологическими режимными параметрами, являющимися воздействиями объекта управления. Стохастический характер управляемой системы требует вероятностного подхода к построению данной математической модели.
Для слабонерегулярной кабельной линии телеграфные уравнения, описывающие распространение волн вдоль линии, могут быть сведены к уравнению Риккати относительно входного коэффициента отражения в произвольном сечении кабеля с координатой х. Решение этого уравнения имеет вид [2]:
1
Г Вх =—12'( Х)е
1-ух
2 20 0
йх, (2)
где 20 - среднее значение волнового сопротивления линии; 2(х) - закон изменения волнового
сопротивления 2 по длине линии х; г'(х) — производная 2(х) по длине х; I - длина линии; у -
коэффициент распространения;
у = а + г( , (3)
где а - коэффициент затухания, ( - коэффициент фазы;
Р = (д^ЬС, (4)
где д - частота, Ь— индуктивность линии на единицу длины, С - емкость линии на единицу длины;
С = ^ = 1Ц-Ц- (Ф/м), (5)
18 1 Б — 1п —
й й
где £ 0 - диэлектрическая проницаемость вакуума, £ г - относительная диэлектрическая проницаемость среды между внутренним и внешним проводниками кабеля, Б - внутренний диаметр внешнего проводника, й - диаметр внутреннего проводника;
Ь = ^ 1пБ = 2 • 10-71пБ (Гн/м), (6)
2л й й
где т - магнитная проницаемость материала проводников.
Емкость и индуктивность линии определяют ее волновое сопротивление на высоких частотах [2]:
60 Б
2 Ч С =Х 7- <7)
Выразив относительную диэлектрическую проницаемость среды между двумя проводниками кабеля £г через параметры изолированного внутреннего проводника, можно получить выражение, отражающее зависимость волнового сопротивления кабеля от его параметров, формируемых в ходе технологического процесса [6]:
2р ‘
іп Б
й
, б і , Б
іп—+—іп—
й єи й
(8)
где Би - диаметр изоляции внутреннего проводника, Єи - относительная диэлектрическая проницаемость изоляции.
Наряду с выражениями (7) и (8), волновое сопротивление может вычисляться по результатам измерения емкости и коэффициента укорочения волны X в соответствии с рекомендацией Международной электротехнической комиссии:
3333
X, (9)
С
где X = С = с^ЬС = ,£, с — скорость света в вакууме, V — фазовая скорость (скорость распро-
V
странения электромагнитной волны в кабеле).
Для симметричного радиочастотного кабеля, содержащего два изолированных проводника в общем экране, волновое сопротивление определяется межосевым расстоянием двух проводников а1, зависящим от диаметров их изоляции, диаметром экрана Ое и эквивалентной диэлектрической проницаемостью среды £ [7]:
120,
2 = —=■ ІП V Є
С 2 2^
2а1 Б; - а{
й + а?
(10)
Выражения (7)—(10) позволяют выразить нерегулярность волнового сопротивления коаксиального кабеля через вариации его параметров П, формируемых в процессе изготовления:
і = 2П +
где
П П о + ^А2/(х) = 2о 2/(x), /=1 /=1 (11)
= -П АП / (х) = К • АП / (х), і 0 (12)
К - чувствительность волнового сопротивления к малым вариациям /-того параметра кабеля.
Детерминированная модель. Обозначим величину коэффициента отражения на произвольной фиксированной частоте/через Гб/ и перепишем выражение (2) в следующем виде:
ГБх/ = 2— Ї 2\х)И(х)е ] 2Ьхйх, 22
(і3)
где
Н (х) = [1(х) -1( х -1 )]е ~2ах. (14)
Величина 2^, входящая в выражение (13), определяет угловую пространственную частоту гармоники функции 2(х), взаимодействующей с гармоническим сигналом частоты / :
2Ь = 2^ = 2р ^ = 2р - = 2яg.
V V Ь
(15)
где V — фазовая скорость; Ь — длина волны взаимодействующей гармоники; g — пространственная частота взаимодействующей гармоники.
Как следует из (13), при а ® 0 и I ® ¥ величина коэффициента отражения на частоте / полностью определяется амплитудой взаимодействующей гармоники с частотой
2 Г
g = -
С учетом (15) выражение (13) можно представить в виде
ГБх/(g) = ?2'(х)И(х)е-2ЩХйх = F{2'(х)И(х)} = 2'С?) *ИС?X
17 г. 1 І7~ І7~
(16)
(17)
0 -¥ 0 0 т.е. величина входного коэффициента отражения определяется преобразованием Фурье от произведения функции 2'(х) и Н (х) или же сверткой их спектров Фурье:
V
2'(я) = ^'(х)},
Н(я) = F{H(х)} = [ Н(х)е~}2pgхdх =-----------1-----(1 - е~{2а+}2щ)1). (18)
2а + 1 2ря
2а + у2ря
Квадрат модуля функции Н (я) имеет вид:
1 о -2а/ „„„ , „-4а/
Н (я)2 = 1 2е 2С0252ря/ +е-------------------------------------------------. (19)
4р я + 4а
I |2
Функция Н(я) при я=0 достигает максимума, равного
А =
' 1 - е ~2а/ Л2
2а
V J
2
(20)
При а/ << 1 максимальное значение А равно / . При удалении точки я от точки я=0 функция
22 Н (я) быстро убывает. Площадь под кривой Н (я) :
Г I I? 1 - е
В = [ Н (я ) 2 dя = —----------------------------------------------------. (21)
Л 1 1 4а
При малых значениях а площадь
-4а/
4а
В »/. (22)
2
Эффективная ширина функции Н (я) :
= В = (1 - е-4а1 )4а2 = а- (1 - е~4а1)
А 4а(1 - е-2а‘)2 (1 - е-2а‘)2 ' ( )
При а/ >> 1 имеем
Ляу » а, (24)
а при а/ << 1 -
Ля у »1 • (25)
2
Зависимости (24) и (25) показывают, что функция Н(я) имеет характер узкого пика, и поэтому свертку (17) можно рассматривать как процедуру сглаживания широкополосного спектра г (х) узкой весовой функцией Н(я) (аналогично сглаживанию периодограммы функцией
спектрального окна при оценивании спектральной плотности случайного процесса). Соответственно выражения (23), (24) и (25) определяют ширину отрезка частот, которому принадлежат
гармонические составляющие функции г(х), обусловливающие величину входного коэффици-
яу
ента отражения на частоте / = —, где я - координата середины указанного отрезка. Эти же
выражения с учетом (16) определяют и характерную ширину пиков частотной рефлектограм-мы. Зависимость (24) совпадает с зависимостью, полученной в [8] для бесконечно длинной линии, а зависимость (25) с учетом (16) с зависимостью, полученной в [9] для линии без затухания. Зависимость (23) справедлива для линии с любым соотношением длины и коэффициента затухания.
На рис. 1 изображена зависимость эффективной полосы частот от суммарного затухания линии. Из него видно, что границей области применения формул (24) и (25) является величина затухания 1 Нп. Для использования графика при конкретных значениях длины линии / необходимо градуировать ось ординат в соответствии с числовым значением указанной точки —. Ширина эффективной полосы частот в герцах может быть получена с использованием зависимости (16). Если не исследуются фазовые характеристики линии, то ее нерегулярность полностью описывается частотной характеристикой модуля входного коэффициента отражения:
Г=|ГВх|. (26)
При малом коэффициенте затухания и большой длине линии величина модуля входного коэффициента отражения на частоте / согласно (13) может быть представлена в виде
где Ї г.(ї)- значение периодограммы функции г'(х) на частоте ї,
П(ї)=у
1
| г'(х)е“12Р *дх
Или, учитывая, что получим
П (ї) = 4р2ї2Я (ї),
Г, = ї-Яї).
/
.-[2і
(27)
(28)
(29)
(30)
Дї, м 1
Из (30) видно, что характер частотной рефлектограммы определяется характером периодограммы функции г (х).
Полученные результаты основываются на теории малых отражений, которая применима и для сверхпроводящих линий передачи.
В связи с тем, что по длине линии ее волновое сопротивление флуктуирует около некоторого среднего значения г0, представим функцию г( х) в виде
1.5
У
г(х) = г0 + г( х) (31)
и определим взаимосвязь периодограммы функции
г(х) с функцией Г(/) при произвольной длине и коэффициенте затухания линии.
Периодограмма П. (я) характеризует распреде-
1_
г
0.5
~т
0 0 5 1 0 1 а1,нп
Р и с. 1. Зависимость эффективной полосы
« „ , •, ч частот от суммарного затухания линии
ление мощности функции 7(х) по различным пространственным частотам. Если рассматривается малый отрезок Дї с центром в точке ї, то, согласно теореме Парсеваля, площадь под периодограммой на отрезке
1.0 1.5
Дї + Дї
ї-Т’ї+т
равна
g 2 1 У 2
Пї = 1 П (ї )дї = У 17ї (х)дх
Дї У 0
(32)
(33)
где ії (х) - составляющая функции і(х) с частотами, принадлежащими отрезку (32). Если
і(х) - функция гармоническая, т.е.
і ї (х) = Дг 8Іп 2п§х,
1
и если У >> —, то амплитуду Дг можно оценить: ї
Дг = Л/2П
(34)
(35)
где Пя определяется как площадь соответствующего пика периодограммы. Оценивание амплитуды периодической нерегулярности по площади пика периодограммы предпочтительнее, чем по его амплитуде, так как амплитуда меняется при изменении разрешающей способности периодограмманализа. Таким образом, возможно использование модели функции г(х) в виде суммы гармонических составляющих, каждая из которых обусловливает соответствующий пик периодограммы. Более удобной представляется модель, которая может быть получена следующим образом. Периодограмма рассматривается как совокупность “фоновой” составляющей и
2
2
пиков, существенно превышающих “фоновый” уровень. Последние интерпретируются как гармонические нерегулярности с амплитудами, вычисляемыми по выражению (35). “Фоновые” составляющие предварительно сглаживаются путем скользящего усреднения, в результате чего
получается медленно меняющаяся (плавная) функция П. (я).
Если на отрезке (35) величина П. (я) меняется незначительно, то
Пя» П-2 (я )Ля. (36)
С учетом того обстоятельства, что свойства линии на частоте / обусловливаются составляющими функции г(х) с частотами, принадлежащими отрезку шириной Ляэ с центром в точке (16), при исследовании характеристик линии должна приниматься во внимание величина
П. = П (я) <001. (37)
где а берется на частоте f.
Влияние на характеристики линии составляющей г я (х), определяемой зависимостью (33) с учетом (37), эквивалентно влиянию г(х), определяемой зависимостью (34), если их мощности равны. Поэтому, подставляя (37) в (35), получим выражение, позволяющее рассчитать амплитуду гармонической нерегулярности, эквивалентной исследуемой нерегулярности с плавным спектром на частоте я:
42 ^ П (я). (38)
Таким образом, при исследовании образца кабеля с произвольной г(х) могут быть использованы известные простые зависимости [5] для случая гармонической нерегулярности, амплитуда которой определяется либо по формуле (38) для случая плавного изменения периодограммы в окрестности рассматриваемой частоты я, либо по формуле (35) для случая узкого пика периодограммы (шириной порядка Ляэ). Так, при плавном изменении периодограммы модуль входного коэффициента отражения
Г/ =
П^). (39)
Для сверхпроводящего кабеля с малым суммарным затуханием выражение (39) сводится к виду (30). Это обстоятельство позволяет сделать вывод, что при исследовании сверхпроводящих кабелей с а/ << 1 можно пользоваться простой зависимостью (30) для нерегулярностей
любого вида, не выделяя из них периодические. Кривая Гf (я) подобна кривой я^П, (я) .
При большом суммарном затухании кабеля периодограмма функции г(х), определяемая на
всей длине У, получается с большей разрешающей способностью
т
Vу,
чем ее обеспечивает час-
тотный рефлектометр (Д/э ~Дїэ). Поэтому при исследовании такого кабеля методом периодо-
грамманализа характеристика П. (ї) должна рассчитываться по функции г(х), оцененной на участке кабеля от места подключения рефлектометра до точки с координатой 1и, где
'■=дї: <40)
Полученная периодограмма также будет подобна частотной рефлектограмме. При необходимости более точного оценивания частотных характеристик кабеля должна использоваться исходная зависимость (2).
Стохастическая модель. Детерминированная модель позволяет использовать результаты периодограмманализа волнового сопротивления и конструктивных параметров кабеля для исследования нерегулярности конкретного образца кабеля. Для определения более точных закономерностей, устанавливающих взаимосвязь вероятностных характеристик нерегулярности, необходимо рассматривать функции г(х) и Г(/) как случайные функции, импульсную и час-
тотную рефлектограммы конкретного образца кабеля - как реализации этих случайных функций. Детерминированные периодические составляющие функции г(х) могут быть выявлены по незатухающим колебаниям корреляционной функции волнового сопротивления или по пику в оценке спектральной плотности, полученной усреднением ансамбля периодограмм волнового сопротивления, и вычленены. В кабелях высокой регулярности периодические нерегулярности удается устранить технологическими мерами, и достижимый уровень качества кабеля определяется уровнем его случайных нерегулярностей.
Закон распределения случайной величины Г является исчерпывающей вероятностной характеристикой при рассмотрении свойств нерегулярной линии на фиксированной частоте. При исследовании линии в полосе частот необходимо знать не только закон распределения Г на каждой частоте, но и многомерные законы распределения и моментные характеристики, так как Г( f) является нестационарной случайной функцией. Однако, учитывая “решетчатый” характер функции Г( f), можно перейти к рассмотрению множества ее независимых отсчетов, разделенных интервалами А/~э:
Частотные зависимости закона распределения Г и интервала Д/, зависящего от длины линии и ее затухания, достаточно полно характеризует вероятностные свойства нерегулярной линии.
Анализ функции Г(/) можно упростить, рассматривая ее характеристики на фиксированной частоте, где ее можно считать случайной величиной, либо в узкой полосе частот, где ее можно считать стационарной случайной функцией.
При изучении свойств кабеля на фиксированной частоте / величина модуля входного коэффициента отражения Г/ может рассматриваться как случайная величина, вероятностные характеристики могут быть оценены по выборке, полученной путем измерения Г/ на достаточно
большом количестве образцов кабелей, изготовленных по одной технологии. С другой стороны, отсчеты функции Г(/), разделенные интервалом Д/, можно считать двумя некоррелированными случайными величинами, имеющими одинаковые законы распределения ввиду малости интервала (41). Поэтому вероятностные характеристики Г/ могут быть оценены и по выборке, полученной путем измерения Г/ на одном отрезке кабеля в узкой полосе частот с центром в точке / [9]. Таким образом, измерения, проведенные на одном образце кабеля, позволяют сделать выводы о вероятностных характеристиках генеральной совокупности образцов кабелей, изготовленных на одном технологическом оборудовании.
Вероятностные характеристики Г(/) полезны, так как, с одной стороны, они описывают достигнутый уровень технологии, с другой стороны, позволяют оценить вероятностные характеристики искажений сигнала, передаваемого по нерегулярной линии [8].
Определим вероятностные характеристики входного коэффициента отражения для бесконечно длинной линии (У = ¥), каковой может считаться и линия конечной длины с суммарным затуханием а1 > 1 Нп (см. рис. 1). Более того, уже при суммарном затухании порядка 5дБ величина измеряемого входного коэффициента отражения кабеля конечной длины отличается не более чем на 10% от соответствующей длины бесконечно длинной линии [4]. Будем рассматривать наиболее важный случай стационарной эргодической функции г(х) с нормальным законом распределения [4, 8].
Возьмем произвольную точку с координатой х и определим, чему равен входной коэффициент отражения в данной точке. Под входным коэффициентом отражения в точке х мы понимаем сумму всех отражений от нерегулярностей правее данной точки. Он определяется выражением
где V - вспомогательная переменная.
Ввиду того, что у - величина комплексная, для фиксированной частоты / функция ГВх (х) получается в результате линейного преобразования г(х), представляет собой комплексную стационарную случайную функцию действительного аргумента х и распределяется по круго-
(42)
вому нормальному закону с нулевым средним. Поэтому стационарная случайная функция ГВх (х) полностью характеризуется своей корреляционной функцией:
Г (т) = M [гВх (х +т )G*e (х)] = ] J Je-2(a+)ve“2(a - j )v'Kz, (т + v - V)dvdV
4z0 0 0
Двухсторонняя спектральная плотность функции ГВх (х) имеет вид
¥ 1
SГД (g) = J Kr (т)е-J 2pgT dt = -T | G (- g )|2 SД z,(g), (43)
-¥ 4 zo
где
G(-g) = Je_2(“+jb)xeJ2pgxdx, (44)
0
21
G(-g) 2 = —---------------(45)
4a 2 + (2b - 2pg )2
Kz.(t), SД z,( g) - корреляционная функция и двухсторонняя спектральная плотность случайной
функции z'(х). Функция G(g) определяется как преобразование Фурье от функции е~2gc.
Необходимо отметить, что для получения вероятностных характеристик ГВх на определенной частоте f необходимо в полученные выражения подставить соответствующие этой частоте значения а и b :
b = — . (46)
v
I |2
Сравнивая выражения (44) и (18), видим, что функция G(-g) может рассматриваться как
I 12
функция H(g) , полученная при l = ¥, смещенная по оси g на 2f / v и имеющая характер узкого пика с эффективной шириной
Dgэ ------J |G(-g)f dg = а (47)
|G(-gA max -”
и площадью, равной (4а)-1.
Величина а мала по сравнению с возможной шириной спектральной плотности волнового сопротивления, поэтому можно считать, что в пределах полосы частот Dg^ с центром в точке 2 f
g0 = — выполняется условие v
Sд,.(g) » SДг,(g0). (48)
Дисперсию функции ГВх (х) можно определить, интегрируя Sp (g):
s2 г =
J SrG" (g)dg = -L J |G(-g)f Sд,,(g)dg »J |G(-g)|:
-¥ 4 Z0 -¥ 4 Z0 -¥
= ^ 5 Д,'(я0) = (я0). (49)
16 г0а 4 г0а
Формула (49) отражает зависимость дисперсии комплексного входного коэффициента отраже-
г я 0 у
ния линии от частоты передаваемого сигнала / = —^~ и позволяет решить задачу определения
диапазона частот, в котором функция ГВх (^ может считаться стационарной. Таким диапазоном можно считать интервал частот, в пределах которого функция 5Д, (я) меняется незначительно. Критерием незначимости изменения функции 5Д, (я) может служить статистическая
2
незначимость соответствующего приращения дисперсии о ГВх .
Дисперсия а гВ1 случайной величины ГВ1, распределенной по круговому нормальному закону, равна удвоенному значению квадрата радиального среднеквадратического отклонения исходного нормального двухмерного распределения. При этом модуль ГВ1 распределяется по закону Рэлея [9, 10]:
F(Г) = Р[Гp < Г] = Рг = l - exp
Г2
а 2 Гв
Дисперсия модуля входного коэффициента отражения задается выражением
4 - p 2
------а 1Г
(4 - p)p g0 16 zOa
SДz (g0).
Отсюда
F[Г] = l - exp Плотность распределения имеет вид
f [Г] = 4 - p Г exp
4 - p Г 2 4 а 2 г
2
а 2 г
4-p Г2
4
а2Г
Математическое ожидание модуля входного коэффициента отражения:
тг =
p
4-p
(50)
(5l)
(52)
(53)
(54)
Величина модуля входного коэффициента отражения, которая не будет превышена с вероятностью Рг:
ГР =
4
-ln-
l
Р(Г)
и 1 г. (55)
'4 - р 1 - рг
Гипотеза о рэлеевском характере закона распределения модуля входного коэффициента отражения была экспериментально проверена и подтверждена при исследовании радиочастотных кабелей [4, 10]. Для проверки этой гипотезы для сверхпроводящих радиочастотных кабелей
исследовался опытный образец № 10 кабеля длиной 166 м на частотах до 1 ГГц. В диапазоне
—2
частот 45 0-5 50 МГц оценка дисперсии оказалась равной 0,39 10 [11].
На рис. 2 приведена гистограмма частотной рефлектограммы кабеля в диапазоне частот 450-550 МГц и теоретическая кривая, построенная по формуле (53) при а2г = 0,39 10—2. Экспериментальные данные не противоречат гипотезе о рэлеевском характере распределения модуля входного коэффициента отражения сверхпроводящего радиочастотного кабеля.
Таким образом, согласно (49), вероятностные характеристики функции Г(/) будут полностью определены, если станет известной спектральная плотность волнового сопротивления.
Определим зависимость, аналогичную (49), для линии произвольной длины с произвольной величиной коэффициента затухания. Для этого снова воспользуемся тем обстоятельством, что на характеристики линии на частоте / влияет часть спектра функции z(х) на отрезке частот с центральной частотой (16) и шириной (23). Тогда из всей дисперсии волнового сопротивления на характеристики линии
„ 2 на частоте / влияет лишь ее часть а g,
определяемая как площадь под кривой
односторонней спектральной плотности Sz (g) над указанным отрезком частот. Ввиду малости
длины этого отрезка частот можно записать:
б
4
2
0
0,05
0,l0
0,l5
0,20
Г
Р и с. 2. Гистограмма модуля входного коэффициента отражения
Г
а2, »Д£э^(я) = ^11 ^(8).
2 *
Здесь а 8 есть дисперсия случайной, но узкополосной (квазигармонической) функции г я (х) с
центральной частотой я . Для математического описания этой функции модель (34) непригодна. Узкополосная случайная функция г я (х) может быть представлена в виде
г я (х) = Дг( х) со8[2л£х + 0( х)], (56)
где Дг( х) и 0(х) медленно меняющиеся (по сравнению с я) случайные функции.
Амплитуда Дг(х) распределяется по закону Рэлея [12], поэтому дисперсия амплитуды
2
а д связана с дисперсией функции (х) следующим образом:
а2 2 4 -Па2 а(4 - р)(1 - е~4а!) е (^ (57)
а Дг = 2_ая = 2(1-е_2а!)2 ^ (Я). (57)
Учитывая, что влияние квазигармонической нерегулярности на характеристики линии эквивалентно влиянию гармонической нерегулярности той же мощности, можно воспользоваться известной зависимостью, связывающей модуль входного коэффициента отражения линии с амплитудой гармонической нерегулярности [4]:
Г = --pg-(1 - е~2Ы ]дг. (58)
4 г оа
Если амплитуда Дг случайна, то величина Г является случайной величиной с дисперсией
2
22
_ 2 р § /л -2а! \ 2 —2 /^о\
аг =------Г“7(1 - е ) а Дг . (59)
16г 0а
Подставляя (57) в (59), получим
2 р 2 (4 - р)1 - е“4а! 2
а г =—--------------------я 2 (я). (60)
32г2 а
2
Дисперсию а г можно выразить и через одностороннюю спектральную плотность функции
аг- = Л 2 я2Sz (я). (62)
г (х):
2 4 - р 1 - е“4а!
а г 2 =------2---------$А я). (61)
128г2 а
Из формул (49) и (51) следует, что для бесконечно длинной линии дисперсия модуля входного коэффициента отражения связана с односторонней спектральной плотностью волнового сопротивления следующим образом:
р (4 - р)
32г^а
Это же выражение получается из (60) при подстановке ! = —, что подтверждает правильность сделанных выкладок
Из (60) и (62) видно, что для линии конечной длины
2 -4а! 2
аг = (1 - е )аг- .
Интересно также сравнить полученные результаты с известными результатами для линии с цепочкой сосредоточенных нерегулярностей, равновероятно распределенных вдоль линии [4, 9, 10], согласно которым наиболее вероятное значение модуля входного коэффициента отражения гт :
2 2 2 тр2 1 - е~4а!
гт2 =--аг2 = —^ =--------------------, (63)
т 4 - р 2 4а!
где т - число нерегулярностей в линии; р2 - дисперсия местных коэффициентов отражения. Учитывая, что — = g , получим
_ 2 4 -р г 2 4 -р 1 - е"4а/ р,2
Г = —г— =—2-------------- ---^^. (64)
2 32 а 2
Приравнивая (64) и (60), получим
— 2 / л ,—\ 1 —4а/ л 1 —4а/ „ 2
р (4-р)1 -е 2о / ч 4-р 1 -е р,
-----------g ^ (g) = —-------------. (65)
32г^ а 32 а 2
Отсюда следует, что зависимость (63) является частным случаем зависимости (60) и справедлива при
^ (g)=Щ- 1, (66)
2р g
т.е. в случае, когда спектральная плотность волнового сопротивления изменяется обратно пропорционально частоте g.
Таким образом, полученные зависимости (50) - (55), (60), (61) позволяют рассчитать все вероятностные характеристики входного коэффициента отражения линии по известной спектральной плотности волнового сопротивления.
Для экспериментальной проверки этих зависимостей была обработана импульсная рефлек-
тограмма опытного образца № 10 и получена оценка спектральной плотности функции г (х), изображенная на рис. 3. Здесь же изображена кривая изменения математического ожидания модуля входного коэффициента отражения, оцененная по частотной рефлектограмме. Теоретическая зависимость между этими функциями имеет вид
тг =
(Г
Для образца № 10 а/ << 1, поэтому из (67) находим
128^2 * а
тг » ' Р/ = 0,08, ма5-ОмЛ (68)
о
Р и с. 3. Статистические характеристики кабеля:
1 - оценка спектральной плотности производной волнового сопротивления; 2 - оценка математического ожидания модуля входного коэффициента отражения
Оцененное по кривым рис. 3 отношение (68) оказалось равным 0,096. Учитывая малую точность частотного и импульсного рефлектометра, а также погрешности, вносимые нелинейностью шкал рефлектометров, можно считать совпадение теоретических и экспериментальных результатов удовлетворительным.
На рис. 4 показаны кривые спектральной плотности производной волнового сопротивления опытного образца № 10 и спектральной плотности производной диаметра изоляции того же ка-
беля. Кривые подобны, что подтверждает предположение о превалирующем влиянии диаметра изоляции на волновое сопротивление данного кабеля. Поэтому спектральная плотность волнового сопротивления кабеля может быть оценена по параметрам его элементов.
Р и с. 4. Статистические характеристики кабеля:
1 - оценка спектральной плотности производной волнового сопротивления;
2 - оценка спектральной плотности производной диаметра изоляции
Таким образом, полученные зависимости (50) - (55), (60), (61) позволяют рассчитать все вероятностные характеристики входного коэффициента отражения на фиксированной частоте / по известному значению спектральной плотности волнового сопротивления на соответствующей частоте g. Кроме того, рассчитав спектральную плотность волнового сопротивления и затем частотную зависимость Гр при вероятности рг = 0,90 — 0,95, можно построить огибающую частотной рефлектограммы кабеля еще на этапе изготовления отдельных элементов кабеля.
Полученные зависимости позволяют, в частности, записать выражение, показывающее, что математическое ожидание модуля входного коэффициента отражения на произвольной частоте / определяется величиной спектральной плотности входного сопротивления на пространствен-
2 /
ной частоте g =-----. С учетом зависимости (1) на плоскости Sz (g) можно построить семейст-
во границ равных отражении , т.е. границ, выход за которые спектральной плотности волнового сопротивления приводит к нарушению допуска на заданные значения КСВн. Поэтому области допустимого расположения кривых спектральных плотностей волнового сопротивления кабелей расположены внутри указанных границ (ниже их).
На рис. 5 показаны границы равных отражений, рассчитанные для магистрального кабеля для систем кабельного телевидения РК 75 - 17 - 13 - С (г0 = 75 Ом; £ =1,5; V = 2108 м/с; С= 67 пФ/м; а = 0,046 дБ/км; КСВнтах<1,35) и соответствующих неравномерности волнового сопротивления, выраженной величинами КСВн, равными 1,1; 1,2; 1,3; 1,4. Так как характеристики кабеля нормируются в полосе частот 30+270 МГц, границы равных отражений с учетом (16) построены в диапазоне пространственных частот от gн = 0,3м —1 до gв = 2,7м —V
Выбрав типовую модель случайной функции z(х) и, выразив ее спектральную плотность через дисперсию
2 „ аz и интервал корреляции Тк, на
2
плоскости (а,т%) также можно выделить области равных отражений.
С учетом связи волнового сопротивления с параметрами кабеля (12) спектральная плотность волнового
Р и с. 5. Границы равных отражений
V
сопротивления может быть выражена через спектральные плотности параметров кабеля, формируемых в ходе управляемых технологических процессов:
^ (я) = £К2Б1, (я) + £К1К]Ь1] (я), (69)
,=1 из=и*}
где Ь у (я) - синфазные компоненты взаимных спектральных плотностей ,-го и у-го параметров кабеля.
Спектральные плотности параметров кабеля, рассматриваемые как функции его длины х, могут быть представлены и как спектральные плотности сигналов датчиков данных параметров, установленных на технологическом оборудовании, регистрируемых как функции времени:
Бп(я) = 2путБп (а), (70)
4рут /
где а =-, Ут - скорость движения кабельного изделия в ходе технологического процесса.
V
На плоскостях спектральной плотности каждого параметра также можно показать границы равных отражений, которые, в соответствии с (69) характеризуют вклад параметров в волновое сопротивление.
Оценив по спектральной плотности волнового сопротивления вероятность непревышения модулем входного коэффициента отражения заданного значения Г по формуле (50), можно определить вероятность того, что из п кабелей заданное значение Г не будет превышено не менее чем в т кабелях:
Рт =Ё Рип , (71)
,=т
где Р,,п = Сп'Рг1 (1 - Рг ) п~',
с , = п!
п ,!(п -,)!
И, наконец, рассматривая кабель в полосе частот Д/, где функцию ГВх (/) можно считать стационарной, и считая значения входного коэффициента отражения, разделенные интервалом
частот (41), независимыми случайными величинами, можно определить вероятность непревышения величиной Г заданного значения в полосе частот Д/ [9, 10]:
Рг= Ргм , (72)
где
N = ^Д- . (73)
Д£ ^
Использование вероятности Рг в формулах (71) позволяет рассчитать вероятность того, что из п кабелей заданное значение Г не будет превышено в заданной полосе частот Д/ не менее чем в т кабелях.
Если же рассматривается более широкая полоса частот, где функцию ГВх (/) нельзя считать стационарной, вероятность Рг должна рассчитываться по формуле
N
Рг=ПРГ , (74)
г=0
где Рг - значение Рг, рассчитанное на частоте / :
/ = /0 +ХД4 (75)
к=0
/0 - нижняя частота рассматриваемой полосы частот; Д/к - значение Д/, рассчитанное по формулам (41) и (23) при подстановке значения коэффициента затухания а, соответствующего /к.
Если рассматриваемая полоса частот достаточно широка, частотно-зависимый коэффициент затухания заметно изменяется, и значения Д/к получаются неодинаковыми: они увеличиваются с ростом частоты.
Для больших значений п формула (71) может быть упрощена. При рассмотрении множества образцов кабелей может быть определено математическое ожидание числа кабелей, в которых не будет превышено заданное значение Г :
М [т] = пРГ. (76)
При рассмотрении одного образца кабеля может быть определено математическое ожидание числа пиков частотной рефлектограммы, имеющих амплитуду больше заданного значения Г :
М[к] = ИРГ . (77)
Полученная модель кабельной линии распространена и на случай математического описания усилительного участка, когда линия собирается из ряда строительных длин. Функция г(х)
линии в данном случае представляется в виде суммы функций г , (х), умноженных на функции временного окна [13]. Но так как интервал корреляции волнового сопротивления много меньше строительной длины кабеля, то модель усилительного участка сходится к исходной модели кабельной линии.
Передача импульсных и цифровых сигналов по линии связи делает необходимым получение аналогичных моделей, когда обобщенными параметрами являются параметры импульсной характеристики линии и ее затухание. На примере сверхпроводящего кабеля для цифровых систем передачи показано, что данные характеристики нерегулярного кабеля также выражаются через характеристики входного коэффициента отражения [14].
Известно, что коэффициент затухания регулярного радиочастотного кабеля определяется двумя слагаемыми. Первое слагаемое зависит от размеров и свойств проводников, а второе - от свойств проводников не зависит и определяется теми свойствами изоляционного материала кабеля, которыми практически невозможно управлять в процессе экструзии [4, 9, 10]. Показано так же, что за счет переотражений волн на случайных нерегулярностях возникает приращение затухания Да!, которое может быть выражено через Гтах [14] и пропорционально дисперсии
местных коэффициентов отражения р ,2 [4].
Учитывая зависимость дисперсии модуля входного коэффициента отражения от спектральной плотности волнового сопротивления (60), данное приращение затухания также выражается через величину спектральной плотности волнового сопротивления на соответствующей пространственной частоте £:
Да! = К —--------/2 Б
а
, V ,
V У
(78)
Здесь К - частотно-зависимый коэффициент, определяемый конкретными особенностями кабеля [4, 14], вид которого желательно эмпирически уточнить для получения адекватной модели.
У симметричных кабелей главное ограничение на полосу частот кабеля накладывается характеристиками взаимного влияния, а не параметрами передачи [1]. Для симметричных кабелей с четверочной скруткой в качестве обобщенного параметра выбирается переходное затухание на ближнем конце А0, которое определяется частотной характеристикой электромагнитной связи между цепями на ближнем конце N12 (/). Частотная характеристика N12 (/) находится решением обобщенных телеграфных уравнений, которое имеет вид (2), где вместо функции г'(х) фигурирует меняющаяся по длине кабеля электромагнитная связь ^2(х) [1]. Учитывая, что вариации N12 (х), благодаря малости вариаций параметров изолированных жил симметричного кабеля, могут быть выражены через их приращения ДП 1 (х) аналогично (12), все полученные при построении стохастической модели коаксиального кабеля результаты могут быть применены и для математического описания формирования симметричного кабеля [6].
Следует отметить, что спектральную плотность /-того параметра кабеля БП (а) можно рассматривать как спектральную плотность динамической ошибки системы автоматической стабилизации /-того параметра. Она определяется спектром возмущающего воздействия , (а) и амплитудно-частотной характеристикой системы А, (а). Поэтому формирование обобщенных параметров нерегулярного кабеля связи в ходе технологического процесса ТП, под управлением и,(/) комплекса систем автоматической стабилизации КСАС, описано построенной математической моделью ММ и может быть представлено схемой, показанной на рис. 6.
Ui (t)
Р и с. 6. Схема формирования обобщенных параметров качества нерегулярных кабелей связи Вектор обобщенных параметров качества нерегулярных кабелей связи может быть представлю в виде Q = (КСВнср;КСВНшах; Да/шах;Ж12шах).
Данная модель позволяет организовать оптимальное управление технологическими процессами производства кабелей связи. Целью управления должна быть выбрана минимизация уровня нерегулярности кабеля путем оптимизации всего КСАС по критерию минимума выбранного главного обобщенного параметра при наложении ограничений на остальные.
Имитационное моделирование управляемого технологического процесса показало, что использование классического критерия максимального быстродействия приводит к существенно худшим частотным характеристикам нерегулярного кабеля, чем критерия, отражающего взаимосвязь обобщенных параметров с параметрами регуляторов [6].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. ГродневИ.И., Шварцман В.О. Теория направляющих систем связи. М.: Связь, 1978. 296 с.
2. Гальперович Д.Я., ПавловА.А., ХренковН.Н. Радиочастотные кабели. М.: Энергоатомиздат, 1990. 256 с.
3. Шварцман В.О. Взаимные влияния в кабелях связи. М.: Связь, 1966. 431 с.
4. ДорезюкН.И., ПоповМ.Ф. Радиочастотные кабели высокой регулярности. М.: Связь, 1979. - 104 с.
5. Чостковский Б.К. Алгоритмизация терминального управления совмещенными технологическими процессами изготовления радиочастотных кабелей // Вестн. Самарск. гос. техн. ун - та,Сер.: Технические науки. Самара, 2002. Вып.14. .С. 33 - 37.
6. Чостковский Б.К. Алгоритмизация и частотная оптимизация управления процессами производств кабелей связи // Вестн. Самарск. гос. техн. ун-та. Серия: Техн. науки, 1998. Вып. 5. С. 28 - 35
7. ЕфимовИ.Е., Останькович Г.А. Радиочастотные линии передачи. М.: Связь, 1977. - 408 с
8. KarbowiakA.E. Investigation ob signal distortion in cables caused by imperfections in cable manufacture // Proc. Jnst. Elec. Eng. 1974. Vol. 121. № 6. Р. 67-73.
9. Mu//en J., Pritchard A. The statistical prediction ob voltage standing wave ratio // IRE Transaction on microwave theory and techniques. 1957. MTT 5. № 2. Р. 127-130.
10. ПерецР.И. Статистические характеристики тракта СВЧ // Антенны. 1974. № 17. С. 104-118.
11. Гальперович Д.Я., Гречков В.И., Коржакова Т.В., Чостковский Б.К. Сверхпроводящие коаксиальные пары для кабелей связи. // Электросвязь. 1990. № 1. С. 38 - 41.
12. СвешниковА.А. Прикладные методы теории случайных функций. М.: Наука. 1968. 464 с.
13. Чостковский Б.К., Тян В.К. Идентификация частотных характеристик кусочно-синтезированного случайного сигнала. // Каталог алгоритмических модулей АСУТП I - II кв. 1984, инв. № 6087, СМОФАП.
14. Попов М.Ф., Чостковский Б.К., Юдашкин А.А. Характеристики нерегулярного сверхпроводящего кабеля для цифровых систем передачи. // Кабельная техника. 1995. Вып. 7 (245). С.11 - 14.
Поступила 23.12.2005 г.