LBM入门小案例matlab格式书写

代码说明：

代码注释清晰，适合初学者学习使用 %给定初始的密度，流速，设定粒子分布函数的初始值 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