湖泊有机物新陈代谢模型问题

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
%微分方程组求解主程序
clc;clear all;clf;close all;
%windows时钟自动计时
: o% [/ k9 r, B# y, V! sT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
% Y5 i$ i- F8 m& S! d* R3 T[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,
9 K4 {  I1 A6 g) wys=(191:end,1)
/ r5 H5 Z' R, q3 d5 \yp=(191:end,2)
& n# M$ S9 G& J- D% c' z& ~yh=(191:end,3)
yr=(191:end,4)
yo=(191:end,5)
$ t' U4 H# @( ]$ @ye=(191:end,6)
# b: j8 H$ D. l% d1 y$ iT2=clock;
6 |$ `8 r* k" t. i8 m; j& bAPI_elapsed_time=T2-T1;
# ]/ w+ ~$ d2 O* m6 o  \if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
$ l2 J# R% S! V& `' e' K* x2 DAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
) G7 s5 L- m+ A. P* I* q$ r5 xAPI_elapsed_time(4)=API_elapsed_time(4)-1;
end
8 Q, k: U* J7 t6 m* M/ o- s$ r+ {* Bstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
9 c( k5 m% Q' O7 IAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
0 s1 e6 y) K2 i9 a0 z" Q  }) l5 ]disp(str);
" h" R; {5 W+ e1 c0 o4 L%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
" {& O5 _, }8 C! A' [+ q- `3 Vydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
4.85*y(3)-4.65*y(4);
2.55*y(2)+6.12*y(3)+1.95*y(4);
2 G. j: w* \1 q# l* P1.10*y(2)+6.90*y(3)+2.70*y(4)];
$ d/ W8 ?4 {8 F; L( n7 S3 {运行结果（我是按书上写出来的程序）
, b$ m+ U* ]2 m* k计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
  T/ }/ a" K) _??? Error using ==> feval
, p8 T  w) X5 U( Z: xUndefined function or method 'lbwfun' for input arguments of type 'double'.
) f# A' `1 P# v1 Y9 [9 s- J4 y! U
. R; ]" H8 n2 |# O: ]( v5 ^7 c7 hError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

7 T: I+ e3 ?3 i  _Error in ==> ode45 at 173
3 {, _  S8 b+ `" y7 D[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...