注意长度(x) = 21,但是你的n = 20。这将导致问题。我建议您初始化有点不同,如下所示。
我不明白为什么你有向量T0和T1。我建议T0和T1是标量。我建议您创建一个数组,称为T T的每一行对应一次。T的每一列对应一个位置在酒吧。
% %初始化
L = 1.0;dx = 0.05;x = 0: dx: L;nx =长度(x);%的空间
达峰时间= 30;dt = 0.05;t = 0: dt:最高温度;元=长度(t);%的时间
α= 0.01;%热扩散率,m ^ 2 / s
T0 = 0;%左温度
T1 = 20;%对termperature
T = 0 (nt、nx);%为温度分配数组
T (: 1) = T0;%的左temp = T0
:T (1) = T0;%的初始temp = T0位置
T(:,结束)= T1;%做出正确的temp = T1
解决这个问题。
为我= 2:nt
为j = 2: nx-1
T (i, j) =(张,j) +α* (dt / dx ^ 2) * (T(张,j + 1) 2 * T(张,j) + T(张,j - 1));
结束
结束
一些结果。
情节(x, T (1:)“- r”,x, T(101年:)“- y”,x, T(201年:)“g”,x, T(301年:)“c”,…
x, T(401年:)“- b”,x, T(501年:)“- m”,x, T(601年:)“- r”);
网格在;标题(温度与位置和时间的)
包含(“位置(m)”);ylabel (“温度(K)”)
传奇(t = 0年代”,“t = 5 s”,“t = 10 s”,“t = 15秒”,“t = 20年代”,' t = 25 s ',“t = 30年代”);
当我跑代码与dt = 0.5年代,该国的解决方案,给结果如T = + / -10 ^ 167后100秒。这种情况如果时间步长太大。所以我减少了时间步长10倍:dt = 0.05 s。然后我得到了合理的结果,如上图所示。
学习如何使用冲浪()做一个漂亮的3 d图与轴的位置,时间,温度。
