用MATLAB 写以下模型的程序1. 初始化：容器中的相对湿度为0%，水滴体积为V0，表面温度为20℃；2. 在每个时间步长中，计算当前时刻的相对湿度RH，并更新容器中的湿度；3. 计算水滴的质量变化量dM/dt，根据质量守恒原理得出新的水滴体积V；4. 根据新的水滴体积V计算液态水表面蒸发的分子数，进而计算水蒸气分压增加的量，从而更新相对湿度RH；5. 根据计算出的dM/dt和上述公式计算水滴表面温度Ts；6. 重复步骤2到5，直到水滴完全蒸发。

以下是可能的 MATLAB 程序，其中使用了欧拉方法进行数值模拟：

% 模型参数
P0 = 101325;    % 初始压强，Pa
T0 = 293.15;    % 初始温度，K
V0 = 1e-9;      % 初始体积，m^3
R = 8.314;      % 气体常数，J/mol/K
M = 18e-3;      % 水分子质量，kg/mol
rho = 997;      % 水的密度，kg/m^3
Cp = 4182;      % 水的比热容，J/kg/K
Dv = 2.26e-5;   % 水的蒸汽扩散系数，m^2/s

% 时间参数
dt = 0.1;       % 时间步长，s
tmax = 60*10;   % 最大模拟时间，s
t = 0:dt:tmax;  % 时间向量

% 初始化
RH = zeros(size(t));    % 相对湿度，%
RH(1) = 0;
V = zeros(size(t));     % 体积，m^3
V(1) = V0;
M = zeros(size(t));     % 质量，kg
M(1) = rho*V0;
Ts = zeros(size(t));    % 表面温度，K
Ts(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;

需要注意的是，该程序中的计算公式和参数均为近似值，实际应用时需要根据具体情况进行调整。