УЧЕТ КОМПЛЕКСНОГО ХАРАКТЕРА
ПАРАМЕТРОВ СХЕМЫЗАМЕЩЕНИЯ И РЕЖИМА [1]
Система уравнений узловых напряжений для цепи переменного тока:
.
| (9.7) |
При решении на ЭВМ системы уравнений узловых напряжений для сети переменного тока, как правило, она приводится к системе действительных уравнений порядка
, где
- число узлов схемы. Для этого представляют матрицы и вектор - столбцы с комплексными элементами в виде сумм матриц и вектор - столбцов с действительными элементами (при этом надо в виде такой суммы представить каждый комплексный элемент и учесть правило сложения матриц):
| (9.8) |
Подставляя (9.8) в (9.7), получим:
| (9.9) |
Уравнение (9.9) переписываем, разделяя действительные и мнимые слагаемые.
.
| (9.10) |
.
| (9.11) |
Иными словами, систему уравнений узловых напряжений для цепи переменного тока можно записать в виде блочного матричного уравнения:
| (9.12) |
Выражение (9.12) является системой действительных уравнений порядка
и содержит
неизвестных действительных и мнимых составляющих узловых напряжений, представленных в форме действительных чисел.
ОПИСАНИЕ РАСЧЕТА УР С ПОМОЩЬЮ ВЕКТОР-ФУНКЦИИ
Нелинейные уравнения установившегося режима в общей форме можно записать в виде системы неявных функций:
| (9.53) |
где
- вектор-функция.
Эти уравнения связывают между собой параметры установившегося режима электрической системы. Часть параметров режима задана (независимые переменные
). Обозначим вектор-столбец независимых переменных при расчете установившегося режима
. Остальные (зависимые переменные
) могут быть найдены из уравнений установившегося режима. Обозначим вектор-столбец зависимых переменных
.
Число зависимых переменных
равно числу уравнений установившегося режима. Это означает, что вектор-функция
и вектор-столбец
имеют одинаковую размерность. В зависимости от постановки задачи и способов задания исходных данных в состав векторов независимых и зависимых переменных
и
могут входить разные параметры режима.
Разделение параметров режима на зависимые и независимые переменные играет важную роль при оптимизации режимов, при определении предельных по статической апериодической устойчивости режимов и при исследовании существования и единственности решения уравнений установившегося режима.
При расчетах установившегося режима вектор независимых переменных задан, то есть
, следовательно, нелинейную систему уравнений (9.53) можно переписать как
| (9.54) |
Число уравнений в этой системе также равно числу зависимых переменных
, то есть равно размерности вектора
. В результате решения уравнений УР можно найти все зависимые переменные
.
МЕТОД НЬЮТОНА
Метод Ньютона для решения систем нелинейных уравнений представляет собой обобщение на многомерный случай метода касательных, применяемого для решения одного нелинейного уравнения.
Идея метода Ньютона состоит в последовательной замене на каждой итерации системы нелинейных уравнений некоторой линейной системой, решение которой дает значения неизвестных, более близкие к решению нелинейной системы, чем исходное приближение. Поясним идею этого метода на примере решения уравнения
| (5.4) | |
Решение уравнения - точка, в которой кривая проходит через нуль (рисунок 5.1):
Рисунок 5.1- Графическая иллюстрация метода Ньютона
|
| |
Зададим начальное приближение к решению уравнения
и вычислим значение функции
. Если точка
достаточно близка к решению, то для его получения целесообразно разложить функцию
в ряд Тейлора в окрестностях точки
:
| (5.5) |
Выражению (5.5) соответствует касательная к функции, проведенная в точке
. Такая касательная показана на рисунке 5.1. Для малых значений
приближение (5.5) хорошо моделирует саму функцию
, поэтому в качестве приближенного решения исходного уравнения целесообразно использовать решение
линеаризованного уравнения (5.5), равное
| (5.6) |
Полученную точку
можно использовать в качестве нового приближения и сделать шаг, в результате которого будет найдено приближение
и т. д. Итерационный процесс получения решения показан на рисунке 5.1.
Аналогично определяется решение для системы нелинейных уравнений. Рекуррентное выражение, представленное в матричной форме записи, имеет вид:
| (5.7) |
где
- матрица Якоби (или якобиан), составленная из частных производных;
- вектор невязок, вычисленный в точке
;
- вектор поправок к
приближению
.
Пример использования метода Ньютона для решения УУН
Для электрической сети, представленной на рисунке 5.2, определить напряжения в узлах, используя метод Ньютона (три итерации).
Для рассматриваемой схемы электрической сети может быть записана система нелинейных УУН в форме баланса токов:
| (5.8) |
| |
| Рисунок 5.2 | Схема электрической сети. |
Для решения методом Ньютона система УУН (5.8) представляется в форме баланса мощностей
| (5.9) |
И приводится к виду
|
Рекуррентное выражение метода Ньютона:
,
| (5.10) |
Где: 1) элементы матрицы Якоби вычисляются по формулам:
,
|
2) вектор невязок вычисляется в точке
по следующим выражениям:
,
|
3)
- вектор поправок к
-му приближению
.
Новые напряжения вычисляются по выражению
| (5.11) |
Для схемы электрической сети, представленной на рисунке 5.2, исходная система (5.9) имеет вид:
|
Итерация 1
Начальное приближение:
кВ
Вектор невязок записывается: 
Элементы матрицы Якоби:
|
Для заданного начального приближения
кВ элементы матрицы Якоби приобретают значение
|
Подставляем все найденные величины в
и получаем систему двух линейных уравнений:
|
Решив ее, находим:
|
По выражению (5.11):
кВ
|
Итерация 2
Вектор невязок:
|
Матрица Якоби:
|
Система двух линейных уравнений:
|
Решив ее, находим:
|
По выражению (5.11):
кВ
|
Итерация 3
Вектор невязок:
|
Матрица Якоби:
|
Система двух линейных уравнений:
|
Решив ее, находим:
|
По выражению (5.11):
кВ
|
Вектор невязок для 4-ой итерации составил бы:
|
| Алгоритмизация задач энергетики Екатеринбург, 1999. |
.
.
.
- точка, в которой кривая
проходит через нуль (рисунок 5.1):
Рисунок 5.1- Графическая иллюстрация метода Ньютона
,
,
,
кВ
кВ
кВ