Double integration using experimental data

1 次查看（过去 30 天）

Manuel 2015-6-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/222786-double-integration-using-experimental-data

编辑： Manuel 2015-6-18

I need to perform a double integration using experimental data, but my integration limits are the same for each integral, in this case the time. After review some old posts I reached to the following code:

   T = 0:0.1:600;
   x = T;
   A = rand(1,length(T)); % my data
   pp_int = spline(T,A );
   DoubleIntegration = integral(@(x)arrayfun(@(T )(integral(@(T ) ppval(pp_int,T ),0,  T  )),x),0,T(end)  );

The code take so long to run, and give huge values. Someone I can indicate whether I am doing something wrong? If this is the case how I can correct ?

Regards