MATLAB00 / MaTLAB_Simvol_LAB
.doc
% Создадим следующий сценарий под именем sdif6:
[x1,x2,x3]=dsolve('Dx1=-x1+2.3','Dx2=x3','Dx3=1.025*x1-0.4*x3-0.8*х2',
'x1(0)=0','x2(0)=0','x3(0)=0')
ezplot(x1,[-1,19]),grid, hold on,ezplot(x2,[-1,19],
title('Неоднородная система 3-го порядка'),
%ezplot(x3,[-1,19])
% Функция ezplot не позволяет строить графики с заданными цветами. Применим fplot. Для этого в функцию fplot следует вставить решения из командного окна MATLAB. Например, решение по первой координате имеет вид
x1 =
exp(-t)*(23/10*exp(t)-23/10)
% Тогда формат записи fplot для графика по х1 будет следующий (с красным цветом):
» fplot('exp(-t)*(23/10*exp(t)-23/10)',[-1,19],'r'),grid
% Через hold on можно добавить еще fplot по второй координате х2 и по третьей х3.
Задание.
Для объекта с передаточной функцией решить соответствующее дифференциальное уравнение при нулевых начальных условиях и построить график решения (переходный процесс).
Для того же объекта с передаточной функцией W(s) записать систему дифференциальных уравнений первого порядка, решить ее и сравнить с решением соответствующего дифференциального уравнения 3-го порядка.
Исследовать переходной процесс по выходным координатам двух колебательных звеньев, соединенных последовательно. Передаточные функции звеньев принять в виде:
-
Решение в символьном виде систем дифференциальных уравнений сравнить с решением численного метода с помощью решателя ode45.