MATLAB构建一个二维扩散场,【Matlab 056期】【物理应用3】Matlab 二维对流扩散温度场...

function A=secdim2

% 本程序可实现两种可选边界条件的对流扩散实时动态模拟

% 采用乘方格式

% 初始温度场为10，速度可选

% 时间差分采用全隐式格式

clc

dt=0.001;

bc=input('南北边界条件：1-第一类，2-第二类: ');

q=0;

ddt=input('迭代次数： ');

L1=0.02;r=1000;k=0.017;sp=0;L2=0.01;

u=input('横向速度： ');

if bc==2

v=0;

else

v=input('纵向速度 ： ');

end

n=100;m=100;

dx=L1/n;dy=L2/m;F1=r*u*dy;F2=r*v*dx;

ap0=r*dx*dy/dt;

D1=k*dy/dx;pe1=F1/D1;D2=k*dx/dy;pe2=F2/D2;

tp0=10*ones(n*m,1);%初场温度

tw=10;te=10; %边界条件

A=sparse(n*m);

b=ap0*tp0;

%二维影响系数采用乘方格式

ae=D1*max(0,(1-0.1*abs(pe1^5)))+max(-F1,0);

aw=D1*max(0,(1-0