MATLAB: Overlap integral of the electric field distribution

MATLABnumerical integrationtrapz function

Hi, I am trying to calculate overlap integral (OI) of electric field distribution, using Matlab, when run the following code:

xc = linspace(0,1,M);
yc = linspace(0,1,M);
A = E1 * E2;
B = (abs(trapz(xc,trapz(yc, abs(A))) ) ).^2;
C = trapz(xc,trapz(yc, abs(E1).^2)) * trapz(xc,trapz(yc, abs(E2).^2 ));
OI = B./C;

where E1 is the electric field distribution of the desired mode and E2 is polarized mode. Both E1 and E2 are matrices of size MxM.

but, I have got this error: Error using trapz (line 57) *LENGTH(X) must equal the length of Y in dim 1.

Best Answer

xc = linspace(0,1,M);
yc = linspace(0,1,M);
A = E1 .* E2;
B = (trapz(xc,trapz(yc, abs(A), 1)) ).^2;
C = trapz(xc,trapz(yc, abs(E1).^2, 1)) * trapz(xc,trapz(yc, abs(E2).^2, 1 ));
OI = B./C;

Best wishes

Torsten.

Related Solutions

MATLAB: If i have two equations with three variables, how to get the ratio (transfer function) using matlab

Here you go:

syms E1 E2 E3 s
Eqn1 = E3 - E1 + E3*s/2 + (E3 - E2)/(4*s);
Eqn2 = (E2 - E3)/(4*s) + E2*s/3 + E2;
E31_sol = solve(Eqn1, E3);
E32_sol = solve(Eqn2, E3);
E2_sol = solve(E31_sol == E32_sol, E2);
E2_by_E1 = E2_sol/E1
E2_by_E1 =
6/(4*s^3 + 20*s^2 + 29*s + 12)

MATLAB: How to make this code like a cycle

JM - perhaps you can just replace the if/else block with

maxIterations = 1000;
atIteration = 0;
while (abs(E1)>E || abs(E2)>E) && atIteration <= maxIterations
    P=4*v2*sin(x1); 
    Q=-4*v2*cos(x1)+(4*(v2^2)); 
    J=[4*v2*cos(x1) 4*sin(x1); 4*v2*sin(x1) 8*v2-4*cos(x1)]; 
    
    A=[b2-P;b1-Q];
    B=inv(J)*A;    
    E1=B(1,:) %(Note that E1 and E2 have a new value)
    E2=B(2,:)
    x1=x1+E1
    v2=v2+E2 %(Note that x1 and v2 have a new value)
    
    atIteration = atIteration + 1;
end

So the while loop would exit when both abs(E1) AND abs(E2) are less than E. (Or when the maximum number of iterations has been reached.)

Best Answer

Related Solutions

MATLAB: If i have two equations with three variables, how to get the ratio (transfer function) using matlab

MATLAB: How to make this code like a cycle

Related Question