有机物新陈代谢模拟程序

有机物新陈代谢模拟程序

下面来介绍一下利用MATLAB来模拟自然湖泊中有机物的新陈代谢的过程，具体如下：

1、首先在MATLAB主界面的编辑器中写入下列代码：

%微分方程组求解主程序
clc;clear all;clf;close all;
%Windows时钟自动计时
T1=clock;%Clock函数返回的值是clock = [year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待……');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,:)
ys=y(191:end,1)
yp=y(191:end,2)
yh=y(191:end,3)
yr=y(191:end,4)
yo=y(191:end,5)
ye=y(191:end,6)
T2=clock;
API_elapsed_time=T2-T1;
if  API_elapsed_time(6)<0API_elapsed_time(6)=API_elapsed_time(6)+60;API_elapsed_time(5)=API_elapsed_time(5)-1;
end
if  API_elapsed_time(5)<0API_elapsed_time(5)=API_elapsed_time(5)+60;API_elapsed_time(4)=API_elapsed_time(4)-1;
end
if  API_elapsed_time(4)<0API_elapsed_time(4)=API_elapsed_time(4)+60;API_elapsed_time(3)=API_elap