非傍轴近似下会聚球面波的matlab代码

非傍轴近似下的球面波可以使用以下代码实现：

% 定义参数
k = 2*pi/lambda;   % 波数
r = linspace(0,R,100);   % 距离
theta = linspace(0,pi,100);   % 极角
phi = linspace(0,2*pi,100);   % 方位角

% 生成网格
[R_grid, Theta_grid, Phi_grid] = meshgrid(r, theta, phi);

% 计算球面波
E = exp(1i*k*R_grid)./R_grid;

% 可视化结果
figure;
isosurface(R_grid.*sin(Theta_grid).*cos(Phi_grid), ...
    R_grid.*sin(Theta_grid).*sin(Phi_grid), ...
    R_grid.*cos(Theta_grid), abs(E).^2);
axis equal;
xlabel('X');
ylabel('Y');
zlabel('Z');

其中，`lambda` 是波长，`R` 是球半径。该代码使用 `meshgrid` 生成网格，然后计算球面波的电场强度 `E`，最后使用 `isosurface` 函数可视化结果。

需要注意的是，非傍轴近似下的球面波通常是高度对称的，因此只需要在第一象限生成网格即可。同时，计算球面波时需要除以距离 `R`，以保证能量守恒。