Отчёт к первой лабораторной работе
Описание файла
Документ из архива "Отчёт к первой лабораторной работе", который расположен в категории "". Всё это находится в предмете "управление техническими системами (утс)" из 8 семестр, которые можно найти в файловом архиве МГТУ им. Н.Э.Баумана. Не смотря на прямую связь этого архива с МГТУ им. Н.Э.Баумана, его также можно найти и в других разделах. Архив можно найти в разделе "лабораторные работы", в предмете "управление техническими системами (утс)" в общих файлах.
Онлайн просмотр документа "Отчёт к первой лабораторной работе"
Текст из документа "Отчёт к первой лабораторной работе"
Отчет о лабораторной работе №1
Целью данной работы было ознакомление с пакетом моделирования MatLAB и освоение основных приемов моделирования систем автоматического управления.
ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ
4.1. Ознакомиться с пакетом прикладных программ MATLAB (см. Приложение А).
4.2. В соответствии с вариантом задания (см. табл.1.1) построить схему моделирования линейной системы автоматического управления, используя уравнение (1.1.-1.3).
4.3. Осуществить моделирование системы при двух видах входных воздействий:
а) u = 1(t). Начальные условия нулевые. При этом график функции будет являться графиком переходной характеристики , а график функции будет являться графиком импульсной переходной функции . На монитор выводить графики сигналов и u(t). Продолжительности интервалов наблюдения выбрать самостоятельно.
б) u = sin(t). Начальные условия нулевые. На монитор выводить графики сигналов и u(t). Продолжительности интервалов наблюдения выбрать самостоятельно.
4.4. Осуществить моделирование свободного движения системы с нулевыми и ненулевыми начальными условиями (см. табл.1.2). Снять выходные характеристики x(t) системы автоматического управления.
Началные условия:
a0=7
a1=3
a2=1
b0=10
Решение:
Текст программы
а) u = 1(t). Начальные условия нулевые.
function dy=ode(t,y)
a0=7;
a1=3;
b=10;
u=1;
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=b*u-a0*y(1)-a1*y(2);
б) u = sin(t). Начальные условия нулевые.
function dy=ode1(t,y)
a0=7;
a1=3;
b=10;
u=sin(t);
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=b*u-a0*y(1)-a1*y(2);
Построение графиков
hold on
[t,y]=ode45('ode',[0 5],[0 0]);
plot(t,1,'*',t,y(:,1))
hPlot = plot(t,1,'*',t,y(:,1));
set( hPlot, 'LineWidth', 2 );
[t,y]=ode45('ode1',[0 5],[0 0]);
plot(t,1,t,y(:,1))
grid on
hold on
[t,y]=ode45('ode',[0 5],[1 0.4]);
plot(t,1,'*',t,y(:,1))
hPlot = plot(t,1,'*',t,y(:,1));
set( hPlot, 'LineWidth', 2 );
[t,y]=ode45('ode1',[0 5],[ 1 0.4]);
plot(t,1,t,y(:,1))
grid on