Blasius方程的数值解matlab程序.docVIP

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