книги / Надежность и диагностика энергетических электромашин
..pdfРис. I . Области значений моментов Л/асин хронных двигателей при изменении диамет ра J) , допустимого перегрева IT и ин тенсивности охлаждения.
тов или |
о твета, |
что в заданном габарите нужный технический пока |
затель |
получить |
н ел ьзя . |
В работе рассматривается один из методов проектирования тя гового асинхронного двигателя с максимально возможным моментом при заданных габаритных ограничениях с использованием непосред ственного жидкостного охлаждения активных элементов.
Для определенного типа электрической машины в интересующем достаточно широком, диапазоне геометрических параметров, перегревов активных частей и интенсивностей охлаждения численным эксперимен том на ЭВМ по разработанному алгоритму расчета, описанному в ли тературе94, было получено поле реализуемых электрических машин. Анализ таких расчетов показал, что зависимость момента, развивае мого электрической машиной, от габаритных размеров при постоянной линейной скорости бочки ротора представляет собой в логарифмичес кой системе координат прямые линии для различных интенсивностей охлаждения и перегревов. На рис.1 представлена зависимость момен та М на валу асинхронного двигателя от его наружного диаметра J)
для шести значений интенсивности охлаждения, охватывающих’ следую-
хАнпилогов Н .Г ., Бандурин В ;В ., Вакуленко К .Н ., Вишникин А.И. Алгоритм проектирования частотно-управляемых короткозамкнутых асинхронных двигателей . - Техн.электродинамика,1981, А 4, с . 49-53.
171
щвй диапазон коэффициентов теплоотдачи на основной поверхности теп
лосъема: 1 - 250 Вт/м2 °С; 2 - |
500; |
3 |
- 750; 2 - |
1000; |
5 |
- |
1500; |
||
6 |
- 2000 Вт/м2 °С. Такой диапазон |
соответствует |
всем |
возможным зна |
|||||
чениям коэффициентов теплоотдачи в системах струйного масляного |
|||||||||
охлаждения электрических машин. Зоны, |
обозначенные на рис.1 номе |
||||||||
рами системы охлаждения, ограничены снизу линией перегрева, |
равно-, |
||||||||
го |
JT = 80°С, а |
оверху - перегревом |
200°С. Зоны I и 2-й |
||||||
систем охлаждения |
ограничены на |
р и с .I |
сплошными линиями, зоны 3 и |
||||||
4 - |
штриховыми линиями, зоны 5 |
и 6 - |
штрихпунктирными линиями. Для |
||||||
электромашин транспортного исполнения |
принятый диапазон по |
перегре |
|||||||
ву охватывает все |
классы применяемой изоляции. В исходных |
данных |
|||||||
задается температура охладителя на |
входе, следовательно, |
по |
вели |
чине максимально допустимой температуры принимается для выбираемо
го варианта величина перегрева. |
|
|
|
|
|
||
|
Для каждой системы охлаждения расчеты выполнялись |
для |
соотно- |
||||
’ шения наружного диаметра электрической машины к |
ее |
длине |
I |
.рав |
|||
ного |
л/1- = |
1 при следующих габаритах: J) - 0 ,3 |
м; |
0 ,4 |
м; |
0,5 3 м; |
|
0,6 |
м; 0,7 |
м и при перегревах ^7 V 8 0 °C ; 120°С; |
160°С; |
200°С. |
|||
|
Выбранные диапазоны указанных параметров в |
приводимом числен |
ном эксперименте являются достаточными для проектирования целого ряда электрических машин такого типа. Графическое представление расчетов наглядно иллюстрирует тенденции изменения момента на ва лу асинхронного двигателя при различных соотношениях исходных па раметров. Поскольку области на р и с.1 , соответствующие каждой из систем охлаждения, перекрываются, то следует, что получение в равном габарите одного и того же момента возможно при различных системах охлаждения■, создавая различные перегревы, или получение равного момента при одной и той же системе охлаждения возможно в разных габаритах, , создавая также различные перегревы .
Если задана температура на входе и максимально допустимая температура изоляция, то выбор вариантов электрической машины це лесообразно проводить по графику, показанному на р и с .2, где для перегрева, равного лТ= 120°С, приведены зависимости момента от
наружного диаметра. Аналогичные зависимости могут быть построены для различных перегревов по данным р и с .1 .
При проектировании новых транспортных электрических машин в заданных габаритах с произвольным выбором типа электротехничес ких материалов и вида системы охлаждения возможны ситуации, когда ограничителями выходного момента являются или нагрев, т . е . недо-
172
статочно эффективная система ох |
МЮ\Н-м |
|
лаждения, или индукция насыщения, |
|
|
стали, т . е . невозможно обеспечить |
|
|
необходимые |
размеры магнитопрово- |
|
да. Очевидно, что проектирование |
|
|
таких электромашин при габаритных |
|
|
ограничениях |
сводится к нахождению |
|
оптимального |
соотношения между |
|
электромагнитными нагрузками я теп |
|
|
ловыми проводимостями, приводящи |
|
|
ми к максимальному моменту на в а |
|
|
лу. Предлагаемый ход проектирова |
|
|
ния с выбором по графику оптималь |
|
|
ного соотношения между габарита |
|
|
ми, энергетическими и тепловыми |
|
|
параметрами |
позволяет оператив |
|
но решать задачу выбора вариан та электрической машины для транс 0,2
портного средства, где компоно
вочные решения могут изменяться
0J
в зависимости от ее выходных па раметров .
На основании полученных зави симостей по задаваемым величине момента и габаритным размерам вы бираются параметры системы охлаж
дения, а затем варьируя этими параметрами,, а также геометрически ми соотношениями и электромагнитными нагрузками добиваются необ ходимого момента в данном габарите. Таким обраэом, без проведения поисковых предварительных расчетов выясняется вопрос возможно ли удовлетворить требованиям предложенному техническому заданию в диапазоне реализуемых основных параметров электрической машины.
Предлагаемый подход к проектированию применим к различным типам электромашин с форсированным охлаждением, при этом умень шаются сроки доводки головных образцов и повышается их эксплуата ционная надежность.
173
Ш: 621.3.012.6.001.24 : 681.3
Г .Л.Баранов
МОДЕЛИРОВАНИЕ ОБЪЕКТОВ С ПЕРЕДАТОЧНОЙ ФУНКЦИЕЙ ' ВТОРОГО ПОРЯДКА
Динамика объектов различной природы описывается дифференци альным уравнением или передаточной функцией. Известно, что переда
точные функции высокого порядка представимы в воде произведения
передаточных функций первого |
и второго |
порядка, |
моделирование |
ко |
||||||||||||||
торых на ЭРМ выполняется с помощью рекуррентных соотношений. |
|
|||||||||||||||||
|
При моделировании непрерывных объектов, |
у |
которых имеется |
пе |
||||||||||||||
редаточная функция второго |
порадка'возможны три случая. |
|
|
|
||||||||||||||
|
В первом случае корни характеристического |
уравнения комплекс |
||||||||||||||||
но-сопряженные |
и передаточная функция объекта имеет вид |
|
|
|
||||||||||||||
|
|
|
|
W,<P) |
Ггр + ГР + (Г+б)г Г2) |
|
|
|
|
|
(1) |
|||||||
|
Во втором случае корни мнимые и сопряженные, |
тогда |
модель |
|
||||||||||||||
объекта ш еет |
вид |
|
|
|
к(Г,р+ f) |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(2) |
|||||
|
|
|
|
|
|
|
Гр + I |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
В третьем случае обычно корни разные действительные |
числа |
и |
|||||||||||||||
передаточная функция объекта эквивалентна передаточной функции |
|
|||||||||||||||||
двух последовательно |
соединенных звеньев первого порядка |
|
|
|
||||||||||||||
|
|
|
|
w , |
\ |
|
« |
ггР + 0 |
|
|
' |
|
|
|
|
О ) |
||
|
|
|
|
|
|
(Гр*1)(Тл р+1) |
|
|
|
|
||||||||
|
В качестве дискретного аналога непрерывной модели объекта в |
|||||||||||||||||
этих случаях целесообразно |
применять рекуррентную |
формулу |
|
|
||||||||||||||
где |
А, В, |
С, J> - постоянные |
eYi-2i’CXi "-DXM ’ |
|
|
|
от |
|
(4) |
|||||||||
коэффициенты, зависящие |
физичес |
|||||||||||||||||
ких параметров, передаточной функции объекта |
и от |
шага |
А , |
опре |
||||||||||||||
деляющего интервал времени между рассматриваемыми |
(7 - /; |
- |
/ - |
|
||||||||||||||
и |
/ -м - |
( / + |
I)-м |
состояниями |
объекта |
в |
переходном |
процессе; |
|
|||||||||
/ / , |
- |
входные |
8начения |
воздействия |
|
соответственно |
в |
|
|
|||||||||
и ( / - 1)-й моменты времени; |
Yh |
|
Y;_,, |
' |
Y}_2 |
|
выходные |
значения |
||||||||||
реркции объекта на входное воздействие |
соответственно |
в |
/', |
U-D, |
||||||||||||||
( i - 2)-й моменты времени. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Дискретная модель вида (4) полностью определяется численны |
|||||||||||||||||
ми значениями коэффициентов А, |
3, |
С, D , |
которые |
зави сят |
от |
пере- |
||||||||||||
174 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
даточной функции объекта моделирования и от споооба описания вход
ного |
воздействия. |
||
|
|
В зависимости от целей моделирования численные значения А . |
|
S, |
С, |
V |
меняются.. |
|
|
Определим численные значения постоянных коэффициентов А , S, |
|
с |
, J) , |
необходимых для дискретного моделирования непрерывного |
объекта, для чего выполним расчет переходного процесса по.методу
предварительного |
интегрирования, примененному в работе / I / |
для |
|||||
объектов первого |
порядка. |
|
|
|
|
||
|
В качестве примера дадим аналитический вывод для объекта с |
||||||
передаточной |
функцией |
|
|
|
|
||
|
|
|
. ___________ * |
|
Y(P) |
|
|
|
|
|
Шр)~ Тгр+2аГр +г |
Х(р) ' |
(5) |
||
где |
р = — |
- оператор дифференцирования. |
|
||||
|
В этом случае дифференциальное уравнение переходного процес |
||||||
са |
на интервале |
/ < t* t. |
имеет |
вид |
|
|
|
|
|
|
|
-Н |
УСв |
= к х С И , |
( 6) . |
|
|
|
|
|
|
|
где x(f) - модель возмущающего воздействия на входе объекта, Y(t> -
реакция объекта на это воздействие.
Если сделать замену переменных г - 1 - t-j f • на расчетном интервале принять линейную модель входного воздействия
|
|
*(*)- *н + |
|
г> |
|
(7) |
||
перейти |
к изображениям по Лапласу /2 7 |
с оператором, то |
получим |
|||||
операторное уравнение |
|
|
|
|
|
|
||
где х |
Y а. |
у \ =——— ■— |
начальные условия или исходные |
|||||
'■> |
О-0) |
/-V (+0) |
п |
|
|
|
|
|
значения соответствующих |
переменных, |
когда |
г |
стремится |
к нулю |
|||
справа. |
|
|
|
объекта y(S) |
|
|
||
Искомое |
изображение |
реакции |
из уравнения (8) |
|||||
найдем |
в ввде |
у (3 )- ^ |
г* |
|
, |
+ |
|
|
+2afS+ , |
- Y(tg) |
|
|
|||||
|
r*S+ 2 <xT |
{ M S * * j+ H t/s -//4 s s) x ,lf |
• . |
|||||
+ |
T2S2 *2аУ5*/ |
r zs e -*■ PetГs +/ |
|
’ |
|
175
Используя соотношение |
j |
f ( s) +-*5 |
f ( t ) d r , |
|
осуществим |
об |
||||||
ратное |
преобразование Лапласа |
/ 2 ,3 / , тогда |
получим требуемые |
ча |
||||||||
стные |
оригиналы |
|
|
|
|
|
|
|
|
|
|
|
|
• S+ZaTS *1 |
|
|
|
|
|
|
|
|
|
|
|
|
r*s+ 2а Г |
|
|
|
|
|
|
|
|
|
|
|
|
|
= в |
Т |
(cOS H it* |
S in |
< v t ) , |
|
|
|
|
|
|
|
r*s**taTS+1 |
|
|
|
|
|
|
|
|
|
|
|
|
s ( r * s 2 + M S + 1) = /- s |
r |
. (a> soir*~ r |
s i n * |
j r ) , |
|
|
|
||||
$2(T*s*i-2ars*f) - = (t-2 a Г)* e r ^ZaTcostvr* ( |
-oiT2^ |
sin d/rj |
|
|||||||||
Здесь |
принято, |
что &= .£ |
/- а 2 . |
|
|
|
|
|
|
|
||
Из операторного уравнении (9) на основе полученных оригиналов |
||||||||||||
после соответствующих преобразований получим рекуррентное |
соотно |
|||||||||||
шение вида (4 ), |
в котором для рассмотренного примера с |
(5) |
при |
|||||||||
расчетном шаге, |
равном |
h «постоянные |
коэффициенты запишутся в |
|||||||||
ввде |
|
|
|
|
|
|
|
|
ah |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Л" [ Ж ' в , п + |
|
|
sin(ah* arcty а>г\е -~ г |
|
|
|
|||||
|
|
ah |
|
|
|
|
|
|
|
|
|
|
|
|
sin ah, |
|
|
|
|
|
|
~T~ | |
|
(Ю) |
|
C ^ ~ ^ ( h - 2 a l)^ \(^ ~ a>Tz '}swti)h+ 2ctrcostnh j e |
|
|
||||||||||
* |
-[(Mr+w |
|
|
|
|
]' |
|
|
||||
Полученные коэффициенты позволяют |
при любом |
//, |
когда вход |
|||||||||
ное воздействие |
на расчетном интервале |
изменяется |
согласно |
(7 ), |
рассчитывать на основе (4) точные значения реакции объекта в диск
ретные моменты времени. Погрешность расчёта переходного процесса
в динамическом объекте появляется только при отклонении закона
изменения входного воздействия на расчетном интервале |
h |
от его |
линейной апроксимации. Таким образом, точность р асч ета |
по данно |
|
му методу зависит не от параметров передаточной функции |
объекта, |
|
а 0Т точности линейной апроксимации входного воздействия |
на рас |
|
четном шаге.* |
|
|
176
. Расчетные модели интегродифференциальных объектов второго порядка
'и/(р ) |
Постоянные |
коэффициенты. А,В, С, |
вычислителя |
для |
данных |
W(p) |
передаточных |
|||||||||||
|
|
|
|
|
|
|
|
|
функций |
|
|
|
|
|
|
|
||
|
|
д= IwTsina^ |
|
|
|
|
s)r>- |
|
arctff ®Г) \ е~АГ* |
В = |
с Л/Тsirr |
|
|
|||||
Щ р *п |
с- ж |
\ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Грг+2Тр+ (7+ДгГ2) |
|
|
|
|
|
|
|
|
|
|
-р/р . |
. . |
2Т-(Гг-Ь)0**?Г*) |
|||||
|
|
„ кТг \ уЛИ-ЛГ+ГГг -/*-йГ,*гГ*)г* 4 ‘Г*(АГ-*»Т7-ЪГ)я' |
||||||||||||||||
|
|
M~hT} |
[ |
|
|
|
r-fa(t+o*r*) |
|
|
е |
*П(«*+/)+ |
(t+a*r*)* |
|
|||||
|
|
|
/+ й>Г2 |
|
] ' |
* |
~ arct9 "Щ-АГЧ- тт, - г * - ц л *г* |
-*а»*9(-*т) |
|
|
||||||||
НтгР+1) |
4= [ Г/ЛSin”,/Г + COS*/Г ] ; |
|
С**11 f a *А)-Ггсо54 - /*,* £ ] ; |
|
|
|
||||||||||||
Т2р2+1 |
B = -T s,nT ’ |
|
|
|
|
|
Г,'A) MS |
|
|
Гу |
in Л /г] |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
ИГтРЧ ) |
. t.\h(T*h) „-Щ |
П Ы ) П-Ш \ |
г. ягз |
Гш - Ь ) - л /% |
|
Г(Г,-Г) - л /г ................,7 ■ |
|
|||||||||||
А~Ыт3- п е |
|
|
щ - п * |
\’ |
|
7 ^ Г |
е |
~сг-г3) е +<ь+П-г-Щ |
|
|||||||||
(Гр+/)($р*Г) |
п |
[ Ъ |
-Щ _ _ £ / _ |
А /п . |
w t o - m |
w |
ш _ Ог&0>+Г3) -Щ |
+<**+»! |
1 |
|||||||||
|
|
в |
[h(r5- r ) e |
|
Щ Т )* |
J ’ 3 |
щ |
1 сггЬ). |
|
|
c r -t) |
|
/ ; j |
|||||
к(.Гьр+1)(Т2р*Г) |
* - |
/ / > [ '♦ > - г * - * " |
] ; |
|
М |
|
|
( |
|
|
|
? |
|
' |
||||
р(Гр*П |
' ~ f ( » * " ) > |
|
' ^ [ ( т |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
A ( fr + n |
. |
А |
|
|
3--Г, |
|
С -к |
(А г/зт - - f f |
)■ |
|
J)*Hft2/ e r + |
) |
|
|||||
Гр* .. |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
— Ш рТГ)--------- |
A - e h/ r |
; |
В*0-, |
|
c^B -\-i-A (rr T)*h\-t |
2 =2 -[(т-А)(ГГ Г)+РА] |
|
|||||||||||
Гр |
f f |
■ ~ ~ 7 Г п |
|
т у - |
|
с , к {? & * -)■ , |
|
|
j > - k ( w / ; - y |
“ |
||||||||
* J 4 |
L - |
|
|
|
|
Выполнив аналогичный расчет и преобразования для всех корней характеристического уравнения объекта моделирования, получим соот ветствующие значения постоянных коэффициентов А, В, С, В , расчетные
выражения для которых сведены в таблицу. Для сравнения в таблице кроме других пяти типов моделей второго порядка приведены две мо дели первого порядка, для которых так же справедливо рекуррентное
соотношение -(4 ). |
|
А. В, |
С, J), |
Полученные расчетные выражения для |
коэффициентов |
||
позволяют, зная аналоговые параметры непрерывных объектов, |
такие |
||
как к - коэффициент усиления и К |
постоянные |
времени, пе |
реходить к точным дискретным моделям вида (4 ) . При изменении физи ческих параметров объекта следует, согласно таблице, пересчитать соответствующие параметры дискретной модели. Расчет коэффициентов /ft в, С7 в осуществляется на этапе настройки модели на задачу исследования.
Расчет по рекуррентному соотношению (4) справедлив для многих конкретных объектов первого и второго порядка, что позволяет реали
зовать эффективный интегродифференциальный вычислитель с адаптив ной параметрической настройкой.
Разработанная модель и метод расчета коэффициентов рекуррент
ного вычислителя нашли эффективное применение в программах струк турного моделирования сложных объектов непосредственно по ' струк турной схеме модели Д 7 , .в специализированных локальных средствах технической диагностики, в цифровых регуляторах.
• I. Баранов Г.Л. Раочет на ЦВМ нелинейных динамических систем непосредственно по структурной охеме модели. - Изв.АН СССР. Энер
гетика и транспорт, » 6, 1969, с . 71-79. |
|
преобраэова- |
|
|
2. Деч Г. Руководство п р акти ческо м у применению |
||
|
3 . Мартыненко В.С. Операционное исчисление. - Киев : Высш. |
||
школа, 1973. - 267 с . |
|
|
|
УДК 621.313.321 |
|
|
|
Гв.П.Коваленко ,| Н.А.Куличихина |
|
|
|
|
СПЛАЙН-МЕТОДЫ В ЗАДАЧАХ ВДЕИТИФИКАЦИИ |
|
|
|
ЭЛЕКТРИЧЕСКИХ МАШИН• |
|
|
• |
Рассматривается задача определения параметров дифференциаль |
||
ного |
оператора Н(р) взаимосвязи между векторами |
входного воздействия |
|
Х(р) |
и реакции объекта Y(p) , Предполагается, |
что для |
искомой функ- |
178
ция Н(р) известна аналитическая форма и поэтому задача сводится к* определению величины коэффициентов ограниченного множества функ ции Н(р) однотипной структуры, но различной сложности (например, дробно-рациональной функции п ервого,'второго и третьего порядков).
Для решения задачи идентификации предложено много различных методов / 2 , 3 , б7, причем их большое разнообразие объясняется некор ректностью задачи, как следствие строгого несоответствия аналити ческих свойств функции ff(p)€ С " и экспериментально определяемых .
динамических характеристик, |
для которых, |
в лучшем случае, справед |
||||
ливо |
утверждение |
об их |
принадлежности к |
классу непрерывных функг |
||
ций, |
т . е . /КО, |
Н ( ъ ) € С ° . |
|
|
||
|
Задача |
определения |
т |
неизвестных коэффициентов функции ^ |
||
сводится к |
решению матричного уравнения |
т -го порядка |
л* -в.
где Л - определитель Гильберта, в процессе.решения которого полу чаются оценки первых т разделенных разностей (производных) экс периментально измеренной динамической характеристики.
Учитывая плохую обусловленность определителя Гильберта / 7 / , попытка определить из (I ) более, пяти-шести неизвестных коэффици ентов операторной функции Н(р) нереалистична. Удовлетворительные
результаты дает практически любой из известных методов идентификации, если определяется не более двух-трех коэффициентов. '
Для получения надежных оценок неизвестных пяти-шести коэффи-: циентов операторной функции Н(р) необходимо не только обеспечить максимально достижимую точность измерения динамических характерис тик, по и сделать все возможное для обеспечения удовлетворитель ной однородности исходной информации и ее качественной фильтрации.
Подходящим для |
этой цели аппаратом являются интерполирующие |
|||
■и аппроксимирующие |
кубические |
сплайн-функции / 1 , 4 / . |
Опыт работы |
|
со сплайн-функциями |
показал, |
что условие |
/ Б / |
|
|
P j'-M /fy/f: Mf/H |
<5, |
(2). |
где М-, - величина второй производной интерполирующей сплайн-фувк- ции в /'-м узле сетки определения динамической характеристики „ яв ляется хорошим критерием для оценки однородности исходной инфар>- мации. Причем нарушение этого условия указывает либо на грубую ошибку измерения динамической характеристики, либо на ошибку перформа'ции исходных данных.
179
При измерении динамическая характеристика обычно задается
в виде простейшей, но наиболее распространенной регулярной (пери одической) дискретной последовательности. Однако известно доста
точно много форм дискретного представления экспериментального ма териала, выгодно отличающихся от периодической дискретизации мощ
ностью подавления шумов измерения. Среди них, по-ввдиыому, наилуч шими являются’ при анализе частотных последовательностей тригономет рические ряда Фурье, а для временных последовательностей - ряды Пуассона /5 7 , при вычислении которых реализуется близкая к опти мальной методика фильтрации исходных данных.
В работе / 3 / показано, что между коэффициентами спектров
Фурье или Пуассона и искомыми коэффициентами операторной функции
НСр) существует простая й однозначная |
взаим освязь. Так, |
в |
частно |
||||
сти, из теоремы о смещении изображения |
непосредственно |
следует,что |
|||||
коэффициенты спектра Пуассона |
отличаются только знаком |
от |
коэффи |
||||
циентов Тейлора - |
Маклорена разложения |
операторной функции |
Я(р) |
||||
в степенной ряд в |
точке |
/>** |
|
|
|
|
|
f K p ) * \ e (Р |
) |
= |
~§~( |
\ |
e * * t* q Ct)c(t = |
|
|
- 'L c -P) t Pk,
где Рк - коэффициенты спектра Пуассона.
Для вычисления коэффициентов спектров Фурье или Пуассона раз работано достаточно много'методов /3 ,7 7 , среди которых для полу чения их неомещеныых оценок целесообразно использовать метод (по идее аналогичной методу Филона), использующий кубическую сплайнфункцию для гладкого восполнения дискретной последовательности динамической характеристики объекта.
Вычисление моментов Пуассона сводится при этом к простой реккурентной процедуре:
■V * |
4 * * и * |
+ |
+ * /'* * > . |
|
. |
\ |
Ы н |
nr |
т-1 , |
Ят s J |
е |
dw, |
где Of, bj. Cj - коэффициенты интерполирующей сплайн-функции в /-м узле сетки определения Ш )
180