错误: 文件: tfwendu3.m 行: 58 列: 1 脚本中的函数定义必须出现在文件的结尾。 请将 "greenhouse_model" 函数定义后面的所有语句都移到第一个局部函数定义前面。

这个错误是因为在 MATLAB 脚本中，函数定义必须出现在文件的结尾，而你将 "greenhouse_model" 函数定义后面的语句放在了函数定义前面，导致了这个错误。要解决这个问题，你需要将 "greenhouse_model" 函数定义后面的语句都移到第一个局部函数定义前面，如下所示：

% 温室参数
A = 10; % 温室面积
H = 2.5; % 温室高度
U = 2.5; % 温室墙体总导热系数
d = 0.2; % 温室墙体厚度
cp = 1005; % 空气比热容
rho = 1.2; % 空气密度
alpha = 0.7; % 温室玻璃透过率
epsilon = 0.9; % 温室玻璃辐射率
sigma = 5.67E-8; % 环境辐射常数
Tout = 10; % 外部温度
RHout = 0.7; % 外部相对湿度
I = 1000; % 光照强度
Tsoil = 15; % 土壤温度
RHsoil = 0.6; % 土壤相对湿度
E = 0.1; % 植物蒸腾率
F = 0.1; % 植物蒸发散热率
tau = 0.8; % 大气透明度
v = 2; % 风速
theta = 30; % 风向
tmax = 86400; % 模拟时间

% 初始条件
T0 = 20; % 初始温度
RH0 = 0.5; % 初始相对湿度

% 时间步长
dt = 60;

% 计算参数
k = U / d;
Qsoil = A * F * (T0 - Tsoil);
Qplant = A * E * (T0 - Tout);
Qrad = A * epsilon * sigma * H * (T0^4 - Tout^4);
Qconv = A * alpha * (T0 - Tout);
Qvent = A * rho * cp * v * (T0 - Tout) / 3600 * 10; % 假设温室有10个通风口

% 建立生长模型
function dTdt = plant_model(T, Qplant)
    E = 0.1; % 植物蒸腾率
    F = 0.1; % 植物蒸发散热率
    Qtrans = A * E * (T - Tout);
    Qevap = A * F * (T - Tsoil);
    Qplant = Qtrans + Qevap;
    dTdt = (Qplant + Qsoil - Qrad - Qconv - Qvent) / (A * rho * cp * H);
end

% 建立模型
function dTdt = greenhouse_model(t, T)
    Qsoil = A * F * (T - Tsoil);
    Qrad = A * epsilon * sigma * H * (T^4 - Tout^4);
    Qconv = A * alpha * (T - Tout);
    Qvent = A * rho * cp * v * (T - Tout) / 3600 * 10; % 假设温室有10个通风口
    dTdt = (Qsoil + Qplant(T, Qplant) - Qrad - Qconv - Qvent) / (A * rho * cp * H);
end

% 数值计算
t = 0:dt:tmax;
T = zeros(size(t));
RH = zeros(size(t));
T(1) = T0;
RH(1) = RH0;
for i = 2:length(t)
    T(i) = T(i-1) + dt * greenhouse_model(t(i-1), T(i-1));
    RH(i) = RH(i-1); % 简化模型，假设相对湿度不变
end

% 输出结果
plot(t/3600, T);
xlabel('时间（小时）');
ylabel('温度（摄氏度）');
title('温室温度随时间变化曲线');

将 "greenhouse_model" 函数定义后面的语句都移到第一个局部函数定义前面之后，就可以避免这个错误了。