clearclcr=0.1;s=0.01;T=0:0.5:20;w=0;f0=(-tanh(r)^0./cosh(r))

clear
clc
r=0.1;
s=0.01;
T=0:0.5:20;
w=0;
f0=(-tanh(r)^0./cosh(r));
f1=(-tanh(r)^1./cosh(r));
B0=(f0./4)*s*sqrt((0+1)*(0+2)).*(1-cos(2.*T/s));
B1=(f1./4)*s*sqrt((0+2)*(0+3)).*(1-cos(2.*T/s));
C0=(f0./2)*s*(0+1)*sin(2.*T/s);
C1=(f1./2)*s*(0+2)*sin(2.*T/s);
w=0.5*(0*(f0^2)+(0+1).*(C0^2)+(2*0+2).*(B0^2)+(0+1).*f0.*f1+C0.*C1.*(0+2)+2*sqrt((0+1)*(0+3)).*B0.*B1);
for n=1:1000;
f(n)=((-tanh(r))^n)./cosh(r);
f(n+1)=((-tanh(r))^(n+1))./cosh(r);
B(n)=(f(n)./4)*s*sqrt((n+1)*(n+2)).*(1-cos(2.*T/s));
B(n+1)=(f(n+1)./4)*s*sqrt((n+2)*(n+3)).*(1-cos(2.*T/s));
C(n)=(f(n)./2)*s*(n+1)*sin(2.*T/s);
C(n+1)=(f(n+1)./2)*s*(n+2)*sin(2.*T/s);
W(n)=0.5*(n.*(f(n)^2)+(n+1).*(C(n)^2)+(2*n+2).*(B(n)^2)+(n+1).*f(n).*f(n+1)+C(n).*C(n+1).*(n+2)+2*sqrt((n+1)*(n+3)).*B(n).*B(n+1));
w=w+W(n)
end
T=0:0.5:20;
plot(T,w)
这个程序的主要目的就时要做出w随T变化的图,w是一个含有T的求和公式,n从0到1000求和 ,0《=T《=20,步长0.5