20、血液透析器动力学与表皮伤口愈合模型解析

血液透析器动力学与表皮伤口愈合模型解析

1. 血液透析器动力学模型

1.1 MOL 例程

在血液透析器动力学模型中，采用了四阶龙格 - 库塔积分方法：

u0=u; t0=t;
[u,t]=rk4(u0,t0,h,nsteps);

这里调用了 ODE 例程 pde_1 进行计算。之后会显示并存储解，用于后续绘图。以下是部分代码实现：

if(ip==1)
    fprintf('%5.2f%10.4f%10.4f\n',t/3600,u(2*n+1),u(n+1));
end
u1plot(iout)=u(2*n+1);
u2plot(iout)=u(n+1);
tplot(iout)=t/3600;

最后，会绘制出终端血液和透析液浓度随时间变化的曲线：

figure(ncase);
plot(tplot,u1plot,'-o');
axis([0 2 -0.1 1.1]);
ylabel('u1R(t),u2(0,t)');xlabel('t');
if(ncase==1) title('ncase = 1; u1 - o; u2 - x'); end
if(ncase==2) title('ncase = 2; u1 - o; u2 - x'); end
hold on
plot(tplot,u2plot,'-x');