用matlab画 ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/pi*sin(n*pi*(R-r+th)/th)*exp(-(n*pi/th).^2*D*tci=1;cl=0.2;D=1.12x10.^(-9);R=0.001;th=0.0004;[r,t]=meshgrid(0.001:0.0001:0.0014,0:0.001:1);for n=1:100000ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/p
来源：学生作业帮助网 编辑：六六作业网 时间：2024/11/17 06:27:28
用matlab画 ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/pi*sin(n*pi*(R-r+th)/th)*exp(-(n*pi/th).^2*D*tci=1;cl=0.2;D=1.12x10.^(-9);R=0.001;th=0.0004;[r,t]=meshgrid(0.001:0.0001:0.0014,0:0.001:1);for n=1:100000ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/p
用matlab画 ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/pi*sin(n*pi*(R-r+th)/th)*exp(-(n*pi/th).^2*D*t
ci=1;
cl=0.2;
D=1.12x10.^(-9);
R=0.001;
th=0.0004;
[r,t]=meshgrid(0.001:0.0001:0.0014,0:0.001:1);
for n=1:100000
ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/pi*sin(n*pi*(R-r+th)/th)*exp(-(n*pi/th).^2*D*t));
end
mesh(ca);
这个程序不对
其实可以这样,只要求出比如当t=0.001,ca关于r的函数图就好,二维的;n是从1到无穷的求和,那个图片中有原函数
方程和希望得到的 图如下
用matlab画 ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/pi*sin(n*pi*(R-r+th)/th)*exp(-(n*pi/th).^2*D*tci=1;cl=0.2;D=1.12x10.^(-9);R=0.001;th=0.0004;[r,t]=meshgrid(0.001:0.0001:0.0014,0:0.001:1);for n=1:100000ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/p
感觉r和t有点问题.其一是r和t的步数.r的步数仅为5,而t的步数为1001,这样不合理.另外,r和t的范围可能存在点问题,不过这个我只是猜测,因为我画的时候发现,t>0.0012后就恒定了.不多说了.把代码和图发给你（我把r的步数改成了41,把t的步数改成了101）：
ci=1;
cl=0.2;
D=1.12*10.^(-9);
R=0.001;
th=0.0004;
rrange=0.001:0.00001:0.0014;
trange=0:0.01:1;
[r,t]=meshgrid(rrange,trange);
temp1=zeros(size(t));
for n=1:200000
temp1=temp1+(-1)^n*2/n/pi*sin(n*pi*(R-r+th)/th).*exp(-(n*pi/th)^2*D*t);
end
ca=cl+(ci-cl)*R./r.*((R-r+th)/th+temp1);
surf(r,t,ca)
