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

度过无风

%微分方程组求解主程序
# I' O- J$ V: M0 ?, g9 N5 Wclc;clear all;clf;close all;
%微分方程组求解主程序
clc;clear all;clf;close all;
8 ^( R3 _# ?& N/ o/ O9 b) {5 S%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
: _/ G& I/ N6 A! w. M7 }+ H5 p[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
1 j) N/ o6 c6 Y1 t) Tt=tt(191:end,
9 ?3 D- i. T: h1 Xys=(191:end,1)
7 w  x9 K/ q+ p9 T1 y. }yp=(191:end,2)
yh=(191:end,3)
4 n7 B3 @; i6 f1 i) S1 x& c9 |yr=(191:end,4)
6 b& h8 Q7 d" w- N" T# f4 Oyo=(191:end,5)
0 |% P6 G3 j6 g3 s! Mye=(191:end,6)
. R5 y0 b! `% y9 [, qT2=clock;
9 h# y# c0 ~3 C. ?. T3 AAPI_elapsed_time=T2-T1;
. J! O" E) T6 [6 c8 xif API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
9 L) z) ~8 c) J. z5 @if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
end
0 a& S! s8 u3 g' y+ W6 nstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
4 {: X8 \1 |6 W8 k2 P6 |0 @3 N%子函数：微分方程组odefile文件并命名为lbwfun.m
( i' d8 S- R: X% |function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
- k) e. L& B  l/ I) ?! R# K0.48*y(2)-17.87*y(3);
3 {+ ^. `" m! S7 D" F; @- |0 y4.85*y(3)-4.65*y(4);
$ _- R9 n" Q  Z( B. `2.55*y(2)+6.12*y(3)+1.95*y(4);
& w" a+ K8 X1 n1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
( _* X: L. Q9 k& \计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
4 h* V& P. j3 w??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.
2 t/ t. p. r# i8 \0 Z5 G4 n
Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
: e; `& C- ~) L7 x
( U7 z, Z* o% R" u3 H" z9 C4 nError in ==> ode45 at 173
1 q& z  E5 D4 p2 i/ a/ \[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...