Анализ систем автоматизированного управления численными методами
Курсовой проект - Компьютеры, программирование
Другие курсовые по предмету Компьютеры, программирование
Тема: Анализ систем автоматизированного управления численными методами
Введение
Бурное развитие новейшей техники и всё большее внедрение современных разделов математики в инженерные исследования неизмеримо повысили требования к математической подготовке инженеров и научных работников, занимающихся прикладными вопросами. В настоящее время, требуется знание многих разделов современной математики и в первую очередь основательное владение методами и приемами вычислительной математики, так как решение почти каждой инженерной задачи должно быть доведено до численного результата. Вычислительная техника наших дней представляет новые мощные средства для фактического выполнения счетной работы. Благодаря этому во многих случаях стало возможным отказаться от приближенной трактовки прикладных вопросов и перейти к решению задач в точной постановке. Разумное использование современной вычислительной техники не мыслимо без применения методов приближенного и численного анализа. Этим и объясняется чрезвычайно возросший интерес к методам вычислительной математики. А именно, ценным является тот метод, который допускает простую реализацию на машинах, а также не вносит в вычислительный процесс значительных погрешностей. Ведь может произойти накопление погрешностей округления, что приводит к неустойчивому состоянию схемы, которая в дальнейшем не пригодна для использования.
1. Техническое задание на проектирование Задана структурная схема следящей системы постоянного тока
Общий коэффициент передачи разомкнутой следящей системы определяется как:
=kиу*kфd*kун*kэму*kдв*kред
Таблица 1.1 Исходные данные для проектирования:
Ткз[с]Тупр[с]Тэм[c]Тя[с]Дискретность KvМетод решения0.2080.0460.1050.0210.05Рунге-Кутта
Необходимо: 1.Найти критический коэффициент (kvкрит) передачи следящей системы используя критерий устойчивости Гурвица и метод исключения Гаусса 2.Найти и графически построить переходные функции замкнутой системы для 2х вариантов: а). kv=1/2kvкрит б). kv=1/4kvкрит используя метод Рунге-Кутта. 3.Вычислить квадратичную интегральную оценку для 2х вариантов по импульсной переходной функции используя метод Симпсона.
2. Расчёт критического коэффициента передачи замкнутой следящей системы
По заданной структурной схеме записываем передаточную функцию разомкнутой следящей системы:
(Tкзp+1)(Тупрp+1)[(Tэмp+1)Tяp+1]
тогда характеристический многочлен замкнутой системы равен:
A(p)=p(Ткзp+1)(Tупрp+1)[(Tэмp+1)Tяp+1]+kv=a5p+a4p+a3p+a2p+a1p+a0
a5=Tкз*Тупр*Тэм*Тя; =Тупр*Tэм*Tя+Tкз*Tэм*Tя+Tкз*Тупр*Tя; =Tэм*Tя+Тупр*Tя+Tкз*Tя+Tкз*Тупр;
a2= Tя+Тупр+Tкз; a1=1; a0= kv;
2.1 Методика расчёта критического коэффициента
Критическим коэффициентом kvкрит называется коэффициент передачи разомкнутой системы находящейся на границе устойчивости В соответствии с критерием устойчивости Гурвица составляем определитель Гурвица. Для заданной системы 5-го порядка этот определитель имеет следующий вид:
Система будет находиться на границе устойчивости, если один из диагональных миноров равен 0, а все остальные - положительны. Процедура поиска критического коэффициента передачи состоит в следующем: задавая, последовательно увеличивая значение kv, начиная от kvнач=0.05 с дискретностью kv=0.05 определяем такое его значение при котором хотя бы один из миноров становится равным 0. Следовательно для каждого значения kv нужно вычислить все пять определителей. Вычисление определителей будем проводить методом исключения Гаусса с выбором главного элемента.
2.2 Метод Гаусса с выбором главного элемента
Метод исключения Гаусса состоит из двух основных этапов:
1)прямой ход - основан на приведении матрицы системы к треугольному виду. Это осуществляется последовательным исключением неизвестных из уравнений системы. 2) обратный ход - с помощью найденной треугольной матрицы последователь но вычисляем искомые неизвестные. Рассмотрим применение метода Гаусса для системы 3-го порядка (прямой ход):
а11х1+а12х2+а13х3=b1, (1)
а21х1+а22х2+а23х3=b2, (2)
а31х1+а32х2+а33х3=b3. (3)
Домножим уравнение (1) на ( - а21/а11 ) и прибавим его к уравнению (2). Затем, умножив уравнение (1) на ( - а31/а11 ) и прибавив результат к уравнению (3), получим равносильную систему уравнений вида
а11х1+а12х2+а13х3=b1, (4) а22х2+а23х3=b2, (5)
а32х2+а33х3=b3, (6)
где
а22 = а22 - (а21/а11)* а12; 2 = b2 - (а21/а11)* b1;
а23 = а23 - (а21/а11)* а13; 3 = b3 - (а31/а11)* b1;
а32 = а32 - (а31/а11)* а12;
а33 = а33 - (а31/а11)* а13;
Домножим уравнение (5) на ( - а32/а22 ) и складываем его с уравнением (6)
а11х1+а12х2+а13х3=b1, (7) а22х2+а23х3=b2, (8) а33х3=b3, (9)
где а33 = а33 - (а32/а22)* а23;3 = b3 - (а32/а22)* b2;
Мы привели матрицу системы к треугольному виду. На этом заканчивается прямой ход метода Гаусса. Обратный ход начинается с решения (9) уравнения системы:
х3 = b3/a33;
Используя это значение, можно найти х2 из (8) уравнения, а затем х1 из (7):
х2 = 1/a22*(b2 - a23x3), х1 = 1/a11*(b1 - a12x2 - a13x3);
Метод Гаусса с выбором главного элемента заключается в том, что перед началом исключения переменных необходимо привести матрицу к такому виду, чтобы максимальный элемент столбца попадал на главную диагональ. Эту процедуру необходимо в?/p>