查看: 228|回复: 0|关注: 0

[未答复] 程序可以运行也可以出图像,但是图像跟我要的有差别,而且数值输出的也都是0,哪位可以帮忙看一下怎么解决吗?

[复制链接]

新手

7 麦片

财富积分


050


1

主题

5

帖子

0

最佳答案
发表于 2019-7-11 21:44:22 | 显示全部楼层 |阅读模式
这是程序
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

            

这是我的差分公式

这是我的差分公式
                     

得出的图片

得出的图片
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表