% 模型参数P0 = 101325; % 初始压强，PaT0 = 293.15; % 初始温度，KV0 = 1e-9; % 初始体积，m^3R = 8.314; % 气体常数，J/mol/KM = 18e-3; % 水分子质量，kg/molrho = 997; % 水的密度，kg/m^3Cp = 4182; % 水的比热容，J/kg/KDv = 2.26e-5; % 水的蒸汽扩散系数，m^2/s% 时间参数dt = 0.1; % 时间步长，stmax = 60*10; % 最大模拟时间，st = 0:dt:tmax; % 时间向量% 初始化RH = zeros(size(t)); % 相对湿度，%RH(1) = 0;V = zeros(size(t)); % 体积，m^3V(1) = V0;M = zeros(size(t)); % 质量，kgM(1) = rho*V0;Ts = zeros(size(t)); % 表面温度，KTs(1) = 293.15; % 初始表面温度为室温% 数值模拟for i = 2:numel(t) % 计算相对湿度 P = P0*exp(-M(i-1)*g/(R*T0*V(i-1))); % 当前压强，Pa RH(i) = 100*P/wsat(Ts(i-1))/P0; % 当前相对湿度，% % 更新湿度和体积 dVdt = 4/3*pi*(3*V(i-1)/(4*pi))^(2/3)*Dv/rho*RH(i)*P0/P; % 体积变化率，m^3/s V(i) = V(i-1) + dt*dVdt; % 更新质量和表面温度 dMdt = -4*pi*rho*Dv*(3*V(i-1)/(4*pi))^(2/3)*RH(i)*P0/P; % 质量变化率，kg/s M(i) = M(i-1) + dt*dMdt; Ts(i) = Ts(i-1) + dt*(dMdt/(4*pi*V(i-1)*Cp*rho) - 1/V(i-1)*dVdt/(4*pi)*Ts(i-1));end% 绘图figure;plot(t/60, Ts-273.15);xlabel('时间（分钟）');ylabel('表面温度（℃）');title('水滴蒸发模型');grid on;优化这段程序