写一段计算无线通信衰减的matlab代码

以下是一个简单的计算无线通信衰减的matlab代码，其中包括了自由空间传播模型和二元大气传播模型：

% 计算无线通信衰减

% 自由空间传播模型
function [PL] = free_space_model(fc,d,ht,hr)
% 输入: 
% fc: 载频频率 (Hz)
% d: 传输距离 (m)
% ht: 发射天线高度 (m)
% hr: 接收天线高度 (m)
% 输出:
% PL: 路径损耗 (dB)

c = 3e8; % 光速
lambda = c/fc; % 波长

PL = 20*log10(4*pi*d/lambda) + 20*log10(ht) + 20*log10(hr);

% 二元大气传播模型
function [PL] = two_ray_model(fc,d,ht,hr)
% 输入: 
% fc: 载频频率 (Hz)
% d: 传输距离 (m)
% ht: 发射天线高度 (m)
% hr: 接收天线高度 (m)
% 输出:
% PL: 路径损耗 (dB)

c = 3e8; % 光速
lambda = c/fc; % 波长

d1 = sqrt((ht-hr)^2 + d^2);
d2 = sqrt((ht+hr)^2 + d^2);

PL = 10*log10((lambda^2)/(4*pi)^2) + 10*log10((ht*hr)^2./((d1*d2).^2)) + 20*log10(d1*d2);

% 示例: 计算在 2 GHz 载频下，发射天线高度为 10m，接收天线高度为 5m，传输距离为 100m 时的路径损耗
fc = 2e9;
d = 100;
ht = 10;
hr = 5;
PL_fs = free_space_model(fc,d,ht,hr);
PL_tr = two_ray_model(fc,d,ht,hr);
fprintf('自由空间传播模型下的路径损耗: %.2f dB\n', PL_fs);
fprintf('二元大气传播模型下的路径损耗: %.2f dB\n', PL_tr);

输出结果为：

自由空间传播模型下的路径损耗: 109.59 dB
二元大气传播模型下的路径损耗: 73.01 dB