9、抗体结合动力学的数值模拟与分析

抗体结合动力学的数值模拟与分析

1. 主要程序概述

在抗体结合动力学的研究中，我们使用了一个主要程序 pde_1_main 来进行数值模拟。这个程序涉及到多个步骤，包括清除先前文件、定义参数、设置初始条件、进行常微分方程（ODE）积分以及输出和绘制结果等。以下是该程序的详细步骤：
1. 清除先前文件并声明全局变量 ：

% Clear previous files
clear all
clc
% Parameters shared with the ODE routine
global zl zu z dz dz2 D kf cbsat kr n cbulk ndss ncall

这部分代码的作用是清除工作区中的所有变量和命令窗口中的内容，并声明一些全局变量，以便在其他函数中可以共享这些变量。
2. 定义参数数值 ：

% Parameter numerical values
D=1.0e-10; kf=1.0e+05; cbulk=4.48e-05;
h=5.0e-05; c0=0; cb0=0;

这里定义了一些在模型方程中出现的参数的数值，这些参数将在后续的计算中使用。
3. 设置接口结合饱和度和接口解绑速率的变化 ：

% Vari