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

книги / Эксергетические расчеты технических систем

..pdf
Скачиваний:
11
Добавлен:
20.11.2023
Размер:
19.71 Mб
Скачать

нечного числа уравнений

- Ц - = 0 (i = 1,2,

k; / = 0, 1, . . . . t).

(9.78)

Решение системы уравнений (9.78) позволяет вычислить значения экстремалей в заданном числе точек и на их основании построить кусочно-линейные функции, аппроксимирующие экстремали функционала (9.72). При этом считается, что в пределах каждого промежутка времени экстремали остаются постоянными, а их изменения происходят скачкообразно на границах участков. Чем больше вы­ бранных участков, тем меньше погрешность полученных результатов, и поэтому при t оо задачу можно решить точно.

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

П = [Zu (е?',

©к11’. Д О

+

Z12 (е2, ек, ATJ +

Z13 (e f\

0<13\ AT™) +

+

Z21 (e3, 0 И> AT5) +

Z22(e?>, AT?') + Z31 (eQoxji,

0 B, ДГВ) +

 

+ 2

32( ^

(32k

 

+ ц ЭЛ V

(Еп (еУ, ©!/’, АТ1») Ат] +

 

AT?')]л ,

т г

 

 

 

 

 

 

 

 

/ ± 1

 

 

 

 

 

+

£

[Е13(еУ, ©if, АТш') Дт] +

£ [Егг( ^ , А?/') Дт] +

 

 

 

/=1

 

 

 

 

 

/=1

 

 

 

+

S

[^ з2 ( С

л . AT?)Ax]\ +

nwi

 

©к71, AT'w) Д т] +

 

 

/=1

 

 

 

 

)

/—1

 

 

 

 

 

 

+

S

 

 

©!/’, ДП71) -

е^1 Дт} +

 

 

 

 

 

/=1

 

 

 

 

 

 

 

+

Б

 

[Е3 (4 0ХЛ, ©в\

АПЛ) -

е?) Дт}

( / = 1 . 2 , . . . , о,

(9.79)

 

/=1

 

 

 

 

 

 

 

 

 

 

где — функциональный множитель Лагранжа;

надстрочный индекс / соот­

ветствует текущему значению номера участка, а индексы 11, 13, 22, 32 — режи­ му максимальной производительности соответствующего (по номеру) элемента оборудования холодильной установки.

Используя выражение функционала (9.79) и соотношения, описываемые зависимостями (9.76) и (9.78), можно получить следующую систему уравнений:

 

 

+

. dZ\\l)

dZ\f

i \ - щ п [цэл (£</‘ +

^

Н - ^ Г Г +

= 0;

 

сэе(11>

d0il3)

/ = 1

 

 

 

к

/

- (гТ е У + Ztf)

0;

 

V

 

/=1

дв?

 

 

 

/=21

dQ(j'

+ Z'i)

= 0;

 

 

 

 

 

 

 

 

 

z\$] +

 

/5 7 (H)

 

л 7(13)

 

[цЭЛ(£$? +

£$) + ц.И§ +

^

т

г

+ 1 ^

| г = 0;

(9.80)

d&lW

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

л 7(22)

 

 

 

^ 0 U » + r f W + « > + * £ r - f c

 

 

 

 

 

 

 

^7(32)

 

 

 

 

 

 

 

•j- (Цэл-^32

+ 4

Л^

+ 4 ? )

+

0 Д 7 2)

=

®

(

/ = 1

. 2

...............О»

 

d&T(J>

 

 

 

 

 

 

 

 

 

 

 

 

г д е

 

 

 

 

 

 

 

 

 

 

 

 

г? = ^ т г [Цэл (М? +

Е% +

nwv?2 + №

 

 

/57(11)

 

/57ОЗ)

 

+

 

- ^ 5 -

+

;

(9.81)

дер

 

 

 

 

 

 

 

 

 

 

 

 

ГФ

=

л

 

 

+

4?) +

 

<97^2’

 

 

ГЗ

(Цэл^М +

 

 

22

 

 

 

 

 

 

 

 

 

 

 

* Г

 

 

Система уравнений (9.80) с учетом выражений (9.81) представляет собой запи­ санное в общем виде решение задачи динамической оптимизации одноступенча­ той одноцелевой холодильной установки, позволяющее получить набор зна­ чений экстремалей, обеспечивающих с учетом сделанных допущений минимум функционала годовых приведенных затрат (9.79).

Решение системы уравнений (9.80) с учетм (9.81) дает наборы оптими­

зирующих

переменных

©Д

АTw\ ©Д ДГД ©Д

AT{J' ( / = 1 , 2,

t),

объемных

расходов

сред

у}#, иД 0Д V{b и

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

этим

объемным расходам мощностей электродвигателей насосов охлаждающей сре­ ды е(/з, промежуточного хладоносителя djl, электродвигателя вентилятора £32

и электродвигателя компрессора

£п, обеспечивающих

минимальное значение

функционала (9.79). Набор значений i//2\ vKJ\

И/’ (j

= 1, 2,

..., t) отражает

закон оптимального управления

рассматриваемой

холодильной

установкой с

учетом использованных допущений. В то же время значения V\$\ vf2\ г432), ИД служат для подбора соответствующего оборудования (насоса охлаждающей воды, насоса хладоносителя, вентилятора воздухоохладителя, компрессора) по их максимальной в течение рассматриваемого периода производительности. Поверхности теплопередающих аппаратов (конденсатора FK, испарителя Fuy воздухоохладителя FBO) согласно принципу их полного использования (9.75) могут определяться по оптимальному значению соответствующего среднелога­

рифмического перепада температур (©[/*, ©Д ©Д для любого рассматриваемого промежутка времени при условии, что тепловой поток и коэффициент тепло­ передачи определены для того же промежутка времени.

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

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

ных значений температур (t = 12), количество этих уравнений увеличивается до тридцати девяти.

Решение системы такого количества уравнений, которое для более сложных схем установок еще более возрастет, даже при наличии современной вычисли­ тельной техники затруднительно. Однако в Ленинградском технологическом институте холодильной промышленности на базе теории динамического про­ граммирования разработана методика решения, позволяющая свести решение /^-мерной задачи к решению tk одномерных задач вида у (х) = 0 [126].

Сущность метода сводится к решению задачи динамической оптимизации в две стадии. На первой стадии при некоторых значениях параметров системы Ft (9.75) за счет выбора оптимального управления ДТц находится локальный минимум; на второй стадии определяется оптимальное значение параметров Fi, обусловливающих оптимальные значения 0 t/, и осуществляется выбор оп­ тимального управления АТ1*/, обеспечивающего минимум рассматриваемого функционала (9.79).

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

На втором этапе на основе методов комбинаторного анализа подбирается реальное, выпускаемое промышленностью оборудование. Подбор осуществля­ ется на основе сопоставления характеристик реального оборудования с харак­ теристиками оборудования, полученными на первом, идеализированном, эта­ пе решения. Возможны три варианта подбора оборудования: полной заводской готовности (конденсаторы, испарители, воздухоохладители), комплектуемое из стандартных модулей (охлаждающие батареи) и с регулируемой производитель­ ностью (компрессоры, насосы, вентиляторы). Критерии для выбора наибо/iee рационального варианта, которыми могут служить, например, максимальная степень приближения по производительности к идеализированному решению первого этапа, минимальная масса оборудования или занимаемая площадь, задаются ЭВМ либо варианты выбираются проектировщиком по комплексу показателей в процессе диалога человек — машина.

На третьем этапе, используя основную программу динамической оптими­ зации, для выбранного оборудования решается первая стадия задачи — опре­ деляется закон оптимального управления подобранным оборудованием. Сопо­ ставление значений функционала, полеченных на третьем и первом этапах, по­ казывает степень удаления от оптимального решения, обусловленную реальной

номенклатурой основного оборудования.

Рассмотренная модель одноцелевой одноступенчатой холодильной установ­ ки принимается в качестве базовой. Моделирование более сложных систем осу­ ществляется путем некоторого дополнения базовой модели.

Моделирование и оптимизация двухступенчатых одноцелевых холодильных установок

При моделировании двухступенчатой холодильной уста­ новки, схема которой приведена на рис. 9.10, в зону / термоэкономической модели (рис. 9.11) дополнительно включается компрессор ступени низкого давления 14 и промежуточный сосуд 15 (см. рис. 9.10). Это приводит к появле­ нию дополнительных оптимизирующих переменных, в качестве которых могут быть выбраны разность температур на холодном конце змеевика промежуточ­ ного сосуда ДГп.с и среднелогарифмическая разность температур в промежу­ точном сосуде 0 П.С. При этом следует учитывать вводимые в первую зону моде­ ли суммарные отчисления от стоимости компрессора ступени низкого давления с электродвигателем z14 и от стоимости промежуточного сосуда z15, а также энергию, необходимую для привода компрессора ступени низкого давле­ ния еи .

При таких уточнениях программа, используемая для моделирования про­ цессов в зоне 1 одноступенчатой установки, позволяет моделировать процессы ступени высокого давления двухступенчатой холодильной установки. Для мо­ делирования процессов ступени низкого давления, на основе изложенных выше принципов, нужно составить аналогичную дополнительную программу.

Рис. 9.10. Функциональная схема двухступенчатой одноцелевой холодильной установки:

14 — компрессор ступени низкого давления с электродвигателем; 15 — промежуточный сосуд (оо. бальные обозначения — те же, ч^ч на рис. 9.7)

% Щ

Ч щ

Рис. 9.11. Термоэкономическая модель двухступенчатой одноцелевой холодильной установки:

®пс — температурный напор в промежуточном сосуде; Д7ПС — разность температур на холодном конце змеевика промежуточного сосуда (остальные обозначения — те же, что на рис. 9.8)

Функционал, определяющий минимум переменной части приведенных за­ трат, для одноцелевой двухступенчатой холодильной установки имеет вид

n = [Z11(<r , e

,K11\

АТ™,

в™, AT™) + Z12(e2, 0 К, ATw, 0П.С| ДГП.С) +

+ Z 13(e f\

eg®,

Д ^ 3’, ДГ'п13',

0п.3с) +

Z 14(e214\ в™, ДТ ™ , 0 <п14с\ ДП!4)) +

+ Zls(e2, 0 К, ДTw, 0 П.С, ДГп.е) +

Z21 (е3, 0И| ATS) + Z22(ef*,

AT™) +

 

 

+

231 (е<зохл,

0В, АТв) +

Z32 (е™хл, ДГ(32’)] тг +

 

 

 

 

+

ЦЭЛ( Х

 

 

 

©к7’, АТ\1\

0 ^ С) Д7(плс) Дт] +

 

 

 

 

 

 

 

1/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

£

 

 

е^,

А71?, ©&, ДГ|/?с) Ат] +

 

 

 

 

 

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

£

[^14 (ей\ 0кЛ,

Д71",

0(„лс, ДТ(Л ) Дт] +

 

 

 

 

 

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

£

[£я (^л, д П л)Дт] +

Е

[£зг ( 4 L ,

 

дГвл)Дт]) +

 

 

 

,=1

 

 

 

 

 

 

 

/=1

 

 

 

 

 

J

 

 

 

 

+

ц „ £

[V12(<^, ©к1, д ^ ',

ей»,

Д 7 ^)Д т] +

 

 

 

 

 

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

£

{ ^ [ £ 2(^л, 0«», A T f ) е{2 ]] Дт] +

 

 

 

 

 

 

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

+

£

{ГЙ» [ £ , ( < „ .

®(в ,

А П ") - 4

л1Ат}

( / = 1 , 2 ......... о.

(9.82)

i = [

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На основе функционала

(9.82)

составляется

следующая

система

урав­

нений:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

}

i l f

 

[цэл (£"

+ £ '1 3 + £ '14} + ц - к'1/2' +

z *2+

г Ц

+

 

 

 

 

 

 

+

t

dZ^

+

a zfl3)13

+

<

4>

 

= 0;

 

 

 

 

 

 

 

 

ае<П)

a i

1131

a e f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

[Цэл (£ ‘n +

£ 13

+

E и )

+

ЦWV {\2 +

2 12 +

Zi§]j -f-

 

 

 

 

 

 

 

/

a z]1,1l 1’

 

dz\\3'z?»

t

a zfl(I144>)

,

_

 

 

 

 

 

 

 

 

+

M

i g n r +

ae«j®

^

a e g 4’

 

— u,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

d

/

(/)!?(/) _i_

7(/)\

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

+

Z2l)

=

 

 

 

 

 

 

 

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

=

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/=1

[Цэл (Еп + Е[>] +

Е\<}) +

и У £ +

+ 2 Й 1 +

dzft1’

d Z \f

dZ\14'

 

++ дЬТ™ + адтй4» ~~ ’

 

 

[Дэл (£# +

Е(й +

Ей) + п У й + У

+ z$ ) +

 

 

 

dZ\\l)

dZfj'

 

dZ\14'

_ 0 .

 

 

 

 

д Д Г ^

дД7<п13)

^ Э Д Г ^

 

 

 

 

 

_

 

 

 

 

 

dZK ’

 

 

 

 

- £ ц г

(Дзл i t f

+ ^

+

ZSP)

+ -

^

=

0;

 

 

 

 

 

Q7(32)

 

 

 

 

 

d ^ jW (Цэл-^32 +

Гз^Ез' +

Z3i’) + д^

(з2)

®

(/ =

1.

2, . . . , / ) , (9.83)

г д е

Bjj

 

 

 

 

 

 

 

 

 

 

Г(2) =

[Дэл (■£'п) +

£ $ +

£ Й )

+ Ц .У &

+ Zi2

+

Z15] +

 

 

 

dZ*,1.11

dZ(,33)

 

dZ^4’

 

 

 

 

 

# - + ^ l r + - ^ r -

 

 

 

r f =

д

 

 

 

Л7*22*dz!21

 

 

 

 

( u „ £ g + / f

+ Z ^ ) +

^ - J l p

( / =

1 , 2 , . . . , / ) .

 

 

 

 

 

 

Ц

 

 

 

Выражения (9.80) и (9.81), относящиеся к зонам II и / / / модели, не пре­ терпевают никаких изменений по сравнению с предыдущим вариантом одноце­ левой одноступенчатой установки, что подтверждает универсальность полу­ чаемых решений.

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

Моделирование и оптимизация холодильных установок с оборотной системой водоснабжения

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

•скольку температура воды TWx, направляемой в конденсатор, зависит от сов­ местной работы градирни и холодильной установки, оптимальные размеры гра­ дирни и режим ее эксплуатации могут быть выбраны только при комплексном рассмотрении градирни и холодильной установки.

Не останавливаясь на моделировании процессов, протекающих в градирне [125J, рассмотрим методику включения градирни в задачу оптимизации холо­ дильной установки. Учитывая, что температура охлажденной воды TWl зависит и от глубины, на которой вода охлаждается АТШУявляющейся одной из оптими­ зирующих переменных первой зоны, водоохлаждающее устройство вклю­ чается в зону / базовой модели.

Схема одноцелевой одноступенчатой холодильной установки с оборотным водоснабжением показана на рис. 9.12, а термоэкономическая модель такой си­ стемы — на рис. 9.13 (аналогично решается и для двухступенчатой системы).

В зону / (см. рис. 9.12) дополнительно вводятся суммарные отчисления от стоимости градирни г1б (см. рис. 9.13), стоимости вентилятора с электродвига­ телем z17 (для вентиляторных градирен) и стоимость свежей воды цс.в vc.B, иду­ щей на компенсацию испарившейся циркулирующей воды и ее механического уноса. Кроме того, к зоне / должна быть подведена дополнительная энергия е17, необходимая для привода вентилятора (для вентиляторных градирен). В ка­ честве дополнительных оптимизирующих переменных выбираются температура охлажденной в градирне воды TWt и удельная гидравлическая нагрузка (плот­ ность орошения) gw.

Рис. 9.12. Функциональная схема одноступенчатой одноцелевой холодильной установки с системой оборотного водоснабжения:

16 — вентиляторная градирня; 17 — вентилятор градирни с электродвигателем (остальные обозна­ чения те же, что на рис. 9.7)

9w Twi вн A7w

ATs

eg Щ

Рис. 9.13. Термоэкономическая модель одноступенчатой одноцелевой холодильной установки с системой оборотного водоснабжения:

gw _ плотность орошения в градирне; Д Tw% — температура охлажденной воды (остальные обозначе­ ния те же, что на рно. 9.8)

Для установки, схема которой изображена на рис. 9.12, а модель — на рис. 9.13, функционал, определяющий переменную часть приведенных затрат,

имеет вид

П = [Zu (eT, А Т?\ Т ? ) + Z12(e2, 0К, АТт Тщ) +

+ Z13{e?\ 0<13\ ДГ"3), T™) + Z21(e3, 0„, ATS) + Z22( e f A 7 f 2’) +

+ Z ai(eQ<JXJI, ©в, АГВ) + Z32( 4 1 . АП32))]тг+

+

Цэл(Б [Д и (4л, 0 к \ АГш’, TSl) А т ] +

S [^ 1з(е2л, 0 ^ , АГ<?, 7 ^ ) Дт] +

 

+ S

АТ?) Ат] + £

[Е32(е?охл, АТ?) Ат]) +

 

/=1

/=1

)

+

£ { М ? ,

0 к \ AT?, g?, T?J][V12(e?, в ? , AT?, g?, Т?)]Ах) +

 

/=1

 

 

+£ {г? [Е2(е?, 9?, АТ?) е?] Ат} +

/=1

+ £ [г? [Е3 (е?охл, в ? , АТ?) - е?] Ат} (/= 1, 2,

t). (9.84)

Выражение (9.84) написано с учетом определения тарифа на оборотную воду по выражению

/ _ “с.вУс.в + W l7 + г!66’ + г<™

Hw --

”l2

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

На основе функционала (9.84) составляется система уравнений:

£ ( - ^ [ « . и * + m

+ ш р v e + z g j ) + , ( - Щ

. + i j g . ) - «

/=1 ' К

/

\

к

к J

2

д0(Л (r?E? + Z?)

=

0;

 

/=1

г

 

 

 

/

0;

 

 

2

=

 

 

 

 

 

 

/=1

 

 

 

 

■(r?E? + Z?)

=

0;

 

/=1 L д в ?

ns

^7(13)

 

 

 

 

dz?>

 

гг P u . ( £ i f + £ 8 ) + Ц « ? + 7 1 Й + — З п г + — i s - _ 0 ;

дАТ<?

 

дАТ?>

<?д г <‘ 3 '

д г [им(£!f + £!Й + ц“Т8 + z|{j + -Ц£- +

аг"3’

Wt

 

О*ш,

 

10,

T 0XJl{ih

охлаждаемый объект,

 

П =

(9.86а)

 

*=1 /=1

Выражение (9.86а) определяет вид термоэкономической модели многоце­ левой холодильной установки. Эта модель может быть представлена рядом параллельных ветвей независимых одноцелевых холодильных установок, каж­ дая из которых предназначена для отвода теплоты от охлаждаемого объекта с фиксированной температурой при условии, что оптимизирующие пере­

менные зоны I(i)y характеризующие работу конденсатора (вк\ ДТ1^) и охлади­

теля воды ТЩ), общие для всех ветвей.

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

m

п= X {[2и(о (£2(!)\ i=i

, Д / ш , ^п.с(0»

п.с{ i ) l w x ) “Г

+ ^ 12 (0 Ь к, Д I w , 1 Wt> (^п.с(0» А * п.с(О) Н

 

+

^2(i) »

6'-'к'fl),> Д ГwГ ,»

Т{™\* »

'0n-'псСО».c(o. А Гп.^сШ/:

 

4

 

0.(Ш)

 

АТ.(.ио

 

ш х

 

й<!1‘L

Ш

 

 

 

/Д141)

 

 

Т.(.Ш)

 

Л Т }( 41)>

\

»

 

I 7

о ( 14^

»

д/р(140

»

T^(l4ft

»

с\(140

 

 

“Г ^14(0 (^2(1) ,

 

 

Д-* w

ill),

™п.с(0> Д* п.с(О)

“Г

 

+

^15(0 (^2(0»

 

 

Д^се»»

^

1» ®п.с(0> Д7*п.с(о) +

 

 

—Ь -2^21(0 (^3(t>, ®И(0»

Д Г 5(о)

+

2 22(0 (бзщ0 ,

Д Т 1? ? 0 )

+

2 з к ;) (бдохл(П,

0 В<£)> ДТ'вСо) “Ь

+

2 32(0 (^ охл (0» ^ ^ ’^ ) ))] х г} +

Дэл X

{ Х

[£*11(0 (£ад* © к \ Т ’Й*!,

 

 

 

 

 

 

 

 

Х =1

1/=1

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

Д Т Й , 0^ 0(0,

Д ^ п ^ Д т И -

S

[^13(0(4),

©(Л

Д Г Й ,

Г Й , 0 ffi*D, ДП^(0)А^] +

 

 

 

/=1

 

 

 

 

 

 

 

 

 

 

+S [^22(0 (взЙ)> ДТ'^Ц)) Ат] +

м

 

+

S

 

[Sum A

 

©к\ Д71Р, ГЙ, 0 п.с(о» А П^(0т] +

 

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ £

з2„-,(<хл,п, АГЙ»)Ах]} +

 

 

+

£

£

 

{ [ L W > ( e & . ,

0 к \

ArSf, Яш

, Г Й , е й с о . А ^сп)] X

 

i=l /=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

[ 1/(120

А

 

0 кУ),

А т {1 \

g i \ Т %

@пх(0,

Д Г Й с(о)] Д т )

+

 

 

+

S

 

S

 

{ ^ Й )

[^*2(0 ( 4 Й ) ,

0 и й ) »

A T ' s P o ) —

^2<о] Дт^} - Ь

 

 

 

 

1=1 /=1

 

 

 

 

 

 

 

 

 

 

 

 

+

£

 

£

 

{Г$, (£эд (< ОХЛ(0, ©йь

АГЙО -

ей,] Дт}

(Э.&7)

 

 

 

(=\

/=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(*

= 1.

2,

. . . .

/п;

/ =

1,

2, . . . ,

Q.

 

Соседние файлы в папке книги