MATLAB 从 dsolve 公式中找到最大值

Question

我正在求解两个微分方程并绘制第二个微分方程的结果。我想找到图中最大点的值。但是，也许是因为我使用的是 sym，我不知道该怎么做。我知道在数学中，您只需找到一阶导数并将其设置为零即可。

% Constants
k1 = 1/10; 
k2 = 1/1.3; 
k3 = 1/1.3; 
k4 = 1/(1/48);
k5 = 1/3.9; 
t = 1:5;
% Initial concentration
A0 = 5891.694;
B0 = 0; 
%%% Mass balance equation 1
syms A(t) t
eqn = diff(A,t) == -k1*A -k2*A -k3*A -k4*A;
cond = A(0) == A0;
A(t) = dsolve(eqn,cond);
% A(t) = exp(-k1*t);
%%% Mass balance equation 2
syms B(t) t
eqn = diff(B,t) == k1*A +k2*A +k3*A +k4*A - k5*B;
cond = B(0) == B0;
B(t) = dsolve(eqn,cond); 

%%% Plot
figure('visible','on');
t = 0:20;
plot(t,B(t))

所以对于上面的图片，我想找到最大点。有时，它不是精确的时间步长，所以我希望 MATLAB 只显示最大点的精确值。

Answer 1

您可以使用 findpeaks 来查找数据中的最大值。为此，首先您需要将您的 ode 解决方案（代数表达式）转换为特定学习间隔的数据列表。

B(t) = dsolve(eqn,cond,'t') ;
t = 0:0.1:20;
z=eval(B(t))

现在使用 findpeaks

[pks, locs] = findpeaks(z)

pks =

5.7310e+03

locs =
 2

MATLAB 从 dsolve 公式中找到最大值

MATLAB Find the max point from a dsolve formula

matlab

max

differential-equations