- 1、本文档共5页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
计算方法数值实验四
姓名:龙斌 班级:地物10-2班 学号
数值实验四
实验课题(一) 数值积分方法的使用和比较
对给定积分
I= (2)I= (3)
分别用下述计算方案求其积分,要求误差不超多 。将计算的结果分别按表6-10的形式排列,比较收敛的速度,并讨论或说明原因。
方案一 复化梯形法 方案二 复化辛普森法
方案三 复化两点高斯-体形法
解:
实验目的:学会编程数值积分,并比较它们的收敛速度。
实验内容:用VB分别编程方案一,方案二,方案三,并计算例题一二三,然后比较它们的收敛的速度,并讨论或说明原因。
实验过程:
(1) 用VB编程复化梯形法,复化辛普森法,复化两点高斯-体形法写例题(1)的程序:
Private Sub Command1_Click()
Dim a As Single, b As Single, n As Integer, h As Double, r
a = Val(Text2.Text)
b = Val(Text1.Text)
r = Val(Text3.Text)
Picture1.Cls
For n = 1 To r
x = a
h = (b - a) / 2 ^ n
X1 = a
X2 = a + h / 2
s = a * a * Exp(a) + b * b * Exp(b)
s1 = 0
s2 = X2 ^ 2 * (Exp(X2))
G = 0
D = 1 / (2 * 3 ^ (1 / 2))
C = 0
For k = 1 To 2 ^ n - 1
x = a + k * h
C = C + (x ^ 2) * (Exp(x))
Tn = 0.5 * h * ((a ^ 2) * Exp(a) + 2 * C + (b ^ 2) * Exp(b))
X1 = X1 + h
X2 = X2 + h
s1 = s1 + X1 * X1 * Exp(X1)
s2 = s2 + X2 ^ 2 * Exp(X2)
Sn = h * (s + 2 * s1 + 4 * s2) / 6
Next k
For i = 0 To 2 ^ n - 1
G = G + (a + (i + 0.5 - D) * h) ^ 2 * Exp(a + (i + 0.5 - D) * h) + (a + (i + 0.5 + D) * h) ^ 2 * Exp(a + (i + 0.5 + D) * h)
Next i
Picture1.Print Tab(0); 2 ^ n; + ; 1; Tab(13); Format(Tn, 0.########); Tab(27); Format(Sn, 0.########);
Picture1.Print Tab(41); Format(h * G / 2, 0.########)
Text3.Text = 2 ^ n
Label5.Caption = 分成的区间数
Next n
End Sub
VB运行的结果如下:
(2) 用VB编程复化梯形法,复化辛普森法,复化两点高斯-体形法写例题(2)的程序:
Private Sub Command3_Click()
Dim a As Single, b As Single, n1, n2, n As Integer, h As Double, r As Double
a = Val(Text2.Text)
b = Val(Text1.Text)
r = Val(Text3.Text)
Picture1.Cls
For n = 1 To r
x = a
h = (b - a) / 2 ^ n
X1 = a
X2 = a + h / 2
s = a * (1 + a * a) ^ 0.5 + b * (1 + b * b) ^ 0.5
s1 = 0
s2 = X2 * (1 + X2 * X2) ^ 0.5
G = 0
D = 1 /
文档评论(0)