Интеграл свертки с дде - PullRequest
Новая
       42

Интеграл свертки с дде

0 голосов
/ 11 июня 2019

Я пытаюсь решить это дифференциальное уравнение, используя dde:

differential equation

У меня проблема с интегральным членом.В приведенном ниже коде tau - это задержка, но я не могу просто указать постоянное значение, поскольку оно также находится в интеграле, который идет от 0 до t.Кто-нибудь знает, как с этим бороться?Спасибо!

function sol = exer_3

   sol = dde23(@exer3f,tau,[0; 0],[0, 10]);
   figure
   plot(sol.x,sol.y)


function v = exer3f(t,y,Z)

    k = 125; m = 5; F = 1; w = 8;
    c=@(t)exp(-t^2);
    ylag = Z(:,1);
    v = zeros(2,1);
    v(1)=y(2);
    v(2) = -(k*y(1) - F*cos(w*t) + integral(@(tau)c(tau).*ylag(1), 0, t,'ArrayValued',true))./m;

Я также разместил этот вопрос на домашней странице MathWorks .

...