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

度过无风

%微分方程组求解主程序
- Y( ]; k- W) B1 u5 }4 {clc;clear all;clf;close all;
9 f: M! m- F; z/ a. G( @5 K%微分方程组求解主程序
clc;clear all;clf;close all;
5 a- a0 m" I; t, j%windows时钟自动计时
1 M2 L; g. F  m9 aT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
2 h7 v; c/ E' G9 z# C[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
7 t3 P% w5 T( f7 K" j# Bt=tt(191:end,
6 y, A& y  ?4 Y( f5 fys=(191:end,1)
  e8 F" r9 P: ?+ i& D8 {% iyp=(191:end,2)
! c9 F/ f0 f' d9 f2 H/ L) yyh=(191:end,3)
3 `% y# ^" x/ a8 R) T& Iyr=(191:end,4)
8 P: T% ]% I- fyo=(191:end,5)
* g1 B" \2 S* q+ a/ u! A4 [1 Uye=(191:end,6)
T2=clock;
5 U: z% V6 g6 H1 U. Q% H7 a: j* YAPI_elapsed_time=T2-T1;
$ j6 ^, c8 Z: xif API_elapsed_time(6)<0
; ^5 f" ~0 l* _( s" K& v' FAPI_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;
4 l& g" ]$ c1 D$ uend
! I" S4 N& O' m# }" j: _, y' O/ P) ystr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
& |" y/ {; K  m! YAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
$ m3 m+ S( I/ |1 zdisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
+ ?% a- @- R) F; ~! ofunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
4 O) W( U7 ^" py(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
; f# k' _. K# t& J: M4.85*y(3)-4.65*y(4);
2.55*y(2)+6.12*y(3)+1.95*y(4);
1 i2 K5 ?# S4 f, c1 i1.10*y(2)+6.90*y(3)+2.70*y(4)];
" o' ^9 m8 a$ Z5 l7 `, ?& g4 H1 S运行结果（我是按书上写出来的程序）
. F) N' h8 D: Q/ N7 ^0 l计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
" g  i2 u' `8 p0 s2 w: D8 T6 i??? Error using ==> feval
" P# S2 t+ o* c( ^/ Q6 CUndefined function or method 'lbwfun' for input arguments of type 'double'.
" ]2 ]/ ^# J6 ^, A( T9 A  D4 V7 R
% i3 }- q" v- ?0 E: fError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

3 j9 S! n5 j- O! v- o6 VError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...