функция управления в Matlab - PullRequest
0 голосов
/ 10 октября 2018

Привет ~ Я учусь с помощью метода Шутига для решения дифференциального уравнения с граничным условием.Проблема в (из учебника Зауэра):

enter image description here

function z=F(s)
a=0;b=1;yb=3;
ydot=@(t,y) [y(2);4*y(1)];
[t,y]=ode45(ydot,[a,b],[1,s]);
z=y(end,1)-yb; % end means last entry of solution y

Мой вопрос о "ydot = @ (t, y) [y (2); 4 * у (1)];».Я знаю, что это дескриптор функции.Но что здесь у (2) и у (1)?Я видел функцию дескриптора с круглыми скобками.Почему у нас здесь квадратные скобки?

1 Ответ

0 голосов
/ 10 октября 2018

Поскольку каждый ODE может быть преобразован в систему ODE первого порядка, почти все решатели ODE требуют преобразования в ODE первого порядка, прежде чем передать правую сторону вашего ODE.

Для вашего второгоODE порядка y '' = 4y set y_1 = y и y_2 = y '.Тогда y_1 '= y' = y_2 и y_2 '= y' '= 4y_1.Теперь вы можете написать:

enter image description here

Теперь ясно, что ваш дескриптор функции находится справа от этого ODE первого порядка.

...