Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Ширман, Я. Д. Разрешение и сжатие сигналов

.pdf
Скачиваний:
16
Добавлен:
22.10.2023
Размер:
13.3 Mб
Скачать

Вводя обозначение т =

I k и новую переменную

интегрирования

t= s (k + /) 772, получаем

 

 

 

 

 

M 1

 

CO

 

р (и f )-

Mtj,

2

2

*!2KkFT einmFr S

u o(*+

 

m=—M+I

k

—со

^

+ r+2mT ] U*Jt— x+2mT ) ei2nFi di I .

( 7 )

Слагаемые (7) с одинаковыми

значениями m образуют геометриче­

скую прогрессию вида е '2я/гЛТ.

При т > 0 значения k = I т

Рис. 2.3.5. Рельеф автокорреляцион­

Рис.

2.3.6. Рельеф автокорреляцион­

ной функции время — частота коло­

ной

функции время — частота пря­

кольного радиоимпульса с линейной

моугольного радиоимпульса с ли­

частотной модуляцией.

 

нейной

частотной модуляцией.

 

меняются от 0 до М т — 1, при т <

0 от \т\ до М — 1. Иначе,

в любом случае, значения /е изменяются

в пределах (\т \ т)/2

до

М — 1 — (I m I +

т)І2. При этом в правой части равенства (7) можно

выделить выражение

 

 

 

 

 

 

М—I—(|иі|+»0/2

 

n fT sin

(M—\m\)FT]

 

0jntnFT

V

 

 

 

к= (\„ц^т)/2

 

 

 

sin (nFT)

 

Окончательно получим

 

 

 

 

 

 

 

1

M—1

sin [jt (М— \m\)FT]

 

р (в

F)

У

 

Мти

т =^ М + 1

 

 

X

 

 

 

 

sin (nFT)

 

X

Uo t + x + mT U*

 

X+ niT

(8)

§ 2.3.3.

299

Рнс. 2.3.7. Рельеф автокоррсляционной функцни время — частота прямоугольной коге­ рентной пачки прямоугольных радиоимпульсов.

Выражение под интегралом представляет собой комплексную дву­ мерную автокорреляционную функцию одиночного радиоимпульса.

Рельеф тела неопределенности,

соответствующий

прямоугольной

форме

последнего,

представлен на

рис.

2.3.7.

 

 

 

В

ш е с т о м примере остановим­

ся

на двухэлементной антенне (Af=2)

с

достаточно большим

разносом эле­

ментов,

принимающей

в течение вре­

мени Т колебания, отраженные от движущейся радиолокационной цели. Вычислим пространственно-времен­ ную автокорреляционную функцию при значениях: угловой координаты

0, О -Ь ft, ее производной 0, б -]- ft. Предполагается, что ожидаемые и истинные значения координат и про­ изводных мало отличаются между собой, скорость радиального движе­ ния цели и соответствующая ей доп­ плеровская частота известны точно. Искомую двумерную функцию р =

= р (ft, ft) найдем из [(3), §2.1.8] н [(14), §2.3.1]. Ее можно получить,

вводя усреднение по времени, полагая ЛІ = 2 и заменяя 0 и 0' соот-

Рис. 2.3.8. Пример автокорреляционной функции угол — угловая скорость пространственно-временного сигнала для системы из двух ненаправленных ' антенн.

ветственно на 0 + 0/ и Ѳ + ft + (0 + ft)/. Пренебрегая величиной второго порядка малости ftft, после преобразований находим

 

 

 

r / 2

ПІ

 

 

 

 

COS

(ft sin 0 -j-ft/ cos0) dt

р ( f t ,

f t ) =

-

 

 

 

L

^ 0

 

или

 

 

T /2

 

 

 

/

 

 

 

 

 

 

 

sin (nlT$ cos Q/2\0)

 

 

 

cos (— •Osin©')

р ( f t ,

f t )

=

(9)

 

 

 

\, 70

j

я/rftcos Ѳ/2Яо

300

§ 2.3.3.

Тело р (Ф, ■&) имеет пик при Ф = 0 и Ь = 0 (рис. 2.3.8). Наилучшие

условия для разрешения по ■& имеют место, когда источник излучения расположен на прямой, соединяющей обе точки приема (cos Ѳ = 1). При М = 2 имеется резко выраженная неоднозначность по {К Тело вдоль оси •0’ имеет многопиковую структуру'. Степень неоднозначности уменьшается с увеличением М.

§ 2.3.4. ПРИМЕРЫ АНАЛИЗА ВОЗМОЖНОСТЕЙ РАЗРЕШЕНИЯ ЦЕЛЕЙ ПО КООРДИНАТЕ И ЕЕ ПЕРВОЙ ПРОИЗВОДНОЙ

Примем, что сигналы от целей являются когерентными. Анализ разрешения по координате и ее производной существенно не отли­ чается от проведенного для разрешения по координате (по одному па­ раметру).

При анализе разрешения по двум таким параметрам как время и частота, что обычно достаточно полно характеризует разрешение по дальности и радиальной скорости в радиолокации, используем, как и ранее, соотношения § 1.1.3. Положим, что комплексная амплитуда мешающего сигнала Ul (/) = U2 (t — т) е -і2лГІ, где U%it) — комплекс­ ная амплитуда полезного сигнала; т — сдвиг сигналов по времени; F — сдвиг по частоте. Решающая функция R (/) определяется по фор­ муле [(31), § 1.1.3], коэффициент использования энергии — по форму­ ле [(55), § 1.1.3]. Величина последнего существенно зависит от коэф­ фициента корреляции сигналов (см. § 1.1.3, 2.3.3).

В качестве п е р в о г о примера рассмотрим оптимальное разре­ шение двух колокольных радиоимпульсов без внутриимпульсной мо­

дуляции длительности т,„

отличающихся по времени запаздывания на

0,5 т п и по частоте на 0,5

(1/т„); длительность отсчитываем на уровне

е—л/4

0,46. Согласно [(1),

§ 2.3.3] при этом значение р « е~ я/ 4 я«

я» 0,46.

Коэффициент использования энергии

при сильном мешаю­

щем сигнале согласно [(55), §

1.1.3] составит k

1 — р2 дя 0,79. Если

бы разрешение проводилось

только по времени запаздывания (даль­

ности) или только по частоте (скорости), то значение р оказалось бы заметно больше р « е~ л/ 8 «я 0,68. Соответственно уменьшился бы коэффициент использования энергии k ля 1 — р2 я» 0,54. Одновремен­ ное разрешение по двум параметрам повышает, таким образом, коэф­ фициент использования энергии.

Оптимальная решающая функция R (/) при сильном мешающем сиг­ нале в рассматриваемом случае разрешения по двум параметрам имеет

вид

 

R(t) = U2( t ) - p U z ( t - x ) e - * * r ‘.

(1)

Эта функция характеризует опорное напряжение при оптимальном разрешении в случае корреляционной обработки.

В качестве в т о р о г о примера рассмотрим оптимальное разре­ шение двух прямоугольных линейно-частотно-модулированных радио­ импульсов. Задаваясь произведением тиД/ = 100, будем считать,

§ 2.3.4.

301

что временной сдвиг радиоимпульсов составляет т = 1,5/Д /= = 0,015/тп. В отсутствие различий по частоте согласно f(4), § 2.3.3] р та та 2/Зя л; 0,21. Мешающий сигнал в данном случае действует по первому боковому лепестку сжатого радиоимпульса. Опти­ мальная процедура разрешения (1) связана с подавлением первого бо­ кового лепестка. Даже в случае сильного мешающего сигнала коэф­ фициент использования энергии оказывается высоким k та 1 — р3 та та 0,96.

Расстройка по частоте может не только улучшить, но и ухудшить разрешение. Пусть F — — 1,5/т„ = —0.015Д/, тогда р та 1 — |т|/т„ = = 0,985. Коэффициент использования энергии k = 1 — pä та 0,03 оказывается низким потому, что эффект сдвига по частоте скомпен­ сировал в данном случае эффект временного сдвига, что связано с осо­ бенностью тела р (т, F) 4M радиоимпульсов. Величина р оказалась такой же, как для прямоугольных радиоимпульсов без частотной мо­ дуляции, разнесенных по времени на т = 0,015 ти.

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

время Т =

0,1 с. Пусть расстояние между антеннами I = 200 м. При­

мем далее,

что в формулах 1(12), (14), § 2.3.1] А0 =

0,1 м, 0о—

90°,

ѵ1 = 0, ѵ2 = 250 м/с, г = 50 км. Разрешение по углу

при 0 = 0

со­

гласно указанным формулам отсутствует. Отсутствует и разрешение по радиальной скорости. Имеется только отличие в скоростях изменения

угловой координаты: O' = 250/50 • 103 = 0,005 рад/с для одного излу­ чателя и Ѳдля другого. Согласно [(9), § 2.3.3] р = sin ^ j у я» 0,21,

так что возможный коэффициент использования энергии при сильном мешающем сигнале k та 1 — р3 я» 0,96.

§ 2.3.5. ОПТИМАЛЬНОЕ ОБНАРУЖЕНИЕ ДВИЖУЩЕЙСЯ ТОЧЕЧНОЙ РАДИОЛОКАЦИОННОЙ ЦЕЛИ В ОБЛАКЕ ПАССИВНЫХ ОТРАЖАТЕЛЕЙ. ПОСТАНОВКА ЗАДАЧИ

Обнаружение объекта в протяженном облаке распределенных отражателей уже рассматривалось в § 1.2.3 как характерный пример, иллюстрирующий методы теории разрешения. В этом параграфе, од­ нако, обсуждалось одно только разрешение по дальности; различие радиальных скоростей цели и мешающих отражателей не учитывалось. Поэтому ставится задача дальнейшего развития проведенного анали­ за с учетом возможной селекции по скорости. Чтобы выявить влияние нестационарного распределения помехи в пространстве, далее рассма­ тривается более широкая, чем в § 1.2.3, совокупность предположе­ ний о структуре облака. Наконец, учитывается обзор пространства, влияющий на возможности разрешения цели и отражателей (рис. 2.3.9).

Постановка задачи заключается в следующем [137].

Исследуется функция /?(/), описывающая закон оптимальной ко­ герентной обработки сигнала, отраженного от радиолокационной

302

§ 2.3.5

точечной цели при импульсном зондировании. Она определяется из интегрального уравнения [(46), § 1.1.3]

/? (0 + 1' ф(*> s)R (s)d s= UA(l),

(1)

где Uа (і) — ожидаемый полезный сигнал; cp (t, s) — нормированная (относительно спектральной плотности шума N 0) автокорреляцион­ ная функция отражений от независимых эле­ ментов облака:

Ф (t,

s) = ± ^

U Ai(t)CrAJ(s)

( /=

1, 2 ,...).

 

 

/

 

 

Мешающий

сигнал от /-го элемента обла­

ка в

приближении [(5), (6), §

2.3.2]

будет:

UAi (1)=ѴЩА (f-£,-fy)l7(f-fy) е - /2я^

Как и полезный

сигнал

UA (/) =

А (t) U (/),

 

 

 

 

он формируется

в

соответствии

с функция­

 

 

 

 

ми: U(t), задаваемой передатчиком, и А (t)

Рис.

2.3.9.

Пояснение

системой обзора по угловой

координате. Ве­

процесса

образования

личина Dj характеризует отношение эффек­

сигнала от цели

и обла­

тивных поверхностей /-го элемента облака и

ка

отражателей

(1р —

цели. Величины £;-,

йу--—запаздывания,

зави­

протяженность элемента

сящие от сдвига /-го элемента облака и цели .

разрешения

по

дально­

 

сти).

 

(по угловой

координате

и

дальности

соот­

 

 

 

 

 

 

ветственно).

Величина

Fj

— допплеровское

смещение

частоты,

При непрерывном распределении элементов облака с относитель-

ной плотностью D (£, й, F) получим

 

ф(/,

S): ^ ä j D ( t , й, Л Л ( * - £ - й ) ( 7 ( * - й ) х

 

X

і4* (s —£—й) О* (s —й) е - / 2лГ ~ s>d ld M F .

(2)

Распределения D (£, й, F) будем задавать произвольно как некото­ рые математические модели облака. При составлении моделей обычно стараются согласовать противоречивые требования приближения мо­ дели к действительности и упрощения математического анализа. На­

ряду с этим постараемся согласовать получаемые

данные

с опубли­

кованными в различных изданиях. В этой связи зададим

ч е т ы р е

основные м о д е л и

облака, предполагая частные распределения по

£, й, F независимыми:

 

 

 

D06(Q8($)y.(F),

(а)

 

P (t,

А >б(й)Ѵ( П

(б)

(3)

й, F) =

(в)

 

ö 06(£)e-°**Y(f),

 

 

Doe-(a#*m*)T (F).

(г)

 

§ 2 . 3 .5

303

Здесь у (F) — спектральная плотность модуляционного множителя помехи, полагаемого стационарным во времени.

П е р в а я м о д е л ь (За) соответствует сосредоточенному рас­ пределению помехи по дальности н угловой координате, когда отра­ жатели расположены в общем с целью разрешаемом объеме. Такую модель будем называть моделью сосредоточенной помехи. Оказывается, что упрощающие предположения, используемые отдельными авторами,

сводят задачу

именно к этой модели.

В т о р а я

м о д е л ь

(36) соответствует сосредоточенному рас­

пределению по дальности,

но неограниченному по угловой координате.

Протяженность помехи по этой координате много больше элемента углового разрешения.

Т р е т ь я м

о д е л ь (Зв)

соответствует сосредоточенному распре­

делению помехи

по угловой

координате и рассредоточенному по

дальности. Протяженность облака может изменяться за счет выбора

параметра а.

Случай а = О или е- “*2 = 1 соответствует безгранич­

ному облаку.

Этот случай назовем н у л е в ы м а с и м п т о т и ч е -

с к и м п р и б л и ж е и и е м. Полагая величину а отличной от нуля, но достаточно малой, так что (аналогично § 1.2.3).

 

 

е - “»3^ 1— ай3,

(4)

приходим

к

п е р в о м у а с и м п т о т и ч е с к о м у

п р и б л и ­

ж е н и ю .

 

Оно применимо, если подынтегральное выражение (2) об­

ращается в нуль при больших Ф, когда соотношение (4) перестает со­ блюдаться. Пользуясь первым приближением, можно установить пре­ делы применимости нулевого приближения. Далее убедимся, что ну­ левое приближение оправдано (как и первое), если полная длитель­ ность зондирующего колебания меньше протяженности облака в еди­ ницах временного запаздывания.

Чаще, однако, встречается случай, когда протяженность облака в единицах временного запаздывания велика по сравнению с длитель­ ностью зондирующего радиоимпульса, но мала по сравнению с дли­ тельностью пачки радиоимпульсов (либо части этой пачки в виде 2, 3 и т. д. импульсов, используемых для компенсации помехи). Тогда от

третьей модели

естественно перейти к ч е т в е р т о й (За), которая

соответствует

большим, но конечным протяженностям облака и по

дальности и по угловой координате'. Имеется в виду протяженность по дальности в единицах времени запаздывания , которая заметно пре-, вышает длительности радиоимпульсов, ио много меньше полупериода следования Т/2. Четвертая модель характеризуется поэтому непере­ крывающимися отражениями радиоимпульсов пачки. Это обычно и имеет место. (Для третьей модели отражения многократно перекры­ ваются, что может наблюдаться при квазинепрерывном излучении). Протяженность облака по угловой координате в единицах времени по­ ворота оси антенны полагается много большей длительности пачки импульсов. В отличие от третьей модели это обеспечивает неодинако­ вое число отражений от цели и облака.

304

§ 2.3.5

С учетом сказанного в дополнение (4) имеем для четвертой модели

e-ߣ’ Ä l—ßts.

(5)

Соотношением (4) пользуемся только при |'&| < 772, при |Ф| ^

772

полагаем

 

e -aö=Ä о.

Являясь простейшими и учитывая нестационарность облака, при­ веденные модели, конечно, не исчерпывают всех особенностей реаль­ ных ситуаций. Так, иногда имеют место детерминированные зависи­ мости у (F) от угловых координат и дальности за счет движения бор­ тового радиолокатора относительно подстилающей поверхности, вы­ сотных и угловых перепадов скоростей ветра, переносящего отража­ тели облака и т. д. На этом остановимся очень кратко в § 2.3.8.

Считаем, что модулирующая функция для всех четырех моделей

U ( t ) = Z U 0( t - t v)

(6)

V

 

состоит из отдельных периодически следующих неперекрывающихся импульсов U0 (ttv), tv — моменты их прихода (ѵ = 0 , ± 1 , ± 2 , ...). Центр тяжести (или начало) импульса U0 (t) условимся относить к мо­ менту t = 0. Для коротких импульсов

А (t) U0 (t - t v ) « A (tv) UQ( t - t v ) = A v U0 (t - tv),

число импульсов в пачке считаем при этом ограниченным.

Задачей анализа является исследование решений уравнения (1)

при моделях облака (За), (36), (Зв), (Зг)

и импульсном зондирова­

нии (6).

 

§ 2.3.6. РАСЧЕТНЫЕ СООТНОШЕНИЯ ДЛЯ РАЗЛИЧНЫХ

МОДЕЛЕЙ ОБЛАКА

Необходимые расчетные соотношения

получим как решения [(1),

§ 2.3.5]. 'Чтобы упростить вычисления,

предварительно

введем н о р ­

м и р о в а н н ы е

о д н о м е р н ы е

а в т о к о р р е л я ц и о н ­

ны е функции:

1)случайных флюктуаций отражений от облака для вспомогатель­ ного случая, когда оно облучается гармоническим колебанием на не­ сущей,

оо

р Д т )= \ у.(Р)е-’2лП dF-,

(1)

оо

2)неслучайной огибающей U0 (t) одиночного радиоимпульса

 

оо

 

Ро(т) = С0

$ U0( t ) U t ( t - x ) d t ,

(2)

 

—со

 

■гдеС0 = 1/230 — нормирующий множитель; Э0 —энергия импульса;

§ 2 .3 .6

305

3) неслучайной огибающей А (() пачки радиоимпульсов

со

Рл|(т) = $ A ( t ) A * ( t~ T ) d t

(3)

— со

 

(нормирующий множитель здесь опущен; имеется в виду выбор масш­ таба А (/), при котором рл (0) = 1);

4) неслучайной огибающей

XJa (/) А (() U0(t tv)

совокуп-

 

V

 

ности радиоимпульсов, образующих пачку,

 

 

со

 

Р£м(т) = С£М

S UA(t)U*A(t— z)dt,

(4)

 

— со

 

где СиА = 1/2Эил\ Зил — энергия пачки радиоимпульсов.

 

Решению уравнения [(1), § 2.3.5] должно предшествовать в ы ч и с-

л е н и е автокорреляционных функций отражений от облака ср (/, s) для различных его моделей по формуле [(2), § 2.3.5].

Подставляя в эту

формулу

выражение [(За), § 2.3.5)], соответст­

вующее п е р в о й м о д е л н ,

и используя (1) и [(6), (7), § 2.3.5], по­

лучаем

 

 

 

 

 

 

 

 

 

 

 

Ф(/, s) =

- ^

2

А%, А* pF(t — s) U0(t — tv) U*0(s~t,).

(5)

 

 

N0

V. I

 

 

 

 

 

 

При коротких радиоимпульсах

в области отличных от нуля их комп­

лексных амплитуд значения медленно меняющейся

нормированной

автокорреляционной

 

функции

р/.- (/ — s) сводятся

к

коэффициентам

межпериодной корреляции

помехи pF (/ѵ — /;) =

Pfvi

в отсутствие

последовательного обзора по угловой координате.

 

 

 

 

Для

в т о р о й

м о д е л и ,

используя (1), (3)

и [(2),

(36),

(6), (7)

§ 2.3.5],

находим

 

 

 

 

 

 

 

 

 

Ф(/, s ) = - ^ 2

;

pF( t ~ s ) p A { t~ s ) U0{t— tv)Ub{s — ii).

(6)

 

« o v ,

 

 

 

 

 

 

 

Вошедшее в выражение (6) произведение pF (/ — s) рд (t s) при ко­ ротких радиоимпульсах заменяется коэффициентом межпериодной корреляции с учетом последовательного углового обзора:

 

РрАѵі = Р р ^ ѵ

Р л ^ ѵ U)-

 

Для т р е т ь е й

м о д е л и , используя (4) и [(2),

(Зв), (6), § 2.3.5],

а также асимптотическое соотношение [(4), § 2.3.5], получаем

Ф

(t, s)

= ф (/ — s) +

афа it, s).

(7)

Здесь ф — s) — функция

разностного

 

аргумента,

характерная для

бесконечно протяженного облака,

 

 

 

 

^ (t) =

^0 ЭиА Рк (т) Рсіа (т);

(8)

 

 

306

§ 2.3.6.

•фа (t, s) — функция, учитывающая поправку на конечную его протя­ женность,

00

% (t, s) = ~ - ^ p F( i~ s )

Л'о

\ UA( t - ü ) U Z { s - e ) W d Q ,

(9)

^

 

со

а— величина, обратно пропорциональная квадрату протяженности облака.

Для ч е т в е р т о й м о д е л и , используя [(2), (Зг), (6), §

2.3.5]

и учитывая малую длительность импульсов, находим

 

cp(t, s) = -jrPF(t— s ) ^

^

е - ( а Ф 3 + Р£3) х

 

V ,

I _

с

 

х А у — 1— й ) А * ( 8 - £ — Ъ)и0у — Ъ— Ъ ) т ( 8 — ^ — Ъ){1М£.

(10)

Напомним, что интервалы между импульсами существенно превы­ шают их длительность. Для заданной пары значений t, s найдется по­

этому не более одной пары чисел ѵ, /, такой,

что разности 11 tv \

и | s — ^l' будут заметно менее полупериода

следования импульсов.

В выражении (10) сохраняется всего одно слагаемое суммы по ѵ, /, име­ ющее вид двойного интеграла. Вычисляя значение последнего, восполь­

зуемся асимптотическими приближениями [(4), (5), §

2.3.5]. После

замены переменных интегрирования получим

 

 

 

 

c p (t,

s) = ~ ^

3 oPF (t — ^ P a V — s ) P o (^ —

?v s

+

* i ) +

 

 

 

-f-aria (t, s) + ßt]ß(*, s).

 

 

 

 

Заменив

tv — tt =

(v — l)T0 = mT0, (m =

0, ± 1 ,

± 2 , ...),

вве­

дем искусственно знак суммы, поясняющий переход

к функции

раз­

ностного аргумента,

 

 

 

 

 

 

Po(t tv — s + tl) = '2l p0{t— s — mT0).

 

 

(Юа)

 

 

«

 

 

 

 

Введение знака суммы ничего не меняет: при небольшой длительности импульса и фиксированных t, s в этой сумме отлично от нуля не бо­ лее чем одно слагаемое суммы. Окончательно имеем для четвертой мо­ дели

ср(£,

s) =

т)(i

s ) ат)а (t, s)-fßtiß(f, s),

(11)

где

 

 

 

 

 

Л М = - ^ £- 5 0р/г(т)рл (т) S po(t—т Г 0),

(12)

 

 

Л'о

tn

 

Ч а ( * . S ) = = — ^ - p F {t — S) pA (t — S ) X

 

 

 

y* /2

 

/Vo

 

xS,

 

 

 

 

 

°(

u 0 {t — tv — #) ü? ( s —

(13)

V , *

 

J

 

 

 

 

-Го/2

 

 

 

§ 2.3.6.

307

 

Щ{і, s)

2Dp Э0рF (t— S) 'Epo(t— s — mT0) x

 

 

 

 

W„

 

 

 

 

rn

 

 

 

 

 

 

X 5 Л ( / - £ ) Л * ( 5 - С К 3 <*£.

 

 

( 1 4 )

Запаздывание ■&огибающей A(t)

в (14) не учитывается.

 

 

 

В соответствии с (5), (6), (7), (11) ядро <p (t, s) уравнения [(1), § 2.3.5]:

в ы р о ж д е н н о е

для

моделей

1,

 

2; р а з н о с т н о е в

н у л е ­

в о м

п р и б л и ж е и и и для моделей 3,

4.

 

 

 

 

В

случае в ы р о ж д е н н о г о

 

 

ядра

(5)

или (6)

решение

[(1),

§ 2.3.5] найдем как линейную комбинацию

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

(15)

 

 

 

 

 

 

 

 

 

 

 

 

с неопределенными

коэффициентами

Н т.

Подставляя

(15)

в

[(1),

§ 2.3.5], умножая на

U0 ( t — /д)

и интегрируя в бесконечных преде­

лах,

находим для

Н т

систему

уравнений,

аналогичную

[(31),

§ 2.2. 11],

 

 

 

 

 

 

 

 

 

 

 

 

ШDnw Нт =

 

(і1==1>

2, ...),

 

 

(16)

где для моделей (1) и (2)

соответственно

 

 

 

 

 

 

 

 

б/пц+ —г- Э0 Ат р _

,

 

 

 

 

 

 

 

Wo

 

 

Fmß'

 

 

 

 

Пяі(Х —'

 

 

 

 

 

 

 

(17)

 

-’ mß Т

2£>о

Эо Pp

 

 

 

 

 

 

 

W„

 

F A m ß '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь öjnn — символ Кронекера (1 при m =

ц и 0 при m Ф р).

 

Решение интегрального уравнения

[(1),

§ 2.3.5] для

р а з н о с т ­

н о г о

ядра с малыми неразностными поправками найдем в виде сум­

мы членов нулевого и первого порядков; для модели 4, например,

 

 

R(t) = R0(t) + aRa(t) + $Rß(t).

 

 

(18)

Подставляя (11) и (18) в [(1), § 2.3.5] и пренебрегая асимптотиче­ ски малыми второго порядка, приравняем порознь в обеих частях ра­ венства три группы величин: нулевого порядка, пропорциональные а или ß. Тогда придем к системе трех независимых интегральных урав­

нений для /? 0 (/), Ra И 7?ß (ty.

R 0{t)+

J

R 0{s)r\{t — s)ds = UA(t),

 

 

— * 03

 

 

 

oo

 

 

R a { t) +

5

Ra(s)l\(t— s)ds=Ua(t),

(19)

 

—00

 

 

 

oo

 

 

# ß ( 0 +

$ Rß(s)r](t—s)ds-=[/ß(t).

 

 

—00

 

 

308

§ 2.3.6.

Соседние файлы в папке книги из ГПНТБ