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

GOS for Iphone / mobile / ТАУ / Книги / ТАУ учебник

.pdf
Скачиваний:
231
Добавлен:
10.05.2015
Размер:
4.33 Mб
Скачать

Опыт проектирования САУ показывает, что для наилучшего приближения желаемого программного движения x0(t) к реальному процессу, протекающему в системе с синтезированными параметрами, коэффициент затухания ρ1 координатных функций целесообразно выбрать

ρ1 = α,

(2.73)

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

Остальные коэффициенты затухания ряда ρm–1 следует выбрать в виде геометрической прогрессии (со знаменателем прогрессии r = 2), т. е.

ρ

q

= ρ rq−1

= ρ 2q−1, q = 1,2,…,m,

(2.74)

 

1

1

 

что обеспечивает меньшее время затухания каждой из m–1 экспонент по сравнению со временем затухания первой координатной функции.

2.5.Синтез импульсных САУ

снесколькими нелинейными элементами

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

Особенности применения метода ортогональных проекций (обобщенного метода Галеркина) в случае синтеза САУ с несколькими нелинейностями заключается в следующем:

– необходимо определять желаемые программные движения x 0(t)

 

 

 

 

 

 

 

 

 

 

 

 

i

(i

=

1,

 

2, ..., r)

x

на

входах всех нелинейных элементов

F x

(t ), x

(t )

, F

 

 

(t ), x

(t )

 

i

i

 

i

 

 

 

 

 

 

 

 

;

 

i

i

 

i

 

– необходимо минимизировать r невязок из условия обеспечения заданных показателей качества работы синтезируемой САУ при ограничениях на грубость и абсолютную устойчивость системы, а также технических ограничениях, наложенных на значения искомых параметров.

Динамика импульсной системы управления с нелинейными элементами r описывается в матричной форме уравнением вида

 

 

*

f

(t ),

(2.75)

Qx + Q x

+ Ry + R y

= Sf (t ) + S

61

где x = x1(t), x2(t),..., xr(t) т; x* = x*1(t), x*2(t),..., x*r(t) т – векторыстолбцы непрерывных и импульсных процессов на входах нелинейных

элементов соответственно; y = y

(t), y

(t),..., y

(t)

т; y* = y*

(t),

y*

(t),..., y*

(t)

 

 

1

 

2

 

 

r

 

1

 

т – векторы столбцы непрерывных и импульсных про-

2

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

цессов на выходах нелинейных элементов; f(t), f*(t) – внешнее воздей-

ствие на входе и выходе импульсного элемента; Q, Q* – диагональные

матрицы порядка r вида

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

0

0

 

0

 

 

Q

0

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

0 Q

 

 

 

 

 

 

 

 

0 Q 0

0

 

 

0

0

 

 

 

Q =

 

2

 

 

; Q; =

 

2

 

 

 

 

 

 

 

0

0 Q

 

0

0

0

Q

 

0

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

Qr

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qr

 

 

 

 

 

 

 

 

 

 

 

 

R, R* – квадратные матрицы порядка r вида;

 

R R R

R

 

 

 

 

 

R11

R12

 

 

 

 

 

 

 

 

 

 

 

 

11

12

13

1r

 

 

 

 

 

R21

R22

R =

R21

R22

R23

R2r

 

 

 

 

=

R R R

R

 

 

;

R

 

 

 

31

32

33

3r

 

 

 

 

R31

R32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rr1

Rr2

Rr3

Rrr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rr1

Rr2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R13

 

R1r

 

 

 

 

 

R23

 

R2r

 

 

 

 

;

R33

 

R3r

 

 

 

Rr3

 

Rrr

S, S* – векторы-столбцы, содержащие r строк, вида

S =

 

S ,

S

 

,

,

S

 

 

T

;

 

=

 

S

 

,

 

,

,

 

 

T

 

 

 

 

 

 

 

 

2

r

 

 

S

 

 

S

2

S

r

.

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Матрицы Q, Q*, R, R*, S, S* являются функциями оператора обобщенного дифференцирования D и в общем случае функциями варьируемых параметров C = ck т, k = 1, 2,..., m.

Подставим вектор желаемого программного движения в уравнение (2.74) и образуем вектор невязки

Ψ(C,t) = Q(C,D)x0 + Q (C,D)x0

+ R(C,D)y0 +

 

(C,D)y

0

 

 

(t),

+R

 

S(C,D) f (t) + S

(C,D) f

где Ψ (C,t) – вектор-столбец невязки, определяемый следующим образом:

62

Ψ(C,t) =

 

ψ

(C,t), ψ

2

(C,t),

... , ψ

r

(C,t)

 

 

 

T .

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

а x0 = x0

(t), x0

(t),...,

x0

(t) т; x0* = x0*

(t), x0*

 

(t),..., x0*

(t) т

– векто-

1

2

 

 

 

r

 

 

 

1

 

 

2

 

 

 

 

r

 

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

нелинейных элементов соответственно; y0 = y01(t), y02(t),..., y0r(t) т; y0* = y0*1(t), y0*2(t),..., y0*r(t) т – векторы-столбцы непрерывных и

импульсных процессов на выходах нелинейных элементов при наличии на их входах желаемых программных движений.

Ортогональность невязки системе координатных функций приводит к следующей системе алгебраических уравнений:

ψi (C,t)ϕ q (t)dt= 0, i= 1,2,…, r.

0

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

r

m

 

2

J = ∑∑

 

 

, i= 1,2,…, r; q= 1,2,…, m. (2.76)

 

ψi (C,t)ϕ q (t)dt

i=1 q=1

0

 

 

Варьируемые параметры оператора управления (регулятора) определяются путем минимизации функционала (2.76) с помощью известных [89, 90] методов поиска экстремума целевой функции. На каждом шаге поиска параметров проверяется ограничение на устойчивость нелинейной импульсной системы по частотному критерию абсолютной устойчивости В. М. Попова [45], представленному в алгебраической форме [89, 90].

Процедура определения вектора желаемых программных движений на входах нелинейных элементов по желаемому программному движению, задаваемому на выходе синтезируемой САУ, была предложена И. А. Орурком [89], поэтому ниже рассмотрим ее примене-

63

ние на примере конкретной структурной схемы САУ с несколькими нелинейными элементами.

Для определения вектора желаемых программных движений записываются уравнения связи между входными координатами нелинейных элементов и желаемым программным движением на выходе системы

x0 = M1z0 (t ) + M2y + M3 f (t ),

(2.77)

где M1 = M1i – векторный оператор преобразований координаты выхода системы управления z0(t) в координаты входов нелинейностей xi0(t); M2 = M2i – векторный оператор преобразований координаты выхода нелинейных элементов yi(t) в координаты входов нелинейностей xi0(t); M3 = M3i – векторный оператор преобразований координаты f(t) в координаты входов нелинейностей xi0(t).

В частных случаях x0 может не зависеть от отдельных составляющих уравнения связи (2.77).

Рекомендуется следующая процедура определения желаемых процессов на входах нелинейных элементов, которая для наглядности иллюстрируется структурой системы управления приведенной на рис. 2.10, где приняты следующие обозначения: Wk (k = 1,..., 5) – передаточные функции линейных звеньев; Fi (i = 1,..., 5) – функциональные операторы нелинейных звеньев yi = Fi(xi, xi ).

f(t)

 

 

 

 

z(t)

W1(p)

F1

W2(ck,p)

F2

W3(p)

F3

 

 

 

 

 

 

 

 

W4(p)

F4

F5

 

 

W5(ck,p)

 

 

Рис. 2.10

При этом возможны следующие случаи:

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

пересчет процессов через нелинейное звено от входа к выходу;

пересчет процессов через нелинейное звено от выхода ко входу;

пересчет процессов через линейные звенья, содержащие варьируемые параметры.

64

В первом случае применяется символический метод подробно рассмотренный [89] и показанный ниже применительно к линейным САУ с несколькими частотами прерывания сигналов.

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

Сложнее обстоит дело в третьем случае, когда необходим обратный пересчет процесса от выхода нелинейного элемента к ее входу. Если по характеристике нелинейного звена можно построить обратную однозначную зависимость, то дальнейший расчет не будет отличаться от второго случая, отмеченного выше. Если же исходная нелинейная характеристика дает неоднозначность в определении обратной характеристики (имеет место при релейных характеристиках и характеристиках гистерезисного типа), то следует использовать итерационную процедуру, которая заключается в следующем. В начале определяется стартовая точка путем обобщенной линеаризации нелинейности секущей или касательной [89], и определяется процесс на входе нелинейного элемента x0 нулевого приближения. Затем по полученному процессу нулевого приближения вычисляются коэффициенты обобщенной линеаризации первого приближения, после чего определяется желаемый процесс x0 на входе нелинейности первого приближения и т. д. Итерационная процедура считается законченной в том случае, когда коэффициенты обобщенной линеаризации вычислены с заданной степенью точности.

Четвертый случай является наиболее сложным, поскольку определить процессы x и y в результате однократного расчета не представляется возможным, так как они должны определяться в каждой точке n-мерного пространства искомых параметров методом последовательных приближений, суть которого заключается в следующем. В начале необходимо задаться вероятными значениями всех или части параметров ck так, чтобы уравнения связи не содержали неизвестных значений. Далее после нахождения процессов на входах нелинейных элементов определяются значения варьируемых параметров. Решение задачи повторяется до тех пор, пока значения параметров i-го приближения не будут отличаться от значений параметров (i–1)-го приближения с требуемой погрешностью.

Теперь в качестве примера, иллюстрирующего предложенный выше подход, рассмотрим систему управления, структурная схема которой 65

приведена на рис. 2.10. Для данной САУ процессы на входах четвертого F4 и пятого F5 нелинейных элементов по процессу z0(t) определяются непосредственно, так как процесс на входе F4 совпадает с z0(t), а процесс на входе F5 определяется символическим методом.

Желаемые процессы y0(t) на выходах F4 и F5 могут быть определены либо с помощью кусочно-линейной аппроксимации, либо с помощью обобщенной линеаризации.

Выходная координата третьего нелинейного элемента F3 совпадает с сигналомz0(t).Следовательно,дляопределенияжелаемогопроцессанавходе данной нелинейности необходимо построить обратную нелинейную характеристику. Если полученная обратная характеристика однозначна, то дальнейшие вычисления не отличаются от рассмотренного выше для F4. В случае неоднозначности обратной нелинейной характеристики необходимо осуществить итерационную процедуру, описанную выше.

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

лишь одного нелинейного звена F5. При движении справа в начале необходимо определять процесс на выходе нелинейного звена F1. Для этого требуется осуществить линеаризацию трех нелинейных элементов

F2, F3, F4 и определять искомый процесс в каждой точке пространства варьируемых параметров. Кроме того, дальнейшее определение желае-

мого сигнала на входе F1 с помощью обратной нелинейной характеристики может дать неоднозначность решения.

Кнелинейному элементу F2 также можно подойти с двух сторон, причем в обоих случаях необходимо предварительно линеаризовать два нелинейных звена. Однако, движение справа предпочтительнее, поскольку в этом случае в уравнение связи не будет входит передаточная функция оператора управления, содержащая искомые параметры.

Таким образом, при движении последовательно от выхода САУ для

двух нелинейных элементов F4 и F5 желаемые процессы на их входах определены точно, а для остальных – приближенно с точностью до погрешности метода обобщенной линеаризации. В то же время при движении от входа системы необходима одновременная линеаризация всех нелинейных элементов, что требует достаточно трудоемкой итерационной процедуры, при этом определение процессов на входах всех нелинейностей должно производиться в каждой точке пространства варьируемых параметров.

66

Глава 3

СИНТЕЗ КУСОЧНО-ЛИНЕЙНЫХ СИСТЕМ УПРАВЛЕНИЯ

САМПЛИТУДНО-ИМПУЛЬСНОЙ МОДУЛЯЦИЕЙ

Вданном разделе проводится детальная проработка общей схемы решения задачи синтеза импульсных САУ, изложенная в предыдущей главе, применительно к нелинейным системам с амплитудно-импульс- ными модуляторами.

Рассматривается получение рекуррентных аналитических выражений, определяющих интегралы Галеркина амплитудно-импульсных систем как с однозначными, так и неоднозначными характеристиками нелинейных элементов, для идеального импульсного элемента и с учетом конечной длительности замыкания ключа.

Полученные соотношения позволили полностью алгебраизировать решение поставленной задачи и построить унифицированный алгоритм синтеза нелинейных амплитудно-импульсных систем произвольно высокого порядка.

3.1. Системы автоматического управления

сидеальным импульсным модулятором

Втеории импульсных систем, как было отмечено выше, импульсный элемент обычно считают идеальным, т. е. сигнал на его выходе может быть представлен в виде (2.1). Рассмотрим решение задачи синтеза параметров САУ, содержащей идеальный АИМ и нелинейный элемент, характеристика которого допускает кусочно-линейную аппроксимацию. Необходимо отметить, что впервые обращение обобщенного метода Галеркина на решение задачи синтеза САУ с идеальным импульсным элементом было рассмотрено в ряде совместных работ В. Ф. Шишлакова и Л. А. Осипова.

Уравнение движения САУ, записанное относительно координаты входа нелинейности (исследуемая координата) x(t), имеет вид (2.17), либо с учетом выражений для полиномов оператора обобщенного дифференцирования D

67

 

 

 

 

n

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

i

k

 

 

 

 

 

 

 

 

 

i

 

 

k

)Di x

 

 

 

 

i

 

k

 

 

 

 

 

 

 

 

 

 

a

(c

)Di x (t) +

 

 

a

(c

(t) +

 

 

b

(c

)Di y (t) +

 

 

i=0

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

+

i

k

)D

y

 

 

 

 

 

i

 

k

)D

 

 

 

 

 

 

 

i

 

k

 

 

f

(t);

 

 

 

b

(c

(t) = e

(c

 

i

f (t) +

 

e

(c )D

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y (t) =

F x (,t)( x)

t

,

 

 

 

y (t)

=

F

 

x

(,t)

(x)

t;

 

(3.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

(

 

t

)

= F x

(

t

)

 

 

, y

 

(

t

)

=

F

 

 

 

(

t

)

 

 

 

 

 

 

 

 

 

 

где

 

, x(t)

 

 

x

x t

 

 

– сигналы на выходе

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нелинейного элемента при непрерывном x(t) и импульсном x*(t) входном сигнале соответственно.

В соответствии с общей схемой решения задачи параметрического синтеза нелинейных САУ методом ортогональных проекций (обобщенным методом Галеркина), изложенной выше, необходимо задаться желаемым программным движением. При синтезе нелинейных САУ n-го порядка (в соответствии с рекомендациями, изложенными во второй главе программное движение целесообразно (в первом приближении) задать в виде

x0 (t) = xу + H eαt cos (β t− ϕ 0 ) 1(t).

(3.2)

Тогда сигнал на выходе идеального АИМ, согласно выражению (2.1), будет

 

0

 

α nT

 

 

 

(t) =

cos(βnT − ϕ 0 ) δ (tnT ).

 

x

 

xу + H

e

(3.3)

 

 

n=0

 

 

 

 

Для определенности задачу синтеза рассмотрим при заданных граничных (начальных (2.13) и конечных (2.14)) условиях и внешнем воздействии вида f(t) = H1(t), которое на выходе импульсного элемента (в соответствии с (2.1)) будет описываться следующим образом:

f (t) = H1 (t) =

(t nT ).

(3.4)

n=0

 

 

Подставляя желаемое программное движение (3.2) в уравнение движения САУ (3.1), образуем невязку, а из условия ортогональности не-

68

вязки координатным функциям (2.19), выбранным в виде (2.72), получаем систему из m линейных алгебраических уравнений

где Aqi

Aqi

Bqi

Bqi

Cqi

Cqi

 

 

 

 

 

 

n

 

 

a (c )A +

n

a

 

(c

)A

 

u

(c

)B +

 

 

 

 

 

 

 

 

 

 

 

+ b

 

 

 

 

 

 

 

 

 

 

 

 

 

qi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

k

 

 

i

 

k

qi

 

i

k

qi

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

+

 

 

 

(c

e

 

(c )C

(c

0,

 

 

 

 

b

)B

 

e

)C =

 

 

 

 

 

 

 

 

i

 

 

 

k

qi

 

 

i

 

 

k

qi

 

i

k

qi

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q =1,2,,m,

 

 

 

 

=

 

 

{x0 (t)} e−ρq tdt,

 

 

 

 

 

 

 

 

 

 

Di

i = 0,1,…,n;

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

{x

0

 

 

 

 

−ρq t

 

 

 

 

 

 

 

 

 

 

D

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t)} e

 

dt,

 

i = 0,1,…,n ;

 

 

 

 

0

 

 

{F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−ρq t

 

 

 

 

 

D

i

 

x

0

 

 

0

 

 

 

 

 

= 0,1,…,u;

 

 

 

 

 

 

(t), x

(t)

} e

dt, i

 

 

0

 

 

{F

 

 

 

 

 

 

 

 

 

 

 

 

 

} e

 

 

 

 

 

 

=

 

 

 

 

 

 

0

 

 

 

0

 

 

 

−ρq t

 

 

 

 

D

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

(t), x

 

(t)

 

 

dt,

i = 0,1,…,u ;

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

{f (t)} e−ρq tdt, i = 0,1,…,v;

 

 

 

 

Di

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

{

 

 

 

 

 

 

}

−ρq t

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

f

 

 

(t)

 

e

 

dt,

i = 0,1,…,v .

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.5)

(3.6)

В работах [90, 187, 201–204] подробно рассмотрена методика вычисления интегралов Галеркина (3.6) импульсных САУ, содержащих идеальный модулятор. Поэтому ниже приводятся лишь окончательные рекуррентные аналитические выражения, определяющие данные интегралы для процессов вида (3.2), (3.3)

 

 

 

 

 

 

e− α t cos(β t− ϕ

 

)

1(t) e−ρqtdt=

 

A

=

Di

x

у

+ H

 

A ρi− 1,

qi

 

 

{

 

 

0

 

}

q q

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.7)

 

 

 

 

 

 

 

i = 1,2, ,n,

 

 

69

 

 

H

 

(α + ρ

q

)cos ϕ

+

βsinϕ

 

ρ

q

 

Aq = xу +

 

 

 

 

0

 

0

 

где

 

 

 

 

 

 

 

 

 

 

,

 

 

(α + ρq )2 + β2

 

 

 

 

 

 

 

 

 

 

 

q = 1, 2,, m.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

αnT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−ρqt

 

 

 

i

 

 

D

i

 

xу

 

 

 

 

 

cos(βnT − ϕ 0 )

δ

(t

 

e

dt=

 

Aqi =

 

 

 

 

+ H e

 

 

 

 

 

nT )

 

 

 

Aqρq ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i = 1,2,…,n

,

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

(

α+ρ

q )

 

 

 

(

 

 

q )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

e

 

T

cos ϕ 0

e

 

α

 

T

cos (βT+ ϕ

 

0 )

 

 

 

 

 

 

 

 

 

 

xу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

где A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−ρqT

 

 

 

 

 

 

 

 

2

(

α+ρ

T

− 2e(

α+ρ

T

 

 

 

 

 

 

 

 

q

 

1

e

 

 

 

 

 

 

 

 

 

 

cosβT + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

q )

 

 

 

q )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q = 1, 2,, m;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cqi = Di {H1(t )}e−ρqt dt = Cqρiq1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i = 1,2,, v, q = 1,2,, m;

 

 

 

 

 

 

 

(3.9)

где C

 

= H;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

D

i

 

 

 

 

 

δ (t

 

 

 

 

 

 

 

−ρqt

 

 

 

i

 

 

 

 

 

 

 

(3.10)

 

 

Cqi

 

 

 

 

nT )

e

dt = Cqρq ,

i = 1, 2,…, v ,

 

 

 

 

 

 

 

 

 

 

0

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где C

=

 

 

H

 

 

 

 

;

 

 

q = 1,2,

 

,

m.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−ρqT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

1 − e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для вычисления интегралов Bqi

и Bqi

– нелинейные функции

F x

 

 

 

 

 

 

 

 

 

,

 

F

 

x

 

 

 

 

 

 

 

 

необходимо представить в виде обоб-

(t ), x(t)

 

 

 

 

(t ), x (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

щенных функций [205]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= F+0 (t)1(t) +

r

 

 

 

 

 

 

 

 

1(t t j ),

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F x

 

 

(t), x (t)

 

 

 

 

 

 

 

 

 

 

(3.11)

 

 

 

 

F+ j (t) F− j (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

70