求助:一个热传导方程的差分解法!谢谢

zhanglei_1

热传导方程转化为差分方程后是:T(i,j,k,n+1)=T(i,j,k,n)+Pd(i,j,k,n)*t(n)/(p*Cp)...<BR>                                                                 V*[T(i-1,j,k,n)-T(i,j,k,n)+T(i+,j,k,n)-T(i,j,k,n)...<BR>                                                                       +T(i,j+1,k,n)-T(i,j,k,n)+T(i,j-1,k,n)-T(i,j,k,n)<BR>                                                                       +T(i,j,k+1,n)-T(i,j,k,n)+T(i,j,k-1,n)-T(i,j,k,n)<BR>                                                                     ]        %t(n)为加热时间<BR>其中T为25*21*10*10矩阵,边界条件为:T(1,:,:,:)=T(25,:,:,:)=T(:,1,:,:)=T(:,21,:,:)=T0=23...<BR>                                                            %加热四方体区域四侧面为常温<BR>                                                           T(:,:,:,1)=T0=28        %加热时间为0温度为常温<BR>          我使用了直接迭代解上面的方程,但结果是错误的,求助,指点!<BR>          谢谢!<BR>