Файл: Алгоритмы - шпоры (Final).doc

З6’ Решение уравнений установившегося режима на основе L-H факторизации матрицы коэффициентов. Вычислительная схема прямой и обратной подста­новки.

1-й этап:

Решается прямой подстановкой

2-й этап:

Решается обратной подстановкой

только везде вместо g -> q

37’ L-H факторизация матрицы коэффициентов системы линейных алгебраиче­ских уравнений (матрицы А). Алгоритм вычисления элементов факторизованной матрицы.

38’ Алгоритмическая и программная реализация L-H факторизации матрицы А.

Subroutine lhfact (A,n,l)

Complex A(l,l),t

Do k=1, n-1

Do j=k+1, n

t=a(k,j)/a(k,k)

do i=k+1, n

a(i,j)=a(i,j)-a(i,k)*t

enddo

enddo

enddo

return

end

Совмещает достоинства методов обращения и Гаусса:

1)Операция факторизации только 1 раз → слабозаполнен. Матрицы L и H.

2)Вычисления существенно проще Гаусса, где на каждой итерации пересчит. матрица; алгоритмич прост.

“-“ ген. узлы только неопорные

Вывод: ↓ объем вычислений и объем требуемой памяти.

Алгоритм см билет 36!

39’/41’ Методы, используемые для расчета УР при записи узловых уравнений в форме баланса токов. / Методы, применяемые для решения комплексного узлового уравнения в форме баланса токов.

1) Метод Зейделя

- прост в алгоритмическом отношении

- естественным образом учитывает слабую заполненность матрицы узловых проводимостей

- при расчете утяжеленных режимов область сходимости невелика

- нечувствителен к начальным приближениям

2) Метод Гаусса

- важным достоинством является высокая скорость решения

- надо хранить в памяти как исходную, так и пересчитанную матрицы

- на каждом шаге прямого хода надо выбирать главный элемент

3) Метод обращения матрицы Y

- обращение матрицы узловых проводимостей проводится 1 раз

- из слабозаполненной матрицы Y получается сильнозаполненная Z

4) LH-факторизация

Совмещает достоинства методов обращения и Гаусса.

- операция факторизации производится 1 раз => слабозаполненные матрицы L и H

- вычисления существенно проще метода Гаусса, где на каждой итерации пересчитывается матрица; алгоритмически прост

- может использоваться только с НЕОПОРНЫМИ ген. узлами.

=> малый объем вычислений и памяти ПК.

40’ Методы расчета режимов, основанные на сочетании методов Зейделя и Гаусса. Достоинства и недостатки.

Матрицу Y разбиваем на блоки, выделяя блок генераторных и опорных узлов.

Решая ур-е (1) методом Зейделя относительно напряжений в ген. оп. Узлах, получим:

Из ур-я (3):

Ур-е (4) решаем на основе метода Гаусса, ищем комплексы напряжений в нагрузочных узлах.

Только метод Зейделя даст возможность учесть опорные ген. узлы.

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

Поскольку брать производные по комплексным величинам нельзя => переходим от записи в комплексном виде к 2м уравнениям с веществ. коэффициентами.

(n-1) ур-й с компл. → (2n-2) с веществ.

Пр. Метод Ньютона

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

43’ Узловое уравнение в форме баланса мощности, записанное в прямоугольной системе координат.

44’ Узловое уравнение состояния эл. сист. в форме баланса S при записи напряжений в полярной, а проводимостей – в прямоугольной системах координат.

45’ Узловое уравнение состояния эл. сист. в форме баланса мощности, записанное в полярной системе координат.

46’ Возможные формы записи нелинейных узловых уравнений установившегося режима для решения их методами, требующими разделения комплексных переменных на действит.сост.

1)узловые уравнения в форме баланса токов в прямоугольной системе координат

2)---//-----//--- в полярной системе координат

3)узловые уравнения в форме баланса мощности в прямоугольной системе координат

4)---//-----//--- в полярной системе координат

47’ Метод Ньютона. Решение узлового уравнения методом Ньютона, записанного в прямоугольной и полярной системе координат.

1)Начальные приближения

2)В точке начального приближения данная система линеаризуется путем разложения в ряд Тейлора и отбрасыванием нелинейных частей.

Т.о. метод Ньютона сводится к многократному решению (на каждой итерации) СЛАУ.

3)

Прямоуг. сист. коорд.

Полярная ;

48’/49’ 3ависимость размерности матрицы Якоби от формы представления генераторных узлов и системы координат, в которой записаны узловые уравнения в форме баланса мощности.

Опорные: k- опорных узлов

- входит в уравнение баланса Q

Вывод: При опорных генераторных узлах кол-во уравнений в полярной форме снижается (при решении уравнения в форме баланса мощности)

50’ Аналитическое выражение элементов матрицы Якоби узловых уравнений в форме баланса мощности, записанных в полярной системе координат.

При решении уравнений установившегося режима, записанных в форме баланса мощности, предпочтительным является решение в полярной системе координат:

50’ Продолжение.

Аналогично,

С билетом тебе не повезло. Улыбайся преподу 

51’ Свойства матрицы Якоби. Свойства метода Ньютона.

Св-ва матрицы Якоби:

1) Слабозаполненная

2) Структурно-симметричная, но числовой симметрии нет:

Структурная симметрия облегчает учет слабой заполненности, поиск и хранение ненулевых элементов.

3) Диагонально-доминирующая: Св-ва совпадают со св-вами м.Y

Св-ва матрицы Ньютона:

1) Сильная чувствительность к начальным приближениям.

2) Квадратичная сходимость

1. При отсутствии активных ограничений:

– Наличие резервов акт. и реакт. мощности, – Большой Кз по ста­ти­­ческой устойчивости, – Отсутствие перегрузок по линиям

2.При наличии активных ограничений:

­– Дефицит Q, – Предел по станциям (вышли на ограничения по генерат. узлам), – Pг, Uг, Pг,Qг

3.Сходимость может ухудшиться при плохой обусловленности матрицы Якоби

(режимы, близкие к пределу по статике; резко неоднородная сеть)

4.Погрешности исходных данных влияют на сходимость метода Ньютона, решение системы уравнений по Гауссу.

5.Хорошо согласуется с методами решения оптимизац-ых задач.

6.Трудоемок в части алгоритмического представления.

52’ Модификации метода Ньютона

Метод Ньютона основан на решении всех уравнений системы на каж­­­­дой итерации, это повышает вычислительную эффективность. Применение этих методов требует перехода от (n-1)-го уравнений с комплексными к (2n-2)-м ур-ям с вещественными коэффициентами и переменными. Это связано с тем, что применение этих методов тре­бует диффер. ур-й по искомым переменным (а про­из­вод­ные по комплексным переменным не определены). Кроме того, для общего случая задания оп. генер. узлов (P, U), иско­мые переменные Q, δ– вещественные. Переход осуществляется на ос­нове записи компл. чисел в прямоуг. или полярной системах ко­ординат.

1) В мет. Н с обращением матрицы Якоби СЛАУ решается на каждой итерации относительно вектора приращения независимых перем. с использованием обратной матр. Якоби:

[ 2]

3) С разделением параметров (блочная диагонализация). Существенно уменьшает объем расчетов на ит., из-за отбрасывания недиаг. блоков матр. Якоби, т.е. полагая ∂P/∂δ=0 и ∂Q/∂δ=0 (см.55’)

4) Модифицированный – матрица Якоби вычисляется только 1 раз.

53’ Модифицированный метод Ньютона. Сущность и область

сходимости.

Основан на том, что если Δx ^(k)→ 0, то и ΔW(x ^(k–1))→ 0, т.е. можно вычислять матрицу Якоби 1 раз.

Эту матрицу можно однократно факторизовать и использовать в процессе итераций в факторизованном виде, что значительно уменьшает объем вычислений на каждой итерации.

Особенность: из-за невысокой сходимости этот метод применяется только для нетяжелых режимов.

54’ Метод Ньютона по параметру. Сущность и область метода.

Введение параметра t (0< t <1) т.е. переход от итерационной формулы к выражению ниже есть переход к мет. Н. по параметру (при t=1 – простой мет. Н)

Данный метод применятся для расчетов режимов близких к преде­лу по статике.

Область сходимости выше чем у обычного метода Ньютона.

Недостаток – большой объем вычислений по итерации: вычисление элементов матрицы Якоби и вектора небалансов, решение СЛАУ

55’ Способы определения параметра t в методе Ньютона по параметру.

.

Если функция не уменьшается, то шаг делится пополам.

На одном шаге:

– уменьшение шага. Многократные вычисления на одном шаге.

Если , то - метод Ньютона,

Если , то .

56’ Метод Ньютона с разделением переменных (блочная диагонализация матрицы Якоби). Сущность и область сходимости метода.

Применяется при незначительной разнице в фазах векторов напряжения в узлах, что соответствует нетяжелым режимам сети. Предполагается, что в виду малой зависимости P от U и Q от δ можно принять и .

Отсюда следует, что итерационный процесс описывается следующими уравнениями:

Решая последовательно эти системы уравнений относительно приращений δΔ и ΔU либо методом Гаусса, либо методом обращения матриц, можно получить значения независимых переменных на следующем шаге:

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

57’ Модификации метода Ньютона, сокращающие объем вычислений на каждой итерации. Скорость и область сходимости этих методов.

К модификациям метода Ньютона, позволяющим сократить объем вычислений на каждом шаге относятся Модифицированный метод Ньютона и метод Ньютона с разделением переменных.

В модифицированном методе Ньютона сокращение числа расчетов достигается за счет устранения необходимости пересчета коэффициентов матрицы Якоби на каждом шаге итерационного процесса. Коэффициенты матрицы Якоби рассчитываются только один раз на первом шаге с нулевыми приближениями, и далее используется эта матрица. За счет подобного упрощения существенно страдает сходимость метода, так как значения коэффициентов матрицы Якоби могут достаточно существенно изменяться при переходе к следующему шагу.

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

58’ Методы расчета установившихся режимов сложных ЭЭС. Их сопоставление.

Абсолютно сходящихся методов нет!

Эффективность метода оценивается по 2 критериям:

1. Количество вычислений на одной итерации;

2. Общее число итераций для достижения заданной точности.

Методы – Зейделя (считает все, что считается), Ньютона + его модификации (сильно чувствительны к выбору начальных приближений) и градиентные методы.

Зейделя – см. 27, 28

Ньютона + модификации – см. 47, 52-54.

Градиентные – см. 67.

59’/60’ Использование метода Гаусса, обращения, L-H факторизации матрицы коэффициентов системы линейных алгебраических уравнений в расчетах режимов методом Ньютона.

При расчетах режимов методом Ньютона на каждом шаге приходится решать СЛАУ, например . Причем метод ее решения может быть выбран из уже рассмотренных:

1. Метод Гаусса – см. 18

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

3. L-H факторизация – см. 35, 36.

61’ Характеристика современных программно-вычислительных комплексов (ПВК), содержащих расчет УР.

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

Основные из них:

А) РАСТР (расчет УР, эквивалентирование, расчет предельных режимов, оптимизация режимов по Q, РПН)

Б) ДАКАР (расчет УР с учетом изменения частоты, возможность учета СХН, потерь на корону, РПН, утяжеление режимов, расчет колебательной статической устойчивости, расчет эл.-мех. п.п., несимметричных и неполнофазных режимов, эквивалентирование)

В) МУСТАНГ (расчет УР с возможностью учета СХН и вставок постоянного тока, утяжеление режимов, расчет эл.-мех. п.п.)

Г) EUROSTAG (расчет УР с учетом изменения частоты, расчет колебательной устойчивости, динамической устойчивости, длительных п.п.)

Д) КОСМОС (расчеты УР по показаниям телеметрии, построение модели режима, утяжеление режима, оптимизация по Q, оптимизация по P для обеспечения функционирования энергорынка)

62’ Основные задачи, решаемые на основе использования ПВК РАСТР.

РАСТР (Учебно- научно-производственное предприятие "УПИ-Энерго Формат БД Metakit)

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

2. Эквивалентирование.

3. Расчет предельных режимов.

4. Оптимизация режимов по реактивной мощности, расчет анцапф трансформаторов с РПН.

63’ Основные задачи, решаемые на основе использования ПВК ДАКАР.

ДАКАР

ОАО "Институт

'Энергосетьпроект", ООО "Элекс"

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

2. Деление энергосистемы в процессе расчета режима на ряд подсистем с разными значениями частоты в каждой из них.

3. Утяжеление режимов и определение пределов мощности.

4. Расчет колебательной статической устойчивости с построе­нием областей устойчивости в координатах стабилизирующих коэффициентов и выбором настроечных параметров АРВ ге­нераторов.

5. Расчет электромеханических переходных процессов с мо­делированием действий любых устройств ПА.

6. Расчет длительных электромеханических переходных ре­жимов.

7. Расчет несимметричных, неполнофазных режимов и токов короткого замыкания.

8. Эквивалентирование режимной схемы.

9. Создание графической схемы сети и коммутационных схем подстанций.

64’ Основные задачи, решаемые на основе использования ПВК МУСТАНГ.

МУСТАНГ

ОДУ Северо-запа­да

1. Расчет установившихся режимов с возможностью учета статических характеристик нагрузки и вставок постоянного тока.

2. Утяжеление и расчет предельного режима.

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

65’ Основные задачи, решаемые на основе использования ПВК EuroStag.

EUROSTAG

"Electricite de France", "Tractebel"

1. Расчет установившихся режимов с учетом изменения частоты и без в режимах off-line и on-line.

2. Расчет колебательной устойчивости в виде экспорта линеаризованной системы уравнений в подпрограмму расчета собственных значений для проведения модального анализа.

3. Расчет динамической устойчивости (электромагнитные и электромеханические переходные процессы).

4. Расчет длительных переходных процессов.

66’ Основные задачи, решаемые на основе использования ПВК Космос.

КОСМОС

Институт электродинамики НАН Украины, Киев

1 Оперативные расчеты режимов ЭЭС на основе телеметрической информации. В состав комплекса входят программы решения следующих задач:

1.1 формирование расчетной схемы;

1.2 оценивание состояния;

1.3 построение модели режима по расширенной схеме, включающей внутренние и внешние ненаблюдаемые фрагменты;

1.4 расчет установившихся и самоустанавливающихся по частоте режимов;

1.5 утяжеление по заданным траекториям;

1.6 оптимизация по реактивной мощности.

2. Оптимизация режимов по активной мощности для обеспечения функционирования энергорынка (оптимизация для ФО-РЭМ).

Источником информации о параметрах текущего или какого-либо ретроспективного режима служат базы данных оперативного управления, создаваемые и обслуживаемые оперативно-информационными комплексами (ОИК).

67’ Градиентный метод расчета уравнений УР. Сущность, область сходимости.

Трудоемкий, медленно сходится. Но: ОЧЕНЬ большая область сходимости. Дает возможность ввести режим в область существования, когда другие методы не позволяют это сделать.

Формируется

Находится минимум :

Т.к. то

grad – движение к max. -> к min – «анти-grad»:

68’ Система уравнений, характеризующая необходимое условие минимума функции в градиентном методе расчета УР.

Формируется

Находится минимум .

Необходимое условие :

Найдя , решим исходную систему уравнений.

(От автора: В принципе, система, про которую спрашивается в вопросе – это то, что в прямоугольнике. Что про это говорить – хрен его знает!..)

69’ Понятие градиента, антиградиента функции . Математическая формулировка итерационного процесса градиентного метода расчета УР.

Находится минимум :

Решение:

grad – движение к max. -> к min – «анти-grad»:

Нахождение градиента:

Итерационный процесс: