Blasius方程的数值解matlab程序(共4页).doc

程序说明该程序用于求解边值问题的非线性Blasius方程的数值解，该方程用来描述通过一块无限大平板的不可压缩的两维稳定流问题。该程序用MATLAB编写，由s.m和q.m两个程序组成，s.m为主程序，采用了打靶法和RKF（龙格库塔费尔伯格）法，在MATLAB中，RKF法选用函数ode45。在用打靶法解题过程中需要先选定两个初始值（即的值）用于叠代计算，在叠代过程中，第一次采用的是一次多项式插值法，以后各次均采用的二次多项式插值法。具体程序执行如下：执行s.m;命令窗口提示：请输入第一个任意初始值:键入 0.3命令窗口提示：请输入第二个任意初始值:键入 0.4程序执行完毕。得到数值解图像需要完整数值解，可在命令窗口执行 t,yans = 0 0 0 0.332050425572608 0.200000000000000 0.006640684785757 0.0664