资 源 简 介
代码注释清晰，适合初学者学习使用
%给定初始的密度，流速，设定粒子分布函数的初始值
rho=ones(nx,ny); %rho=1
u=zeros(nx,ny); %x方向速度分量
v=zeros(nx,ny); %y方向速度分量
f=zeros(nx,ny,9); %nx*ny个节点，每个节点有9个方向
for m=1:9 %粒子分布函数f(:,:,m)表示同一页的数值，f1=1/9 f2=1/9...1/36, 4/9
f(:,:,m)=w(m)*rho.*(1 ...
+3*(e(1,m)*u+e(2,m)*v)...
+9/2*(e(1,m)*u+e(2,m)*v).^2 ...
-3/2*(u.^2+v.^2));
end
ftemp=zeros(nx,ny,9);
flag=ones(nx,ny); %nx*ny的单位矩阵Flag，用于区分上下边界的位置
flag(:,1)=0; %flag矩阵的第一列赋值为0
flag(:,ny)=0; %矩阵的ny列赋值为0
