下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业程序说明该程序用于求解边值问题的非线性Blasius方程的数值解,该方程用来描述通过一块无限大平板的不可压缩的两维稳定流问题。该程序用MATLAB编写,由s.m和q.m两个程序组成,s.m为主程序,采用了打靶法和RKF(龙格库塔费尔伯格)法,在MATLAB中,RKF法选用函数ode45。在用打靶法解题过程中需要先选定两个初始值(即的值)用于叠代计算,在叠代过程中,第一次采用的是一次多项式插值法,以后各次均采用的二次多项式插值法。具体程序执行如下:执行s.m;命令窗口提示
2、:请输入第一个任意初始值:键入 0.3命令窗口提示:请输入第二个任意初始值:键入 0.4程序执行完毕。得到数值解图像需要完整数值解,可在命令窗口执行 t,yans = 0 0 0 0.2608 0.0000 0.5757 0.0735 0.1647 0.0000 0.1103 0.8084 0.2616 0.0000 0.8740 0.4027 0.3863 0.0000 0.7842 0.7182 0.5304 1.0000 0.4709 0.1296 0.9412 1.0000 0.8780 0.8191 0.8660 1.0000 0.5193 0.0615 0.1080 1.0000
3、0.1678 0.1145 0.0281 1.0000 0.9120 0.3704 0.6929 2.0000 0.0154 0.4094 0.6562 2.0000 0.9520 0.0408 0.5080 2.0000 0.1476 0.3094 0.9961 2.0000 1.2861 0.8408 0.8860 2.0000 1.1495 0.8625 0.5294 3.0000 1.4558 0.6905 0.6103 3.0000 1.5089 0.5320 0.7636 3.0000 1.3586 0.5514 0.3319 3.0000 1.1937 0.2812 0.1445
4、 3.0000 2.1031 0.5617 0.2537 4.0000 2.5077 0.2905 0.6456 4.0000 2.6955 0.4841 0.9064 4.0000 2.6070 0.8678 0.4526 4.0000 2.5013 0.1711 0.2806 4.0000 3.6712 0.2908 0.1081 5.0000 3.2969 0.2639 0.9442 5.9999 3.1607 0.6837 0.2965 5.0000 3.0485 0.0022 0.6624 5.0000 3.3193 0.9695 0.2210 5.0000 4.9796 0.801
5、7 0.9394 6.0000 4.3525 0.1846 0.1230 6.9999 4.2378 0.3418 0.7374 6.0000 4.0177 0.5153 0.6568 6.0000 4.8227 0.2873 0.2299 6.0000 5.0512 0.2870 0.9775 7.0000 5.5719 0.1772 0.5724 7.0000 5.1113 0.8357 0.1169 7.0000 5.0127 0.7110 0.8928 7.0000 5.6138 0.6780 0.7492 7.0000 6.4299 0.8038 0.3282 8.0000 6.59
6、44 0.7916 0.3429 8.9999 6.5597 0.9281 0.0147 8.0000 6.5744 0.5243 0.3786 8.0000 6.9892 1.3378 0.2564 8.0001 7.4974 1.1604 0.5283 9.0000 7.9148 1.0094 0.9078查看叠代过程,可在命令窗口执行 YY = 0.0000 0.0000 0.6614 0.2608 xx = 0.9033 1.2684 1.5041 1.0094%Y表示时的值,x表示(实际为9)时对应下的值。从中可以看到趋近过程。数值解和精确解的比较数值解精确解数值解精确解数值解精确解
7、000000.26080.3320610.47090.165770.12960.329790.94120.3230120.01540.650030.40940.629770.65620.2667531.45581.396820.69050.846050.61030.1613642.50772.305760.29050.955520.64560.0642353.29693.283290.26390.991150.94420.0159164.35254.279640.18460.998680.12300.0024075.57195.279260.17720.999220.57240.0002286
8、.59446.279230.79161.000000.34290.000018.87.49747.079231.16041.000000.52830.00000从上表对比可知,数量解和精确解的结果几乎一致。该程序最大的优点就是可以对输入不同的初始值,来观察叠代次数和对结果的影响。例如:当输入初始值为0.1和1时,与前面的比较见下表:0.00000.00000.90330.75830.00001.00001.26842.79450.66140.11931.50411.28040.26080.74381.00941.21210.18550.72260.62901.4372很明显可以看出叠代次数不一
9、样,第二次比第一次多计算两次,但对结果的影响几乎可以忽略不计。说明叠代初值对叠代次数的影响非常大,如选初值为3和40,则需要计算172次,但仍然能够将精确计算到0.5133。s.m:clear,clfY(1)=input(请输入第一个任意初始值:)t,y1=ode45(q,0:0.2:9,0,0,Y(1);x(1)=y1(length(t),2);Y(2)=input(请输入第二个任意初始值:)t,y2=ode45(q,0:0.2:9,0,0,Y(2);x(2)=y2(length(t),2);Y(3)=polyval(polyfit(x(1) x(2),Y(1) Y(2),1),1);t,y3=ode45(q,0:0.2:9,0,0,Y(3);x(3)=y3(length(t),2);n=0;while 1n=n+1;Y(n+3)=polyval(polyfit(x(n) x(n+1) x(n+2),Y(n) Y(n+1) Y(n+2),2),1);t,y=ode45(q,0:0.2:9,0,0,Y(n+3);x(n+3)=y(length(t),2);if abs(x(n+3)-1)10e-7 break;endendplot(t,
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 财务管理实习报告14篇
- 4.2 行使权利 注意界限(课件) 2024-2025学年高中政治统编版选择性必修二法律与生活
- 2024-2025学年小学数学五年级下册浙教版教学设计合集
- 小学数学苏教版二年级下《口算两位数加两位数》教学设计
- 足球脚背外侧运球 教案-2023-2024学年高一上学期体育与健康人教版必修第一册
- 比例尺(教案)-2023-2024学年六年级下册数学人教版
- 篮球基本功训练 教学设计-2023-2024学年高一上学期体育与健康人教版必修第一册
- 黑龙江省绥化市青冈一中2024-2025学年高考绝密冲刺卷:化学试题试卷含解析
- 黑龙江省大庆十中2025届高三统一质量检测试题化学试题试卷含解析
- 第1章 第2课时 有理数2023-2024学年七年级上册数学高效课堂教学设计(人教版)
- 2024秋国开学习网《形势与政策》形考任务专题测验1-5答案
- 部编人教版小学四年级上册道德与法治全册教案
- 2024-2025学年八年级上学期英语第一次月考模拟考试卷01(人教版)
- 2024年北京移动校园招聘(高频重点提升专题训练)共500题附带答案详解
- 小学主题班会:爱眼护眼-保护视力课件
- 2024年湖北省高考真题政治试卷及答案
- 统编版语文二年级上册第四单元 我是家乡代言人单元任务群整体公开课一等奖创新教学设计
- 2024年全国乡村医生考试复习题库及答案(共420题)
- 账户变更补充合同范本
- 幼儿园大班语言文学活动《小鱼睡在哪里》课件
- DL-T 1071-2023 电力大件运输规范
评论
0/150
提交评论