Convolution integral with matrices

1 次查看（过去 30 天）

显示更早的评论

Nicolás Faedo 2018-3-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/391246-convolution-integral-with-matrices

编辑： Walter Roberson 2019-4-30

在 MATLAB Online 中打开

Hey everybody! My question is quite simple I guess:

I would like to perform the following convolution integral (written in latex form for the sake of clarity):

, where A is 2x2, B is 2x1 and u is scalar.

I would like to compute this WITHOUT symbolic functions.

So far, this is what I've done:

T = 0:0.01:20 % Assume 20 seconds of simulation.
% Compute the values of expm(A) and B*u for each time step ---------
  for i=1:length(T)
    tt=t(i);
    E(:,:,i) = expm(A*tt);
    U(:,i) = B*u(tt);
  end
 nu = size(A,2); % Assume that A is 2x2.
 result2 = 0;
 for j=1:nu
    AA=squeeze(A(:,j,:));
    result2 = result2 + median(diff(T))*conv2(AA,B(j,:),'full'); 
 end