23、聚合物基质中药物分布的MOL解法详解

聚合物基质中药物分布的MOL解法详解

1. 变量与参数概述

在研究聚合物基质（PM）中的药物分布时，涉及多个变量和参数，如下表所示：
| 变量与参数 | 解释 | 代表单位 |
| — | — | — |
| u1, u2 | 水和药物的浓度（归一化） | - |
| r | PM中的径向位置 | cm |
| z | PM中的轴向位置 | cm |
| t | 时间 | s |
| r0 | PM的半径 | cm |
| zL | PM的长度 | cm |
| u10, u20 | u1和u2的初始值（归一化） | - |
| u1e, u2e | u1和u2的外部值（归一化） | - |
| Du1, Du2 | 扩散系数 | cm²/s |
| ku1, ku2 | 传质系数 | cm/s |

这里假设PM的扩散系数和尺寸是恒定的，但实际上，PM的孔隙率和尺寸可能随时间变化。不过，为了保持模型的基本性和适度复杂性，暂未考虑这些因素，重点在于介绍MOL（Method of Lines）解法。

2. 主程序分析

以下是用于求解方程（7.2） - （7.7）的MOL解法的主程序：

%
% Dynamic analysis of drug delivery
%
% Clear previous files
clear all
clc
%
% Global area
%
% Global area
global
nr
nz
dr
dz
drs
dzs...
r