Алгоритм оценки стойкости радиоответчика к воздействию МЭМП приведен на рис. ПЗ.

ФОРМИРОВАНИЕ ФОРТРАН-ПРОГРАММЫ

Постановка задачи. Оценка параметров помех в приемном тракте (см. рис. П1), вызванных прохождением части МЭМП через АФУ.

Начальные сведения о РЭС и воздействующей МЭМП. Те же, что и в предыдущей задаче.

Используемые математические выражения:

= - j Sn (со) cos [со t- ф (со)] d со; (П2>

Wf = [Sn ( )]2 da; (ПЗ) to oil

Ua(t)=.E (t)h/2; (П4)

Для импульсных МЭМП, имеющих сравнительно низкочастотный спектр, источниками которых являются молнии, ЛЭП в режиме КЗ, и рабочего диапазона частот радноответчика выше 100 МГц справедливы следующие выражения:

р 2(Ti-т,)Д/£

Еп max =-;-;-. U Щ

(Ti -Та) Ejriax [sin (соа /) -Sin (щ t)] р /А 2 (Ti - Та) gmaa: sin (Асо f/2)

Eon (О =--; cnsj

Tl 2 нес

(-i-2)gL.A/(l2C + A/)

12яЧоТ?т(2/ 2(Ti-Г2)£>паж

12яЧо?П24с-АЛ)

12 нес -En.:

(ПИ)

В выражениях (П1) ... (ПИ) кроме общепринятых в книге использованы следующие обозначения: £п(). Ел max - напряженность электрической составляющей поля МЭМП, приведенная к полосе пропускания приемного тракта Дш=2яД/, и ее максимальное значение [при выполнении условия (П1) влиянием МЭМП на работу радиоответчика можно пренебречь], В/м; Eon{t), £ом(0-соответственно огибающие Eu{t) и £00(0. В/м (см. рис. 2.8); Wf -плотность потока энергии помехи, приведенная к полосе пропускания приемного тракта, Дж/м; й!а -энергия помехи на входе приемника, приведенная к полосе его пропускания, Дж; Ua {t) - напряжение помехи на входе пр:иемника, приведенное к полосе его пропускания. В; taan - время ухудшения, прекращения радиосвязи за счет прохождения части энергии МЭМП через АФУ, с; Т - абсолютная

Расчет промеж-уточных величин:

Расчет /Гптаз: Г

<0 =0 >о

Расчет: En,(/,P,W,

Печать: М7,/Гп тах,пм,пм/.

Расчет промежс/точных величин : ТП, ТП1

г- GOTO 17

-\Р0тУ=7,30

Лрисвоение начальных значений Y1, YZ. Расчет М,Н,К7

Расчет [Sn при ui=tOf и и>=Ы2

---\/7йв1 = 7.

7<1

Расчет функции \

величин ff,iVnp с точностью 7/S4

Печать W-f,W

------qy

I GOTO 75

ПО 8 KlM8/i9,M10

\ Расчет TI

I-----075J=7,30

Присвоение начальных значений YZ, Y4. Расчет М,Н,К4

Расчет SnM-cos (<yf - f((4} при ь>=ы. и ы=ы2 \

®

г>с----W9K31.K4

Расчет S CM)-ws(o)f - f (и ) \ с шааом Н /

и величин пГгУ, Uf, (t) , Pf, Cf)

Печать t, e(t), (f), Pf, (f)

i----

GOTO 8

Puc. П4. Структурная схема подпрограммы АНТ

-\ -

\ \расчет:С (tXBov (t),Eon (t),C/ Ш,РаШ

Рис. 175. Структурная схема подпрограммы ЦИКЛ

температура антенны. К; - коэффициент, учитывающий собственные шумы приемника; г]-коэффициент полезного действия АФУ; ыи Мг - крайние частоты полосы пропускания приемника, рад/с; Мнес = 2nfнес - несущая частота полезного сигнала, рад/с; 5п(са) и фп(а>)-модуль и аргумент спектральной плотностч МЭМП; Ti = l/ai, с; Т2=1/а2, с; п.доп - максимально допустимая амплитуда напряженности электрического поля помехи на входе АФУ в рабочем диапазоне частот приемника, не влияющая на его работу, В/м.

Структурные схемы подпрограмм АНТ и ЦИКЛ, а также подпрограмм-функций YF н №1, используемых для расчета Е.тах, En(,t), W, Ua(t), Wa, £оп(0, £ом(0 и зап, лриведены на рис. П4... П7, их покарточные распечатки на с. 248. Обращение к подпрограммам ЦИКЛ, YF и YF\ производится из подпрограммы АНТ. Обращение к подпрограмме АНТ имеет вид: CALL АНТ (ti, Т2, £m<ix, So, Т, ifaec. Л/, Г], N, йвф.а, п.доп. Ml, М2, МЗ, М4, М5, Мб, М7, S4, л, k, TBI, DI, М8, М9, МЮ), где т xs, /нес А/, k, D1 описываются как величины с двойной точностью; М1=/ при использовании подпрограммы АНТ в цикле Z)01/=l, N (N - целое число), в противном случае М =1; М2... М10 - величины, задающие шаг квантования времени t при выводе результатов расчета на печать (в частности, величины М2... М10 могут быть приняты равными М2=М4=М5=М7=М8=М10=1, М3=20, М6=М9=100); S4 -величина, задающая точность вычисления интегралов, входящих в выражения (П2) и (ПЗ), в частности при 54=100 интегралы вычисляются с погрешностью 17о; ТБ1=9-10*, Dl = 108 - часто встречающиеся в вычислениях константы.

ПОСТРОЕНИЕ И ИССЛЕДОВАНИЕ

МОДЕЛИ ИЗМЕРИТЕЛЯ ВРЕМЕНИ ЗАДЕРЖКИ

Постановка задачи. Исследование зависимости точности измерения от параметров измерителя и внешних дестабилизирующих факторов.

Начальные сведения. Измерение времени задержки осуществляется в соответствии с методом моментов [60], т. е. задержка определяется между фронтом импульса на входе исследуемой цепн н центром тяжести задержанного импульса на ее выходе (рис. П8). Структурно-функциональная схема измерителя приведена на рис. П9.

FUNCTION YFfX)

\CDMM0N/PZ/P7,P3J=7,P8,P3

Расчет [Spfe;] при различных ы

Рис. П6. Структурная схема подпро-.граммытфункцин YF

FUNCTION YF7(X)

Рис. П7. Структурная схема подпрограммы-функции YF1