Программирование и основы алгоритмизации

По структурной схеме системы автоматического управления составлена система дифференциальных уравнений, описывающих её функционирование. Разработана программа на языке Pascal для решения системы методом Эйлера и построен график переходных процессов.
Краткое сожержание материала:

3

Федеральное агентство образования Российской Федерации

Тульский государственный университет

Кафедра «Системы автоматического управления»

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

к курсовой работе по дисциплине

«Программирование и основы алгоритмизации».

Разработал студент группы 120161

Матях Р.И.

Проверил доцент кафедры САУ:

Илюхин А. С.

Тула 2007 г.

ЗАДАНИЕ

По заданной структурной схеме САУ составить систему ДУ, описывающих её функционирование.

Разработать программу на языке Pascal для решения этой системы методом Эйлера.

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

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

Общий вид функциональной схемы САУ

3

Обозначения:

СУ - суммирующее устройство;

УУП - устройство усиления и преобразования;

ИУ - исполнительное устройство;

ОУ - объект управления;

УОС - устройство обратной связи;

g - управляющее воздействие;

Э - сигнал рассогласования;

u - сигнал на выходе УУП;

б - координата исполнительного органа исполнительного устройства;

y - регулируемая величина;

z - сигнал обратной связи.

Варианты УУП:

1. Пропорциональный регулятор

3

2. Пропорциональный интегральный регулятор

3

3. Пропорциональный дифференциальный регулятор

3

4.Интегрально-дифференциальный регулятор

3

5. Пропорциональный интегрально-дифференциальный регулятор

3

Варианты схем исполнительных механизмов:

1. Двигатель постоянной скорости

3

2. Автоматический привод

3

3. Автоматический привод

3

Варианты объектов управления:

W(p)=K4/((T4*p+1)*(T5*p+1)*(T6*p+1))

W(p)=K4/(T4²*p²+2*з*T4*p+1)

W(p)=K4/(p*(T4*p+1)*(T5*p+1))

W(p)=K4/(p*(T4²*p²+2*з*T4*p+1))

W(p)=K4/((T4*p+1) (T5²*p²+2*з*T5*p+1))

W(p)=K4/((T4*p+1)*(T5*p+1))

W(p)=K4/(p*(T4*p+1))

Передаточная функция:

W(p)=K6

1. Решение дифференциальных уравнений методом Эйлера

1.1 Описание метода

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

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

Пусть задано дифференциальное уравнение первого порядка с начальным условием (задача Коши)

dy/dx=F(x,y), (1)

У(x₀) = Уо (2)

где F(x,y) - заданная функция двух переменных x и у, x₀, у₀ - известные числа. Требуется определить функцию у = у(x) при х=>x₀. Уравнение (1) можно рассматривать как задание кривой через ее производную в координатной плоскости x, у, поскольку известно как вычислить производную в каждой точке этой кривой через ее координаты. В общем случае уравнению (1) удовлетворяет целое семейство кривых; начальное условие (2) позволяет выбрать из этого семейства одну определенную кривую, которая проходит через заданную точку x₀, y₀.

Для численного решения (1), (2) заменим область непрерывного изменения аргумента х дискретным множеством точек, т.е. введем сетку. Положим, что величина х изменяется от значения х=x₀ до значения х=b. Тогда, рассматривая равномерную сетку, получаем узловые точки x₀, x₁,… x_k,…b, находящиеся на расстоянии h друг от друга, т.е.

x_k+1 - x_k= h, k =0,1... , (3)

где h - шаг сетки. Соответствующие значения функции будем обозначать у_k, т.е.

y_k=y*(x_k).

Здесь у* (х) - функция, которая является приближенным решением (1), (2). Для получения численного решения, дифференциальное уравнение (1) заменяется уравнениями относительно значений функций у*(x) в узловых точках. Эти уравнения называются разностными. Простейшее разностное уравнение для (1) имеет вид

(y_k+1-y_k)/h=F(x_k, y_k) k=0,l ,... , (4)

Уравнение (4) следует из (1), если производную dy/dx приближенно представить через значения функции у(x) в соседних узлах.

Соотношения (2.12.4) можно записать в виде

Y_k+1=y_k+h*F(x_k,y_k). (5)

Тогда, учитывая (2), с помощью формулы (5) можно последовательно определить значения у₁,y₂,.… Этот метод приближенного решения (1), (2) называется методом Эйлера. Геометрическая интерпретация этой схемы дана на рис.1, где изображено поле интегральных кривых. Использование только первого члена формулы Тейлора означает движение не по интегральной кривой, а по касательной к ней. На каждом шаге мы заново находим касательную; следовательно, траектория движения будет ломаной линией. Из-за этого метод Эйлера иногда называют методом ломаных.

Рис.1

Доказывается, что если шаг сетки h стремится к нулю, то приближенное решение, определяемое (5), стремится к точному решению (1), (2), т.е. имеется факт сходимости приближенного решения к точному при h 0. Однако в условиях реальных вычислений на компьютере при конечном шаге целесообразно знать насколько полученное приближенное решение близко к точному. В качестве меры отклонения (нормы ошибки) часто используют величину

y= max| y(x_k) - y_k|. (6)

Здесь у(х) -- строгое решение (1), (2), y_k -- приближенное значение...