Использование математического процессора MathCad в программировании

Модель цепи в пространстве состояний, специфика её построения. Аналитическое решение систем линейных дифференциальных уравнений, использование матричной экспоненты в MathCad. Реакция цепи на периодическую последовательность прямоугольных импульсов.

Рубрика Программирование, компьютеры и кибернетика
Вид учебное пособие
Язык русский
Дата добавления 10.12.2014
Размер файла 1,6 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

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

Указание. Для построения рисунка можно пользоваться любым редактором растровых изображений, например Microsoft Paint. Для получения приемлемого качества изображения при небольшом расходе памяти можно рекомендовать к использованию 16-ти уровневый монохромный формат bmp. матричный цепь импульс mathcad

2. Построение системы уравнений Кирхгофа

Составить систему уравнений Кирхгофа, записав ее в текстовом регионе рабочего листа MathCad. Уравнения следует снабдить комментариями, поясняющими, к какому узлу или контуру данное уравнение относится.

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

3. Формирование регионов с определениями основных функций

3.1. Формирование функций расчета матричных элементов модели цепи в пространстве состояний (матриц A, b и K).

3.2. Формирование функций расчета матричной экспоненты.

3.3. Формирование функций расчета реакции цепи на включение единичного источника ЭДС.

3.4. Формирование функций расчета реакции цепи на одиночный прямоугольный импульс.

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

3.6. Формирование функций расчета установившегося процесса в цепи при воздействии периодической последовательности прямоугольных импульсов.

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

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

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

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

4.1. Задать значения параметров для реактивных элементов
L = 25 мГн, С = 1 мкФ.

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

4.3. Варьируя численные значения сопротивлений R1, R2,… добиться того, чтобы собственные числа матрицы A стали комплексно - сопряженными. При этом переходный процесс в цепи станет колебательным с показателем в пределах 0.1 - 0.2.

4.4. Построить графики изменения во времени составляющих матричной экспоненты Ф(t).

4.5. Продолжая вариацию величин сопротивлений, добиться того, чтобы на каком - либо из графиков, например, для элемента Ф(t)0,0 , просматривалось 2 - 3 периода колебаний.

4.6. Определить численное значение периода собственных колебаний цепи T.

5. Исследование отклика цепи на включение источника ЭДС единичной амплитуды.

5.1. Построить графики переменных состояния iL(t), Uc(t) и переменных наблюдения.

Указание. Здесь и далее при построении графиков в качестве t использовать ранжированную переменную, изменяющуюся в необходимых пределах с шагом 10-6 сек.

6. Исследование отклика цепи на прямоугольный импульс

6.1. Построить графики переменных состояния iL(t), Uc(t) и переменных наблюдения при воздействии прямоугольного импульса единичной амплитуды при длительностях импульса
ф = 0.1T, ф = 0.5T и ф = T. Здесь и далее Т - период собственных колебаний цепи.

Указания:

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

· диапазон изменения времени на графиках принять равным 5Т.

7. Исследование динамики цепи при включении внешнего воздействия воздействии в виде периодической последовательности импульсов единичной амплитуды скважности
г = 0.5 с периодом следования Tp = T.

7.1. Построить графики изменений переменных состояния iL(t) или Uc(t).

7.2. Построить графики наблюдаемых переменных

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

Указание. На графиках, получаемых в этом пункте необходимо отразить импульсную последовательность источника ЭДС.

8. Исследование установившегося процесса в цепи при воздействии периодической последовательности импульсов

8.1. Построить графики переменных состояния iL(t), Uc(t) в установившемся режиме при воздействии периодической последовательности прямоугольных импульсов единичной амплитуды с периодом Т со скважностями г = 0.25, г = 0.5, и г = 0.75

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

8.3. Для установившегося процесса Uc(t) при г=0.5 и Tp=T по графику определить максимальное и минимальное значение на периоде.(Воспользоваться средством трассировки графика)

Указание. Графики процессов в установившемся режиме следует строить на временном интервале [0, T] с использованием в качестве t ранжированной переменной с шагом 10-6 сек.

Масштабы всех графиков осям ординат (осям Y) должны быть согласованы с фактическим диапазоном изменения изображаемых процессов с тем, чтобы обеспечить наиболее детальное представление воспроизводимых на них колебаний.

9. Оформить пояснительную записку в виде файла MathCad с комментариями (см. «Приложение В»).

Указание. Если разработка проекта производилась с использованием MathCad версии более старшей, чем MathCad 11, то перед отсылкой на рецензию его следует перезаписать как файл в формате MathCad 11

Приложение А

Вариант

схема

наблюдаемые

переменные

Вариант

схема

наблюдаемые

переменные

Y0

Y1

Y0

Y1

1

1

U1

iC

26

7

iC

U2

2

1

UL

i4

27

7

UL

U1

3

1

U3

iE

28

7

i3

iL

4

1

U2

i4

29

8

U1

i2

5

2

U4

UL

30

8

U2

iC

6

2

U3

i5

31

8

U3

U4

7

2

U2

iC

32

8

UL

i3

8

2

U5

i3

33

9

U1

iC

9

3

iC

U3

34

9

U2

i4

10

3

U2

UL

35

9

U3

i1

11

3

iE

U4

36

9

U5

UL

12

3

U1

R2

37

10

U2

iC

13

4

i1

U3

38

10

U3

iE

14

4

iC

U5

39

10

U1

UL

15

4

UL

U2

40

10

U4

i3

16

4

U4

i5

41

11

U1

iC

17

5

i1

U4

42

11

U2

i5

18

5

iC

U3

43

11

U3

UL

19

5

UL

U2

44

11

U4

i1

20

5

U5

i3

45

12

U1

i3

21

6

UL

U4

46

12

U2

i1

22

6

U2

i3

47

12

U3

iC

23

6

iC

U1

48

12

U4

UL

24

6

U3

U4

49

13

U1

i2

25

7

U1

i2

50

13

U2

iE

51

13

U3

iC

76

19

U4

U5

52

13

U4

UL

77

20

U1

UL

53

14

U1

i3

78

20

U2

iC

54

14

U2

iC

79

20

U3

i1

55

14

U3

i2

80

20

U4

i2

56

14

U4

UL

81

21

U1

i3

57

15

U1

iE

82

21

U2

iC

58

15

U2

iC

83

21

U3

i5

59

15

U3

i4

84

21

U4

UL

60

15

U4

UL

85

22

U1

i3

61

16

U1

i3

86

22

U2

i4

62

16

U2

i4

87

22

U3

iC

63

16

U3

i5

88

22

U4

i5

64

16

U4

UL

89

23

U1

i2

65

17

UL

i2

90

23

U2

iC

66

17

U1

i3

91

23

U3

i5

67

17

U2

iC

92

23

U4

UL

68

17

U3

iE

93

24

U1

iC

69

18

UL

i2

94

24

U2

i5

70

18

U1

iC

95

24

U3

UL

71

18

U2

i3

96

24

U4

i1

72

18

U3

i4

97

25

U1

iC

73

19

U1

i5

98

25

U2

iE

74

19

U2

UL

99

25

U3

i1

75

19

U3

iC

100

25

UL

i2

Приложение Б

Приложение В

Резонансная цепь с периодическим воздействием

Курсовой проект по дисциплине "Программирование"

Вариант №78

Выполнил Петров Семен Васильевич

г.Нерюнгри

Код: nr001psv

Пароль 12345678

1. Схема цепи с разметкой

Наблюдаемые переменные: UR4 и iC.

2. Система уравнений Кирхгофа

2.1.Система уравнений в исходном виде:

1) iE - iL - i2 = 0 для узла 1

2) i2 - iС - i3 = 0 для узла 2

3) i3 + iL - i4 = 0 для узла 3

4) i2R2 + Uc = E для контура I

5) i3R3 + I4R4 - Uc = 0 для контура II

6) UL + iLR1 - i3R3 - i2R2 = 0 для контура III

2.2 Запись системы уравнений Кирхгофа в матричной форме

3. Формирование регионов с определениями основных фунций

3.1 Операторы определений функций для расчета матричных элементов модели цепи в пространстве состояний

3.1.1 Получение матрицы системы уравнений Кирхгофа (1.5):

3.1.2 Получение матрицы левой части системы уравнений (1.7):

3.1.3 Получение матрицы правой части системы уравнений (1.7):

Проверка выражения для определителя матрицы D0:

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

3.1.4 Получение матрицы правой части системы уравнений (1.9):

3.1.5 Получение матрицы A и вектора b системы дифференциальных уравнений (1.11):

3.1.6 Получение матрицы наблюдения:

Формирование 1-й строки матрицы наблюдения:

(для напряжения на R4)

Формирование 2-й строки матрицы наблюдения:

(для тока через емкость)

Компоновка К1 и К2 в единую матрицу наблюдения К

3.2 Формирование функций расчета матричной экспоненты

3.2.1 Определение скалярной функций "физически реализуемой"

экспоненты

3.2.2 Определение функций расчета собственных чисел

матрицы A

3.2.3 Определение функций расчета собственных векторов

матрицы A

3.2.4 Определение функции "физически реализуемой" матричной

экспоненты в собственном базисе матрицы А, диагональная

форма (2.17)

3.2.5 Получение функции расчета матричной экспоненты в

исходном базисе

3.3 Определение функции расчета реакции цепи на включение единичного источника ЭДС

3.3.1 Определение функции скачка единичной амплитуды

3.3.2 Определение функции для расчета установившегося значения

вектора переменных состояния

3.3.3 Определение функции для расчета значений вектора

переменных состояния в отклике на включение ЭДС

3.3.4 Определение функции для расчета значений вектора

наблюдаемых переменных

3.4 Определение функции расчета реакции цепи на одиночный импульс единичной амплитуды длительности

3.4.1 Определение функции получения прямоугольного импульса единичной амплитуды длительностью

3.4.2 Построение функции расчета вектора переменных состояния:

3.4.3 Построение функции расчета вектора наблюдаемых переменных:

3.5 Определение функций расчета переходного процесса в цепи при подаче периодической последовательности прямоугольных импульсов единичной амплитуды с периодом следования T и скважностью

3.5.1 Определение функции получения последовательности прямоугольных

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

3.5.2 Построение функции расчета вектора переменных состояния:

3.5.3 Построение функции расчета вектора наблюдаемых переменных:

3.6 Определение функций расчета переходного процесса в цепи при подаче периодической последовательности прямоугольных импульсов единичной амплитуды с периодом следования T и скважностью

3.6.1 Построение функции расчета вектора начальных условий

для переменных состояния в установившемся режиме:

3.6.2 Построение функции расчета вектора переменных состояния:

3.6.3 Построение функции расчета вектора наблюдаемых переменных:

3.7 Построение функций, используемых при выборе величин сопротивлений

3.7.1 Построение функции для расчета показателя колебательности переходного процесса

3.7.2 Построение функции для расчета периода колебаний переходного процесса Т0

4. Задание численных величин параметров

4.1 Задание величин индуктивности и емкости

4.2 Выбор величин сопротивлений

Выбор сопротивлений, входящих в контур обмена:

Выбор сопротивлений, не входящих в контур обмена:

Расчет собственных чисел матрицы А

Собственные числа получились комплексно-сопряженными

Расчет показателя колебательности

Расчет периода колебаний переходного процесса

Получаем графики компонентов матричной экспоненты на интервале времени [0, 3T]

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

5. Получение реакции цепи на включение единичного источника ЭДС

- создаем ранжированную переменную времени,

изменяющуюся в интервале [0;3Т] с шагом 1мкс

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

5.1 Графики переменных состояния

Ток через индуктивность

Размещено на Allbest.ru


Подобные документы

  • Использование ранжированных переменных в программном пакете Mathcad. Создание матриц без использования шаблонов матриц, описание операторов для работы с векторами и матрицами. Решение систем линейных и нелинейных уравнений с помощью функций Mathcad.

    контрольная работа [964,6 K], добавлен 06.03.2011

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

    курсовая работа [240,5 K], добавлен 18.12.2011

  • Возможности математического пакета MathCad в среде Windows 98 для использования матричной алгебры и решения системы линейных алгебраических уравнений. Методы решения систем линейных алгебраических уравнений. Сравнение метода Гаусса с методом MathCad.

    практическая работа [62,6 K], добавлен 05.12.2009

  • Решение линейных дифференциальных уравнений численными и символьными методами в рамках пакета компьютерной математики MathCAD. Сравнения результов решений и применение их при исследовании функционирования автоматических систем и электрических агрегатов.

    контрольная работа [51,5 K], добавлен 07.05.2009

  • Назначение и состав системы MathCAD. Основные объекты входного языка и языка реализации. Характеристика элементов интерфейса пользователя, настройка состава панелей инструментов. Задачи линейной алгебры и решение дифференциальных уравнений в MathCAD.

    курс лекций [1,6 M], добавлен 13.11.2010

  • Схема и основные параметры элементов цепи. Вывод системы дифференциальных уравнений. Реализация алгоритма на языке программирования высокого уровня Pascal. Решение дифференциальных уравнений в пакете MathCAD. Решение интерполяции в пакете Excel.

    курсовая работа [375,4 K], добавлен 06.01.2011

  • Решение системы дифференциальных уравнений переходных процессов в RLC-цепи численным методом. Анализ графиков в Excel. Расчет переходного процесса в математическом пакете MathCad по точным формулам. Разработка программы на языке программирования Pascal.

    курсовая работа [777,3 K], добавлен 22.10.2012

  • Суть метода Рунге-Кутта и его свойства. Решение дифференциальных уравнений первого порядка. Вычислительный блок Given/Odesolve. Встроенные функции rkfixed, Rkadapt, Bulstoer. Решения линейных алгебраических уравнений в среде MathCad и Microsoft Excel.

    курсовая работа [1,1 M], добавлен 02.06.2014

  • Решение однородных дифференциальных уравнений в MathCad. Расчет значений функций напряжения на конденсаторе и тока в цепи второго порядка в свободном режиме при отсутствии гармонического воздействия с использованием системы MathCAD. Графики этих функций.

    курсовая работа [705,0 K], добавлен 21.01.2011

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

    отчет по практике [1,5 M], добавлен 11.09.2014

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.