在上面function程序加上从fengdian.m程序调取求出的t1

好的，我们可以在 `fsec` 函数中添加一个参数 `t1`，用来接收从 `fengdian` 函数中计算得到的 `t1`。具体实现如下：

function f = fsec(toff, t1)
    % 根据toff计算出M和C
    M = ...
    C = ...
    % 根据公式计算fsec
    fsec = 1 - M * exp(A) * cos(pi - B0 - C);
    % 计算适应度函数的值
    f = fsec;
end

然后，在 `fengdian` 函数中调用 `fsec` 函数时，需要将 `t1` 作为参数传递给 `fsec` 函数。修改后的 `fengdian` 函数代码如下：

function toff_opt = fengdian(PL, t1)
    % 定义适应度函数
    function f = fsec(toff, t1)
        % 根据toff计算出M和C
        M = ...
        C = ...
        % 根据公式计算fsec
        fsec = 1 - M * exp(A) * cos(pi - B0 - C);
        % 计算适应度函数的值
        f = fsec;
    end

    % 粒子群算法的代码
    % ... 省略粒子群算法的代码 ...

    % 输出最优解
    toff_opt = swarm(global_best_index,:);
    fprintf('最优的toff为：%f\n', toff_opt);
    fprintf('此时fsec的最小值为：%f\n', global_best);
end

这样，我们就可以在 `fengdian` 函数中调用 `fsec` 函数，并传递 `t1` 参数。