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

程序说明 该程序用于求解边值问题的非线性Blasius方程的数值解，该方程用来描述通过一块无限大 平板的不可压缩的两维稳定流问题。 2f + #=0 rj = QJ = 0j=0^ = ooj = \ 该程序用MATLAB编写，由s.m和q.m两个程序组成，s.m为主程序，采用了打靶法和RKF （龙格一库塔一费尔伯格）法，在MATLAB中，RKF法选用函数ode45。在用打靶法解题 过程中需要先选定两个初始值（即/”的值）用于壳代计算，在叠代过程屮，第一次采用的 是一次多项式插值法，以后各次均采用的二次多项式插值法。具体程序执行如下： 执行s.m; 命令窗口提不： 请输入第一个任意初始值： 键入 0.3 I 命令窗口提示： 请输入第二个任意初始值： 键入 0.4 I 程序执行完毕。 得到数值解图像 Blasius方程数值解 ?Ityl 0 0332050425572608 0.200000000000000 0.006640684785757 0.066406419000735 0.331976930281647 0.4000000000()0000 0.026559272491103 0」32761397348084 0.331462940642616 ().6()()0000000()0000 0.8000000(X)0000001.0000()00000()0000 1.2