热传导方程在工程仿真中就像老朋友的聚会——大家总得找个方式坐下来聊聊温度怎么传。今天咱们用MATLAB折腾几个一维和二维的导热算例，顺便聊聊不同计算格式的脾气

【MATLAB各种计算格式的一维/二维导热算例】 1. 一维下的全显格式全隐格式半隐格式的导热算例 2.二维下的全显格式全隐格式隐式格式ADI迭代隐式格式点迭代隐式格式高斯-赛德点迭代GS的导热算例 3.可以自己设定导热几何距离进口温度等边界一维战场显式与隐式的博弈先看个一维金属杆传热的例子。假设这根杆子长1米左端100度右端常温25度。全显格式就像急性子抄起代码就冲% 参数设置 L 1; dx 0.01; x 0:dx:L; alpha 1e-5; dt 0.5*dx^2/alpha; % 稳定性条件 T 25*ones(size(x)); T(1) 100; for n 1:500 T_old T; for i 2:length(x)-1 T(i) T_old(i) alpha*dt/dx^2*(T_old(i1)-2*T_old(i)T_old(i-1)); end plot(x,T); pause(0.01); end这段代码里的dt必须满足稳定性条件就像走钢丝时不能迈太大步。显式虽然计算快但时间步长被卡得死死的这时候全隐格式就带着矩阵解法来了A gallery(tridiag,length(x),1,-2,1); A(1,1) 1; A(end,end) 1; % 边界处理 for n 1:50 T (eye(size(A)) - alpha*dt/dx^2*A)\T; plot(x,T); pause(0.1); end隐式方法每次都要解线性方程组计算量飙升但时间步长可以放开手脚。至于半隐格式像是个折中的和事佬把部分项隐式处理这里就不展开代码了。二维战场迭代的艺术【MATLAB各种计算格式的一维/二维导热算例】 1. 一维下的全显格式全隐格式半隐格式的导热算例 2.二维下的全显格式全隐格式隐式格式ADI迭代隐式格式点迭代隐式格式高斯-赛德点迭代GS的导热算例 3.可以自己设定导热几何距离进口温度等边界升级到二维平板传热边界条件变成四边固定温度。全显格式这时候就像在雷区跳舞[X,Y] meshgrid(0:dx:L); T 25*ones(size(X)); T(:,1) 100; % 左边界热源 for n 1:1000 T_old T; for i 2:size(X,1)-1 for j 2:size(X,2)-1 T(i,j) T_old(i,j) alpha*dt/dx^2*(... T_old(i1,j) T_old(i-1,j) ... T_old(i,j1) T_old(i,j-1) - 4*T_old(i,j)); end end surf(X,Y,T); zlim([25 100]); pause(0.01); end二维显式的稳定性条件更苛刻这时候ADI交替方向隐式方法带着它的分治策略登场for n 1:100 % X方向隐式 A gallery(tridiag,size(X,1),1,-2,1); for j2:size(X,2)-1 T(:,j) (eye(size(A)) - alpha*dt/(2*dx^2)*A) \ T(:,j); end % Y方向隐式 for i2:size(X,1)-1 T(i,:) (eye(size(A)) - alpha*dt/(2*dx^2)*A) \ T(i,:); end surf(X,Y,T); pause(0.1); end这种交替解法把二维问题拆成两个一维隐式计算计算量友好。要是想更省事高斯-赛德尔迭代就像个会抄近路的老司机for n 1:500 for i 2:size(X,1)-1 for j 2:size(X,2)-1 T(i,j) (T(i1,j) T(i-1,j) ... T(i,j1) T(i,j-1))/4; end end if mod(n,50)0 surf(X,Y,T); pause(0.5); end end边算边更新最新值收敛速度比雅可比迭代快一倍适合对实时性有要求的场合。参数自由定制指南在代码开头区域可以随意修改L 1; % 区域尺寸 dx 0.02; % 空间步长 alpha 2e-5; % 热扩散系数 boundary [100, 25, 25, 25]; % 四边温度[左,右,上,下]改这些参数就像调鸡尾酒不同的配比会呈现完全不同的温度场演变效果。比如把右边界的25改成50立马能看到热流方向的改变。从显式的横冲直撞到隐式的稳扎稳打再到迭代方法的灵活应变每种算法都有它的生存哲学。下次遇到传热问题不妨像选工具一样根据精度需求和计算资源挑个合适的解法。毕竟在数值计算的世界里没有最好的算法只有最适合场景的选择。