Главная » Просмотр файлов » Диссертация

Диссертация (1145311), страница 10

Файл №1145311 Диссертация (Применение алгебраических методов для анализа сложных систем) 10 страницаДиссертация (1145311) страница 102019-06-29СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 10)

Однако метод, представленный в главе I, может быть применен также к системам с бо́льшим количеством параметров, а также для болееширокого класса областей Ω.73Замечание 7. Здесь хотелось бы отметить следующее. Все приведенные методы таковы, что можно было бы вычислить число корней исходного полинома в заданной алгебраической области комплексной плоскости как функцию параметров крана и построить графики, определяющие области устойчивости в пространстве параметров. Однако имеющихся вычислительныхсредств для этих целей, к сожалению, недостаточно. Поэтому вычисленияпроведены только для конкретных значений параметров. Решение данной задачи в общем виде приведено далее.В работе [109] приводится также другое решение той же задачи, котороеиспользует параметры Маркова и результаты Кронекера.

Это решение такжеалгебраическое, результат можно получить, проделав конечное число алгебраических операций над коэффициентами полиномов.Сначала, как и ранее, найдем число вещественных корней полинома f (x),удовлетворяющих неравенству g(x, 0) > 0. Сделать это можно изложеннымвыше способом.Рассмотрим теперь задачу отделения комплексных корней. Так же, как иранее, сводим нашу задачу к исследованию системы уравнений (1.51)Φ(x, Y ) = 0, Ψ(x, Y ) = 0.(Здесь использованы те же обозначения, что и выше.)Исключим Y из этой системы уравнений.

РассмотримX (x) = RY (Φ, Ψ).Разложим дробь iΨ(x, Y )/Φ(x, Y ) в ряд по отрицательным степеням Y :iΨ(x, Y )µ0 (x)µ1 (x)µ2 (x)= iµ−1 (x) + i− i 2 + i 3 − ...Φ(x, Y )YYYОчевидно, что функции iµ−1 (x), iµ0 (x), iµ1 (x), . . . являются параметрами Маркова полинома f (x + iy), если рассматривать его как полином от y с коэффициентами, зависящими от x.74По следствию 6 к теореме 14 результант X (x) может быть выражен черезпараметры Маркова:Теорема 28.X (x) = b2n−εM[n/2] ,0где1b0 = f (n) (x), ε = 0n!1b0 =f (n−1) (x),(n − 1)!если n четно,ε = 1 если n нечетно.Несложно также показать, что deg X = n(n − 1)/2.Далее будем считать, что предположение 1 раздела 1.2 справедливо.Тогда мы можем представить вторую компоненту решения (1.51) в видевещественной рациональной функции первой:Bj = r(αj ), j = 1, N .Используя следствие 5 к теореме 14, мы можем выразить функцию r черезпараметры Маркова:Теорема 29.Bj = r(αj ) = −µ0. .

. µ[n/2]−2µ1. . . µ[n/2]−1...... ...µ[n/2]−3 . . . µ2[n/2]−4µ[n/2]−1 . . . µ2[n/2]−2M[n/2]−1−a1.a0Пример 4. Выразим вторую компоненту решения (1.51) через первую дляполиномаf (z) = a0 z 6 + a1 z 5 + a2 z 4 + a3 z 3 + a4 z 2 + a5 z + a6 .75Мы можем записатьf (x + iy) = b0 (x)y 6 + b2 (x)y 4 + b4 (x)y 2 + b6 (x)+iy[b1 (x)y 4 + b3 (x)y 2 + b5 (x)]= −a0 y 6 + (15a0 x2 + 5a1 x + a2 )y 4 + (−15a0 x4 − 10a1 x3−6a2 x2 − 3a3 x − a4 )y 2 + a0 x6 + a1 x5 + a2 x4 + a3 x3 + a4 x2 + a5 x + a6+iy[(6a0 x + a1 )y 4 + (−20a0 x3 − 10a1 x2 − 4a2 x − a3 )y 2+6a0 x5 + 5a1 x4 + 4a2 x3 + 3a3 x2 + 2a4 x + a5 ].Вычислим параметры Маркова f (x + iy) (считая его полиномом от y скоэффициентами, зависящими от x): iµ−1 (x), iµ0 (x), iµ1 (x), .

. .Рассмотрим систему уравненийΦ(x, Y ) = 0, Ψ2 (x, Y ) = 0.Результант этой системы по исключении Y b0 (x) b2 (x) b4 (x) 0b0 (x) b2 (x)RY = 00b1 (x) 0b1 (x) b3 (x) b1 (x) b3 (x) b5 (x)b6 (x)0 b4 (x) b6 (x) b3 (x) b5 (x) b5 (x)0 00 µ0 µ 1 µ 2 5 3= b0 i µ1 µ2 µ3 = b50 M3 . µ2 µ 3 µ 4 Первый субрезультант b0 b2 b4 b0 b2 b4 µµ01 = b30 M2 ;R1 = 0 b1 b3 = − b1 b3 b5 = −b30 µµ12 b1 b3 b5 0 b1 b3 76 1 b0 b2 b6 00 b0 b2 b6 2R̃1 = 0 b1 b5 = −i µ0 −µ1 µ2 0 b0 b4 0 µ0 −µ1 0 0 b2 b1 b3 0 b0 b2 0 −b0 i b1 b3 −µ3 = −b20 b2 M2 + b30 µ0 µ3 − b30 µ2 µ1 0 b 1 µ2 µµ01= −b20 b2 M2 + b30 µ2 µ 3 Таким образом,r(x) = −b2+b0 µ 0 µ1 µ 2 µ3 M2=−b2−b0 iµ0 iµ1 iµ2 iµ3 M2.Таким образом, используя изложенные выше результаты и теорему 26, мыможем определить, сколько комплексных корней полинома f (z) удовлетворяютусловию g(x, y) > 0.

В данном методе, в отличие от раздела 1.2, используютсярезультаты Кронекера и параметры Маркова вместо теории исключения.Разумеется, с помощью приведенных выше результатов можно решить поставленную задачу и в более общем виде, а именно:Определить число корней полинома f (z), удовлетворяющих условиям:g1 (x, y) > 0, g2 (x, y) > 0, . . .

, gk (x, y) > 0.Замечание 8. Хотя мы решили задачу в более общем случае, чем случай Рауса-Гурвица, необходимо помнить, что если наша область находится в левойполуплоскости комплексной плоскости, условия Рауса-Гурвица должны выполняться. Имеется связь между этими двумя видами условий:Теорема 30. (Раус, [10]). Полином f (z) (a0 > 0) является устойчивым втом и только в том случае, когда выполнены следующие условия:77а) все коэффициенты f (z) положительны;б) если старший коэффициент X (x) положителен, то и все остальныекоэффициенты X (x) положительны.Пример 5.

(Решаем пример 2 с помощью метода, описанного в настоящемпараграфе.) Для полиномаf (z) = z 4 − z 3 − 2z 2 + 6z − 4найти число корней, удовлетворяющих условиюg(x, y) = 2y 2 − x > 0.Решение. Число вещественных корней f (z) равно 2 (пример из раздела 1.2.1).Найдем сначалаnrr {f (z) = 0|g(x, 0) > 0} = 1.Таким образом, один из двух вещественных корней f (z) удовлетворяет неравенству g(x, 0) > 0 (см.

раздел 1.2.1).Для того чтобы определить, удовлетворяют ли нашему неравенству комплексные корни, рассмотрим систему (1.51):Φ1 (x, Y ) = Y 2 + (−6x2 + 3x + 2)Y + x4 − x3 − 2x2 + 6x − 4 = 0(1.66)32Φ2 (x, Y ) = (−4x + 1)Y + 4x − 3x − 4x + 6 = 0.Вычислим параметры Маркова:µ0 = −4x + 1, µ1 = 20x3 − 15x2 − x − 4,µ2 = −116x5 + 145x4 − 6x3 + 17x2 − 36x − 4.Отсюда результантX (x) = −4(16x6 − 24x5 − 4x4 + 14x3 + 10x2 − 7x − 5)(D(X ) ≈ 2.01852 × 1013 6= 0).78Решения (1.66) связаны формулой6f 0 (x)Y = (3) .f (x)Теперь найдем число комплексных корней f (z) таких, что g(x, Y ) > 0:nrr {X = 0|G(x, r(x)) > 0} = 1.Итак, по теореме 26, число комплексных корней f (z), удовлетворяющих данному неравенству, равно 2.Таким образом, три корня полинома f (z) удовлетворяют условию2y 2 − x > 0.Устойчивость полиномовТеперь рассмотрим множествоf (z) = z n + a1 z n−1 + .

. . + an (a1 , a2 , . . . , an ) ⊂ Rn(1.67)всех полиномов от одной переменной как множество точек в n-мерном пространстве коэффициентов этих полиномов.Справедлива следующая теорема [90]:Теорема 31. Области множества (1.67), соответствующие полиномам содинаковым числом вещественных корней, разделяются дискриминантной поверхностьюD(a1 , . .

. , an ) = Dz (z n + a1 z n−1 + . . . + an ) = 0.Здесь дискриминант вычисляется относительно переменной z.Нас будет интересовать устойчивость семейства вещественных полиномов.Для одного полинома необходимые и достаточные условия устойчивости даеттеорема Гурвица [7]:79Теорема 32. Вещественный полином f (z) (a0 6= 0, n ∈ N) устойчив (т. е.все его корни лежат в левой полуплоскости комплексной плоскости) тогда итолько тогда, когда выполняются следующие неравенстваH1 > 0, H2 > 0, H3 > 0, .

. . , Hn > 0.ЗдесьHj = a1 a3 a5 . . .a0 a2 a4 . . .0a1 a3 . . .0a0 a2 . . ....aj ,(1.68)(ak = 0 при k > n).j×jОпределение 14. Определители Hk (k = 1, 2, . . . , n) называются определителями Гурвица.Рассмотрим последний определитель из условия (1.68). Оказывается, чтос точностью до знака этот определитель равен результанту четной и нечетнойчастей полинома f (z), умноженному на an :Hn = an R a1 Y b(n−1)/2c + a3 Y b(n−1)/2c−1 + . . . + a2b(n−1)/2c+1 ,bn/2cbn/2c−1a0 Y+ a2 Y+ .

. . + a2bn/2c .Обращение этого определителя в нуль при выполнении условий H1 > 0, H2 >0, . . . , Hn−2 > 0 эквивалентно тому, что полином f (z) имеет корень λ = 0, либо он имеет пару чисто мнимых корней λ1,2 = ±iβ, β ∈ R. Это свойство Hnимеет важное значение при решении проблемы устойчивости семейства полиномов (1.1)p = f (z, ν1 , ν2 , . . . , νk ).В самом деле, при изменении вектора параметров (ν1 , . .

. , νk ) внутри параллелепипедаB = {(ν1 , ν2 , . . . , νk )|ν 1 ≤ ν1 ≤ ν 1 , ν 2 ≤ ν2 ≤ ν 2 , . . . , ν k ≤ νk ≤ ν k }80условие устойчивости для соответствующего полинома f может быть нарушенотолько при обращении в нуль Hn .В дальнейшем мы будем также использовать следующее необходимое идостаточное условие устойчивости семейства полиномов [123]:Теорема 33. Все полиномы семейства (1.1) устойчивы тогда и только тогда, когда выполняются следующие условия:(i) Существует точка (ν10 , .

Характеристики

Тип файла
PDF-файл
Размер
1,69 Mb
Высшее учебное заведение

Список файлов диссертации

Применение алгебраических методов для анализа сложных систем
Свежие статьи
Популярно сейчас
А знаете ли Вы, что из года в год задания практически не меняются? Математика, преподаваемая в учебных заведениях, никак не менялась минимум 30 лет. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6376
Авторов
на СтудИзбе
309
Средний доход
с одного платного файла
Обучение Подробнее