When simulating the case (“Diffuse Scattering Film for Automotive Display”) in the case base, I found that some of the script commands as as follows. What is the meaning of ” norm_calc = exp( -(Theta\/beam_width)^2 ); “<\/p>\n
norm_calc = exp( -(Theta\/beam_width)^2 ); # Where Theta_out is the matrix of reflection angles,and Theta = meshgridx(theta_out,phi_out); beam_width=2<\/p>\n
beam_norm = integrate(norm_calc*sin(Theta*pi\/180),1:2,theta_out*pi\/180,phi_out*pi\/180);<\/p>\n
beam_norm = 1\/beam_norm;<\/p>\n
BTDF(fc,1,1,:,:) = pinch(BTDF(fc,1,1,:,:)) + beam_norm*Ttotal(fc)*order_T(n,m,fc)*exp( -(ThetaDiff\/beam_width)^2 );<\/p>\n
The complete code is as follows:<\/p>\n
theta_out = linspace(0, 90, n_theta_out);<\/p>\n
phi_out = linspace(0, 360, n_phi_out);<\/p>\n
<\/p>\n
monitor_nameT = “T”;<\/p>\n
monitor_nameR = “R”;<\/p>\n
<\/p>\n
f = getdata(monitor_nameT,”f”);<\/p>\n
lambda = c\/f;<\/p>\n
<\/p>\n
# Get refractive indices<\/p>\n
dim = getdata(“R”, “dimension”);<\/p>\n
if ( (dim==2) & (%polarization angle% == 90) ) {<\/p>\n
index_R = real(getdata(“index_R”, “index_z”));<\/p>\n
index_T = real(getdata(“index_T”, “index_z”));<\/p>\n
} else {<\/p>\n
index_R = real(getdata(“index_R”, “index_x”));<\/p>\n
index_T = real(getdata(“index_T”, “index_x”));<\/p>\n
}<\/p>\n
# calculate Theta and Phi grids<\/p>\n
Theta = meshgridx(theta_out,phi_out);<\/p>\n
Phi = meshgridy(theta_out,phi_out);<\/p>\n
Ux = sin(Theta*pi\/180)*cos(Phi*pi\/180);<\/p>\n
Uy = sin(Theta*pi\/180)*sin(Phi*pi\/180);<\/p>\n
Uz = cos(Theta*pi\/180);<\/p>\n
norm_calc = exp( -(Theta\/beam_width)^2 );<\/p>\n
beam_norm = integrate(norm_calc*sin(Theta*pi\/180),1:2,theta_out*pi\/180,phi_out*pi\/180);<\/p>\n
beam_norm = 1\/beam_norm;<\/p>\n
<\/p>\n
Ttotal = abs(transmission(monitor_nameT));<\/p>\n
if (dist>0) {<\/p>\n
Rtotal = 1 – abs(transmission(monitor_nameR));<\/p>\n
} else {<\/p>\n
Rtotal = abs(transmission(monitor_nameR));<\/p>\n
}<\/p>\n
<\/p>\n
# T<\/p>\n
BTDF = matrix(length(lambda),1,1,length(theta_out),length(phi_out));<\/p>\n
Ttest = matrix(length(lambda));<\/p>\n
order_T = grating(monitor_nameT, 1:length(f), index_T, direction);<\/p>\n
order_T_ux = gratingu1(monitor_nameT, 1:length(f), index_T, direction);<\/p>\n
order_T_uy = gratingu2(monitor_nameT, 1:length(f), index_T, direction);<\/p>\n
for(fc = 1:length(f)) {<\/p>\n
for(n=1:length(order_T_ux(:,fc))) {<\/p>\n
for(m=1:length(order_T_uy(:,fc))) {<\/p>\n
order_uz = sqrt( 1 – order_T_ux(n,fc)^2 – order_T_uy(m,fc)^2 );<\/p>\n
if( (order_T(n,m,fc)!= 0) & (imag(order_uz) == 0.) ) {<\/p>\n
Tux = cos(-phi_in*pi\/180)*order_T_ux(n, fc) – sin(-phi_in*pi\/180)*order_T_uy(m, fc);<\/p>\n
Tuy = sin(-phi_in*pi\/180)*order_T_ux(n, fc) + cos(-phi_in*pi\/180)*order_T_uy(m, fc);<\/p>\n
ThetaDiff = real(acos( Ux*Tux + Uy*Tuy + Uz*order_uz ) * 180\/pi);<\/p>\n
BTDF(fc,1,1,:,:) = pinch(BTDF(fc,1,1,:,:)) + beam_norm*Ttotal(fc)*order_T(n,m,fc)*exp( -(ThetaDiff\/beam_width)^2 );<\/p>\n
}<\/p>\n
}<\/p>\n
}<\/p>\n
}<\/p>\n
Ttest = real(integrate(pinch(BTDF)*sin(Theta*pi\/180),1:2, theta_out*pi\/180, phi_out*pi\/180));<\/p>\n","protected":false},"template":"","class_list":["post-372589","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t