值积分算法的程序;
2、对基本题,分别取不同步长hba
,试比较计算结果(如
1020
等),
并比较其结果;
4、对应用题,用给定精度ε,试用1用逐次分半梯形法。2用逐次分半辛普生
法,并确定最佳步长。
逐次梯形法的主函数程序:
fu
ctio
jife
umzhuci_tixi
gabtolif
argi
3
eps10e3e
ddelt1
1hbaThzhuci_tixi
g_fazhuci_tixi
g_fb2whiledelt3tol
hh2
fT0Ts0fori1
xah2i1sszhuci_tixi
g_fxe
dTT02hs
2
deltTT0e
djife
T
um
e
d
执行程序后的结果:
f实验分析:
逐次分半的积分算法具有很高的精度,对于复杂的工程实践问题具有很高的应用价值。同时,在利用复合梯形公式、复合Simpso
公式以及Romberg算法等计算数值积分时,精度已经很高,并且随着步长越小,积分精度越高。另外,当给定一个精度的时候,我们还可以利用自适应步长法,确定所需要的最佳步长和结果。
同样在对比中可见:simpso
公式的收敛速度比梯形公式的收敛速度快。
f实验四线性方程组的直接解法
一、问题提出给出下列几个不同类型的线性方程组,请用适当算法计算其解。1、设线性方程组
4231210000x15
8
6
536
5
0
1
0
0
x2
12
42213
0
2
1
51
23
1011
39
14
x3x4
3
2
42
616
7
33
2
3
x5
3
86857172635x646
0
2
1
34
2
5
3
0
1
x7
13
161011917342122x838
4
6
27139
2
0
12
4
x9
19
00183248631x1021
参考解:x1101203112T
2、设对称正定阵系数阵线方程组
42402400x10
2
2
12
1
3
2
0
x2
6
41141
0
2
1
6
2181
8122
344
5310
6
3
3
x3x4x5
20239
433441114x622
0
2
53101
14
2
x7
15
006334219x845
参考解:x11021102T
3,三对角形线性方程组
f4100000000x17
1410
0
0
0
0
0
0
x2
5
014100
0
01410
00
00
00
00
x3x4
13
2
0
0
01410
0
0
0
x5
6
0000141000x612
0
0
0
0
01410
0
x7
14
0000001410x84
0
0
0
0
0
0
0
1
4
1
x9
5
0000000014x105
参考解:x2130123011T
要求:1、对上述三个方程组分别利用Gauss顺序消去法与Gauss列主元消去法;平方根法与改进平方根法;追赶法求解(选择其一);2、应用结构程序设计编出通用程序;3、比较计算结果,分析数值解误差的原因;
4、r