| 查看: 268 | 回复: 2 | |||
| 当前主题已经存档。 | |||
lxx244lxx木虫 (小有名气)
虫
|
[交流]
【求助】谁帮我看下matlab的程序。。。谢谢了
|
||
|
n=5+1;m=100; dr=0.01; dz=0.0585; Ramda=0.45; G=2500;Cp=2.18; dH=140e3;rh0=1440;a=0.05; u0c0=0.069/(pi*a^2); Cp1=1.0e3; F=130/3600; a1=Ramda/(G*Cp*1000)*3600; b1=dH*rh0/G/Cp; a2=0.000427; b2=rh0/u0c0; C1=2*pi*a*Ramda/(F*Cp1); C2=4*dr/(dz*C1); M=a1*dz/dr^2; M1=a2*dz/dr^2; C3=M*C2/2*(1+0.5/(n-1))-(1+M); C4=M*C2/2*(1+0.5/(n-1))-(1-M); T(1:n-1)=873; T(n)=893; x(1:n)=0; X0=[T;x]; Tresult(1,  =X0(1,; xresult(1,=X0(2,;FGrc(T,x); for j=2:m X=fsolve(@TxEquations,X0); T=X(1, ; x=X(2,; Tresult(j,=T, xresult(j,=xxa(j, =sum(2*xresult(j,2)+4*xresult(j,3)+6*xresult(i,4)+8*xresult(j,5)+5*xresult(j,6))/25if xa(j)>0.45, break, end FGrc(T,x) end z=dz*[0:j-1];r=dr*[0:n-1]; z=z',L=spline(xa,z,0.45) plot(z,xa), xlabel('z(m)'), ylabel('x_a_v') figure, surf(r,z,Tresult), xlabel('r(m)'), ylabel('z(m)'), zlabel('T(K)') figure, surf(r,z,xresutt), xlabel('r(m)'), ylabel('z(m)'), zlabel('x') function f=ReactionRate(T,x) k=0.027*exp(0.021*(T-773)); f=15100*exp(-11000/T)*((1-x)/(11+x)-1.2*x^2/k/(11+x)^2); function f=FGrc(T,x); global F G rc n dz M M1 b1 b2 C2 C3 C4 rc=ReactionRate(T,x); F(1)=((1-2*M)*T(1)+2*M*T(2)-b1*dz/2*rc(1))/(1+2*M); G(1)=(1-2*M1)*x(1)+2*M1*x(2)+b2*dz/2*rc(1))/(1+2*M1); i=(2:5); var1=1-0.5/(i-1); var2=1+0.5/(i-1); F(i)=(M/2*(varl*T(i-1)+var2*T(i+1))+(1-M)*T(i)-b1*dz/2*rc(i))/(M+1); G(i)=(M1/2*(var1*x(i-1)+var2*x(i+1))+(1-M1)*x(i)+b2*dz/2*rc(i))/(M1+1) F(n)=(-M*T(n-1)+C4*T(n)+b1*dz/2*rc(n)/C3; G(n)=(M1*x(n-1)+(1-M1)*x(n)+b2*dz/2*rc(n))/(1+M1); function f=TxEquations(X) global n F G rc dz M M1 b1 b2 C2 C3 T=X(l, , x=X(2,fT(1)=F(1)+(2*M*T(2)-b1/2*dz*rc(1))/(l+2*M)-T(1); fx(1)=G(1)+(2*M1*x(2)+b2/2*dz*rc(1))/(1+2*M1)-x(1); for i=2:n-1 var1(i)=(1-0.5/(i-1)); var2(i)=(1+0.5/(i-1)); fT(i)=F(i)+(M/2*(var1(i)*T(i-1)+var2(i)*T(i+1))-b1/2*dz*rc(i))/(M+1)-T(i); fx(i)=G(i)+(M1/2*(var1(i)*x(i-1)+var2(i)*x(i+1))+b2/2*dz*rc(i))/(M1+1)-x(i);% (21) end fT(n)=F(n)+(-M*T(n-1)+b1/2*dz*rc(n))/C3-T(n); fx(n)=G(n)+(M1*x(n-1)+b2/2*dz*rc(n))/(1+M1)-x(n); f=[fT;fx]; 不知道为什么老是说 Undefined command/function 'FGrc'，，，请教各位大神了。。。 [ Last edited by nono2009 on 2009-12-8 at 11:24 ] |
回复此楼» 猜你喜欢
数一英一271专硕（085401）求调剂，可跨
已经有8人回复
-
295求调剂
已经有4人回复
-
0856求调剂
已经有5人回复
-
材料学硕333求调剂
已经有10人回复
-
299求调剂
已经有6人回复
-
一志愿郑州大学，080500学硕，总分317分求调剂
已经有6人回复
-
289求调剂
已经有12人回复
-
070305高分子化学与物理 304分求调剂
已经有6人回复
-
0703化学调剂，求导师收
已经有10人回复
-
各位老师好，我的一志愿为北京科技大学085601材料专硕
已经有5人回复
dic213
木虫 (著名写手)
第三军团儿童团团长
- 应助: 0 (幼儿园)
- 金币: 1777.4
- 散金: 5
- 帖子: 2476
- 在线: 37.1小时
- 虫号: 42931
- 注册: 2004-04-03
- 性别: GG
- 专业: 计算数学
2楼2009-12-08 14:45:45
3楼2009-12-11 18:31:36
6