Duhamel's Integral in MATLAB

I am trying to solve the equation of motion numerically by using Duhamel's Integral:

How do I code this with MATLAB?

Answer 1

[Torsten]

udotdot_g=@(x)...;
h=@(x)...;
u=@(t)-integral(@(tau)udotdot_g(tau).*h(t-tau),0,t);

Best wishes

Torsten.

Answer 2

[Daniele Casagrande]

Hi, I dealt with this topic when I studied the dynamic of structures. You could solve your problem by using the " conv" function of Matlab. When you start it, you can start testing with a simple function, such as sin(t) or cos(t) and then you can try an arbitrary time function, for example a seismic or wind signal. Search with the Matlab help for the " conv" function. An advice uses the ' same' type convolution of this function.

Ing. Daniele Casagrande, Italy