где Y — вероятность интерференциального подавления сигнала [6]; X o — интенсивность помехи; т -длительность импульса сигнала. Принимая за пуассоновское распределение вероятностей появления
импульсов помехи в интервале совпадения т с = 2т , равном разрешающей способности дешифратора, можно записать
Р01 = [1 - exp(-Xотc)]n ;р1г = [! - exp(-Xoтc)]n-1.
На рисунке представлены результаты расчета вероятности стирания информации при n=2 и N= 128, K=2 (сплошная), N=512, K=4 (штриховая) N=512, K=8 (штрихпункгирная).
0 1500 3000 4500 6000
Вероятности расчета стирания информации
Как следует из приведенных зависимостей, вероятность стирания информации существенно зависит от плотности помехи, значности используемых кодов и от объема передаваемой информации. Приведенные расчеты показывают, что для практической
реализации передачи принятой информации в пункт совместной обработки в синхронных сетях рассредоточенных систем первичной локации достаточно использования двухзначных кодов.
Необходимо отметить, что благодаря передаче всего блока принятых сигналов в одном цикле связи на пункте совместной обработки возможно осуществление межпериодной обработки принимаемых сигналов, что позволяет существенно повысить достоверность передачи принятых сигналов. Литература: 1.Ямпольский Э.М., ТолстовГ.Б. Использование многочастотных сигналов для одновременной связи и локации// Методы и устройства обработки сигналов в радиотехнических системах. 1988. С.48-52. 2. Обод И.И. Синхронные сети получения и обработки информации рассредоточенных систем первичной и вторичной локации. Рукопись деп. в УкрИНТЭИ. N192-Y197 от 25.02.97.
3. А.с. 1510555 СССР. Устройство вычисления ВКФ. Обод И.И. // Открытия. Изобретения. 1989. №46. С.184.
4. А.с. 1552199 СССР. Устройство вычисления ВКФ. Обод И.И. // Открытия. Изобретения. 1990. N 11. С. 173.
5. А.с. 1616379 СССР. Устройство вычисления ВКФ. Обод И.И. // Открытия. Изобретения. 1990, N 42. С. 158.
6. Ярлыков М. С., Черняков М.В. Оптимизация асинхронных адресных систем радиосвязи. М.: Связь, 1979. 215с.
Поступила в редколлегию 23.03.98
Обод Иван Иванович, канд. техн. наук, старший научный сотрудник, докторант кафедры № 34 ХВУ. Научные интересы: получение, обработка и передача информации в синхронных сетях разнесенных систем первичной и вторичной локации. Адрес: 310166, Украина, Харьков, ул. Коломенская, 27, кв.1.
УДК 537.87; 621.371
ПРЕОБРАЗОВАНИЕ ПЛОСКОЙ ВОЛНЫ ПОВТОРЯЮЩИМИСЯ ИМПУЛЬСНЫМИ ИЗМЕНЕНИЯМИ ПАРАМЕТРОВ БЕЗГРАНИЧНОЙ СРЕДЫ
НЕРУХА.Г., РЫБИН О.Н, ЩЕРБАТКО И.В.
Рассмотрено преобразование плоской монохроматической электромагнитной волны серией повторяющихся изменений параметров среды во временной области. Изменение состоит в модуляции диэлектрической проницаемости и проводимости по закону прямоугольного импульса. Получены выражения, описывающие эволюцию электрического поля на всей положительной временной полуоси.
Распространение электромагнитных волн в средах с изменяющимися во времени параметрами описывается уравнениями Максвелла с переменными коэффициентами. Исследованию такого распространения посвящен ряд работ [3-9]. Как теоретический, так и практический интерес представляет повторяющееся во времени изменение параметров [9]. Практическая реализация такого изменения характеризуется его ограниченностью во времени. Исследуем преобразование плоской гармонической волны
повторяющимися импульсными изменением диэлектрической проницаемости є(t) и проводимости o(t) безграничной среды
N
e(t) = є0 + (є і - є0) £ {0(t - (k - 1)T) -k=1
N
-0(t - ті - (k - 1)T)}, CT(t) = cti £ {(t - (k - 1)T) -
k=1
-0(t - ті - (k - 1)T)},
которые представляют собой скачкообразные одновременные изменения диэлектрической проницаемости и проводимости до значений є і и сті. Новые значения держатся в течение времени т і и затем скачком возвращаются к исходным значениям є = є 0, ст = 0 . Длительность периода равна Т, а всего таких периодов — N.
Из эволюционного подхода [7] для поперечной составляющей вектора электрического поля следует интегральное уравнение Вольтерра второго рода
ТО ТО
E (t ,x) = £0(t,x) + J dt'Jdx'K (t, t' ,x,x')E (t',x'). (1)
0 0
Здесь E 0(t, x) — первичное поле, существующее
в среде до начала серий; K (t,t',x,x') — ядро интегрального уравнения
28
РИ, 1998, № 1
K(t,t',x,x') = -(і/a2(t'))S(vo(t - t') - |x - x'|) X xjc(t') + |(1 - a 2(t'))д|,
где a(t) = yjz0/s(t); Vq = c/; a(t) = 2na(t)/є:; c
— скорость света в вакууме; 5(t) — дельта-функция Дирака. Уравнение (1) удобно решать методом резольвенты [2], согласно которому
ТО ТО
E(t, x) = Eо (t, x) + J dt' J dx' R(t, t', x, x' )Eq (t', x'). (2)
о о
Здесь R(t,t',x,x') — резольвентазадачи, явное выражение для которой легко получить для постоянных значений параметров путем решения интегрального уравнения
R (t,t ',x,x') = K (t,t',x,x') +
ТО ТО
+ J dt'' J dx'K (t, t", x,x")R (t", t', x",x'). (3)
0 0
Для решения данной задачи разобьем временную ось на интервалы, в которых параметры среды остаются постоянными. Тогда из [3] следует записать, что составляющая вектора электрического поля на n-м временном интервале разбиения рекуррентными соотношениями
Пусть e0(t, x) = el(rot-kx), где k = ю/vo . Тогда,
подставляя формулы (5) и (6) в выражения (4), c помощью метода математической индукции получаем выражения для электрического поля на произвольном периоде изменения параметров среды:
E(t, x) = (C+ + D+)e-CTitel(Qt - kx) +
+|c- + D-|e-CTlte-l(^t+kx), Q2 =tJa2ю2 - сі2, (7)
(k - 1)T < t <т1 + (k - 1)T, k = 2,..,N.
Далее имеем
E(t, x) = Ake1(rat - kx) + Bke-1(rat+kx),
т1 + (k - 1)T < t < kT, k = 2,..,N. (8)
Здесь
C± = a2Ak-1 (a2 - 1)ю - 2ст1(ст^+ iQ) e(k-1)(a1 +i(Q+ro))T k k 1 2iQ(i(ro + Q)±o1) ’
D± = a2Bk-1 (a2 -йю-2с1(с1 + iQ) e(k-1)(CT1 +1(Q±ro))T k k 1 2iQ(1(ro ±Q) ±o1) '
Ak =--
e-(с 1 + 1ю)Т1
~i\Q
с1 + (1- a 2)i^-
t To
E n(t,x) = J dt' J dx' R n(t, t',x,x')Fn(t',x') +
tn-1 -TO
+Fn(t, x),
N-1 11 to
Fn(t, x) = X J dt' J dx' K 1(t, t', x, x')E1(t', x') +
1=1 t i-1 -to
(4)
+E 0(t, x).
В этих соотношениях как структура резольвенты, так и структура ядра инвариантны относительно периода изменения параметров. Из уравнения (3) для скачка диэлектрической проницаемости и проводимости в момент времени t = 0 следует
K i(t,t',x,x') = K(t,t',x,x'), 1 = 1,...,N,
(1 - 1)T < t < т1 + (1 - 1)T;
K 1(t,t',x,x') = 0, т1 + (1 - 1)T < t < 1T, 1 = 1,...,N.
(5)
1 I 2 З2 _ д I
R 1(t, t', x, x') = --!—j(1 - a2)^- + 2ct^- x
a+ 1TO
X
J
dp e
2v1
(P-C1)(t -1')-1
St2
St I
V1
4
P 2 -C12
. 2ni
a- 1TO
Vp 2 -C12
1 = 1,...,N, a>CT1, (1 - 1)T < t <T1 + (1 - 1)T; R 1(t,t',x,x') = 0,т1 + (1 - 1)T < t < 1T, 1 = 1,...,N,
(6)
где a
= -\/є0/ є1 , с1 = 2лст^ є1, v1 = сД/єЬ
Re -Jp2 - CT12 > 0.
XJAk_1( a1 - i(Q + ra) eiQT1 aj + i(Q-m) є-Юх/1 +
| V01 - 1(Q-ro) a1 + 1(Q + ro) j
+2iBk-1e-21ro(k-1)T sin(Qx1)},
Bk =-
e-(a1 - r 24 1®
-----------ja1 - (1 - a2) —
21Q I 2
C - 1(Q - ®) iQr1 C1 + 1(Q + ®) -iqt 1
xjBkJ ^----:----- e
k ^ a1 - 1(Q + ®)
a 1 + 1(Q - ®)
+21A
1Ak-1e2i®(k 1)T sin(Qr1)},
На первом периоде возбуждения среды поле определяется отличными от (7) выражениями
E(t, x) = C+e-a1te1(Qt - kx) + +C-e-C1te-1(Qt +kx), 0 < t < T1.
(9)
E(t, x) = A1e1(rat-kx) + B1e-1(rat +kx), т1 < t < T. (10)
rj ~± c 2n± (a2 - 1)® - 2o1(CT1 - lQ)
Здесь Cr = E0a2,C--—- 4—-. Также 1 1 2iQ(i(® + Q) ± a1)
e -(a1 + 1ro)T1 J- ,, xlra
A1 =--------2lQ------iC1 +(1 - a 2)T
X Ia1 ~ 1(Q + ®) eiQt1 a1 + i(Q~®) e-iQl1
|a1 - 1(Q-ro) a1 + 1(Q + ro)
Кроме того,
e-(a1 -1ю)т1 J- 2\ 1® I ■ /0 \
B1 =-------q------jc1 - (1 - a2) — \ sin(QT1).
РИ, 1998, № 1
29
Как видно из (7) — (10), поле на любом временном интервале постоянства параметров среды представляет собой суперпозицию двух плоских волн: прямой (первые слагаемые в формулах (7) — (10), которая распространяется в направлении первичной, и обратной (вторые слагаемые в формулах (7) — (10), которая распространяется в противоположном направлении. В случае возбужденного состояния пространства
(o(t) = ст 1, e(t) = є і) частоты этих волн отличаются, а волновые числа совпадают c соответствующими параметрами первичной волны, амплитуды убывают во времени по экспоненциальному закону вследствие наличия проводимости в пространстве. В период
невозбужденного состояния (ст = 0, є = є о) фазовые характеристики прямой и обратной волн совпадают с соответствующими параметрами первичной волны. Если ввести приведенные амплитуды
ak = AkeikT, bk = Bke-ikT, уравнения (7) — (10) можно переписать более компактно:
ak+1 bk+1
" f - h
v-h
*Л
f* J
ak
bk
ai = f*e iT,bi = -he
iT
(11)
Здесь введены параметры, не зависящие от номера периода:
f = (e -(s+і)Ті + iT/2q){2qcos(qTi) + i(1 + a2)sin(qxi)}, h = (e-(s-i)Ti -iT/2q){2s - i(1 - a2)} sin(qxi),
ю
Из рекуррентного соотношения (11) следует [1]
ak+11= Rk (a 1І,
Фк+1J vbiJ
(12)
где R k
uk-1(Re(f))
*Л
f* J
Uk - 2(Re(f))I.
Здесь uk (a) — многочлен Чебышева порядка k ;
I — единичная матрица. Из (12) нетрудно получить амплитуды прямой и обратной волн в случае невозбужденного состояния среды:
Ak+1 = e-ikT{(f2 + |h|2e2iTJUk-i(Re(f)) --Uk - 2(Re(f))f}, k = 1,.., N.
Bk+1 = -heikT {(f + f*e2iT )uk-1(Re(f)) -
-e2iT Uk - 2(Re(f))}, k = 1,.., N. (14)
Дополнительные аналитические исследования формул (13) и (14) показали, что амплитуда прямой
волны Ak по модулю всегда больше амплитуды
обратной волны Bk.
Численный анализ формулы (13) показал резонансный характер амплитуды прямой волны в зависимости от изменения длительности возбуждающих импульсов. На рисунке подан график изменения
волны на (k +1) -м периоде изменения параметров
среды к приведенной амплитуде a k прямой волны на k -м периоде (ось ординат) в зависимости от изменения длительности возбуждающих импульсов t1 (ось абсцисс) с использованием безразмерных параметров a, s = стр/ю ,t1 =rax1,t = raT .
Заключение. В работе показано, что возбуждение безграничной среды синхронными сериями одинаковых импульсов диэлектрической проницаемости и проводимости приводит к трансформации плоской электромагнитной волны, в результате которой поле на любом временном интервале представимо в виде суперпозиции прямой и обратной плоских волн. При этом в случае невозбужденного состояния среды фазовые характеристики этих волн и первичной волны совпадают. Также показано, что амплитуда напряженности электрического поля имеет резонансный характер в зависимости от периода изменения параметров среды.
Литература: 1. Борн М, Вольф Э. Основы оптики. М.: Наука, 1973. 719 с. 2. Владимиров В. С. Уравнения математической физики. М.: Наука, 1976. 527 с. 3. Nerukh A.G., ScherbatkoI.V., Nerukh D.A. Using evolutionary recursion to solve an electromagnetic problem with time-varying parameters // Microwave and Optical Technology Letters. 1997. Vol. 14, N 1. P. 31-36. 4. Аверков С.И., Болдин В.П. Волны в не диспергирующих неоднородных средах // Изв. вузов Радиофизика. 1980. T.23, N 9. C. 1060-1066. 5. Борисов В.В. Неустановившиеся электромагнитные волны. Л.: Изд-во ЛГУ, 1987. 240 с. 6. Борисов В.В. Трансформация электромагнитного поля при изменении проводимости среды во времени // Геомагнетизм и аэрономия. 1989. T.29, N 5. C. 730-737. 7. Нерух А.Г., Хижняк Н.А. Современные проблемы нестационарной макроскопической электродинамики. X.: НПО Тест-Радио. 1991. 280 c. 8. Нерух А.Г., Шаворыкина И.Ю. Расщепление электромагнитного импульса при скачке проводимости ограниченной cреды // ЖТФ. 1992. T.62, № 5. C. 108-118. 9. Harfoush FA, Taflove A. Scattering of electromagneticwaves by a materuial half-space with a time-varying conductivity // IEEE Trans. on Antennas and Propag. 1991. Vol.39, N 7. P. 898-906.
Поступила в редколлегию 12.03.98
отношения приведенной амплитуды ak+1 прямой 30
РИ, 1998, № 1
Нерух Александр Георгиевич, д-р физ.-мат. наук, профессор, зав. кафедрой высшей математики ХТУРЭ. Научные интересы: нестационарная электродинамика. Адрес: 310726, Украина, Харьков, пр. Ленина, 14, тел. (0572) 40-93-72.
Рыбин Олег Николаевич, аспирант кафедры высшей математики ХТУРЭ. Научные интересы: распростране-
ние электромагнитных волн в нестационарных средах. Адрес: 310726, Украина, Харьков, пр. Ленина, 14, тел. (0572) 40-93-72.
Щербатко Игорь Владимирович, доцент кафедры высшей математики ХТУРЭ. Научные интересы: нестационарная электродинамика. Адрес: 310726, Украина, Харьков, пр. Ленина, 14, тел. (0572) 40-93-72.
УДК 537.87; 621.371
ПРЕОБРАЗОВАНИЕ ЭЛЕКТРОМАГНИТНОГО ИМПУЛЬСА ВРЕМЕННЫМ ВОЗБУЖДЕНИЕМ СРЕДЫ В ПОЛУПРОСТРАНСТВЕ
Если E o(t, x) — первичное поле в невозмущенной среде, то электрическое поле в нестационарной области описывается интегральным уравнением Воль-терра второго рода [1,2,6]:
ГО ГО
E(t, x) = E 0(t, x) + J dt' J dxK (t, t', x, x')E(t', x'), (1)
0 0
РЫБИН O.H., CAXHEHKO H.K.
Рассмотрено преобразование электромагнитного импульса, вызванное скачкообразным изменением во времени параметров ограниченной среды. Изменение представляет собой прямоугольные импульсы диэлектрической проницаемости и проводимости произвольной длительности и амплитуды. Получены и проанализированы выражения для электрического поля во всем пространстве.
1. Введение
Распространение электромагнитных волн в нестационарных средах описывается уравнениями Максвелла с зависящими от времени параметрами. Как правило, такие задачи приходится решать численными методами [1,2]. Однако основные особенности эволюции преобразования электромагнитных сигналов нестационарными средами можно получить аналитическим путем [3-9].
Наиболее простым с точки зрения аналитического исследования является случай трансформации поля при скачкообразном изменении диэлектрической проницаемости или проводимости среды. Однако, как правило, рассматривается изменение только одного параметра. В настоящей работе исследуется воздействие на электромагнитный импульс синхронного импульсного изменения диэлектрической
проницаемости e(t) и проводимости c(t) среды, которое происходит в области полупространства x > о . Изменение во времени этих параметров описывается формулами
є() = єo[0(-t) + 0(t - т) + ei[0(t) - 0(t - т)]
CT(t) = ai[0(t) - 0(t - т),
где єо — диэлектрическая проницаемость невозмущенной среды; єі и сі — диэлектрическая проницаемость и проводимость среды в полупространстве x > о на временном интервале t є (0, т); 0(t) — единичная функция Хевисайда. Среда в полупространстве x < 0 все время остается непроводящей и имеет диэлектрическую проницаемость є0.
Предположим, что первичное поле падает по нормали на границу области полупространства x > 0 , имеет только составляющую, не зависящую от поперечных координат.
Здесь K (t,t',x,x') — ядро интегрального уравнения, K(t,t',x,x') = -(s(v0(t -1') - |x-x'|)/a2(t')j x
xjc(t') + j(1 - a2(t'))-dtj, (2)
где a(t) = д/є0/є(0; v0 = C;c(t) = 2nc(tVє1; c
— скорость света в вакууме; 8(t) — дельта-функция Дирака.
В области x < 0 поле определяется с помощью того же соотношения (1), которое в этом случае представляет собой квадратурную формулу, выражающую внешнее поле через внутреннее.
Решение уравнения (1) в области x > 0 записывается через резольвенту R(t,t',x,x') :
ГОГО
E(t, x) = E 0(t, x) + J dt' J dx' R (t, t', x, x')E o(t', x'), (3)
00
которая может быть найдена из уравнения R (t,t', x, x') = K (t, t', x, x') +
ГО ГО
+ Jdt'' Jdx'K(t,t",x,x")R(t",t',x",x'). (4)
0 0
В случае, если єД) = є і, c(t) = сі, резольвента представляется с помощью обратного преобразования Лапласа
- а+ іго dp
R (t, t', x, x') = е-сіт J ^nepT[Sl(p) + S2(p)]. (5)
а- іго
Здесь
Дг2 =-2
S1(P)=^p 2-сГ [p-а - -41 e-~vrv p 2 -Ci2
2v0 (p +сі a 2
S2(p) =
Wp2 - сі2 (
2v0
і p-сі
a у p + сі
-+Ф =-2
V1 v
p -сі
e
РИ, 1998, № 1
31