用matlab编写代码，请使用指针的方法来解决这个问题，计算任意输入信号波长，输入收发端高度和坐标、任意数量的刃峰高度和坐标，输出采用Deygout方法计算出的总绕射损耗： 以下给出步骤提示：1、计算各峰的v参数 2、比较v参数大小确定第一主峰； 3、计算从发射端到接收端的主峰的损耗L1（d1,d2,h）,注意h为相对等效高度 4、以第一主峰作为接收端，发射端不变，计算这两者之间的次级主峰的损耗，按照3步骤计算次级主峰的损耗L2，再以次级主峰作为接受或发射端计算再次一级刃峰的损耗 5、以第一主峰作为发射端重复4步骤 6、将所有刃峰损耗计算完后累加，得到总损耗，注意每个刃峰只需要算一次即可，不要重复计算 菲涅尔—基尔霍夫绕射参数：v=h根号下[2（d1+d2）/(lambdad1d2)] 其中h表示刃峰等效高度，lambda表示波长 损耗计算公式为：L=20*log10（0.225/v），只需要考虑x坐标即可

好的，以下是使用 MATLAB 编写的指针方法来解决这个问题的代码：

function total_loss = deygout_method(lambda, tx_height, tx_pos, rx_height, rx_pos, obstacles)
% lambda: 波长
% tx_height: 发射端高度
% tx_pos: 发射端坐标
% rx_height: 接收端高度
% rx_pos: 接收端坐标
% obstacles: 一个 n x 3 的矩阵，每行表示一个障碍物，第一列为障碍物高度，第二列为障碍物 x 坐标，第三列为障碍物 y 坐标

% 计算各峰的 v 参数
n = size(obstacles, 1);
v_values = zeros(n, 1);
for i = 1:n
    obstacle = obstacles(i, :);
    d1 = abs(obstacle(2) - tx_pos(1));
    d2 = abs(obstacle(2) - rx_pos(1));
    v_values(i) = obstacle(1) * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
end

% 比较 v 参数大小确定第一主峰
[~, idx] = max(v_values);
first_peak = obstacles(idx, :);

% 计算从发射端到接收端的主峰的损耗 L1
d1 = abs(first_peak(2) - tx_pos(1));
d2 = abs(first_peak(2) - rx_pos(1));
h = min(tx_height, rx_height);
v = h * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
L1 = 20 * log10(0.225 / v);

% 计算次级主峰的损耗
obstacles(idx, :) = [];
second_peak = [];
second_loss = 0;
while ~isempty(obstacles)
    % 以第一主峰作为接收端，发射端不变，计算这两者之间的次级主峰的损耗
    n = size(obstacles, 1);
    v_values = zeros(n, 1);
    for i = 1:n
        obstacle = obstacles(i, :);
        d1 = abs(obstacle(2) - tx_pos(1));
        d2 = abs(obstacle(2) - first_peak(2));
        v_values(i) = obstacle(1) * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
    end
    [~, idx] = max(v_values);
    second_peak = obstacles(idx, :);
    obstacles(idx, :) = [];
    
    % 计算次级主峰的损耗 L2，再以次级主峰作为接受或发射端计算再次一级刃峰的损耗
    d1 = abs(second_peak(2) - tx_pos(1));
    d2 = abs(second_peak(2) - first_peak(2));
    h = min(tx_height, second_peak(1));
    v = h * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
    loss = 20 * log10(0.225 / v);
    if ~isempty(obstacles)
        d1 = abs(first_peak(2) - tx_pos(1));
        d2 = abs(first_peak(2) - second_peak(2));
        h = min(tx_height, second_peak(1));
        v = h * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
        loss = loss + 20 * log10(0.225 / v);
    else
        d1 = abs(second_peak(2) - tx_pos(1));
        d2 = abs(second_peak(2) - rx_pos(1));
        h = min(tx_height, rx_height);
        v = h * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
        loss = loss + 20 * log10(0.225 / v);
    end
    second_loss = second_loss + loss;
end

% 以第一主峰作为发射端重复计算
obstacles = [first_peak; obstacles];
while ~isempty(obstacles)
    % 以第一主峰作为发射端，计算这两者之间的次级主峰的损耗
    n = size(obstacles, 1);
    v_values = zeros(n, 1);
    for i = 1:n
        obstacle = obstacles(i, :);
        d1 = abs(obstacle(2) - first_peak(2));
        d2 = abs(obstacle(2) - rx_pos(1));
        v_values(i) = obstacle(1) * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
    end
    [~, idx] = max(v_values);
    second_peak = obstacles(idx, :);
    obstacles(idx, :) = [];
    
    % 计算次级主峰的损耗 L2，再以次级主峰作为接受或发射端计算再次一级刃峰的损耗
    d1 = abs(first_peak(2) - tx_pos(1));
    d2 = abs(first_peak(2) - second_peak(2));
    h = min(rx_height, second_peak(1));
    v = h * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
    loss = 20 * log10(0.225 / v);
    if ~isempty(obstacles)
        d1 = abs(second_peak(2) - first_peak(2));
        d2 = abs(second_peak(2) - rx_pos(1));
        h = min(rx_height, second_peak(1));
        v = h * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
        loss = loss + 20 * log10(0.225 / v);
    else
        d1 = abs(first_peak(2) - tx_pos(1));
        d2 = abs(first_peak(2) - rx_pos(1));
        h = min(tx_height, rx_height);
        v = h * sqrt(2 * (d1 + d2) / (lambda * d1 * d2));
        loss = loss + 20 * log10(0.225 / v);
    end
    second_loss = second_loss + loss;
end

% 将所有刃峰损耗计算完后累加，得到总损耗
total_loss = L1 + second_loss;
end

需要注意的是，在 MATLAB 中使用指针需要使用引用传参的方式，因此这里使用了类似于指针的方式，将 `obstacles` 参数作为一个输入输出参数，每次计算完一个次级主峰后将其从 `obstacles` 中删除。