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

度过无风

%微分方程组求解主程序
% y+ K& }% h  M! Qclc;clear all;clf;close all;
+ V7 [: k- i" T/ @% V%微分方程组求解主程序
clc;clear all;clf;close all;
) L, B: r9 v, E( w0 @%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
. N' B& ^  a1 G5 B3 M6 Rdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
: |4 n2 l& ~& i# {% D[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,
1 O, L9 N6 F/ Oys=(191:end,1)
yp=(191:end,2)
/ b1 C( j- R5 y( ?& myh=(191:end,3)
yr=(191:end,4)
3 J1 f1 b& q+ L/ S  [% Dyo=(191:end,5)
* Q* r) Q( v% N0 t' X9 Tye=(191:end,6)
T2=clock;
3 l& c# z& T/ R6 Y( d+ UAPI_elapsed_time=T2-T1;
( J& n) c7 Y- p* }! U  l# \if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
1 G: m- y# N5 ?7 U, `if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
% Q5 I+ t' _. _7 n% {: Rend
6 w5 I8 r& X. U$ xstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
. @3 q+ N+ p  kfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
3 G- m6 p* ^. M  z; K# Ly(1)-4.03*y(2);
5 f' L( O! w9 L9 M( `. k0.48*y(2)-17.87*y(3);
4.85*y(3)-4.65*y(4);
9 \2 ~7 q- U: x# [3 K2.55*y(2)+6.12*y(3)+1.95*y(4);
6 H. w. f0 }0 G1 C4 o% U. l# j  g1.10*y(2)+6.90*y(3)+2.70*y(4)];
5 a( ?) m: J$ P6 L% k1 ]4 L运行结果（我是按书上写出来的程序）
6 `) l1 n. j8 @# J计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
8 k2 E- A" Z1 \! l& b1 c??? Error using ==> feval
: R: O5 {" o+ S7 EUndefined function or method 'lbwfun' for input arguments of type 'double'.
- |& K" ?) x& H" F: W
  |; i3 V+ @, Z  ]( D5 sError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
8 S3 Z6 q3 V, e& N[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...