请将 "calculate_k" 函数定义后面的所有语句都移到第一个局部函数定义 前面。

当然，以下是按照您的要求调整后的代码：

% 定义常数和变量
r_values = 2.5 - 3.5;
t_values = 0.4 - 0.6;
K_0 = 3 * exp(-4); % K_0 的值
ratio = 35/33;

% 函数c的计算
function c_val = calculate_c(s)
    E = 1; % 假设弹性模量E为1，实际应用中应替换为正确的值
    b = 1; % 假设泊松比b为1，同样应替换为正确的值
    c_val = (3 / (2 * E * b)) * ((2 * (s + 1)^3 * (3 * s^2 - 4 * s - 2) / (2 * s + 1)^(2.5)) * atan(2 * s + 1) + ...
        s * (3 * s^4 + 6 * s^3 + 11 * s^2 + 8 * s + 2) / ((s + 1) * (2 * s + 1)^2 + pi));
end

% 函数k的计算，用于寻找与线性模型的交点
function [s_intersection, t_intersection] = find_intersection_at_i(i)
    linear_k = K_0 * (ratio)^(i - 1);
    
    for t = t_values
        s_guess = (linear_k * t) / 2; % 初始猜测s值
        found = false;
        while abs(calculate_k(s_guess) - linear_k) > 0.01 && ~found
            s_guess = s_guess - fzero(@(s) calculate_k(s) - linear_k, s_guess); % 寻找零点
            found = true;
        end
        
        if found
            s_intersection = s_guess; % 存储交点对应的时间s
            t_intersection(i) = t; % 存储交点对应的时间t
            return; % 找到一个交点后立即返回
        end
    end
    
    warning('No intersection found for i=%d', i); % 如果没有找到交点，则发出警告
end

% 批量查找所有交点并存储结果
t_intersections = zeros(1, 10); % 初始化结果向量
for i = 1:10
    [s_intersection, t_intersection(i)] = find_intersection_at_i(i);
end

% 输出结果
t_intersection % 显示10个交点对应的t值

现在，`find_intersection_at_i`函数会依次计算每个i对应的交点，并在找到第一个交点后立即返回。这样，函数定义更紧凑，易于理解和管理。