内蒙古农业大学数值计算方法作业.docxVIP

下载本文档

30
0
约1.08万字
约 35页
2017-01-30 发布于湖北
举报
版权申诉

内蒙古农业大学数值计算方法作业.docx

1、本文档共35页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

内蒙古农业大学数值计算方法作业

数值计算方法学院: 水利与土木建筑工程学院学号: 姓名: 班级: 任课教师: 贾德彬目录1.非线性方程求根31.1迭代法31.2牛顿法31.3弦截法41.4二分法42.插值52.1线性插值52.2二次插值52.3拉格朗日插值62.4分段线性插值72.5分段二次插值73.数值积分83.1复化矩形积分法83.2复化梯形积分法93.3复化辛甫生积分法93.4变步长梯形法104.线性方程组数值解法104.1约当消去法104.2高斯消去法114.3三角分解法124.4雅可比迭代法135.常微分方程数值解法145.1显式欧拉公式145.2欧拉公式的预测校正系统145.3两步欧拉公式155.4改进欧拉公式165.5四阶龙格库塔法171.非线性方程求根1.1迭代法1.1.1程序代码Private Sub Command1_Click()Dim x1!, x0!, ep!x0 = InputBox(请输入初始值x0)ep = InputBox(请输入误差限ep)f = 0While f = 0x1 = (Exp(x0) - x0) / 2If (x1 - x0) ep ThenPrint x1f = 1Elsex0 = x1End IfWendEnd Sub1.1.2例题求方程在x=0.5附近的根（误差限为0.0000001）1.1.3计算结果1.2牛顿法1.2.1程序代码Private Sub Command1_Click()Dim x0!, x1!, ep!, c!x0 = InputBox(请输入初始值x0;)ep = InputBox(请输入误差限ep;)c = InputBox(请输入常数c;)f = 0While f = 0x1 = (x0 + c / x0) / 2If Abs(x1 - x0) ep ThenPrint x1f = 1Elsex0 = x1End IfWendEnd Sub1.2.2例题用牛顿法求15的开方数的近似值（误差限为0初值为3）1.2.3计算结果1.3弦截法1.3.1程序代码Private Sub Command1_Click()Dim x0!, x1!, ep!, x2!, g0!, g1!x0 = InputBox(请输入初值x0)x1 = InputBox(请输入初值x1)ep = InputBox(请输入误差限ep)f = 0While f = 0g0 = x0 * Exp(x0) - 1g1 = x1 * Exp(x1) - 1x2 = x1 - (g1 / (g1 - g0)) * (x1 - x0)If Abs(x2 - x1) ep ThenPrint x2f = 1Elsex0 = x1x1 = x2End IfWendEnd Sub1.3.2例题求xex-5=0在(1,2)内的根，误差上限为0.00011.3.3计算结果1.4二分法1.4.1程序代码Private Sub Command1_Click()Dim x!, x1!, x2!, ep!, f!, f1!, a!x1 = InputBox(请输入初始坐标x1)x2 = InputBox(请输入初始坐标x2)ep = InputBox(请输入误差限ep, 误差, 0.000001)While Abs(x1 - x2) epx = (x1 + x2) / 2f = (Exp(-(x / 4))) * (2 - x) - 1f1 = (Exp(-(x1 / 4))) * (2 - x1) - 1If f * f1 0 Thenx1 = xElseIf f = 0 ThenPrint xElseIf f * f1 0 Thenx2 = xEnd IfWendPrint xEnd Sub1.4.2例题求f(x)=e(3-x)-1=0在[0,50]内的实根1.4.3运行结果2.插值2.1线性插值2.1.1程序代码Private Sub Command1_Click()X0 = Val(InputBox(请输入插值结点X0:))Y0 = Val(InputBox(请输入插值结点Y0:))X1 = Val(InputBox(请输入插值结点X1:))Y1 = Val(InputBox(请输入插值结点Y1:))f = 0While f = 0x = Val(InputBox(请输入插值点X:))L0 = (x - X1) / (X0 - X1)L1 = (x - X0) / (X1 - X0)y = L0 * Y0 + L1 * Y1Print x=; x, y=; yf = InputBox(是否继续（0/1）:)WendEnd Sub2.1.2例题已知（36,6），（64,8）求49的开方数2