та отенонИаны оьн
МЕ 
НаиНоароНе оЕ
оке 
"
ме 
в5 
вех. 
‘оивео
о. 
200-00 — 
9.10000
. 
361:09 —| 019735
2 
149.30 —| 0:2465
с 
450.90 — | 
0:22545
4 
35005 
| | 0:225И
5 
35095 
| 
о22548
* Использоваи. метод переменных состояния.
нии остается постоянной. Из результатов сравнения 
(табл. 4)
вндно, что за 2 суток, почти достигаются _ равновесные значения
при ТТ 
== 2000 КЕ5!. 
= 022548 
и СРНВЕР = 
168,55. 
При
"Т\Т 
= 10000 эти 
значения ‘соответственно составляют 0.0312.
м 387.92,
В настоящее время разработано несколько 
нтерационных мето-
дов, однако здесь мы рассхотрим только два из 
них В методиче-
ской дитературе можно найти более детальное описание методов
‘численного анализа. Общнй момент для всех итерационных мето-
‚дов состоит в том, что требуется найти такое значение нензвест-
ной переменной (в данном случае ВЁЗ!.), при которой сопостав
ляемые скорости равны между собой. При этом_ неявно допуска:
тся, что шаг ннтегрирования является достаточно бодьшим, чтобы
устамовилось равновесне. В этом случае скорость образования
ассимилятов СРНЕЕР должна быть равна скорости 
их потребде-
ния. СОК--МА!МТ. В графической форме решение этой. задачи
своднтея к нахождению точки 
пересечения двух кривых 
(см.
рис. 18).
'В С5МР для решения этой задачи предусмотрен иеявный цикл.
1МР..Этот метод основан на повторяющихся подстановках Зна-
чений переменной. Первое эначение 
ВЕ5!. задается приближенно,
рассчитывается — соответствующая — интенсивность _ фотосинтеза,
затем определяется_ ур 
‘резервов, при. котором общая ско:
рость потребления ассимилятов равна 
ранее рассчитанной. нитен-
сивности фотосннтеза; эта. процедура повторяется до тех пор, пока
ме будет достигнута достаточная сходимость. Первый недостаток
этого метода состонт в том, что необходимо обращать связь между
скоростью н нензвестной 
переменной: ВЕЗ!. должна быть напи-
сана или как функция общей скорости. потребления. ассимилятов,
ли как функция 
фотоснитеза. Если в программе нспользованы
функции АЕСЕМ, 
такое обрашение становится невозможным. Сле-
довательно, 
уровень резервов ВЕЗ!. нельзя записать в виде функ
ции фотосинтеза: 
эту переменную следуст аписать в программе
в виде функиии 
СОК, придав ВЕЗ!. явный вн в гиперболическох
86 
:
`уравнении. В результате неявный цикл 1МРГ. появляется в сле-
‘дующей имитационной программе:
"ПТОЕ 
ВЕЗЕВУЕ ГЕУЕ!. МИТНАМРЫСТТ ГООР.
1NoСОМ 
ВЕб1 — 023
МАМТ 
= 0,015*Т\УТ
СРНОТ 
= СРН5Т* (1. — ЕХР (—0.
ТАГ 
== АМИМ! (М5Н/500., 5)
м5Н 
= ОЛУТМТ
ВЕ5! — 1МРЬ(ВЕЗ, 0.0001, КЕЗ)
ВЕРЕ == АЕСЕМ (ВЕОЕТ, NoЕ5!.)
СРНЕЕР 
= ОРНОТКЕРЕ
СОВ == СРНВЕР 
— МАТМТ
ВЕб!| == СОК"КЕЕЗ51/(0.1°Т\Т 
— СОК)
"А0)
неявный цикл.
РАКАМ ОРНУТ == 400.
РАКАМ КВЕ5!. = 0.1, Т\УТ 
= (2000., 10 000.)
ЕОМСТОМ ВЕОЕТ = 0. 
1, 0.2, ., 0.25, 0., 0.5, 0.
МЕВ РИМТИМ == 
1, РЕОЕ! 
= 1
МЕТНОР ВЕСТ
РЕ1МТ ВЕЗ!,, СРНВЕР
| ЕМр
5тоР
‚ ЕМОЛОВ
этой программе отсутствует интеграл КЕ, вместо которого рас-
итывается КЕЗ!, в результате неявного цикла, первая строка
го имест вид
‚ ВЕ5!. 
— 1МРА(КЕЗЫ, 0.0001, КЕЗ1Л),
ВЕ1 
— начальная 
_оценка; 
`0.0001 
— критернй _ сходимости;
Ы. —- имя переменной, которая замыкает неявный цикл. Этот
'ыкающий 
оператор 
фактически 
является 
гиперболическим:
‘внением, описывающим. зависнмость СОК от ВЕЗГ,, в котором.
ой форме представлен уровень 
резервов КЕЗГ.. Три оператора
н цикла заменяют формулу для переменной состояния.
“