Okay, so there seem to be a lot of things going on (i.e. wrong) here, too much for someone else to solve for you. So a few points to get you started:
- When you calculate E3, it appears you have taken all three E-field directions and smashed them into the same equation. You probably only want the y hat part of (2.21), namely ( z' here is "z prime")
% if z > z'
E3_y = (ky.^2-k^2).*sin(kx*x).*cos(ky*y) .* sin(beta*z ).*sin(beta*(d-z'));
% if z < z'
E3_y = (ky.^2-k^2).*sin(kx*x).*cos(ky*y) .* sin(beta*z').*cos(beta*(d-z ));
- In your for loops, you'll want to use indexing, for example:
for i_x = 1:length(x)
for i_z = 1:length(z)
% Equations use x(i_x) and z(i_z), for example,
E4(i_x,i_z) = x(i_x).*z(i_z);
end
end
- m and n in the photo start at 0 and go to infinity. Why did you choose 4:100? (Note: You might have issues if the equation doesn't converge quickly.)
- m and n have separate indexing, so you need to separate them somehow, by indexing for example: (There are other, probably better, ways to do this though.)
for i_kx=1:length(kx)
% f(x,z,kx,ky) is the part in the summation
E_whatever(i_kx) = sum(f(x(i_x),z(i_z),kx(i_kx),ky;)
end
- Add more points to x and z (make it about 100, at least) to get a better surface.
- You'll eventually want to use meshgrid instead of your current linspace of x and z, and surf or mesh instead of your current plot3.
After you try all that, come back with your latest problems and maybe someone can get you a better answer. Good luck!
