- •1. Матричная алгебра и теория графов в электроэнергетике
- •1.1. Матричная алгебра
- •1.1.1. Классификация матриц
- •1.1.2. Определитель
- •1.1.3. Миноры и алгебраические дополнения
- •1.1.4. Действия с матрицами
- •1.1.5. Транспонированная матрица
- •1.1.6. Обратная матрица
- •1.1.7. Алгоритм вычисления обратной матрицы
- •Основные свойства обратной матрицы
- •Возведение матрицы в степень
- •1.1.8. Нормы матрицы
- •1.1.9. Ранг матрицы
- •1.2. Теория графов в электроэнергетике
- •1.2.1. Некоторые сведения об электрических системах
- •1.2.2. Геометрический образ электрической сети
- •1.2.3. Уравнения законов Ома и Кирхгофа в матричной форме
- •1.2.4. “Прямой ” расчет токораспределения в электрической сети
- •1.2.5. Метод узловых напряжений
- •1.2.6. Метод контурных токов
- •Контрольные вопросы к разделу 1
- •2. Решение систем линейных и нелинейных алгебраических уравнений
- •2.1. Методы решения систем линейных уравнений
- •2.1.1. Метод обратной матрицы
- •2.1.2. Метод определителей
- •2.1.3. Метод Гаусса
- •2.1.4. Вычисление определителя методом Гаусса
- •2.1.5. Вычисление обратной матрицы методом Гаусса
- •2.1.6. Трангуляция матрицы
- •2.1.7. Метод Жордана-Гаусса
- •2.1.8. Метод простой итерации
- •2.1.9. Метод Зейделя
- •2.2. Методы решения систем нелинейных уравнений
- •2.2.1. Понятие о системах нелинейных уравнений и методах их решения
- •Решение нелинейного уравнения методом Ньютона
- •Решение системы нелинейных уравнений методом Ньютона
- •Контрольные вопросы к разделу 2
- •3. Расчет установившихся режимов электрических сетей
- •3.1. Инженерные методы расчета установившихся режимов электрических сетей
- •3.1.1. Расчет установившегося режима разомкнутой электрической сети методом «в два этапа»
- •3.1.2. Расчет установившегося режима кольцевой электрической сети методом в два этапа
- •3.2. Матричные методы расчета установившихся режимов электрических сетей
- •Линейные уун в формуле баланса токов в узлах.
- •3.2.1. Линейные уун в формуле баланса токов в узлах
- •3.2.2. Определение параметров установившегося режима электрической сети по известным значениям
- •3.2.3. Учет поперечной проводимости при использовании уун
- •3.2.4. Нелинейные уун в форме балансов токов в узлах
- •3.2.5. Нелинейные уун в форме балансов мощностей в узлах
- •Контрольные вопросы к разделу 3
- •4. Математическое программирование в электроэнергетике
- •4.1. Математическая модель задачи оптимизации
- •4.2. Характеристика задач оптимизации
- •4.3. Оптимальное распределение потоков мощности в электрической сети
- •4.4. Градиентный метод
- •4.5. Метод наискорейшего спуска
- •4.6. Учет ограничений в форме неравенств
- •Контрольные вопросы к разделу 4
- •Метод наискорейшего спуска.
- •Список литературы
4.4. Градиентный метод
Для поиска экстремума непрерывной дифференциальной нелинейной функции широкое применение имеет градиентный метод. Существуют различные его модификации. Рассмотрим суть градиентного метода.
Г
радиент
– это вектор, указывающий направление
в пространстве переменных, в котором
функция возрастает быстрее всего. Пусть
в точке 1 функция имеет значение F1.
Предположим, что мы нашли градиент grad
F1,
в котором функция возрастает быстрее
всего, значит в обратном направлении
она быстрее всего убывает. Пройдем
немного в обратном направлении и придем
во вторую точку, в которой F3
> F1
. Если затем во второй точке определить
новый градиент и опять сделать шаг в
обратном направлении, то пойдем в третью
точку, в которой F3
> F2.
Двигаясь таким образом мы можем попасть
в минимум целевой функции.
Итак, под градиентными методами понимают методы, в которых направления движения к точке оптимума функции определяются направлением градиента этой функции. Использование градиентных методов для решения задач нелинейного программирования в общем случае позволяет найти точку локального экстремума. Глобальный минимум можно найти для выпускных функций и выпуклых ОДР.
Рассмотрим
простейший случай, когда имеет функцию
от двух переменных:
– это выпуклая поверхность, которая
имеет минимальное значение
.
Пересечем ее двумя плоскостями F1 = C1 и F2 = C2 и получим в пересечении замкнутые линии равного значения целевой функции, они называются линиями одного уравнения или изолиниями. Спроектируем их на плоскость X10X2. При любых значениях х1 и х2 целевая функция на изолинии имеет одно и то же значение.
Г
радиент
всегда изолинии в точке его определения.
grad
F
= G
- вектор,
компоненты которого:
;
;
←
интерпространство
При n > 2 мы имеем многомерное пространство. Для иллюстрации будем представлять двумерное пространство.
Зададимся
начальными значениями
и
найдем значение
и
в точке, определяемой этими координатами.
Затем нам нужно двигаться по прямой в
обратном направлении
.
Для этого воспользуемся параметрическими
условиями прямой, проходящей через
начальную точку
,
т.е. определять координаты новой точки
I
и
,
которая находится ближе к F*.
Из высшей математики параметрические уравнения прямой:
,
где а и b – компоненты вектора направления,
t - параметр или шаг движения по прямой.
Вместо а и b подставляем компоненты вектора направления G c обратным знаком и получаем:
.
При t = 0 имеем начальную точку, при t > 0 приходим в точку I с координатами и .
Важнейшая задача, которую нужно решать - это выбор шага t. Если выберем маленький шаг, то будем медленно приближаться к минимуму, не рискуя проскочить его.
1
-2-3-4-5-6
– траектория движения к оптимуму в
случае малых шагов. Если функция вытянута
и много переменных, то будим делать
много шагов, и процесс будет медленно
сходиться.
Надо применять более “смелую“ стратегию, которой обладает метод наискорейшего спуска.
