Blasius方程的数值解matlab程序.pdf
版权申诉
135 浏览量
2022-11-11
20:42:50
上传
评论
收藏 120KB PDF 举报
程序说明
该程序用于求解边值问题的非线性 Blasius 方程的数值解，该方程用来描述通过一块无限大
平板的不可压缩的两维稳定流问题。
2 f
'''
ff
'
0
0, f 0, f 0;
, f 1
' '
该程序用 MATLAB 编写，由 s.m 和 q.m 两个程序组成，s.m 为主程序，采用了打靶法和 RKF
（龙格—库塔—费尔伯格）法，在 MATLAB 中，RKF 法选用函数 ode45。在用打靶法解题
过程中需要先选定两个初始值（即
f
的值）用于叠代计算，在叠代过程中，第一次采用的
是一次多项式插值法，以后各次均采用的二次多项式插值法。具体程序执行如下：
执行 s.m;
命令窗口提示：
请输入第一个任意初始值:
''
键入 0.3↓
命令窗口提示：
请输入第二个任意初始值:
键入 0.4↓
程序执行完毕。
得到数值解图像
Blasius方 程 数 值 解
2
f
df/η
1.5
d
2
f/dη
2
1
f
/
d
f
/
d
f
2
0.5
0
-0.5
0 1 2 3
η
4 5 6
需要完整数值解，可在命令窗口执行
>> [t,y]
ans =
0 0 0 0.332050425572608
0.200000000000000 0.006640684785757 0.066406419000735 0.331976930281647
评论0
最新资源