嗨,肖恩,
代码中的问题在于变量的计算
β
在第二个循环。目前,计算
β
外的循环,在循环迭代及其价值保持不变。为了解决这个问题,你需要计算
β
内循环,对应的当前值
问
。这是修改后的代码:
clc
关闭所有
清晰的所有
λ= 332.8 e-9;%在真空波长的光
k = 2 *π/λ;
= 3 e;%衍射圆孔半径
Io = 1.0;%相对强度
R = 1;%的距离屏幕从光圈
half_angle = 13.3;%的光锥半角
center_angle = 0;%中心光锥角
Y = (-0.55 e-0:1e-3:0.9e-0);
Z = Y;%屏幕的坐标
我(1:长度(Y), 1:长度(Z)) = 0;
为i = 1:长度(Y)
为j = 1:长度(Z)
q = (Y(我)^ 2 + Z (j)。^ 2)。^ 0.5;
β= k * R * q /;%计算β内循环
我(I, j) = Io。* ((besselj(1(β))/(β+ 1 e-12)) ^ 2);
结束
结束
%计算分数锥的权力
θ=量化(Y, Z);
cone_indices =(θ> center_angle-half_angle) &(θ< center_angle + half_angle);
power_fraction =总和(sum(我(cone_indices))) /笔(金额(I));
%计算远场使用弗劳恩霍夫近似
N = 1000;%在远场的分
θ= linspace(-π/ 2π/ 2 N);%在远场角坐标
E_theta = 0 (1, N);
为i = 1:长度(θ)
q = k * *罪(θ(i));%孔中心距离远场
β= k * R * q /;%计算β内循环
E_theta (i) = (Io /(2 *π* R * q)) * ((2 * besselj(β)/(β+ 1 e-12)) * exp(我* k * q) 1) / (q);
结束
%计算功率在远场锥
θ= abs(θ)* 180 /π;
cone_indices =(θ> center_angle-half_angle) &(θ< center_angle + half_angle);
power_fraction_farfield =总和(abs (E_theta (cone_indices)) ^ 2) /笔(abs (E_theta)。^ 2);
disp ([“权力分数内锥:”num2str (power_fraction)]);
通过计算
β
在循环中,确保它对应的当前值
问
为每个迭代。这应该提供正确的分数在锥和远场。
