проведен корректно, минимальное значение найдено1. Иначе, следует уточнить начальное приближение (см. предыдущий пункт, пп е). Т.к. вектор x состоит из двух векторов – Pg и U, нужно составить новый вектор x0 из векторов Pg0 и U0.
9.Вычислите сумму потерь мощности в электрической сети.
10.Проигравшие студенты по найденным узловым напряжениям вычисляют токи в ветвях электрической сети.
11.Ввод ограничений пропускной способности сети. Чтобы отыграться, проигравшие студенты выбирают ветвь электрической сети и устанавливают ограничение на величину тока в ней так, чтобы на
следующем такте помешать победителю выиграть. Это линейное ограничение следует выразить через узловые напряжения2 и внести в программу в раздел ограничений.
12.Повторите п. 5. Проведите расчет. Проверьте и запишите результаты расчета.
5.Вопросы к защите
1.Объясните результаты расчета. Почему вы проиграли или выиграли?
2.В чем сложность расчета узловых уравнений в терминах мощностей?
3.Как ограничения пропускной способности сети повлияли на результаты расчета?
4.Имеются ли в схеме особые разрезы?
Литература
1.Демирчян К. С., Нейман Л. Р., Коровкин Н. В. Теоретические основы электротехники: Учебник для вузов. 5-у изд. Т. 2. – СПб.: Питер, 2009.
2.Информационные технологии электротехники: учебное пособие / под ред. П.А. Бутырина. – М.: Издательский дом МЭИ, 2007.
3.Сборник задач по теоретическим основам электротехники. Под ред. Чл.-корр. РАН П.А. Бутырина. – М.: Издательский дом МЭИ, Том 2. – 2012.
4.В.Е.Шмелев "Partial Differential Equations Toolbox. Инструментарий решения дифференциальных уравнений в частных производных"
(http://MATLAB.exponenta.ru/pde/book1/index.php)
5.В.Е.Шмелев "Примеры применения PDE Toolbox" (http://MATLAB.exponenta.ru/pde/book5/index.php)
6.И.В.Черных. "SimPowerSystems: Моделирование электротехнических устройств и систем в Simulink" (http://matlab.exponenta.ru/simpower/book1/index.php)
1Если оно равно нулю, это означает, что достигнуто максимальное значение итераций. Его можно увеличить в параметре 'MaxFunEvals'. Такое решение также можно использовать, если оно проходит проверку.
2Учтите, что в MATLAB отсчет индексов в векторах идет с 1, в Mathcad c 0.
81