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