# 7 麦片

050

1主题 0最佳答案

 这是程序 clear all clc h=0.01; NTU1=1.9807; NTU2=1.9807; NTU3=0.2957; NTU4=2.8605; NTU5=0.7923; NTU6=4.7785; NTU7=4.7785; NTU8=289.4532; NTU9=72.1932; NTU10=313.5230; NTU11=741.4073; a=735.9435; b=179656.8108; c=0.0012; d=-0.0107; x=0:h:1;y=0:h:1; m=length(x);n=length(y);M=round(n/3); tp=zeros(n,m);tl=zeros(n,m);tw=zeros(n,m);ts=zeros(n,m); wp=zeros(n,m);wl=zeros(n,m);wtw=zeros(n,m);ws=zeros(n,m); ql=zeros(n,m);tpl=zeros(n,m); tp(:,m)=34;tl(:,m)=26;tw(n,:)=20;ts(n,:)=25;wp(:,m)=22;wl(:,m)=8.35; wtw(n,:)=10;ws(n,:)=9.88;ql(:,m)=0.4;tpl(:,m)=20; E=c*det(tpl)+d; m0=-7.8582;m1=1.8399;m2=-11.7811;m3=22.6705;m4=-15.9393;m5=1.7752; n0=0.28;n1=4.30;n2=0.60;n3=0.21;n4=5.10;n5=0.49;n6=0.362;n7=-4.75;n8=-0.40;n9=0.03; W=det(tl)+273.15; V=1-W/647.3; P=22090*exp((m0*V+m1*V^1.5+m2*V^3+m3*V^3.5+m4*V^4+m5*V^7.5)/(1-V)); A=2-(1+(det(tw)/n0)^n1)^n2; B=(1+(det(tw)/n3)^n4)^n5-1; Y=1-((1+det(tw)/n6)^n7)^n8-n9/exp((det(tw)-0.1)^2/0.005); Pa=Y*(A+B*W/647.3)*P; C=0.622*Pa/(101.325-Pa);     for i=2:n-1         for j=2:m-1           tp(i+1,j)=tp(i-1,j)-2*h*NTU1*(tp(i,j)-tl(i,j));           wp(i+1,j)=wp(i-1,j)-2*h*NTU2*(wp(i,j)-wl(i,j));           tl(i+1,j)=tl(i-1,j)-2*h*NTU3*(tp(i,j)-tl(i,j))+2*h*a*(wp(i,j)-wl(i,j))-2*h*NTU4*(tl(i,j)-tpl(i,j));           ql(i+1,j)=ql(i-1,j)-2*h*NTU5*(wp(i,j)-wl(i,j));           ts(i+1,j)=ts(i-1,j)-2*h*NTU6*(ts(i,j)-tw(i,j));           ws(i+1,j)=ws(i-1,j)-2*h*NTU7*(wtw(i,j)-ws(i,j));           tw(i+1,j)=tw(i-1,j)-2*h*NTU8*(ts(i,j)-tw(i,j))-2*h*b*(wtw(i,j)-ws(i,j))+2*h*NTU9*(tw(i,j)-tpl(i,j));           tpl(i+1,j)=4*tpl(i,j)-tpl(i-1,j)-tpl(i,j-1)+h^2*NTU10*(tw(i,j)-tpl(i,j))+NTU11*(tw(i,j)-tpl(i,j))-tpl(i,j+1);         end     end [X,Y]=meshgrid(x,y); figure(1)%)%显示t=(100-1)*tao时刻X-Y平面个温度（3D） surf(X,Y,tp(:,:)) xlabel('x'),ylabel('y'),zlabel('z'),shading('interp') colormap(jet),colorbar axis([0 1 0 0.5 0 1]) box on grid on              这是我的差分公式

### 相关帖子

 您需要登录后才可以回帖 登录 | 注册 本版积分规则 回帖后跳转到最后一页