How to stop my iteration Tn(j,i)=Tn+1(j,i)

Question

esat gulhan 2020-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/583733-how-to-stop-my-iteration-tn-j-i-tn-1-j-i

编辑： Alan Stevens 2020-8-24

clear clc
dx=0.5
nx=uint32(5/dx+1)
ny=uint32(5/dx+1)
[X Y]=meshgrid(linspace(0,5,nx),linspace(0,5,ny))
Tint=0
T=Tint*ones(ny,nx)
Tleft=100
Tright=0
Ttop=0
Tbottom=0
 
T(:,1)=Tleft
T(:,end)=Tright
T(1,:)=Ttop
T(end,:)= Tbottom
k=3
itr=200
 for k=1:itr
     for i=2:nx-1
        for j=2:ny-1
            T(j,i)=((T(j-1,i)+T(j+1,i)+T(j,i-1)+T(j,i+1))/4)            
        end   
    end
 end
 
[Xq Yq]= meshgrid(linspace(0,5,nx*10),linspace(0,5,ny*10));
colormap jet
Vq=interp2(X,Y,T,Xq,Yq,'cubic',0);

This is my code. My iteration 200 but i want my code stop when Tn(j,i)=Tn+1(j,i) . But i dont know how to break for end...

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Alan Stevens 2020-8-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/583733-how-to-stop-my-iteration-tn-j-i-tn-1-j-i#answer_484256

编辑：Alan Stevens 2020-8-24

在 MATLAB Online 中打开

Replace your for loop by:

tol = 10^-6; % Or whatever is appropriate
err = 1;
 while err>tol
     Told = T;
     for i=2:nx-1
        for j=2:ny-1
            T(j,i)=((T(j-1,i)+T(j+1,i)+T(j,i-1)+T(j,i+1))/4);            
        end   
     end
    err = max(max(abs(T - Told)));
 end