Ирина Эланс

Автор который поможет с любыми образовательными и учебными заданиями

Математические задачи электроэнергетики. 2

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

ФЕДЕРАЛЬНОЕ АГЕНСТВО ПО ОБРАЗОВАНИЮ

ИРКУТСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

Кафедра электроснабжения и электротехники

Допускаю к защите

Руководитель_____________

_________________________

И.О. Фамилия

Вариант №24

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

К курсовой работе по дисциплине

Математические задачи электроэнергетики

Выполнил: студент группы ЭП-10-1______________________

Нормконтроль _____________________Свеженцева О.В.

Иркутск 2012

СОДЕРЖАНИЕ

ВВЕДЕНИЕ 3

1. Математическое моделирование задач электроэнергетики с помощью аппарата линейной алгебры и теории графов 4

1.1 Представление синусоидального тока комплексными величинами 4

1.2 Матричная алгебра 4

1.3 Определитель матрицы и его свойства 4

1.4 Вычисление обратной матрицы 5

2. Расчет установившихся режимов электрических систем 7

2.1.Схема замещения электрической сети как связный граф 7

2.2.Первая и вторая матрицы инциденций 8

3. Методы решения линейных алгебраических уравнений 11

3.1. Решение систем линейных алгебраических уравнений методом Гаусса 11

3.2. Решение систем линейных алгебраических уравнений в MATLAB 11

4. Приближенные методы решения нелинейных алгебраических уравнений. 16

4.1.Метод деления отрезка пополам 16

4.2.Метод Ньютона………………………………………………………………………………………………..

4.3.Метод простой итерации.................................................................................

5. Применение вероятностно - статистических методов в задачах электроснабжения 19

5.1.Основные определения...................................................................................

5.2.Прогнозирование уровня электропотребления на промышленном предприятии............................................................................................................

5.3.Вычисление числовых характеристик случайных величин в системе MATLAB....................................................................................................................

6.Расчетная часть 1 задания.....................................................................................

7.Расчетная часть 2 задания......................................................................................

8.Расчетная часть 3 задания......................................................................................

9.Расчетная часть 4 задания......................................................................................

10.Расчетная часть 5 задания......................................................................................

ЗАКЛЮЧЕНИЕ 23

СПИСОК ИСПОЛЬЗУЕМОЙ ЛИТЕРАТУРЫ 24

Задание 1. Математическое моделирование задач электроэнергетики с помощью аппарата линейной алгебры и теории графов.

Краткие теоретические сведения

1.1 Представление синусоидального тока комплексными величинами

Любое вещественное число можно изобразить графически на числовой оси Ох. Понятие вещественного числа можно обобщить, если ввести в рассмотрение число z, образованное парой вещественных чисел х,у, взятых в определенном порядке. Такое число z = (x,y) называется комплексным. Вещественные числа х,у составляют соответственно вещественную и мнимую части комплексного числа г. Часто используются обозначения

x=Rez y = lmz

Комплексные числа также можно изображать графически. Это изображение будет двумерным на плоскости, образованной двумя взаимно перпендикулярными осями Ох и Оу. Комплексное число на плоскости хОу представляется точкой М{х,у); эту точку также называют изображением комплексного числа и обратно, пару чисел (х,у), образующих комплексное число z , называют аффиксом точки М .

Любое комплексное число можно представить в одной из трех форм.

Алгебраической
Тригонометрической
Показательной

z = x + jy

z = |z|cosα + j[z|sinα

z = |z|e^jα

Где |z| = - модуль комплексного числа

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

Если аргумент α явлеятся линейной функцией времени t , т.е. α =, то

И графическое представление комплексной функции аналогично представлению синусоидального тока вращающимся вектором.

Комплексная функция , у которой модуль и аргумент равны соответственно амплитуде и аргументу синусоидального тока, называется комплексным мгновенным синусоидальным током.

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

Два комплексных числа z = (x,y) и z' = (х',у') считаются равными, если совпадают изображающие их точки. Это означает, что равенство z и z' имеет место в том, и только в том случае, когда

х=х', y=у'

Т.е. другими словами два комплексных числа равны, когда равны их действительные и комплексные части.

Для алгебраической формы представления комплексных чисел справедливо; при сложении двух комплексных чисел складываются отдельно их действительные и комплексные части.

(х, у) + (х', у') = (х + х',у + у').

Умножение двух комплексных чисел следует производить как умножение двух алгебраических двучленов, приводя подобные при нулевой и первой степени числа j и помня, что j² =-1.

(х, у) * (х', у') = (хх'- уу,ху' + х'у)

Если число z = (х,у), то число z'=(х-у) называется комплексно сопряженным к числу z .

Вычитание и деление определяются как операции обратные операциям сложения и умножения, деление на 0 для комплексного числа не определено.

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

1.2 Матричная алгебра

Матрицей размера (mхn) называется прямоугольная таблица:

составленная из элементов a_ij и содержащая m строк и n столбцов.

Положение элементов в таблице определяется двойным индексом ij, первый означает номер строки, второй номер столбца на пересечении которых стоит данный элемент. Запись группы величин в виде матрицы не предусматривает каких-либо действий над ними. Это лишь одна из форм упорядоченной записи в виде условной таблицы.

Если в матрице А строки сделать столбцами, а столбцы строками, то получается транспонированная матрица А^Т.

Квадратной матрицей называется матрица, в которой число строк совпадает с числом столбцов. Если элементы в квадратной матрице располагаются симметрично относительно главной диагонали, то такая матрица называется симметричной.

Диагональной матрицей называется матрица, в которой все элементы, кроме стоящих на главной диагонали, равны 0.

Единичная матрица, это диагональная матрица, у которой на главной диагонали стоят 1.

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

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

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

Умножение матрицы на число. Для того чтобы умножить матрицу А на число , необходимо каждый элемент этой матрицы умножить на число .

Умножение матриц. Умножение матриц в алгебре матриц не коммутативно. Для того, чтобы произведение матриц существовало, необходимо чтобы число столбцов первой матрицы равнялось числу строк второй матрицы. Если матрица А имеет размерность (и хот), а матрица В размерность (mхn), то матрица С=А-В имеет размерность (nхm). В качестве элементов расположенных на пересечении i -той строки и j -го столбца матрицы произведения, принимают суммы попарных произведений, расположенных на одинаковых местах указанных строк матрицы множимого и столбцов матрицы- множителя.

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

1.3 Определитель матрицы и его свойства

Основной числовой характеристикой квадратной матрицы является ее определитель. Рассмотрим квадратную матрицу второго порядка:

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

Рассмотрим теперь квадратную матрицу третьего порядка:

Определителем третьего порядка называется число, вычисленное по следующему правилу:

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

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

Основные свойства определителей матрицы

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

Минором матрицы называется определитель, полученный вычеркиванием из квадратной матрицы одинакового числа столбцов и строк.

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

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

Таким образом

где соответствующий минор порядка.

Вычисление определителя матрицы путем разложения по элементам строки или столбца

Определитель матрицы равен сумме произведений элементов какой- либо строки (какого- либо столбца) матрицы на соответствующие алгебраические дополнения элементов этой строки (этого столбца). При вычислении определителя матрицы таким способом следует руководствоваться следующим правилом: выбирать отроку или столбец с наибольшим числом нулевых элементов. Этот прием позволяет значительно сократить объем вычислений.

1.4 Вычисление обратной матрицы

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

Квадратные матрицы одинакового порядка, произведение которых дает единичную матрицу Е, называются взаимообратными или обратными. Обозначается обратная матрица А^-1 и для нее справедливо

А*А^-1 =Е.

Вычислить обратную матрицу можно только для такой матрицы А, для которой

Классический алгоритм вычисления обратной матрицы

Записывают матрицу А^T, транспонированную к матрице А.
Заменяют каждый элемент матрицы А^Т определителем, полученным в результате вычеркивания строки и столбца, на пересечении которых расположен данный элемент.
Этот определитель сопровождают знаком плюс, если сумма индексов элемента четная, и знаком минус - в противном случае.

4. Делят полученную матрицу на определитель матрицы А.

1.5 Решение задач линейной алгебры в системе MATLAB

При создании матриц в системе MATLAB символы пробел и запятая используются для отделения элементов внутри строки в матрице, символ точка с запятой отделяет строки в матрице.

При создании матриц необходимо следить за равенством длин строк, ее образующих.

В MATLAB операция транспонирования матрицы выполняется с помощью либо оператора «.'», либо функции transpose()

Элементарными матричными преобразованиями являются:

перестановка местами двух строк матрицы,
умножение всех элементов строки матрицы на число, отличное от нуля,
прибавление ко всем элементам строки матрицы соответствующих элементов другой строки, умноженных на одно и тоже число.

Система MATLAB позволяет выполнять поэлементное умножение матриц.

При выполнении поэлементного умножения размерности матриц должны быть одинаковыми.

При поэлементном умножении матриц умножаются значения соответствующих элементов этих матриц и записываются в результирующую матрицу.

Для нахождения определителя (детерминанта) и ранга матриц в MATLAB имеются следующие функции:

det(X) — возвращает определитель квадратной матрицы X. Если X содержит только целые элементы, то результат — тоже целое число. Использование условия det(X)=0 как теста на вырожденность матрицы действительно только для матрицы малого порядка с целыми элементами.

Если является квадратной матрицей, то обратной по отношению к называется матрица, которая при умножении на (как слева, так и справа) дает единичную матрицу:

Для того чтобы квадратная матрица имела обратную матрицу необходимо и достаточно, чтобы она была невырожденной.

Обратная матрица находится с помощью функции .

Расчетная часть

1.6 Расчет индивидуального задания варианта 24

Расчетная часть должна содержать:

вычисление определителя матрицы двумя аналитическими способами;
вычисление определителя в системе MATLAB;
вычисление обратной матрицы классическим способом;
вычисление обратной матрицы по системе MATLAB;

Дано:

1.6.1 Вычисление определителя классическим способом

1.6.2 Вычисление определителя методом разложения по элементам строки и столбца

Расчет проводим по разложению 1-ого столбца.

1.6.3 Вычисление определителя в системе MATLAB

A=[5,2,3;3,1,3;4,5,1];

det(A)

ans =

-19.0000

1.6.4 Вычисление обратной матрицы классическим способом

Записываем матрицу , транспонированную к матрице .

Заменяем каждый элемент матрицы определителем, полученным в результате вычеркивания строки и столбца, на пересечении которых расположен данный элемент.

Этот определитель сопровождаем знаком плюс, если сумма индексов элемента четная, и знаком минус – в противном случае.

Делим полученную матрицу на определитель матрицы .

1.6.5 Вычисление обратной матрицы в системе MATLAB

A=[5,2,3;3,1,3;4,5,1];

inv(A)

ans =

0.7368 -0.6842 -0.1579

-0.4737 0.3684 0.3158

-0.5789 0.8947 0.0526

Задание 2. Расчет установившихся режимов электрических систем

Краткие теоретические сведения

2.1 Схема замещения электрической сети как связный граф

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

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

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

Ветвью называется участок электрической цепи, в которой в любой момент времени ток имеет одно и то же значение.

Узлом называется место соединения двух или большего числа ветвей. Одна из ветвей, соединяющихся в узле, может быть источником тока.

Контуром называется любой замкнутый путь, проходящий но нескольким ветвям. Если схема электрической цепи не содержит контуров, то она называется разомкнутой.

Любая электрическая схема состоит из некоторого числа элементов: линий электропередач, трансформаторов, источников питания, потребителей электрической энергии и т.д. Для проведения расчетов электрическую схему в начале представляют схемой замещения, а затем переходят к направленному графу электрической сети.

Основные определения теории графов

Рис.1 Граф электрической сети

Граф можно представить, если представить множество точек на плоскости X={1,2,3,4,5}, называемых вершинами графа и множества направленных отрезков U = { u₁,u₂,u₃,u₄,u₅,u₆), соединяющих все или несколько вершин и называемых дугами. Таким образом любой граф G можно определить как пару множеств X, UG = (X, U).

Путем в графе G называется такая последовательность дуг = (u_l,u₂---u_k), в которой конец каждой предыдущей дуги совпадает с началом следующей. Длиной пути называется число равное числу дуг на этом пути.

Иногда граф рассматривают без учета ориентации его дуг, в этом случае граф называют неориентированным. У неориентированного графа понятие дуга заменяется на понятие ребро. Цепью называют последовательность ребер.

Говорят, что граф связен, если любые две его вершины можно соединить цепью. Циклом называется конечная цепь, у которой начальная и конечная вершины совпадают.

Деревом называется конечный, связный, неориентированный граф, не имеющий циклов. Если в дереве n вершин, в нем всегда (n-1) ребро.

Все элементы схем замещения делятся на активные и пассивные.

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

Пассивные элементы схем замещения: сопротивления и проводимости создают пути для протекания электрического тока. Пассивные элементы обычно разделяются на поперечные и продольные.

Поперечные пассивные элементы - это ветви, включенные между узлами схемы и нейтралью.

К продольным пассивным элементам относятся ветви, соединяющие все узлы, кроме узла с напряжением равным нулю.

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

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

Источники электрической энергии могут быть представлены в виде источника напряжения с ЭДС и внутренним сопротивлением , либо в виде источника задающего тока с внутренним сопротивлением равным бесконечности.

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

В практике расчета электрических сетей используются схемы замещения с источниками тока.

В схеме замещения электрической системы трансформаторы подстанций и электростанций, ЛЭП, реакторы и другие элементы электрической сети представляются в виде сопротивлений.

Для двухобмоточных силовых трансформаторов, автотрансформаторов обычно применяется Г-образная схема замещения, параметры которой определяются по каталожным данным соответствующего трансформатора. Параметры поперечной ветви в Г-образной схеме замещения моделируют потери мощности в стали.

2.2 Первая и вторая матрицы инциденций

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

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

Здесь {u₁,u₂,…u_m} ветви соответствующего графа схемы замещения, {1,2,…n} узлы соответствующего графа схемы замещения. Для графа представленного на рис.1 первая матрица инциденций будет выглядеть следующим образом:

Условием правильности заполнения матрицы является наличие в каждом столбце только одной 1 со знаком «+» и только одной 1 со знаком «-». Это обусловлено тем, что любая ветвь всегда исходит из одного узла и всегда входит в один узел. Таким образом, сумма всех элементов в первой матрице инциденций равна нулю.

Если в матрице выделить строку, соответствующую некоторому узлу, принятому за балансирующий узел, то матрица без последней строки называется матрицей соединений без балансирующего узла и обозначается М.

Вторая матрица инциденций называется также матрицей контуров и обозначается N. Она связывает ветви и независимые контуры соответствующего графа схемы замещения. Для составления матрицы N нужно определить число независимых контуров схемы. Это число определяется по формуле к=т-n+1, где к - число независимых контуров, т - число ветвей, n- число узлов.

Строки матрицы N соответствуют независимым контурам схемы, столбцы ветвям. Элементы матрицы N определяются по следующим правилам:

2.3 Матричная форма записи уравнений состояния электрической сети

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

Исходными данными о нагрузках потребителей обычно служат значения мощностей S_н1=P_н1+ jQ_н1, которые принимаются постоянными, или как функции напряжения узлов. Исходными данными об источниках питания служат постоянные значения активных и реактивных мощностей или значения активной мощности и напряжения в точках их подключения.

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

Обобщенное уравнение состояния

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

A*I = F (1)

Матричная форма записи уравнения, где А матрица параметров схемы замещения, где I - вектор- столбец токов в ветвях, т - число ветвей в схеме замещения, F- вектор-столбец исходных параметров режима.

Уравнение (1) объединяет два матричных уравнения.

Математические задачи электроэнергетики. 2