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

度过无风

%微分方程组求解主程序
, _7 b$ w' {0 Y6 Z: kclc;clear all;clf;close all;
%微分方程组求解主程序
3 h6 x& D7 z+ aclc;clear all;clf;close all;
$ g5 D) H: C$ J4 P/ c%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
8 F1 Y# V' S4 U5 N3 a7 a3 X& `disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
3 r8 F2 a2 ^: L/ p& n! h, Y8 l[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
5 J/ y( @# m- Q/ v) Xt=tt(191:end,
ys=(191:end,1)
2 z' j" u# J( G$ p. d8 syp=(191:end,2)
' `6 i3 n' p' q" s5 ?yh=(191:end,3)
5 O# {, U" T: j# g4 Fyr=(191:end,4)
1 x8 A9 l9 d$ y; ayo=(191:end,5)
( L, g$ _9 ^- |: I+ zye=(191:end,6)
T2=clock;
API_elapsed_time=T2-T1;
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
- c' R4 d1 Q% G& N& Y' Tif API_elapsed_time(4)<0
3 L! O" G1 O3 eAPI_elapsed_time(4)=API_elapsed_time(4)+60;
/ f3 Q1 j% `4 e9 M0 S5 l2 v$ |API_elapsed_time(4)=API_elapsed_time(4)-1;
end
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
) U% W- q8 E. _/ v$ Q1 s1 VAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
6 @9 e0 J$ m- R; a4 Wydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
; d; M' {, g  q6 P4.85*y(3)-4.65*y(4);
2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
9 X3 ]( V5 s1 v- V8 l运行结果（我是按书上写出来的程序）
; Q1 Z9 s+ [' N" d6 y3 ?. V* o计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.
6 E5 L) r, _" Z  J
Error in ==> odearguments at 110
. g4 a2 Z2 W1 _/ y, v7 y2 F2 H5 Bf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
' C: k: r: Q( l6 ?# X- t
Error in ==> ode45 at 173
* O$ j) i7 x7 Z5 D# w& _" U[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...