用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=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的代码
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的代码
组合恒等式的证明:C(r,r)+C(r+1,r)+C(r+2,r)+…+C(n,r)=C(n+1,r+1) C(n,1)
证明C(r+1,n)+ 2C(r,n)+C(r-1,n) = C(r+1,n+2)
matlab求反函数T=(K*m*p*U^2*R*((n0-n)/n0))/(2*pi*f*(R^2+((n0-n)/n
n=3r. 又有[(3r)(3r-1)(3r-2)……(2r+1)]/(r!)2^r=60 怎么解得?
当2=4)时,证明C(n,r)=C(n-2,r-2)+2C(n-2,r-1)+C(n-2,r)
设A为n阶(n≥2)方阵,证明r(A*)= n ,r(A)=n r(A*)= 1,r(A)=n-1 r(A*)= 0,r
C(0,n)+2C(1,n)+3C(2,n)+...+(r+1)C(r,n)+...+(n+1)C(n,n)=___(n
设A为n阶矩阵,证明r(A^n)=r(A^(n+1))
如果 r满足 r + 1/r 是整数,那么证明 r^n + (1/r^n) 也是整数
C(m,n+r+1)=C(m,n)C(0,r)+C(m-1,n-1)C(1,r+1)+...+C(0,n-m)C(m,r
线性代数 证明题1.设A,B,C,D都是n阶矩阵,r(CA+DB)=n (1)证明:r( A )( B )=n (A,B
已知等比数列an,m=(sn-r,2^n-1),N=(2.1)且向量m//n求常数r