在matlab上实现，有墙体遮蔽情况下，辐射源为立方体体源在房间内形成的辐射场强度分布，用surf函数绘制热图示例代码

以下是一个MATLAB画有墙体遮蔽情况下，辐射源为立方体体源在房间内形成的辐射场强度分布的示例代码，使用surf函数绘制热图：

% 定义房间的尺寸
Lx = 5; Ly = 5; Lz = 3;

% 定义立方体的尺寸和位置
ax = 0.2; ay = 0.2; az = 0.2; % 立方体的尺寸
x0 = 2; y0 = 2; z0 = 1; % 立方体的位置

% 定义计算点的坐标
[X,Y,Z] = meshgrid(0:0.1:Lx,0:0.1:Ly,0:0.1:Lz);

% 计算到立方体各表面的距离
dx = min([abs(X-x0-ax/2);abs(X-x0+ax/2)]);
dy = min([abs(Y-y0-ay/2);abs(Y-y0+ay/2)]);
dz = min([abs(Z-z0-az/2);abs(Z-z0+az/2)]);

% 计算到立方体的距离
d = sqrt(dx.^2+dy.^2+dz.^2);

% 计算辐射场强度分布
E = 1./(4*pi*8.85e-12*d.^2);

% 考虑墙体遮蔽情况
E(X==0 | X==Lx | Y==0 | Y==Ly | Z==0 | Z==Lz) = 0;

% 绘制辐射场强度分布
figure;
surf(X(:,:,1),Y(:,:,1),Z(:,:,1),E(:,:,1)); % 绘制热图
shading interp; % 设置颜色插值方式
xlabel('X'); ylabel('Y'); zlabel('Z'); % 设置坐标轴标签
title('辐射源为立方体的辐射场强度分布热图'); % 设置图像标题
colorbar; % 显示颜色条

在这个示例代码中，我们首先定义了房间的尺寸和立方体的尺寸和位置，然后定义了计算点的坐标，并计算了每个计算点到立方体各表面的距离和到立方体的距离。接着，我们计算了辐射场强度分布，考虑了墙体遮蔽情况，并使用surf函数绘制了热图，并设置了颜色插值方式、坐标轴标签、图像标题和颜色条。