convolution integral with dde

8 次查看（过去 30 天）

Angie 2019-6-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/466372-convolution-integral-with-dde

评论： Satheesh oe 2019-12-20

Hi guys, I am trying to solve this differential equation using dde. I have a problem with the integral term. In the code shown below, tau is the delay but I cannot just specify a constant value because it is also in the integral which goes from 0 to t. Does anywone know how to deal with it? Thank you!

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;