图解法程序:
t=0:0.1:20;
g=0.1;
r=2;
Q(t)=(8-g*t).*(80+r*t)-4*t-640;
Plot(t,Q);
运行结果:
由图可知t=10时,函数取最大值,最大值约为20。
代数法程序:
syms t
r=2;g=0.1;
Q=(8-g*t)*(80+r*t)-4*t-640;
dQ=diff(Q,'t')
t=solve(dQ,t)
t=eval(t)
Q=eval(Q)
运行结果:
方法二与方法一比较:
与方法一都是在t=10时取得最大值,最大值Q=20。