MATLAB答案

为什么结束错误在我的代码中不起作用？？

1视图（过去30天）

埃萨特·古尔汉 2020年8月26日

评论：埃萨特·古尔汉2020年8月26日

接受的答复： Rik

                    清楚的clc
                   
                    L=0.02；H=40；dx=0.01；dy=dx；dt=15
                   
                    nx=uint32（L/dx+1）；ny=uint32（H/dy+1）；rx=dt/dx^2；ry=dt/dy^2；
                   
                    k=28；alfa=12.5*10^-6；g=5*10^6；Ts=30；h=45
                   
                    [xy]=meshgrid（linspace（0，L，nx），linspace（0，H，ny））；色调=200；T=色调*个（纽约，nx）；
                   
                    Tleft=0；Tright=色调；Ttop=色调；t底部=色调；
                   
                    T（：，1）=T英尺；T（：，end）=Tright；T（1，：）=Ttop；T（end，：）=Tbottom；Fo=alfa*dt*（dx^2）；
                   
                    时间=0；n=时间/dt；
                   
                    误差=1
                   
                    tol=10^-1
                   
                    虽然厄尔>托尔
                   
                    n=n+1
                   
                    Tn=T；
                   
                    对于i=2:nx-1
                   
                    对于j=2:ny-1
                   
                    T（j，i）=Tn（j，i）+（Tn（j，i+1）+Tn（j，i-1）+Tn（j+1，i）+Tn（j+1，i）-4*Tn（j，i））+g*dx^2/k）*
                   
                    T（j，end）=Fo*（（g*dx^2）/k-（2*h*（Tn（j，end）-Ts）*dx）/k-2*Tn（j，end）+2*Tn（j，end-1）+Tn（j，end）/Fo）
                   
                    终止
                   
                    终止
                   
                    误差=最大值（最大值（abs（Tn-T）））
                   
                    终止