Технический Университет «мами»




Скачать 46,21 Kb.
НазваниеТехнический Университет «мами»
Дата публикации16.09.2013
Размер46,21 Kb.
ТипОтчет
pochit.ru > Математика > Отчет
Федеральное агентство по образованию РФ Московский Государственный Технический Университет

«МАМИ»
Кафедра: «Автоматика и процессы управления»

Факультет: «АиУ»

Дисциплина: «Прогр. МП контр.»


Отчёт по лабораторным работам


Группа: 5-УИк-7

Студент: Филоненко А. В.

Преподаватель: Мурачёв Е. Г.

Вариант 10

Москва, 2008 г.

Исходный код программы.
Dim x As Double 'действительная часть переходной функции

Dim y As Double 'мнимая часть числителя переходной функции

Dim z As Double 'знаменатель переходной функции

Dim a, b, c, d As Double 'вспомогательные параметры

Dim omega As Double 'частота

Dim alert_naikvist As Boolean 'проверка по критерию Найквиста

Dim alert_mikhailov As Boolean 'проверка по критерию Михайлова

Dim mikhailov As Integer 'степень многочлена, используемого для построения АФЧХ

Dim max_mikhailov As Integer 'число корней по АФЧХ на годографе

Dim e As Double

Function vis_cls() 'функция очистки экрана и легенды графика

LineRed.Visible = False

LineGreen.Visible = False

LineBlue.Visible = False

LblRed.Caption = ""

LblBlue.Caption = ""

LblGreen.Caption = ""

End Function

Private Sub Form_Load()

vis_cls

status.Caption = ("Приложение запущено")

End Sub

Private Sub BtnAFCH_Click() 'построение АФЧХ

pic.Cls

vis_cls

pic.Scale (-12, 12)-(12, -12)

pic.Line (-12, 0)-(12, 0)

pic.Line (0, 12)-(0, -12)

For n = -12 To 12 Step 3

pic.PSet (n, 0)

pic.Print n

pic.PSet (0, n)

pic.Print n

Next n

For omega = 0 To 0.25 Step 0.00001

a = 110 - 10000 * (omega ^ 2) 'определяем вспомогательные параметры

b = 300 * omega + 10000 * (omega ^ 3) 'определяем вспомогательные параметры

c = 21 - 2120 * (omega ^ 2) 'определяем вспомогательные параметры

d = 52 * (omega) + 100 * (omega ^ 3) 'определяем вспомогательные параметры

z = 441 - 86336 * (omega ^ 2) + 4504800 * (omega ^ 4) + 10000 * (omega ^ 6)

x = (a * c - b * d) / z

y = (a * d + b * c) / z

pic.PSet (x, y), vbRed

Next omega

LineRed.Visible = True

LblRed.Caption = " АФЧХ"

status.Caption = ("Построение графика АФЧХ закончено")

status.ForeColor = &H8000&

End Sub

Private Sub BtnReIm_Click() 'построение графиков действительной и мнимой составляющих переходной функции

pic.Cls

vis_cls

LineBlue.Visible = True

LblBlue.Caption = " Действительная"

LineGreen.Visible = True

LblGreen.Caption = " Мнимая"

pic.Scale (-50, 10)-(500, -110)

pic.Line (-10, 0)-(500, 0)

pic.Line (0, 10)-(0, -110)

For n = 0 To 500 Step 50

pic.PSet (n, 0)

pic.Print n

Next n

For n = 0 To -110 Step -110

pic.Line (-30, n)-(0, n)

pic.PSet (-30, n)

pic.Print n

Next n

For omega = 0 To 999 Step 0.01

a = 110 - 10000 * (omega ^ 2) 'определяем вспомогательные параметры

b = 300 * omega + 10000 * (omega ^ 3) 'определяем вспомогательные параметры

c = 21 - 2120 * (omega ^ 2) 'определяем вспомогательные параметры

d = 52 * (omega) + 100 * (omega ^ 3) 'определяем вспомогательные параметры

z = 441 - 86336 * (omega ^ 2) + 4504800 * (omega ^ 4) + 10000 * (omega ^ 6)

x = (a * c - b * d) / z

y = (a * d + b * c) / z

pic.PSet (omega, x)

pic.PSet (omega, y), &H8000&

Next omega

status.Caption = ("Построение графиков действительной и мнимой составляющих переходной функции закончено")

status.ForeColor = &H8000&

End Sub

Private Sub BtnACHFCH_Click() 'построение АЧХ и ФЧХ

pic.Cls

vis_cls

LineBlue.Visible = True

LblBlue.Caption = " АЧХ"

LineGreen.Visible = True

LblGreen.Caption = " ФЧХ"

pic.Scale (-0.3, 15)-(3, -2)

pic.Line (0, 0)-(3, 0)

For n = 0 To 3 Step 0.3

pic.PSet (n, 0)

pic.Print n

Next n

For omega = 0 To 3 Step 0.0001

a = 110 - 10000 * (omega ^ 2) 'определяем вспомогательные параметры

b = 300 * omega + 10000 * (omega ^ 3) 'определяем вспомогательные параметры

c = 21 - 2120 * (omega ^ 2) 'определяем вспомогательные параметры

d = 52 * (omega) + 100 * (omega ^ 3) 'определяем вспомогательные параметры

z = 441 - 86336 * (omega ^ 2) + 4504800 * (omega ^ 4) + 10000 * (omega ^ 6)

x = (a * c - b * d) / z

y = (a * d + b * c) / z

pic.PSet (omega, Sqr((x ^ 2) + (y ^ 2))), vbBlue

pic.PSet (omega, (-Atn(y / x))), &H8000&

Next omega

status.Caption = ("Построение графиков построение АЧХ и ФЧХ закончено")

status.ForeColor = &H8000&

End Sub

Private Sub BtnLACH_Click() 'построение ЛАЧХ

pic.Cls

vis_cls

LineRed.Visible = True

LblRed.Caption = " ЛАЧХ"

pic.Scale (-5, 10)-(100, -5)

pic.Line (0, 0)-(100, 0)

For n = 0 To 100 Step 10

pic.PSet (n, 0)

pic.Print n

Next n

For omega = 0 To 100 Step 0.001

a = 110 - 10000 * (omega ^ 2) 'определяем вспомогательные параметры

b = 300 * omega + 10000 * (omega ^ 3) 'определяем вспомогательные параметры

c = 21 - 2120 * (omega ^ 2) 'определяем вспомогательные параметры

d = 52 * (omega) + 100 * (omega ^ 3) 'определяем вспомогательные параметры

z = 441 - 86336 * (omega ^ 2) + 4504800 * (omega ^ 4) + 10000 * (omega ^ 6)

x = (a * c - b * d) / z

y = (a * d + b * c) / z

pic.PSet (omega, Log(Sqr((x ^ 2) + (y ^ 2)))), vbRed

Next omega

status.Caption = ("Построение графика ЛАЧХ закончено")

status.ForeColor = &H8000&

End Sub

Private Sub BtnPerehod_Click() 'построение графика переходного процесса

pic.Cls

vis_cls

LineRed.Visible = True

LblRed.Caption = " Переходный процесс"

pic.Scale (-1, 7)-(700, 4)

pic.Line (-1, 4.5)-(700, 4.5)

For n = 0 To 700 Step 100

pic.PSet (n, 4.5)

pic.Print n

Next n

e = 2.718 'задаём значение экспоненты

For t = 0 To 700 Step 0.1

y = 5.3 * (e ^ (-2.1 * t)) + 5.2 - 0.57 * (e ^ (-0.01 * t)) * Cos(0.099 * t) - 0.066 * (e ^ (-0.01 * t)) * Sin(0.099 * t)

pic.PSet (t, y), vbRed

Next t

status.Caption = ("Построение графика переходного процесса закончено")

status.ForeColor = &H8000&

End Sub

Private Sub BtnStable_Click() 'проверка на устойчивость

BtnStable.Caption = "Идёт проверка..."

max_mikhailov = 6 'задаем вручную степень уравнения, описывающего АФЧХ

mikhailov = 0 'устанавливаем счетчик корней в нуль

For omega = 0 To 10000 Step 0.01

a = 110 - 10000 * (omega ^ 2) 'определяем вспомогательные параметры

b = 300 * omega + 10000 * (omega ^ 3) 'определяем вспомогательные параметры

c = 21 - 2120 * (omega ^ 2) 'определяем вспомогательные параметры

d = 52 * (omega) + 100 * (omega ^ 3) 'определяем вспомогательные параметры

z = 441 - 86336 * (omega ^ 2) + 4504800 * (omega ^ 4) + 10000 * (omega ^ 6)

x = (a * c - b * d) / z

y = (a * d + b * c) / z

If x = -1 Then alert_naikvist = True

If x = 0 And y = 0 Then alert_mikhailov = True 'при обращении функции в нуль сработает 'критерий Михайлова

If x = 0 Or y = 0 Then mikhailov = mikhailov + 1 'отсчет действительных и мнимых корней '(квадрантов годографа АФЧХ)

Next omega

'критерий Найквиста: если АФЧХ системы прошла через (-1;j0), то система устойчива,

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

If alert_naikvist = True Then MsgBox "По критерию Найквиста устойчивость системы определить нельзя.", , "Критерий Найквиста" Else MsgBox "Система устойчива по критерию Найквиста.", , "Критерий Найквиста"

'критерий Михайлова: если АФЧХ системы не проходит через точку (0;0) и число 'квадрантов,

'которые пересекает годограф равно числу корней характеристического уравнения, 'система устойчива,

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

If mikhailov = max_mikhailov And alert_mikhailov = False Then MsgBox "Система устойчива по критерию Михайлова.", , "Критерий Михайлова" Else MsgBox "По критерию Михайлова нельзя судить об устойчивости системы.", , "Критерий Михайлова"

status.Caption = ("Проверка на устойчивость закончена")

status.ForeColor = &H8000&

BtnStable.Caption = "Проверка на устойчивость"

End Sub

^ Интерфейс программы. Полученные графики.
Пользовательский интерфейс.



Построение АФЧХ фрагмента системы.


Построение действительной и мнимой характеристик передаточной функции.



Построение амплитудно-частотной и фазово-частотной характеристик фрагмента системы.




^ Построение логарифмической амплитудно-частотной характеристики фрагмента системы.




Построение переходного процесса фрагмента системы.


Проверка на устойчивость




Обратное преобразование Лапласа (с использованием ПО MathCad):




Похожие:

Технический Университет «мами» iconТехнический университет
Московский государственный инженерно-физический институт (технический университет)
Технический Университет «мами» iconЧенстоховский технический университет
Ченстоховский технический университет (по польски- политехника Ченстоховска) существует с 1949 года
Технический Университет «мами» iconУльяновский государственный технический университет
Ульяновский государственный технический университет в период с января по апрель 2012 года проводит интеллектуальный конкурс «Формула...
Технический Университет «мами» iconКрасноярский государственный технический университет сибирский государственный...
Методы нейроинформатики / Под ред. А. Н. Горбаня; отв за выпуск М. Г. Доррер. Кгту, Красноярск, 1998. 205 с
Технический Университет «мами» iconВолгоградский государственный технический университет (вгту) Автотракторный факультет

Технический Университет «мами» iconВысшего и профессионального образования липецкий государственный технический университет

Технический Университет «мами» iconУчебное пособие москва 2006 федеральное агентство по образованию...
А. С., Мурашова В. Г., Семенова Т. Н., Смирнов В. В., Тихонова С. В. Социология: Учебное пособие / Государственное образовательное...
Технический Университет «мами» iconРейтинг вузов России
Московский государственный технический университет им. H. Э. Баумана Москва [472691]
Технический Университет «мами» iconГубенко Станислав Николаевич
Харьковский Государственный Технический Университет Радиоэлектроники 2000-2006 г
Технический Университет «мами» iconНациональный технический университет
Распределение номеров теоретических и контрольных вопросов, тестов и задач по
Вы можете разместить ссылку на наш сайт:
Школьные материалы


При копировании материала укажите ссылку © 2019
контакты
pochit.ru
Главная страница