code de chevauchement intégrale dans matlab

Question

J'ai des données recueillies à partir de la simulation, Ces données représentent la distribution du champ électrique.

Est-ce une bonne façon de calculer le ovetlap integral (OI), où E1 est la distribution de champ électrique du mode souhaité et E2 est en mode polarisé. Les deux E1 et E2 sont des matrices de taille MxM.

xc = linspace(0,1,M); 
yc = linspace(0,1,M); 

A = desiredMode .* PolzMode; 

B = (trapz(xc,trapz(yc, abs(A)))).^2; 

C = trapz(xc,trapz(yc, abs(desiredMode).^2)) * trapz(xc,trapz(yc, abs(PolzMode).^2)); 

OI = B./C; 

Answer 1

C'est:

OI=sum(sum(E1.*conj(E2)*dx*dy))... 
    /sum(sum(E1.*conj(E1)*dx*dx))/sum(sum(E2.*conj(E2)*dy*dy)) 

• Puis-je remplacer le sum par trapz?
• Quelles sont les valeurs de dx, dy? Pourquoi ai-je utilisé la fonction conj au numérateur?
• Pourquoi ai-je utilisé la fonction conj au dénominateur?