Вычислительные аспекты

1.6 Вычислительные аспекты

Специфические особенности ЭЭС и МОНО играют решающую роль в рациональной организации вычислительного процесса.

Используемые при оценке состояния ЭЭС матрицы – матрица узловых проводимостей, матрица частных производных, матрица коэффициентов системы линейных алгебраических уравнений

 (1.23)

содержат незначительное число ненулевых элементов, т.е. являются разреженными: значительного сокращения времени счета и существенной экономии используемого объема оперативной памяти ЭВМ можно добиться, если хранить ненулевые элементы и оперировать с ними.

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

Для решения систем линейных алгебраических уравнений вида (1.23)

  (1.24)

( , )

используется метод Гаусса или его модификации. В методе Гаусса система уравнений (1.24) решается в два хода – прямой и обратный. При прямом ходе матрица коэффициентов приводится к верхней треугольной форме. Для этого к системе (1.24) с t неизвестными применяется (t -1) – шаговый процесс исключения неизвестных. В результате на (t -1) – ом шаге будет получена треугольная система:

(1.25)

Обратный ход метода Гаусса состоит в последовательном вычислении неизвестных из (1.25), начиная с последнего уравнения.

Рассмотренные преобразования удобно реализовать в матричном виде. Если обозначить матрицу коэффициентов (1.25)

(1.26)

и ввести матрицу преобразований на r – том шаге

(1.27)

то

. (1.28)

Операция обращения матрицы преобразования (1.27) равносильна инвертированию недиагональных элементов, а произведение нижних треугольных матриц дает такую же матрицу, поэтому

 (1.29)

где

 (1.30)

Выражение (1.29) – т. н. LU – разложение матрицы А в виде произведения нижней треугольной матрицы L и верхней треугольной матрицы U.

Замена z=Uh показывает, что h можно получить, решая треугольные системы:

Lz=b (1.31)

Uh=z (1.32)

Выражение (1.31) – матричная запись заключительной части прямого хода метода Гаусса (пересчета свободных членов), а (1.32) – матричная запись обратного хода. Для симметричной матрицы

где D – диагональная матрица с элементами

 i=1,2…., t,

разложение

(1.33)

называется  – разложением.

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

Одно из основных достоинств МОНО состоит в том, что гауссово исключение не требует выбора главных элементов для поддержания численной устойчивости. Это означает, что матрицу коэффициентов можно переупорядочить, не заботясь о численной устойчивости, причем до начала численного решения: выбирается такая последовательность исключения неизвестных, которая приводит к появлению минимального числа ненулевых элементов. Еще одна важная особенность такого выбора исключаемой переменной состоит в симметричном переупорядоченииматрицы коэффициентов – имеет место симметричная перестановка строк и столбцов [1].

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

Т.к. матрица коэффициентов симметрична, достаточно пересчитывать и хранить только ее верхнюю треугольную часть. Если для каждой строки имеется список столбцов с ненулевыми элементами, то он полностью определяет, в каких строках элементы каких столбцов пересчитываются. Для удобства поиска в этом списке индексы столбцов желательно располагать в порядке возрастания. Например, если на r‑ом шаге в r‑ой строке ненулевые элементы находятся в столбцах r, s, q, то пересчитываются коэффициенты в s‑ой (в столбцах s и q) и в q‑ой (в столбце q) строках.

В матрице частных производных каждому i‑му узлу соответствует два столбца 2i‑1, 2i, а в матрице коэффициентов А – блочная матрица второго порядка:

.

Измерению ветви (i, j) соответствуют четыре ненулевые блочные матрицы: A_ii,A_ij, A_ji, A_jj(i < j). .

Для каждого такого блока местоположение всех четырех элементов однозначно определяется номером строки и номером столбца блока, что равносильно указанию места установки измерительного датчика. Аналогично, номер узла и список смежных с ним узлов определяют блоки ненулевых элементов для измерения в r‑ом узле.

Элементы матрицы коэффициентов хранятся блоками по строкам. Для каждого блока ненулевых элементов номер столбца указывается в массиве «индексы столбцов». Местоположение первого блока каждой строки задается в массиве «указатель индексов строк» [1]. Блочное представление дает существенную экономию памяти как при хранении, так и при формировании системы уравнений. В действительности кодировка расстановки ТИ непосредственно определяет местоположение блоков ненулевых элементов в схеме хранения, следовательно, отпадает необходимость запоминания промежуточных результатов (матрицы частных производных).

Т.о., все необходимые предпосылки для рациональной организации вычислительного процесса гарантируются МОНО. Наиболее трудоемкая часть расчетов должна выполняться вне реального времени на подготовительном этапе. К ним относятся:

1. Формирование структуры первоначального заполнения матрицы коэффициентов. Она (структура) однозначно определяется расстановкой ТИ и топологией электрической сети; формируется с учетом всех ТИ для типовой схемы электрических соединений, в которой все объекты, оснащенные устройствами телесигнализации (ТС), считаются включенными. Текущее состояние объектов, не оснащенных устройствами ТС, отражается в исходной схеме электрических соединений. На подготовительном этапе резервируется место для всех возможных ненулевых элементов. Сформированная таким образом структура заполнения и, следовательно, схема ее хранения может использоваться при оценивании состояния ЭЭС в реальном времени с любым составом ТИ и при любых производимых в сети коммутациях, не приводящих к появлению новых узлов: отключение части ТИ и (или) ветвей отражается только на числовых значениях элементов матрицы.

2. Упорядочение – определение последовательности исключения неизвестных. Это равносильно перенумерации узлов расчетной схемы. Предпочтительней такая последовательность исключения неизвестных, которая приводит к появлению минимального числа новых ненулевых элементов. Наиболее часто используются два алгоритма динамического упорядочения. В первом из них на каждом шаге метода Гаусса исключается неизвестная, соответствующая строка которой содержит минимальное число ненулевых элементов (если таких неизвестных несколько, то выбор произволен), во втором – неизвестная, исключение которой приводит к появлению наименьшего числа новых ненулевых элементов. Оба алгоритма дают достаточно близкие результаты, но первый алгоритм динамического упорядочения предпочтительней, так как проще, требует меньших затрат времени и памяти.

3. Имитация исключения Гаусса с резервированием места под новые ненулевые элементы и формирование схемы хранения матрицы коэффициентов. Структура первоначального заполнения (верхняя треугольная часть) запоминается блоками по строкам. С учетом установленной последовательности исключения неизвестных резервируется место для новых ненулевых элементов.

4. Формирование матрицы узловых проводимостей. Для ускорения процесса поиска нужного элемента ненулевые элементы этой матрицы, несмотря на симметричность, хранятся полностью. Схема хранения аналогична рассмотренной ранее.

В реальном времени выполняются вычисления, необходимые собственно для оценивания ЭЭС:

1. Ввод текущих ТИ.

2. Формирование системы уравнений (1.24).

3. Решение системы уравнений (1.24).

4. Реализация одного шага итерационного процесса (1.19).

5. Проверка критерия окончания счета . Если условие не выполняется, перейти к п. 2.

6. Расчет потокораспределения по результатам оценивания.

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

а) хранение и обработка только верхней треугольной части матрицы коэффициентов системы уравнений (1.24);

б) блочное хранение ненулевых элементов;

в) использование кодировки расстановки ТИ, непосредственно определяющей местоположение блоков ненулевых элементов;

г) вычисление в неявном виде матрицы частных производных без запоминания промежуточных результатов;

д) использование статичной схемы хранения.

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

Похожие работы

Совершенствование систем электроснабжения подземных потребителей шахт. Расчет схемы электроснабжения ЦПП до участка и выбор фазокомпенсирующих устройств

Скачать

179075

32

127

... (от передвижения источников загрязнения) 1180,48 Всего за год: 211845,25 10. Совершенствование системы электроснабжения подземных потребителей шахты Расчет схемы электроснабжения ЦПП до участка и выбор фазокомпенсирующих устройств Основными задачами эксплуатации современных систем электроснабжения горных предприятий являются правильное определение электриче­ ...

Анализ режимов работы электрических сетей ОАО "ММК им. Ильича" и разработка адаптивной системы управления режимами электропотребления

Скачать

118979

22

26

... luc – программа используется для разложения матрицы на треугольные сомножители; rluc – программа, которая отвечает за решение системы уравнений. 4. Разработка адаптивной системы управления режимами электропотребления 4.1 Функции автоматизированной системы Сбор, накопление и передача информации, характеризующей режим электропотребления комбината (информация о нагрузках). Сбор, накопление ...

Расчет, анализ и оптимизация режимов и потерь электроэнергии в предприятии "КАТЭКэлектросеть"

Скачать

138956

15

3

... Еловка ТМН-2500/35 ±6×1,5% Ужурсовхоз ТМН-4000/35 ±6×1,5% 2. Характеристика задачи расчета, анализа и оптимизации режимов РЭС 110-35 кВ по напряжению, реактивной мощности и коэффициентам трансформации Питающие электрические сети напряжением 110 кВ, ...

Потери электроэнергии в распределительных электрических сетях

Скачать

67860

2

2

... линиям относят линии, для которых верхняя граница интервала неопределенности потерь превышает установленную норму (например, 5%). 3. Программы расчета потерь электроэнергии в распределительных электрических сетях 3.1 Необходимость расчета технических потерь электроэнергии В настоящее время во многих энергосистемах России потери в сетях растут даже при уменьшении энергопотребления. При ...