代码可以修改和匹配附件的pdf吗

1次(最近30天)

显示旧的评论

MINATI智利的 2021年8月12日

评论道: MINATI智利的2021年8月19日

1.WALAILAK(合作伙伴). pdf

                         可以有人修改代码来运行和匹配的图在附件pdf和也想画一些冲浪图
                        
                         % u_t (x, t) = u_ {xx} (x, t) + Gr * t + Gc * C - Q * u, T_t (x, t) =(1 /公关)* T_ {xx}识别(x, t) -φ* t, C_t (x, t) = (1 / Sc) * C_ {xx} (x, t) - Kc * C,
                        
                         % u(x,0) =0,T(x,0) =0,C(x,0) =0, x, T <=0
                        
                         % u(0,t) = a,t (0,t) = t, C(0,t) = exp(t)，
                        
                         % u(Inf,t) = 0, t (Inf,t) = 0, C(Inf,t) = 0。
                        
                         Gr = 4;Gc = 4;φ= 0.2;公关= 0.71;Sc = 0.22;Kc = 1;M = 2;Kp = 0.5;Q = M + (1/Kp);
                        
                         xl = 0;xr = 1;J = 100;dx = (xr-xl)/ J;tf = 0.1;元= 50;dt = tf / Nt;μ= dt / (dx) ^ 2;
                        
                         考虑以下3行，取J = 10;
                        
                         %如果> 0.5%确保dt满足稳定性条件
                        
                         % error('mu应该< 0.5!')
                        
                         %结束
                        
                         X = xl: dx: xr;
                        
                         % f = 0 (2,J);g = 0 (2, J);h = 0 (2, J);
                        
                         f = 0;g = 0;h = 0;u = 0 (J + 1, Nt);v = 0 (J + 1, Nt);w = 0 (J + 1, Nt);
                        
                         为n = 1: Nt
                        
                         t = n * dt;
                        
                         %边界条件在左右两侧
                        
                         一个= 0.2;gl = [* t;t;exp (t)];gr = [0;0;0);
                        
                         如果n = = 1第一次步长%
                        
                         为j = 2: j%内部节点
                        
                         % u (j n) = f (j) +μ* (f (j + 1) - 2 * (j) + f (j - 1) - Gr * g (j) + Gc * h (j) - Q * f (j));
                        
                         % v (j n) = g (j) +(μ/ Pr) * (g (j + 1) - 2 * g (j) + g (j - 1) -φ* g (j));
                        
                         % w (j n) = h (j) +(μ/ Sc) * (h (j + 1) - 2 * h (j) + h (j - 1) - Kc * h (j));
                        
                         U (j,n) = exp(t);v(j,n) = 0w (j n) = 0;
                        
                         结束
                        
                         u (1, n) = gl (1);v (1, n) = gl (2);w (1, n) = gl (3);%左端点
                        
                         u (J + 1, n) = gr (1);v (J + 1, n) = gr (2);w (J + 1, n) = gr (3);%为右端点
                        
                         其他的
                        
                         为j = 2: j%内部节点
                        
                         u (j n) = (j, n - 1) +μ* (u (j + 1, n - 1) - 2 * u (j, n - 1) + u (j - 1, n - 1) + Gr * v (j, n - 1) + Gc * w (j, n - 1) - Q * u (j, n - 1));
                        
                         v (j n) = (j, n - 1) +(μ/ Pr) * (v (j + 1, n - 1) - 2 * v (j, n - 1) + v (j - 1, n - 1) -φ* (j, n - 1));
                        
                         w (j n) = w (j, n - 1) +(μ/ Sc) * (w (j + 1, n - 1) - 2 * w (j, n - 1) + w (j - 1, n - 1) - Kc * w (j, n - 1));
                        
                         结束
                        
                         结束
                        
                         结束
                        
                         %绘制结果
                        
                         tt = dt: dt: Nt*dt;
                        
                         图(1)中,冲浪(x, tt, u ');持有在,包含(“x”), ylabel (“t”), zlabel (“u”)、标题(一维抛物型方程的数值解）
                        
                         图(2),图(x, u (: 1));持有在,包含(“x”), ylabel (“u”）