Пример использования метода Ньютона для решения УУН




УЧЕТ КОМПЛЕКСНОГО ХАРАКТЕРА

ПАРАМЕТРОВ СХЕМЫЗАМЕЩЕНИЯ И РЕЖИМА [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.


Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2016-04-27 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: