Plot(t,x),grid,title('dif32'),legend('X1','X2','X3')

Задание к примерам 1-5: произвести графический вывод координат каждой в отдельности

Пример 6. Интегрирование систем линейных дифференциальных уравнений в матричном виде, т.е.

(6.1)

где — вектор переменных состояния размера n´1, — вектор входного воздействия размера rx1, — матрицы чисел размера n´n, n´r, соответственно.

% Матрицы и берем из примера 2

% Создаем М-функцию под именем syst36.m

function ds36=syst36(t,x);

A=[-3,0,0;1,-2,0;0,4,-1]; % Матрица А размера 3´3

B=[10 0 0;0 0 0;0 0 0]; % Матрица B размера 3´3

u=[1;1;1]; % Входное воздействие размера 3´1

ds36=A*x+B*u; % Описание правой части матричного дифференциального уравнения (6.1)

% Создаем М-сценарий под именем ddd45_36

%Сценарий решения примера 3.6

T=[0 9];

x0=[0;0;0];

[t,x]=ode45('syst36',T,x0);

Поиск по сайту: