книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfСовместное распределение вероятностей последова тельностей дискретных переменных її, как для шума, так и для сигнала с шумом, при преобразовании нестацио нарного процесса определяется из формулы (1.7). Коэф фициент правдоподобия при проверке двух простых гипо тез с учетом соотношения (1.7) принимает вид
t  | 
	П її Р№  | 
	(*>  | 
Л ( ц ) = p ( " l 6 , ) =  | 
	* = 0 a = 1  | 
	(2.188)  | 
п п 'bw<*>
Логарифмируя (2.188), получаем
шл й=gs (А) щ  | 
	=s s < *> ^ <*  | 
А = 0 а = !  | 
	А = 0 а = 1  | 
  | 
	(2.189)  | 
Поскольку значения л ; а , а = 1 , г ,  | 
	для каждого /г =  | 
= 0,/г—1 образуют полную группу несовместимых собы
тий, вторая  | 
	сумма в (2.189) содержит для данного  | 
|||||||
/г =  | 
	0, п—  | 
	1  | 
	всегда  | 
	одно  | 
	слагаемое,  | 
	соответствующее  | 
||
попаданию  | 
	случайной  | 
	величины щ в интервал  | 
	квантова  | 
|||||
ния  | 
	Д« .  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	а  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	Структурная схема  | 
	алгоритма (2.189)  | 
	представлена  | 
|||||
на  | 
	рис. 2.36.  | 
	Дискретные  | 
	переменные  | 
	йи  | 
	i = 0,  | 
	1, 2, ... ,  | 
||
с выхода преобразователя записываются в /п-канальный
регистр сдвига (Рг.  | 
	с)  | 
	с числом разрядов  | 
	п.  | 
	Каждое  | 
||||
выборочное  | 
	значение  | 
	нг- параллельным кодом  | 
	вводится  | 
|||||
в регистр адреса постоянного запоминающего  | 
	устройства  | 
|||||||
( П З У ) , хранящего  | 
	весовые  | 
	коэффициенты  | 
	Wa  | 
	(k).  | 
||||
В каждом из п ПЗУ  | 
	хранится  | 
	по 2т  | 
	различных  | 
	ве  | 
||||
сов Wi(k),  | 
	W^ik),...,  | 
	W.r{k). С помощью  | 
	параллельного  | 
|||||
сумматора образуется сумма весовых коэффициентов, которая сравнивается с порогом в цифровом компарато ре ( Ц К ) . Такая схема используется при обнаружении сигнала с неизвестным временем прихода.
Если последовательности дискретных переменных ІЇі для обеих гипотез соответствуют стационарным случайISO
u(t)  | 
	П  | 
	и. (к)  | 
	
  | 
	Рг.с  | 
	
  | 
  | 
	
  | 
	
  | 
|||
  | 
	
  | 
	п-1  | 
	п-2  | 
	0  | 
|
  | 
	
  | 
	1  | 
|||
  | 
	
  | 
	Сдвиг  | 
	Щп-2) .  | 
	
  | 
	
  | 
ЦК
Рис. 2.36. Структурная схема алгоритма (2.189).
НЬІМІ процессам  | 
	«(•/), то в соответствии с  | 
	(1.9)  | 
	
  | 
||||||||
  | 
	іплй =sА .1 п -^=Е X J «•  | 
	( 2 Л 9 0 )  | 
|||||||||
Структурная  | 
	схема,  | 
	
  | 
	соответствующая  | 
	
  | 
	алгоритму  | 
||||||
(2.190) , содержит г счетчиков, определяющих  | 
	
  | 
	Яа ; г ПЗУ,  | 
|||||||||
хранящих веса  | 
	Wа; г  | 
	умножителей  | 
	(XaWa)  | 
	и  | 
	
  | 
	параллель  | 
|||||
ный сумматор с г входами.  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||||
Более удобна при реализации другая,  | 
	эквивалентная  | 
||||||||||
интерпретация  | 
	алгоритма  | 
	(2.190)  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
  | 
	
  | 
	Ш Л ( « ) = 5  | 
	£ •  | 
	da(k)Wa.  | 
	
  | 
	
  | 
	
  | 
	(2.191)  | 
|||
На рис. 2.37 приведена структурная схема  | 
	
  | 
	алгоритма  | 
|||||||||
(2.191) . Дискретные переменные гТг- с выхода  | 
	преобра  | 
||||||||||
зователя (П) поступают  | 
	в регистр  | 
	числа  | 
	ПЗУ, храняще  | 
||||||||
го весовые коэффициенты  | 
	Wa,  | 
	ia=l,  | 
	г. С помощью на  | 
||||||||
капливающего  | 
	сумматора  | 
	образуется  | 
	сумма  | 
	(2.191), ко  | 
|||||||
торая  | 
	сравнивается  | 
	с  | 
	порогом  | 
	с  | 
	либо в  | 
	
  | 
	цифровом  | 
||||
компараторе (ЦК) , либо  | 
	в сумматоре, если  | 
	в него был  | 
|||||||||
введен  | 
	обратный код числа с.  | 
	При этом  | 
	отрицательные  | 
||||||||
весовые коэффициенты удобно хранить в ПЗУ в допол нительном коде.
При малых интервалах квантования Д'«<^; (0,5ч-0,7)а, где а—среднеквадратическое значение помехи или сиг-
181
alt)  | 
	ЦК  | 
ПЗУ  | 
Рис. 2.37. Структурная схема алгоритма (2.191).
нала с помехой, из (2.46) получим
  | 
	р а  | 
	=  | 
	j" ю,  | 
	(ut) dui  | 
	^  | 
	wl  | 
	( « * J Лив , а =  | 
	ГТгТ  | 
	
  | 
|
где « * а £  | 
	Диа  | 
	— средняя  | 
	точка  | 
	интервала  | 
	Дма . Тогда ве  | 
|||||
совые  | 
	коэффициенты будут  | 
	равны  | 
	
  | 
	
  | 
	
  | 
|||||
  | 
	
  | 
	W (k) =  | 
	ln-  | 
	777- =  | 
	In  | 
	.- 5,  | 
	v - •  | 
	(2 192\  | 
||
Алгоритм  | 
	('2.189)  | 
	с учетом  | 
	(2.192)  | 
	соответствует  | 
	(2.55).  | 
|||||
В  | 
	задачах  | 
	обнаружения  | 
	некогерентных сигналов на  | 
|||||||
фоне шумов при независимых отсчетах процессов наи большее распространение получило бинарное квантова
ние. Распределения  | 
	вероятностей и алгоритмы  | 
	обработки  | 
||||
при бинарном  | 
	квантовании  | 
	получаются  | 
	по  | 
	формулам  | 
||
(2.188) — (2.192)  | 
	при подстановке  | 
	a='0,I,  | 
	г=2.  | 
	Алгоритм,  | 
||
соответствующий (2.191), имеет вид  | 
	
  | 
	
  | 
||||
In Л (и) = 2  | 
	К (k) W0  | 
	(k) +  | 
	dx (k) Wx (k)}.  | 
	
  | 
||
Следует учесть, что теперь dB(k)=l  | 
	—dt{k).  | 
	Тогда  | 
||||
In Л (и) = §  | 
	* {W0 (k) + rf, (k) [W, (k) -  | 
	Wa (k)\) =  | 
||||
- q + S ' - w  | 
	?™тЛ  | 
	
  | 
	° q + S " - w y w -  | 
|||
fe=o ft=o  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	(2.193)  | 
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
где величина
Q = S V ( * ) = const
ті ее следует отнести к пороговому уровню.
Реализация алгоритма (2.193) аналогична реализа ции,182 представленной на рис. 2.36, но вместо /и-каналь-
його регистра сдвига используется одноканальний и каж дый из весовых коэффициентов W(k) принимает одно фиксированное значение. По сравнению с алгоритмом (2.68) весовая функция теперь определяется не произ вольно, а наилучшим образом, хотя, как показывает со ответствующий анализ, эффективность подобных обнару жителей слабо зависит от весовой функции. Можно по казать i[95], что при обнаружении слабых некогерентных сигналов W{k)^G\, k = 0, п—1, т. е. весовая функция равна квадрату диаграммы направленности по мощно сти. При изменении 'интенсивности сигналов весовая функция будет также изменяться, хотя с практической точки зрения это не существенно.
Таким образом, для независимых выборочных значе
ний сигнала и шума  | 
	рассматриваемый метод синтеза да  | 
||
ет те же  | 
	алгоритмы  | 
	обнаружения, что и  | 
	метод синтеза  | 
цифровых  | 
	эквивалентов по аналоговому  | 
	прототипу (см.  | 
|
§ 2.3).  | 
	
  | 
	
  | 
	
  | 
Анализ эффективности некогерентных обнаружителей бинарно-квантованных сигналов проведен в § '2.3. Сле дует указать, что алгоритм (2.193) соответствует опти мальному обнаружению целей при отсутствии флуктуа ции и при' шумоподобных флуктуациях. В случае друж ных флюктуации пачки эхо-сигналов коэффициент правдоподобия определяется в соответствии с (2.56). В результате получается алгоритм, реализующий макси мально правдоподобное обнаружение (2.8) с весовой функцией, зависящей от интенсивности отраженного сиг нала. Однако такое усложнение системы обработки, как отмечалось в п. 2.3.3, практически не оправдано.
Перейдем теперь к вопросам обнаружения некоге рентных импульсных сигналов методом последователь ного анализа. Преимущества последовательной процеду ры Вальда по сравнению с процедурой Неймана — Пир сона хорошо известны [61]. Однако реализация алгоритма последовательного анализа в, многоканальных многоце левых Р Л С стала возможной при использовании фазиро ванных антенных решеток, программного обзора про странства и применении ЭЦВМ для обработки информа ции.
Последовательные процедуры в многоканальных си стемах можно разделить на три группы:
I . С игнорированием разрешающей способности.
183
2.С однократными пересечениями (независимым при нятием решения в каждом дальномерном канале и с по глощающими порогами). При этом в каждом канале вы числяются парциальные коэффициенты правдоподобия, значения которых сравниваются с двумя порогами. После достижения одного из них выносится соответствующее решение 'и испытания в данном канале прекращаются до окончания процедур в остальных каналах.
3.С многократными пересечениями порогов (иепо-
глощающис пороги). Испытания в каждом канале про должаются до тех .пор, пока во всех каналах значения парциальных коэффициентов правдоподобия не выйдут из зоны «неопределенности». Это системы с вынужден ным продолжением испытаний в каналах.
Алгоритм процедур первой группы [96] связан с обра зованием коэффициента правдоподобия по формулам (2.3), (2.4), где параметр и соответствует времени за паздывания отраженного сигнала Т. Практическое ис пользование подобных алгоритмов сильно ограничено требованиями наличия только одной цели и независи мостью дальномерных каналов. Кроме того, такая систе ма не обладает разрешением по дальности.
Процедуры второй группы обладают разрешающей способностью и могут использоваться при любом числе целей на одном направлении антенного луча. Ограничим ся рассмотрением процедур второй группы только с би нарным квантованием некогерентных сигналов. Лога рифм коэффициента правдоподобия в каждом дально
мерном канале  | 
	с учетом  | 
	(2.189)  | 
	можно представить  | 
|||
в виде  | 
	
  | 
	
  | 
	
  | 
	4^f"  | 
||
<«>=S \  | 
	+ (1 -  | 
	d (k)) In  | 
||||
In Л , (и) ) ) \  | 
	dС (k) In ^  | 
	
  | 
	
  | 
|||
fc=l L  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	ln \ ~ P s »  | 
	n,  | 
||
PN  | 
	Л—РЫ  | 
	1 —PN  | 
||||
  | 
||||||
  | 
	
  | 
	
  | 
||||
(2.194) где Ai — число единиц в выборке объема п. На каждом
—>
шаге 1пЛп(й) сравнивается с двумя порогами (2.19). Путем преобразования (2.194) можно перейти к весо вому варианту последовательной процедуры [61] с поро гами, заданными величинами
В  | 
	' = 0 '  | 
	A - = S =  | 
	. . K l - X w - ^ , ! •  | 
	< 2 ' 1 9 В >  | 
184  | 
	
  | 
	
  | 
	
  | 
	
  | 
с которыми на каждом шаге сравнивается величина
. . , , ~  | 
	— i n В  | 
	.  | 
	^ [ p s N ( l - p N ) / p N ( \ - p S N ) ] - i  | 
	
  | 
|||
^ • ( ^ - 1 ^ + — i n m - P s N w - p N ) ] —  | 
	х  | 
||||||
п  | 
	
  | 
	
  | 
	п  | 
	
  | 
	п  | 
	
  | 
	
  | 
Х ' £  | 
	d(k)-  | 
	J ] ( l - d ( * ) ) = r 0 + 2][d(*)r -  | 
	
  | 
||||
ft=l  | 
	
  | 
	/6=1  | 
	
  | 
	ft=l  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	-  | 
	(1 -d.{k))s\.  | 
	(2.196)  | 
||
Первое  | 
	слагаемое  | 
	(2.196)  | 
	интерпретируется  | 
	как на  | 
|||
чальный вес процедуры го, а  | 
	множители перед  | 
	сумма  | 
|||||
ми — как веса  | 
	г и s,  | 
	добавляемые при каждом  | 
	появле  | 
||||
нии единицы или нуля  | 
	соответственно.  | 
	
  | 
	
  | 
||||
Процедуры второй группы неоптимальны, поскольку  | 
|||||||
информация,  | 
	поступающая с  | 
	момента окончания  | 
	испы  | 
||||
тания в отдельном  | 
	канале до завершения испытаний во  | 
||||||
всех М каналах, не используется. Процедуры с много кратными пересечениями используют эту информацию,
так как окончание испытаний происходит  | 
	одновременно  | 
|
во всех  | 
	каналах. Среднее число выборок для систем вто  | 
|
рой и  | 
	третьей групп будет определяться  | 
	наибольшим  | 
средним числом выборок, требуемых для элемента с мак
симальным временем  | 
	анализа. Теоретически  | 
	процедуры  | 
|
с многократными пересечениями имеют некоторый вы  | 
|||
игрыш по сравнению с  | 
	процедурами  | 
	второй  | 
	группы, од  | 
нако, как показывают  | 
	расчеты, этот  | 
	выигрыш очень мал.  | 
|
Для анализа эффективности последовательных про
цедур при бинарном квантовании воспользуемся  | 
	весьма  | 
|||||||||||
простым и наглядным  | 
	методом  | 
	[97] . Последовательная  | 
||||||||||
процедура от шага к шагу представляет собой  | 
	случай  | 
|||||||||||
ное -блуждание, описываемое цепью Маркова.  | 
	Каждому  | 
|||||||||||
состоянию  | 
	Ах=Аы  | 
	такой  | 
	цепи  | 
	ставится в  | 
	соответствие  | 
|||||||
комбинация  | 
	из k единиц и / нулей,  | 
	которая  | 
	получается  | 
|||||||||
после  | 
	п-го  | 
	испытания  | 
	
  | 
	данного  | 
	элемента  | 
	разрешения  | 
||||||
(n = k-\-t).  | 
	Комбинации,  | 
	содержащие одинаковое  | 
	число  | 
|||||||||
единиц и пулей, но с различным  | 
	их чередованием,  | 
	отно  | 
||||||||||
сятся к одному состоянию. Вероятность попадания  | 
	в со  | 
|||||||||||
стояние  | 
	Аы  | 
	
  | 
	равна  | 
	р(Ам)  | 
	=С(п,  | 
	k)ph{\—p)n-k.  | 
	
  | 
	
  | 
	
  | 
	Число  | 
||
С {її, k)  | 
	соответствует  | 
	возможным  | 
	способам  | 
	образова  | 
||||||||
ния комбинации из к единиц и п—к нулей.
Случайные блуждания, интерпретирующие последо вательную процедуру, представим на плоскости, где по оси ординат откладываем шаги испытаний, а по оси
185
абсцисс — число нулей (рис. 2.38). Каждая точка на та кой дискретной плоскости представляет собой состояние цепи Маркова, и каждому состоянию ставится в соответ ствие вероятность попадания в это состояние, зависящая только от вероятностен в двух точках, находящихся на предыдущей строке — непосредственно под определяемой и слева от нее. Если цепь имеет поглощающие состояния,
то они выделяются на ри сунке значками — кру жочками н квадратиками. В дальнейшем для про стоты опускаем обозначе ния состояний и записы ваем только соответствую щие им вероятности. По глощающие состояния со ответствуют выходу за пороги • последовательной процедуры или достиже нию их.
Плоскость блужданий делится на три области: принятия решения о на личии сигнала, отклоне
О1 2 3 * ния решения и продолже
ния испытаний. Эти обла сти можно найти, исполь зуя метод проб и ошибок, причем необходимо, что бы при заданных вероят ностях F и D и отношении
сигнал/шум <72 среднее число выборок многоканальной процедуры при отсутствии сигнала HN было •минимальным. Считаем, что решения принимаются отдельно в каждом канале и что сигнал от цели присутствует в данном ка нале все время анализа пли не присутствует совсем. При поиске оптимальных областей конкретная структура обнаружителя не задается, 'Поскольку в результате ана лиза можно получить либо усеченную процедуру, либо процедуру с постоянными или переменными порогами, либо весовой вариант испытаний. Однако при использо вании неусеченного последовательного анализа стремле ние к.уменьшению EN ведет к увеличению среднего числа выбцрок при нал'ичи сигнала HSN, которое может стать
186
таким, что условие наличия цели в данном канале за время анализа не будет выполняться. При этом необхо димо учитывать вероятности перехода сигнала из канала в канал, что существенно усложняет как анализ, так и реализацию последовательного обнаружителя.
Рассчитаем вероятности F и D для примера, пред ставляющего усеченную на пятом шаге последовательную процедуру (рис. 2.38). Вероятность ложного обнаруже ния в элементе разрешения равна
Отсюда по заданной величине F можно найти вероят ность превышения порога квантования м0 шумовым вы бросом pN. По известной величине «о и вероятности PN из таблиц или графиков {98] находится вероятность psx [см. (2.47), (2.42), (2.43)]. Вероятности правильного об наружения и пропуска цели определяются как
D = РІК + 3P3SN  | 
	4SN +  | 
	<?SN '  | 
	
  | 
M = \ - D = q*SN+2pSN  | 
	q3SN+5plN  | 
	g3Sf/ •  | 
	(2.198)  | 
Вероятности пересечения или достижения порогов при отсутствии сигнала от цели за і шагов рг определяются следующим образом:
pl  | 
	=  | 
	0,p2 = irN,pl  | 
	= pa +  | 
	p3lt.  | 
	
  | 
р, =  | 
	р3  | 
	+ Зр3, qN +  | 
	2pNq3N ,P i =  | 
	1.  | 
	(2.199)  | 
Вероятности пересечения порогов за і шагов при отсут ствии сигнала в М элементах разрешения рм, І вычис ляются по формуле
Рм.і=Р"-  | 
	i = l - 2 5 -  | 
	(2-20°)  | 
Вероятность пересечения порогов на і-м шаге
•*М.І = РМ,І-РМ.І-І-  | 
	С 2 - 2 0 1 )  | 
Среднее число испытаний для шумовых элементов раз решения .
.(2.202)
І = І
187
Алгоритм работы устройства состоит• в накоплении первых трех единиц на пяти позициях, при этом выдается сигнал обнаружения; сброс осуществляется по появле нию на первых двух позициях двух нулей подряд пли по накоплению трех нулей па пяти позициях.
2.6.3. Синтез цифровых обнаружителей и анализ их эффективности при зависимых отсчетах [164—166]. За висимость между отсчетами квантуемых процессов соот ветствует ситуации, когда период временной дискретиза ции меньше времени корреляции процессов: Г д < т л . При чины существования соотношения Tn<Xh таковы.
1.Период дискретизации Гд выбирается не из соот ношения Гд^І/2/max, а из других соображений, например из условия однозначной дальнометрии в импульсной ра диолокации: Гд = Г п = '2 /?щах/С.
2.Часто с целью уменьшения ошибок восстановления днскретизированной информации (см. п. 4.2.2) период Г д
приходится  | 
	брать в несколько раз  | 
	меньше, чем l/2/mas-  | 
|
3. Ширина спектра помехи AfP  | 
	может  | 
	быть неизвест  | 
|
на и Ді/р max</max, а период дискретизации  | 
	выбран из со  | 
||
отношения  | 
	Г д = 1 / 2 / ш а х .  | 
	
  | 
	
  | 
4. При  | 
	проверке двух гипотез,  | 
	соответствующих рас  | 
|
пределениям с неперекрывающимися спектрами, выберем Г д = l/2/max, где /max—максимальная частота в спектре процесса, расположенного в области более высоких ча стот. Тогда для процесса, спектр которого находится в области более низких частот, отсчеты будут статисти чески зависимыми.
В соответствии с п. 1.7.1 последовательность дискрет
ных  | 
	переменных Гц,  | 
	1 = 1 , 2, ... , принимающих значения  | 
||
ха,  | 
	<х = 1, г, можно  | 
	аппроксимировать регулярными це  | 
||
пями Маркова  | 
	с г состояниями и связностью, определяе  | 
|||
мой  | 
	из (1.81)  | 
	или (1.83). В общем случае связность це  | 
||
пи Маркова будет различной для обеих гипотез Я 0  | 
	и # i .  | 
|||
Ясно, что при  | 
	синтезе оптимального обнаружителя  | 
	сле  | 
||
дует учитывать наибольшую связность vm ax, полагая, что обе гипотезы соответствуют распределениям со связ ностью Vmax-
Для выявления общих закономерностей и удобства сопоставления с обнаружителем при независимых отсче тах (2.189) рассмотрим синтез обнаружителя, когда по
следовательности  | 
	дискретных переменных, соответствую  | 
||||
щие  | 
	гипотезам Н0  | 
	и Ни  | 
	представляют  | 
	собой  | 
	односвяз-  | 
ные  | 
	(v = l ) регулярные  | 
	неоднородные  | 
	цепи  | 
	Маркова  | 
|
188  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
с г состояниями. Коэффициент правдоподобия для рас пределении типа (1.10) имеет вид
II—1 г
п
Л ( и ) = I '
/г=1  | 
	В=1  | 
	(2.203)  | 
|
И — 1  | 
	Г  | 
||
  | 
/»(2|8,)  | 
	п  | 
	
  | 
  | 
	/г=1 р=1  | 
|
  | 
	о = 1  | 
В результате логарифмирования (2.203) получаем алго ритм оптимальной обработки
In Л (и) = J] ^  | 
	(0) + 5 '  | 
	X с/ар (/г)  | 
	(/г), (2.204)  | 
а = 1  | 
	ft=l  | 
	а,В  | 
	
  | 
где  | 
	
  | 
	
  | 
	
  | 
WJ0)=^(P0aSNIP0«N)-  | 
	К?  | 
	{k)F=b(pa%SN(k)/Pe9N{k)).  | 
|
Ha рис. 2.39 представлена структурная схема алго ритма (2.204). Ее отличие от схемы рис. 2.36, соответст вующей алгоритму (2.189), состоит в том, что т:перь весовые коэффициенты Wa? (k) зависят от значений ха и
в два смежных момента времени. Исключение состав
ляет вес  | 
	Wa[0), зависящий только  | 
	от значения  | 
	ха. Соот-  | 
|
u(t)  | 
	й(к)  | 
	
  | 
	
  | 
	
  | 
  | 
	n-1  | 
	n-2 P Z ' C  | 
	1  | 
	О  | 
  | 
	Сдвиг  | 
	; . .  | 
	. і ' \  | 
	ШО)  | 
  | 
	
  | 
	
  | 
	_J  | 
	
  | 
  | 
	
  | 
	тт  | 
	т т  | 
	
  | 
  | 
	
  | 
	z  | 
	
  | 
	ЦК  | 
Рис. 2.39. Структурная схема алгоритма (2.204).
ветственно ПЗУ, хранящие весовые коэффициенты, долж ны иметь 2т-разрядный регистр адреса и хранить в об щем случае 2 2 т различных весовых коэффициентов. В остальном обе схемы аналогичны.
189
