autocorrelation in Matlab
Dear all,
I am trying to plot the energy spectrum from the autocorrelation function. For this aim I am using the velocity fluctuation of a point of the domain varying through the time. the expression of the autocorrelation from Pope is: <u'(t)u'(t+s)>/<u'>^2, so in case s=0 the autocorrelation is 1. The code that I have written in Matlab is: u_primequad=zeros(n,1); %autocorr=0; for i=1:n-1 u_primequad(i,1)=((u_prime(i)*u_prime(i))); end u_primequad=mean(u_primequad); f=zeros(n,1); autocorr=zeros(n,1); for j=0:n-1 for i=1:n-j autocorr(i,1)=(u_prime(i)*u_prime(i+j)); end autocorr1=mean(autocorr); f(j+1,1)=autocorr1/u_primequad; end Someone can say me if are there some errors? Thank you in advance |
Hi,
I think you missed the correct forum. This forum is about OpenFOAM and not Matlab (: |
I am doing DNS of pipe flow, and I am wondering how to do autocorrelation using Matlab? Thank you.
|
Quote:
I m in a great need. Will be very helpful if you share. |
Quote:
|
Quote:
|
All times are GMT -4. The time now is 13:55. |