Matlab-решение дифференциального уравнения третьего порядка

y''' + 41y'' + 360y' + 900y = 600x' + 1200x;
y(0)= 2 ; y'(0)= 1 ; y''(0) = -0.05

Как я могу решить это уравнение с помощью функции ODE45?

Я попытался это:

==>
function dydt=f(t,y)

dydt = [y(2) ; y(3) ; -41*y(3)-360*y(2)- 900*y(1)]

==>
clear all;

timerange=[0 1.4]; %seconds
initialvalues=[2 1 -0.05];


[t,y]=ode45(@dydt, timerange, initialvalues)

plot(t,y(:,1));

но мне нужно поставить часть X в уравнение - я не знаю, как...

6
differential-equations matlab ode

1 ответов

чтобы использовать ODE45 (или аналогичный), вам нужно преобразовать оду третьего порядка в систему оду первого порядка.

сделать так, пусть

y0 = y
y1 = y0'
y2 = y1'
y3 = y2'

затем 

y0' = y1
y1' = y2
y2' = y3

и

y3' = y''' = -41*y2-360y1-900y0+600dx+1200x

теперь вы можете использовать ODE45 для интеграции системы, вложив функцию, где x(t) и dx (t) доступны.

function test()

    % some random x function 
    x = @(t) exp(-t);
    dx = @(t) -exp(-t);

    % integrate numerically
    [T, Y] = ode45(@linearized, [0 1], [2 1 -0.05 0]);

    % plot the result
    plot(T, Y(:,1))

    % linearized ode
    function dy = linearized(t,y)
       dy = zeros(4,1);
       dy(1) = y(2);
       dy(2) = y(3);
       dy(3) = y(4);
       dy(4) = -41*y(3)-360*y(2)-900*y(1)+600*dx(t)+1200*x(t);
    end
end

plot

3