大家好!我是算Third harmonic generation的。当我在算到radiation pattern的时候,老是不对。有算过的虫友能否给个提示?感激不尽!!!
主要matlab code:
n1=1.5; % refractive index
NA=1.4; % 数值孔径
almax=asin(NA/n1);
%thetaR=linspace(0,almax,N);
phiR=linspace(0,2*pi,N);
U=yy; %yy是 THG power
for n=1:length(thetar)
thetar(n, =thetar(1, ;
end
% expanding phi to span entire space
phir=phir';
for mm=1:length(phir)
phir(:,mm)=phir(:,1);
end
% expanding U to span entire space
for k=1:length(yy)
U(k, =yy(1, ;
end
% converting to spherical coordinates
[x,y,z]=sph2cart(phir,thetar,U); %%2pi在前,pi在后面,最后r %%%%
%plotting routine
surf(x,y,z)
title('Radition Pattern (D=1.0lamda)')
xlabel('x')
ylabel('y')
zlabel('z') |