代码如下:
%给参数赋值,模拟一个实际场景
t=0:1:40;alpha=0.01;l=3;S=1.5*10^5;r=10^3;beta=2*10^3;d=0.6;mu=10;gama=600;
s=100;a=10^5;
t1=(1/(alpha*l))*log((((S*(r+beta*d)+r*mu*l-gama*s*beta*d)/(a*d*beta)))); %计算达到运力峰值的时间
syms x; %定义变量
f = (a*beta*d*exp(alpha*x*l)-r*mu*l+gama*s*beta*d)./(r + beta*d)-S;%定义函数
d1=((a*beta*d*exp(alpha*t*l)-r*mu*l+gama*s*beta*d)./(r + beta*d)).*(t<t1)+...
((a*beta*d*exp(alpha*t*l)-r*mu*l+gama*s*beta*d)./(r + beta*d)+double(int(f,x,t1,t))).*(t>=t1); %计算平台需求量
报错见下,为int函数求解定积分错了:
所欲求解的积分如下表达式如下: