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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
0 I1 y" ]) R9 u* l- d* D( d3 E' a%微分方程组求解主程序
" H4 N7 `. @& Pclc;clear all;clf;close all;
+ o0 v( s" M. R+ P* r$ n8 a%windows时钟自动计时
- A& Z. _/ W; x7 s; D9 TT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
/ G9 ^8 p5 [# a1 |+ Q' Fdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
; R4 j- @+ k5 W( j3 D1 i6 V[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
; T8 X1 A7 Y8 q8 b& E: I9 Lt=tt(191:end,
0 g& y" B% Y* t' a. T; jys=(191:end,1)
& Q1 l/ W1 V. j1 g# }$ gyp=(191:end,2)
, X6 A* m# ~: Zyh=(191:end,3)
yr=(191:end,4)
- ?9 E) |) n; p, v- gyo=(191:end,5)
ye=(191:end,6)
, G  N; p' @' `5 \8 W  aT2=clock;
/ _6 f. U" o7 X3 V0 T0 _. _4 yAPI_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
1 {3 R: E+ v: ~' R% M# p: ^  [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;
2 ?( u. s5 g0 U* [API_elapsed_time(4)=API_elapsed_time(4)-1;
end
0 t; ?+ w2 a7 [6 Y; @( ~+ f- ^0 n9 Xstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
" L% c- [4 L; F' ~* a7 VAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
5 d/ X; y9 ~! P5 v. D" Sdisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
) T( e' K( _! P( |1 n) U- f: kfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
. u# t- v+ b, r1 W) ^y(1)-4.03*y(2);
  v# s" {6 g, F* U0.48*y(2)-17.87*y(3);
4.85*y(3)-4.65*y(4);
+ z1 ^0 e( d2 w% G# M2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
  p5 h3 L" C. q' b9 X/ [运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
3 }* E0 v! U# FUndefined function or method 'lbwfun' for input arguments of type 'double'.

Error in ==> odearguments at 110
( A. B& c' Q) F/ E  j/ _6 S# E- Yf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...