книги из ГПНТБ / Ширман, Я. Д. Разрешение и сжатие сигналов
.pdfВводя обозначение т = |
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. |
|
|
||
|
В |
ш е с т о м примере остановим |
|||
ся |
на двухэлементной антенне (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 (t— tv), 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. |