АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомДругоеЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция

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);


1 | 2 | 3 | 4 | 5 | 6 | 7 |

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



Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Студалл.Орг (0.002 сек.)