基于点核积分法的复杂辐射场的示例代码

以下是基于点核积分法的复杂辐射场示例代码：

% 定义辐射源位置和方向
source1_position = [-2, 0];
source1_direction = [0, 1];
source2_position = [2, 0];
source2_direction = [0, -1];

% 定义积分网格
x = linspace(-5, 5, 501);
y = linspace(-5, 5, 501);
[X, Y] = meshgrid(x, y);
dx = x(2) - x(1);
dy = y(2) - y(1);

% 计算每个网格上的辐射强度
I1 = zeros(size(X));
I2 = zeros(size(X));
for i = 1:numel(x)
    for j = 1:numel(y)
        r1 = [X(i,j) - source1_position(1), Y(i,j) - source1_position(2)];
        r2 = [X(i,j) - source2_position(1), Y(i,j) - source2_position(2)];
        I1(i,j) = dot(source1_direction, r1) / norm(r1)^2;
        I2(i,j) = dot(source2_direction, r2) / norm(r2)^2;
    end
end

% 绘制辐射场图像
I = I1 + I2;
imagesc(x, y, I);
axis equal;
colorbar;

在这个示例代码中，我们同样定义了两个辐射源source1和source2，位置分别为(-2, 0)和(2, 0)，辐射方向分别为[0, 1]和[0, -1]。定义了一个501×501的积分网格，用于计算每个网格上的辐射强度。使用了点核积分法来计算每个网格上的辐射强度，其中dot函数用于计算向量点积，norm函数用于计算向量模长。

最后使用imagesc函数将辐射场绘制成图像，其中x和y分别为网格的横纵坐标，I为每个网格上的辐射强度。可以看到，这个辐射场同样比前面的示例要复杂得多，但是通过点核积分法，我们可以计算出每个网格上的辐射强度并绘制成图像。