作业帮 > 综合 > 作业

用matlab画 ca=cl+(ci-cl)*R/r*((R-r)/th+1+(-1).^n*2/n/pi*sin(n*

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/05/20 21:36:12
用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,打错了)
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的代码