写出 MATLAB代码 用来计算并绘制路径损耗与距离的关系图，采用的是Walfish-Ikegami模型对于非直射(NLOS)环境。

在MATLAB中，你可以使用以下代码来计算并绘制路径损耗与距离的关系图，基于Walfish-Ikegami模型处理非直射环境。首先，我们需要定义一些基本参数，然后计算路径损耗，最后用`plot`函数生成图形。这里是示例代码：

% 参数设置
fc = 2.4e9;           % 中心频率, GHz
h_r = 30;            % 发射机(地面站)高度, m
h_bs = 40;           % 基站高度, m
h_ms = 1.5;          % 移动站高度, m
theta = 30;          % 视角, 度
gamma = 3.2;         % 均匀阴影因子
sigma = 4.2;         % 均方差
d = linspace(1, 10, 100); % 距离范围, km

% 将度转换为弧度
theta_rad = deg2rad(theta);

% Walfish-Ikegami模型计算
alpha = sqrt(pi*gamma^2/2);
L0 = 32.4 + 20*log10(fc);        % 自由空间损耗
Ls = L0 + 20*log10((h_r+h_bs+h_ms)/sqrt(2*h_r*h_bs));   % 直射损耗
Lshad = 10*log10(sigma^2*sin(theta_rad)^2) + alpha*log10(1 + sin(theta_rad)^2);  % 阴影衰落
LnonLoS = Ls + Lshad;

% 绘制路径损耗与距离的关系图
figure;
plot(d, LnonLoS, '-o', 'LineWidth', 1.5);
xlabel('Distance (km)');
ylabel('Path Loss (dB)');
title('Path Loss vs Distance (Walfish-Ikegami NLOS)');
grid on;

% 查看详细信息
disp(['Mean Path Loss: ', num2str(mean(LnonLoS)) ' dB']);
disp(['Max Path Loss: ', num2str(max(LnonLoS)) ' dB at ', num2str(find(LnonLoS == max(LnonLoS))) ' km']);

%