用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=0.000001;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/pi*sin
来源：学生作业帮助网 编辑：六六作业网 时间：2024/11/17 08:52:25
用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=0.000001;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/pi*sin
用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=0.000001;
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到无穷的求和,那个图片中有原函数；（那个T其实是t,打错了）
用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=0.000001;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/pi*sin
ci=1;
cl=0.2;
D=1.12*10.^(-9);
R=0.001;
th=0.0004;
r=0.001:0.00001:0.0014;
t=0.01;
temp1=zeros(size(r));
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);
plot(r,ca)
只针对一个变量r的代码