Метод средней точки для уравнения четвертого порядка - PullRequest
Новая
       63

Метод средней точки для уравнения четвертого порядка

0 голосов
/ 23 апреля 2020

Я пытаюсь решить ОДУ четвертого порядка, используя метод средней точки.

Система уравнений: 

system = @(t,z)[z(2);

    (c*z(2) + (k/a)*(exp(a*(z(1) + l*sin(z(3)))) - 1 + exp(a*(z(1) - ...
    l*sin(z(3)))) - 1) ...
    + 0.2*m*v*sin(w*t))/-m;

    z(4);

    -c/m*z(4) - 3/m/l*cos(z(3)) * k/a*(exp(a*(z(1) + l*sin(z(3))))...
    - exp(a*(z(1) - l*sin(z(3)))))];

или в виде математической формулы .

константы: 

m = 2500;
k = 1000;
a = 0.2;
l = 6;
c = 25;
w = 1.2*pi;
v = 50;

метод должен повторяться до t = 1000.

Кто-нибудь знает о функции, которая может решить эту проблему, используя метод средней точки в Matlab или Python?

...