Double numerical integration of an Experimental data (no function)
3 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Hi all,
I have a rectangular domain ( in x and y) and have 64 points equally distributed data points around this rectangular domain (lets call the data points as A). I have the numerical value of each point on the grid. (sizes: x = 1 x 64, y = 1 x 64, A = 8 x 8)
B = double_integration(A) dx dy.
I tried the already given solution , result = trapz(x,trapz(y,mat,1),2) , but did not get any result. (Error was, LENGTH(X) must equal the length of Y in dim 1.)
Any ideas?
Many thanks
2 Commenti
John D'Errico
il 10 Ago 2016
Are the points scattered irregularly around the domain? Or are they on a regular lattice of some sort? Like an 8x8 grid?
Risposte (1)
Star Strider
il 10 Ago 2016
I would use trapz on each dimension of ‘A’:
int_1 = trapz(A);
Result = trapz(int_1)
5 Commenti
Star Strider
il 10 Ago 2016
The trapz function works very well with matrices. It considers each column of a matrix as individual vectors, and integrates them, creating a vector. The second trapz call integrates the vector. I do not know what ‘x’ and ‘y’ are with respect to ‘A’, but you can easily transpose ‘A’ to integrate with respect to each variable in order, if necessary.
Vedere anche
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!