数值分析
蕤蕊
这个作者很懒,什么都没留下…
展开
-
数值积分:龙贝格求积
通过步长的变化-折半和线性外推法,可以将梯形公式、辛普森公式以及柯特斯公式联系起来,使得精度迅速提高,这就是龙贝格算法的思想,也称为理查森线性外推算法。从计算过程可以得出,只计算了4次的复化梯形求积,但是通过外推技巧,其计算结果从只有 位有效数字提升到 位有效数字,其计算精度有了大幅提升,而计算量增加不多。就能将粗糙的梯形值Tn逐步加工成精度较高的辛普森值Sn、柯特斯值Cn和龙贝格值Rn。原创 2023-02-24 10:17:18 · 2528 阅读 · 0 评论 -
积分的数值解法
本实验由1个小实验组成,在进行积分计算时,若需要精度,可取ε=1e-8.积分的计算结果保留六位有效数字。分别取5个节点和8个节点,用牛顿-柯特斯公式计算积分。原创 2023-02-24 10:02:26 · 147 阅读 · 0 评论 -
非线性方程的数值解法
用Newton法求方程的一个实根近似解,已知该实根在3附近。分别取初值x0=[0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6],用牛顿法计算近似根,统计对应的迭代次数,并画出初值与迭代次数的关系图-柱状图。1输入初值,最大迭代次数Nmax,和精度要求Eps,设置i:=0;4 若i=Nmax,则停止运算,否则i:=i+1,转向步骤2。2 计算X(K+1)=X(K)-f(Xk)/f(Xk)3 若|X(K+1)-X(K)原创 2023-02-24 09:57:47 · 119 阅读 · 0 评论 -
线性方程组的数值解法
分别计算Jacobi迭代法和Gauss–Saidel迭代法其迭代矩阵的谱半径,再用Gauss–Saidel迭代法求解下列线性方程组(3-1),用误差估计e≤10−6和最大迭代次数Nmax=500来控制迭代次数。对于给定的方程组可以构造各种迭代公式,但并非全部收敛。,反复不断地使用迭代式,逐步逼近方程组的精确解,直到满足精度要求为止,这种方法称为迭代法。线性方程组 Ax=b 有惟一解,经过变换构造出一个等价同解方程组。常用的迭代法有两种:雅克比迭代和Gauss–Saidel迭代。原创 2023-02-24 09:45:04 · 198 阅读 · 0 评论 -
插值方法。
我们称这样的问题为插值问题;并称φ(x)为 f (x)的插值函数;f (x)为被插函数,x0 , x1, x2, …, xn是插值节(基)点。以下为两种插值方法的特点比较,当插值节点相同时,无论何种方法得到的插值多项式本质上是相同的,即其插值曲线是重合的。上按等距11个节点,画出牛顿插值函数曲线,分析当节点数较多时,构造的多项式与原函数的逼近关系。求一个简单函数φ(x)作为f(x)的近似表达式,以满足。原创 2023-02-24 09:37:05 · 114 阅读 · 0 评论 -
误差与误差分析
下图为实验不同计算方法的部分计算结果对比图,中间一行为真值,第一行为正推法,第三行为逆推法。(n=15,16,17,18,19,20),将结果填入下表中,分析误差,并画出当n=2、5、8时,该积分对应的几何图形,使用matlab脚本文件实现。从图中可以看出随着n的增加,其积分数值必然是递减的,趋近于0,但一定大于0。⑵.新建脚本文件或函数文件:在新建菜单项中点击下拉箭头,选择脚本或函数,即可实现。⑷函数文件在编写完成后,需要在命令行窗口输入命令,回车确认后,程序即可运行。注意:函数名称与函数文件一致。原创 2023-02-24 09:27:12 · 203 阅读 · 0 评论