Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теоретические основы электротехники. Ч.2. Анализ нелинейных магнитных цепей и расчёт переходных процессов в линейных электрических цепях. Учебное посо

.pdf
Скачиваний:
0
Добавлен:
11.08.2026
Размер:
550 Кб
Скачать

ЗАКЛЮЧЕНИЕ

Учебное пособие направлено на освоение численных методов расчёта переходных процессов с применением современного программного обеспечения. В частности рассматриваются две среды – Mathcad и R.

Сегодня среда Mathcad является достаточно известной и распространённой. Большинство пользователей предпочитает версию Mathcad 15, которой в данный момент без малого 10 лет. Как альтернатива этой среде предлагается язык программирования R. В примерах видно, что между программами в Mathcad и программами в R есть определённое сходство. Таким образом, освоение R не представляется достаточно сложной задачей, а его интенсивное развитие и свободное распространение делают его наиболее перспективным для изучения.

В пособии рассмотрены численные методы решения. Здесь, как и впервой части пособия, поднимается вопрос эффективности методов с точки зрения применения математического ПО. Как правило, наиболее распространённым методом расчёта переходных процессов является операторный метод, так как он позволяет уйти от дифференциальных уравнений к алгебраическим, что значительно облегчает расчёт. Однако его использование совместно с математическим ПО ограниченно из-за необходимости символьных вычислений. По этой причине в задании к работе № 3 рассматривается метод переменных состояния, который не только позволяет рассчитывать цепи с нелинейными элементами, но и позволяет эффективно использовать математическое ПО.

31

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1.Теоретические основы электротехники: [в 3 т.]: учеб. для вузов: Т.1 / К.С. Демирчан [и др.]. – 4- е изд., доп. для самостоят. изучения курса; Гриф МО.

СПб.: Питер, 2006. – 462 с.

2.Теоретические основы электротехники: [в 3 т.]: учеб. для вузов: Т.2 / К.С. Демирчан [и др.]. – 4- е изд., доп. для самостоят. изучения курса; Гриф МО.

СПб.: Питер, 2006. – 575 с.

3.Бессонов, Л.А. Теоретические основы электротехники. Электрические цепи: учеб. для вузов / Л.А. Бессонов. – 11- е изд., испр. и доп. – М.: Гардарики,

2007. – 701 с.

4.Сборник задач по теоретическим основам электротехники: учеб. пособие для вузов / Л.А. Бессонов [и др.]; под ред. Л.А. Бессонова. – 4- е изд., перераб. и испр.; Гриф МО.– М.: Высшая школа, 2003. – 528 с.

5.Бычков, Ю.А. Основы теории электрических цепей: учеб. для вузов / Ю.А. Бычков, В.М. Золотницкий, Э.П. Чернышев. – 3- е изд., стер. – СПб.: Лань,

2004. – 464 с.

6. Шакурский, В.К. Теоретические основы электротехники. В 2 ч. Ч. 2. Аналитические и численные методы анализа переходных процессов в линейных, нелинейных и параметрических электрических цепях: учебное пособие / В.К. Шакурский. – Тольятти: Тольят. гос. ун-т, 2013. – 211 с.

32

ПРИЛОЖЕНИЯ

Приложение 1

Построение простых графиков в среде R

Команда plot()

Общий вид команды plot():

plot(x, y, …)

Синтаксис: plot(x, y)

где x – массив значений аргумента; y – массив соответствующих значений функции; … – набор дополнительных команд.

Пример:

Определим зависимость и построим график:

x=1:100

y=x^2

plot(x,y)

В результате получим следующий график (рис. П.1.1). На рис. П.1.1 дан пример использования функции plot().

Рис. П.1.1

Дополнительные команды включают в себя:

type – переменная определяющая тип линии. Следующие значения типов линий доступны:

“p” – точки (points);

“l” – линия (line);

“b” – совместно прерывистая линия и точки (both);

33

“c” – линия от “b” с исключением точек (прерывистая линия);

“o” – непрерывная линия с точками (overplotted);

“h” – гистограмма (histogram);

“s” – ступени (steps);

“S” – ступени c ранним переходом (steps);

“n” – без начертания (no plotting).

Синтаксис: plot(x, y, type = “o”)

Пример На рис. П.1.2 показан пример изменения типа линии в команде

plot().

x=1:100

y=x^2

plot(x,y, type="l")

Рис. П.1.2

main – основной заголовок графика. sub – нижняя подпись под графиком. xlab – наименование оси абсцисс. ylab – наименование оси ординат.

asp – задание пропорции координат y/x.

Синтаксис: plot(x, y, main= “Graph1”, sub=“Example”, xlab= “x” , ylab=

“y”, asp = 1)

Пример. На рис. П.1.3 приведен пример использования подписей графика.

x=1:100

y=x^2

34

plot(x,y, main="Квадратичная зависимость", sub="График функции y=x^2",

xlab="x", ylab="x^2", type="l")

Рис. П.1.3

xlim – пределы шкалы абсцисс. ylim – пределы шкалы ординат.

Синтаксис: plot(x, y, xlim=c(0, 15), ylim=c(0, 5))

Пример. На рис. П.1.4 приведен пример выделения области графика заданием пределов координатных осей.

x=1:100

y=x^2

plot(x,y, xlim=c(40, 60), ylim=c(1000,4000), type="l")

Рис. П.1.4

35

axes – булева переменная – контролирует отображение осей (значения

TRUE или FALSE).

ann – булева переменная – контролирует отображение наименования осей

(значения TRUE или FALSE).

Синтаксис: plot(x, y, xlab= “x”, ylab= “y”, axes = FALSE, ann = TRUE)

Пример. На рис. П.1.5 приведен пример отключенного отображения осей координат.

x=1:100

y=x^2

plot(x,y, xlim=c(60,80), ylim=c(1000,4000), type="l")

Рис. П.1.5

log – переменная переводящая одну или две оси в логарифмическую систему. Параметры:

“x” – переводит ось абсцисс логарифмической

“y” – делает ось ординат логарифмической

“xy” – делает обе оси абсцисс и ординат логарифмическими. Пример. На рис. П.1.6 приведен пример использования логарифмической

шкалы координат.

x=1:100

y=x^2

plot(x,y,log="x", type="l")

36

Рис. П.1.6

pch – переменная задаёт тип маркера. При этом тип линии должен подразумевать использование маркера (например, “b” или “o”) Количество возможных маркеров огромно. Рассмотрим основные:

Стандартные маркеры включают следующие значения.

Синтаксис: plot(x, y, pch =2)

Пример. Пример изменения типа маркера приведен на рис. П.1.7.

x=1:100

y=x^2

plot(x,y, type="p", xlim=c(40,60), ylim=c(1000,4000), pch=2)

37

Рис. П.1.7

Печатные символы.

Синтаксис: plot(x, y, pch = “p”)

Пример. Пример использования символов в качестве маркеров рис. П.1.8.

x=1:100

y=x^2

plot(x,y, type="p", xlim=c(40,60), ylim=c(1000,4000), pch="p")

Рис. П.1.8

Расширенный набор символов. При формировании графика можно задать другой тип шрифта (переменная font) и использовать символы другого шрифта. Например, при значении шрифта font = 5 можно использовать следующие символы:

38

Чтобы узнать необходимый численный код маркера, необходимо y- координату из приведенного выше рисунка умножить на 10 и к получившемуся значению прибавить соответствующую x-координату. Например, маркеру в виде «сердечка» соответствует код 169 (= 16 х 10 + 9)

Синтаксис: plot(x, y, pch = 169, font = 5)

Пример. Использование расширенного набора символов в качестве маркеров показано на рис. П.1.9.

x=1:100

y=x^2

plot(x,y, type="p", xlim=c(40,60), ylim=c(1000,4000), pch = 169, font = 5)

Рис. П.1.9

39

cex – размер маркера. По умолчанию cex = 1.

col – данный аргумент позволяет настраивать цвет компонент графика. Цвет может быть задан своим наименованием (например, “red”, “green”, “blue” и т.д.) или численным кодом. Перечислим основные варианты окраса компонент графика:

col = “red” – окрашивает график и маркеры в красный цвет;

col.main = “blue” – задаёт синий цвет заголовка графика;

col.lab=”green” – задаёт зелёный цвет наименования осей графика;

bg = “yellow” – для стандартных маркеров 21–25 задаёт цвет заливки.

lwd – ширина линии (по умолчанию lwd =1).

lty – тип линии. Может принимать следующие значения:

1 – сплошная линия;

2 – пунктирная линия;

3 – короткий пунктир;

4 – штрих-пунктирная линия;

5 – длинный пунктир;

6 – штрих-пунктир с удлинённым штрихом.

bty – тип рамки графика. Может принимать следующие значения:

“O” – охватывает график вокруг;

“L” – охватывает график слева и снизу;

“7” – охватывает график справа и сверху;

“C” – охватывает график слева, снизу и сверху;

“U” – охватывает график слева, снизу и справа;

“[” – охватывает график слева, снизу и сверху.

Команды lines() и points()

Данные команды позволяют добавить на уже созданный график, с помощью команды plot(), дополнительных зависимостей в виде линий (команда lines) или точек (команда points). Аргументы данных команд аналогичны аргументам команды plot().

Пример. Пример начертания нескольких зависимостей на одном графике приведен на рис. П.1.10.

x=1:100

y=x^2; y2=x^3; y3=x^4 plot(x,y, type="l") lines(x,y2,lwd=2) points(x,y3, type="b")

40

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]