Files

copied

Scanning the repository...

Last update 6 years 9 months by JosefMatous
FilesMATLAB optimizationlevitation
..
objectiveGrad.m
optFunc.m
p.m
px.m
pxx.m
pxxx.m
pxy.m
pxyy.m
pxz.m
pxzz.m
py.m
pyxx.m
pyy.m
pyyy.m
pyz.m
pyzz.m
pz.m
pzxx.m
pzyy.m
pzz.m
pzzz.m
solverscript.m
pxy.m
function pxy = pxy(P0,k,r,tx,ty,pos) x = pos(1); y = pos(2); z = pos(3); xr = x - tx; yr = y - ty; d = sqrt(xr.^2 + yr.^2 + z.^2); epsilon = 1e-6; % small error sinTheta = sqrt(xr.^2 + yr.^2)./d + epsilon; pxy = (2.*exp(d.*i.*k).*P0.*xr.*yr.*(2.*k.*r.*sinTheta.*(-1 + sinTheta.^2).*(-2 + (-1 + d.*i.*k).*sinTheta.^2).*besselj(0, k.*r.*sinTheta) + ... (-8 + sinTheta.^2.*(4 + sinTheta.^2 + k.*((-d.^2).*i.^2.*k.*sinTheta.^2 - d.*i.*(-4 + sinTheta.^2) + k.*r.^2.*(-1 + sinTheta.^2).^2))).*besselj(1, k.*r.*sinTheta)))./(d.^5.*k.*r.*sinTheta.^5); end
Report a bug