使用MATLAB求解二维顶盖驱动流问题的详细代码和说明

以下是一个使用MATLAB求解二维顶盖驱动流问题的详细代码和说明。顶盖驱动流是计算流体力学（CFD）中的经典问题，涉及在封闭方腔内由移动顶盖引起的流体运动。本代码采用涡量-流函数法求解不可压缩Navier-Stokes方程。

%% 二维顶盖驱动流问题 - MATLAB求解
% 作者: MATLAB助手
% 日期: 2025年6月18日
% 描述: 本程序使用涡量-流函数法求解二维顶盖驱动流问题
% 参考文献: Ghia, U., Ghia, K.N., & Shin, C.T. (1982). J. Comput. Phys.%% 参数设置
Re = 1000;          % 雷诺数
nx =