Материал: Lektsia_-_metod_modalnogo_formirovania

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

располагаться внутри круга радиуса R с центром в начале координат.

Преобразование Ds в Dq посредством q f (s) является конформным.

На основе проведённых рассуждений для того, чтобы найти условия принадлежности спектра si , i 1, n области Ds, требуется указать условия нахождения спектра qi ,i 1, n матрицы B в круге.

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

"РКН – ГБ" требуется найти такую область DP

Pf , для которой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Spec B Dq , p j DP Pf Р,

j 1, k;

k 1,

(2.43)

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где B функционально – преобразованная посредством оператора L

 

 

) .

 

 

 

 

 

 

 

 

 

 

матрица

B

L(

A

 

 

 

 

 

 

 

 

 

 

 

 

 

S и её спектр si DS , то необходимо и достаточно,

 

Если матрица

A

 

A

 

 

было по

чтобы наибольшее по величине собственное число матрицы B

модулю меньше величины R. Тогда для того, чтобы si

DS , необходимо и

достаточно, чтобы удовлетворялось условие

 

 

 

Lim L

 

k 0 ,

 

 

 

A

(2.44)

k

 

 

 

где 0 – нулевая матрица.

 

 

 

 

 

 

 

 

 

 

 

 

В качестве функционала выберем спектральный радиус матрицы B :

 

 

 

 

 

 

 

 

J Rq max

qi

,i 1, n ,

(2.45)

 

 

i

 

 

 

 

 

 

 

 

 

где qi - собственные числа ФП матрицы

 

.

 

 

 

B

 

 

 

Рисунок 2.21 – Единичный круг на плоскости q

 

 

 

 

Таким образом, если

 

 

 

 

max

 

qi

 

1,

 

(2.46)

 

 

 

i

 

 

 

 

 

 

 

 

 

 

то оценка выполняется относительно единичного

круга

 

с центром

 

 

 

в начале координат (рисунок 2.21). Отсюда следует, что Spec

B

Dq

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

и значит Spec

A

DS .

 

 

 

 

 

 

 

 

S

 

 

 

 

Аналитическая структура q f (s) может

быть

различной и

отображать границы области заданного расположения спектра si на области,

ограниченные кругами или алгебраическими кривыми высших порядков,

расположенными внутри круга. Исходим из того, что близкие по динамическим характеристикам системы ("РКН АС" с различными ГЧ),

отвечающие определённым требованиям, будут иметь близкие по форме спектры. Поэтому область Ds принимается для всех однотипных систем одинаковой, а отличия будут только в размерах и положении этой области на комплексной плоскости s. При исследовании влияния степени быстродействия на динамическую совместимость "РКН ГЧ АС"

использовалось функциональное преобразование левой полуплоскости s (с

учётом смещения относительно мнимой оси Im) на внутренность единичного круга с центром в начале координат комплексной плоскости q.

При этом мнимая ось, рассматриваемая как окружность бесконечного радиуса, переходила в единичную окружность. При учёте остальных критериев качества на основе спектрального анализа, а также с учётом того,

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

Рисунок 2.22 – Наложение областей Dp на плоскость проектных параметров p1, p2 для различных участков полёта

На основе высказанных соображений выбрана алгебраическая кривая четвёртого порядка – овал Кассини (рисунок 2.20).

Поиск области Dp производится путём построения её границы Lp. В

принципе возможно построение некоторой объёмной области Dp в

пространстве трёх проектных параметров, однако это является малоприемлемым с практической точки зрения, тем более, что для данной конкретной задачи области Dp в пространстве выбранных проектных параметров строятся в нескольких точках траектории (старт, район

максимальных скоростных напоров, окончание работы первой ступени и т.д) (рисунок 2.22) с последующим наложением областей.

Поэтому граница Lp строится на плоскости двух выбранных проектных параметров из множества Pf, а остальные параметры либо изменяются по определённому алгоритму, либо остаются постоянными. При этом k-мерное пространство проектных параметров рассекается k(k-1)/2 плоскими сечениями, в каждом из которых строится своя граница Lp области Dp. Из комплекса рассматриваемых задач (динамическая совместимость систем

"РКН КГЧ АС") каждый раз выбирается конкретная задача (влияние параметров РКН, ГЧ или АС на совместимость системы), и поэтому реальное число областей Dp

c

 

k(k 1)

 

 

Dpi

 

, c k, с 1,

(2.47)

 

i 1

2

 

 

 

 

 

 

где c – количество доминирующих проектных параметров системы

"РКН – КГЧ – АС", определяющих области Dp в конкретной задаче динамической совместимости.

Так как граница Lp является линией равного значения функционала

(2.45, 2.46), то она может быть выделена путем расчёта максимального собственного числа ФП-матрицы B . Для этой цели используется одна из приближённых оценок спектрального радиуса:

для верхнего приближения

R Bm m

 

 

 

 

m

 

, m 1,2,...;

 

 

 

B

 

(2.48)

 

для нижнего приближения

 

 

 

 

 

 

 

 

R Hm

 

1

 

 

 

 

 

 

m .

 

spB

(2.49)

 

 

 

 

n

 

Здесь m – индекс итерации; n – порядок матрицы B . Норма в (2.48)

может быть вычислена по любому из следующих соотношений: l1 – норма

 

 

 

 

 

n

 

 

 

 

 

bij

;

B

 

 

 

 

l1

 

 

 

 

 

 

 

i, j 1

 

 

 

 

 

евклидова норма

 

 

 

 

 

 

n

 

2

12

 

 

 

 

 

 

 

 

 

 

 

;

B

 

e

 

bij

 

 

 

 

 

 

 

i, j 1

 

 

 

l - норма

 

 

 

 

 

max

 

bij

 

.

B

 

l

 

 

 

 

 

 

1 i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j n

Алгоритм вычисления области Dp сводится к определению её границы

Lp и состоит из следующих этапов:

 

 

1.

Определение глобального минимума функционала (2.45) на

множестве Pf с

использованием

алгоритма

случайного поиска с

направляющим косинусом.

 

 

2.

Проведение

сечения Pf (из

k k 1

возможных (2.47)) через

 

 

 

2

 

полученную точку глобального минимума функционала (2.45).

3.Организация сканирующих направлений из полученной точки и вычисление вдоль этих направлений оценок функционала (2.45) для верхнего

инижнего приближений по соответствующим формулам (2.48), (2.49).

4.Определение граничных точек области DР вдоль сканирующих направлений.

5.Апроксимация границы Lp по полученным граничным точкам алгебраическими кривыми заданного порядка.

Расчёт внутренней точки области Dp допустимых значений проектных параметров в выбранном сечении сводится к решению задачи нелинейного программирования для функционала (2.45) с учётом ограничений (2.38). Для отыскания глобального экстремума (2.45) предлагается использовать метод случайного поиска с направляющим косинусом. Метод применим как для случая многоэкстремальных задач, так и для случая, когда функционал (2.45)

не всюду дифференцируем, особенно в точке экстремума. Он может быть также применён для определения экстремума (2.45) на границе области Dp.

Источник: https://studfile.net/preview/15724095/