MATLAB根据复杂的公式u(z,t)编写程序绘图怎么搞

Question

有一个关于两个变量的函数，计算式复杂，u=c01+q-(c01+q-qexp(-rTv))exp(求和2/Msin(Mz/H)exp(-MMTv)Tm(t))
其中M=(2m-1)*pi/2,m=1,2,3,4等正整数;c01=40;q是一个关于t的循环荷载；详见程序

%%定义循环荷载
T=182;       %循环周期，单位：d
a1=0.3;      %梯形循环荷载加载系数
b1=1.4;      %梯形循环荷载加载系数
t0=130;      %梯形循环荷载时间参数，单位：d
N=1:4;       %循环周期，整数，可取1,2,3,4
%%循环荷载表达式，第一周期
t1=0:0.1:39;
q1=100/39*t1;
t2=39:0.1:91;
q2=0*t2+100;
t3=91:0.1:130;
q3=-100/39*(t3-130);
t4=130:0.1:182;
q4=0*t4;
t=[t1 t2 t3 t4];
q=[q1 q2 q3 q4];
plot(t,q);
axis([0 190 0 110]);
syms t;
dq=diff(q,t);
%%定义相关参数
Z=cv/H^2;
m=1:1000;
M=(2.*m-1).*2*pi;
Bm=cv.*M.*M./H^2;
R(t)=-dq./(c01+q)-(r.*q.*Z.*exp(-r.*t.*Z))./(c01+q-q.*exp(-r.*Z.*t));
ft=R(t).*exp(M.*M.*Z.*t);
syms t
Tm(t)=int(ft,t,0,t)
m=1:1000;
x=0
for i=1:m
    fx=2./M(i).*sin(M(i).*z/H).*exp(-M(i).*M(i).*Z.*t).*Tm(i);
    x=fx+x;
end
u=c01+q-(c01+q-q.*exp(-r.*Z.*t)).*exp(x);
%%绘制深度z=1.0m和z=3.0m处孔压随时间变化曲线图
Tv=0.01:0.1:1;
z=1.0;
plot(Tv,u);

Answer 1

是这个是什么意思呀？我没有看明白 嗯