一维稳态导热.doc

#include<>#include<>#include<>#defineN500#10#C5#100####f300#w400voidmain(){ doublec[N],l[N],r[N],b[N],p[N],q[N],y[N],t[N];inti,j;doubleh;h=L/N;b[0]=2**dy*dz/h+2*TC*(dy+dz)*h*Tf;for(i=1;i<=N-2;i++)b[i]=2*TC*(dy+dz)*h*Tf;b[N-1]=(2*TC*(dy+dz)*h+*dy*dz/(CC+T*h/2))*Tf;     //初始化右端源项c[0]=*dy*dz/h+2*TC*(dy+dz)*h;for(i=1;i<=N-2;i++)c[i]=*dy*dz/h+2*TC*(dy+dz)*h;c[N-1]=*dy*dz/h+2*TC*(dy+dz)*h-*CC*dy*dz/h/(CC+T*h/2);//初始化Tp对应系数for(i=0;i<=N-1;i++)l[i]=-CC*dy*dz/h;for(i=0;i<=N-1;i++)r[i]=-CC*dy*dz/h;           //初始化Tp左右温度TW、TE对应系数p[0]=c[0];for(i=0;i<=N-2;i++){q[i]=r[i]/p[i];p[i+1]=c[i+1]-l[i+1]*q[i];}y[0]=b[0]/p[0];for(j=1;j<=N-1;j++)y[j]=(b[j]-l[j]*y[j-1])/p[j];t[N-1]=y[N-1];for(i=N-2;i>=0;i--)t[i]=y[i]-q[i]*t[i+1]; for(i=0;i<=N-1;i++)printf("T%d