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

度过无风

%微分方程组求解主程序
/ v' P7 v5 d- Rclc;clear all;clf;close all;
%微分方程组求解主程序
$ m  Y9 l( X4 }clc;clear all;clf;close all;
3 X. g- P7 A7 t+ W! X%windows时钟自动计时
4 |% I- H: ~( C; X3 v3 [' i* kT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
& t  X" ?+ ~% W4 Y- n, Ydisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
- M2 t5 C) l! U, K[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
$ M* U$ C1 i. ht=tt(191:end,
ys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
yr=(191:end,4)
0 H6 T2 A7 {- N1 E, Zyo=(191:end,5)
; {  B. u% S- Jye=(191:end,6)
/ I! T* m! x$ i) dT2=clock;
API_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
- _1 M/ m! Z) e0 ^8 {1 ~- `# yAPI_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
! ]1 n! J8 K4 z6 ^+ M+ Bif API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
% _. T- D/ N6 j$ [1 B4 K& A- WAPI_elapsed_time(4)=API_elapsed_time(4)-1;
end
* w! J1 G1 x& }/ k4 W) M4 C: estr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
" H4 M: Z. e3 ^" Gfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
7 q' k3 g0 q( y9 f; s6 L0.48*y(2)-17.87*y(3);
% V/ w9 F, S. k7 J6 S* r* z4.85*y(3)-4.65*y(4);
- a1 \/ Q# g- @- {) z2.55*y(2)+6.12*y(3)+1.95*y(4);
6 y  Y2 T0 ?( t4 {( w1.10*y(2)+6.90*y(3)+2.70*y(4)];
( w" V: Z; ]- R: I运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
( \6 I- ^" x3 k??? Error using ==> feval
/ ~, r% H0 Y, L, l9 V9 p. E1 Q% {Undefined function or method 'lbwfun' for input arguments of type 'double'.
! Q; u6 U1 }5 h( }% Z6 s4 `  s
Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
, ]" G5 Y$ z/ E. b* Y, }( z
Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...