How can I measure the interferometric autocorrelation of an electric field
10 views (last 30 days)
Show older comments
I am trying to evaluate the interferometric autocorrelation of an electric pulse. This means to evaluate
IAC(tau) = \int( (E(t) +E(t+tau))^2)^2 dt
where E is my electric field. I am trying to do that straightforward just by integrating my function E and its copy shifted in time. Has anyone an idea of how to that in a clever manner?
What I do is here:
n = length(E);
n_tau = 2*n; % 2* length vector
E_tau_app = zeros(1,n_tau); % temporary variable
E_long = [zeros(1,round(n/2)), E, zeros(1,round(n/2))];
IAC = zeros(n,1); % Interferometric Autocorrelation
E_tau = E_tau_app; % The Delay tau is the independent variable
E_tau(i:n+i-1) = E;
IAC(i) = sum(((E_tau(2:n_tau-1)+E_long(2:n_tau-1)).^2).^2)*dt + ...
0.5*((E_tau(1)+E_long(1))^4+(E_tau(n_tau)+E_long(n_tau))^4)*dt; % integral with trapezoidal method
Thank you in advance!
Find more on Pulse and Transition Metrics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!