建立数学模型描述恒定温度20摄氏度下从水滴落入密闭的边长为10厘米的正方体容器开始到最后完全蒸发各个时刻水滴周围空气中的湿度变化规律，并给出matlab代码

在水滴落入容器时，水滴周围的空气是干燥的，因此初始时刻空气中的水分含量为0。随着水滴蒸发，容器内的湿度会不断增加，进而影响水蒸气压和蒸发速率。因此，需要同时考虑水滴体积和周围空气中的湿度变化。

假设水滴的初始体积为$V_0$，蒸发速率为$M$，则水滴体积的变化率为$\frac{dV}{dt}=-\frac{4\pi r^2 k(P_s-P)}{\rho}$，其中$\rho$为水的密度，$r$为水滴半径，$P_s$为水的饱和蒸气压，$P$为当前空气中的水蒸气压。空气中的湿度变化率可以表示为$\frac{dP}{dt}=\frac{M}{V}-\frac{P}{\tau}$，其中$\tau$为空气中水蒸气的扩散时间常数。

综上所述，可以得到如下的微分方程组：

$$
\begin{cases}
\frac{dV}{dt}=-\frac{4\pi r^2 k(P_s-P)}{\rho} \\
\frac{dP}{dt}=\frac{M}{V}-\frac{P}{\tau}
\end{cases}
$$

其中，$k$为蒸发系数，$P_s$为水的饱和蒸气压，$P$为当前空气中的水蒸气压，$\rho$为水的密度，$r$为水滴半径，$\tau$为空气中水蒸气的扩散时间常数。

根据上述微分方程组，我们可以使用MATLAB进行数值求解。具体代码如下（假设容器边长为10厘米，初始时刻水滴体积为1毫升）：

% 模型参数
r = 0.28; % 水滴半径，单位：厘米
rho = 1000; % 水的密度，单位：千克/立方米
T = 20; % 恒定温度，单位：摄氏度
V0 = 1; % 初始时刻水滴体积，单位：毫升
P0 = 0; % 初始时刻容器内的水分含量，单位：毫克
P = P0 / (18 * 10^3); % 初始时刻空气中的水蒸气压，单位：mol/L
P_s = exp(20.386-(5132/(T+273.15))); % 水的饱和蒸气压，单位：帕斯卡
k = 2.72e-5; % 蒸发系数，单位：mol/(cm^2 s)
tau = 100; % 水蒸气扩散时间常数，单位：秒

% 数值求解
tspan = [0, 100]; % 时间范围，单位：秒
options = odeset('RelTol',1e-6,'AbsTol',1e-8); % 求解选项
y0 = [V0; P]; % 初始时刻状态
[time, Y] = ode45(@(t, Y) func(t, Y, r, rho, P_s, k, tau), tspan, y0, options);

% 绘图
figure;
subplot(2, 1, 1);
plot(time, Y(:, 1), 'LineWidth', 2);
xlabel('时间（秒）');
ylabel('水滴体积（毫升）');
title('水滴蒸发过程中体积的变化');
grid on;
subplot(2, 1, 2);
plot(time, Y(:, 2) * 18 * 10^3, 'LineWidth', 2);
xlabel('时间（秒）');
ylabel('空气中的水分含量（毫克）');
title('水滴蒸发过程中空气中的湿度变化');
grid on;

% 方程
function dydt = func(t, y, r, rho, P_s, k, tau)
    V = y(1); % 水滴体积
    P = y(2); % 空气中的水蒸气压
    A = 4 * pi * r^2; % 水滴表面积
    M = k * A * (P_s - P); % 蒸发速率
    dVdt = -M * 10^(-3) / rho; % 体积变化率
    dPdt = M / V - P / tau; % 湿度变化率
    dydt = [dVdt; dPdt];
end

运行以上代码，就可以得到水滴蒸发过程中水滴体积和周围空气中的湿度变化图像。需要注意的是，以上模型和代码仅供参考，实际情况可能会更加复杂，需要根据具体情况进行调整。