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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
4 S. G+ o4 d/ e7 a  }+ _, X%微分方程组求解主程序
* x! h7 E. S3 n/ X/ _3 vclc;clear all;clf;close all;
* `4 |6 N0 e( |* G  c: R  E%windows时钟自动计时
3 U# B" L$ m, q; i! p! j! tT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
6 {0 |, u4 Q$ h0 r8 E" J2 f/ pdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,
ys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
* X: f( Z& X* F" ^yr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
2 y  r& ]: ^/ E: J7 o5 i' ZT2=clock;
# N+ ~4 T( X* G6 d, M' zAPI_elapsed_time=T2-T1;
: ]4 Z! j) f8 gif API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
) d+ i: Y8 m% A6 Tend
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
/ k% v3 x1 L# C0 T  tAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
6 T: s0 q2 P  \7 @- ~8 h& N$ \. Ddisp(str);
  I, b8 ^+ N' y) m1 s! k- @%子函数：微分方程组odefile文件并命名为lbwfun.m
: L( L9 w2 S5 Nfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
; v* X1 t! `3 Z8 _( c  ny(1)-4.03*y(2);
6 l+ }- q0 L/ `; g+ F0.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);
0 Z0 e" _- w% b, Q  S/ F1.10*y(2)+6.90*y(3)+2.70*y(4)];
5 P8 {: l: e- r, A0 M1 ~7 l+ E/ f运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
6 k4 `0 q9 m) E/ [" |' LUndefined function or method 'lbwfun' for input arguments of type 'double'.
3 `' I4 g/ h$ T! h+ ?5 g9 Y
9 d, X& ^; q) R+ X9 o, c+ X( H. pError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
9 v: }8 l* D- F% v$ n
3 N( ~& P) B# Y0 oError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...