考研必备数值分析

发布 2022-06-08 06:18:28 阅读 4936

第二章多项式插值 (习题)

1. 利用lagrange插值公式求下列各离散函数的插值多项式(结果要简化):

解(2):方法一。 由 lagrange 插值公式。

可得: 方法二。 令。

由,, 定a,b (称之为待定系数法) □

2. 设是以为节点的次多项式插值问题的基函数。

(1)证明。

(2)证明。

证明(1) 由于

故: ,当时有: ,

也即为的插值多项式,由唯一性,有:

证明(2):

方法一。 用数学归纳法。

方法二。 利用newton插值多项式。

差商表:f(x) 一阶二阶n阶差商。

代入式有:

为次代数多项式,由插值多项式的唯一性:

有。8.设 ,证明。

而且。证明

利用数学归纳法,可证:

因为。故有:

15.设,求在区间上的分段线性插值函数,并估计误差,取等距节点,且。

解 ,设 ,则:

误差估计:第四章数值积分方法与数值微分 (习题)

1.直接验证梯形公式(1.2)与中矩形公式(1.3)具有1次代数精度,而辛甫生公式(1.4)则具有3次代数精度。

解。梯形公式:.

矩形公式: .

以上两求积公式以代入公式两边,结果相等,而以代入公式两边,其结果不相等。故梯形公式的代数精度等于1.

simpson公式:

容易验证:以分别代入simpson公式两边,结果相等,而以代入simpson公式两边,其结果不相等,故simpson求积公式的代数精度为3

2.设在[0,1]上连续,在[0,1]上可积,证明:用复化梯形公式计算的误差形式为。

当。其中是复化梯形和,为积分区间的分划节点。

证明: 当。又 于是:

4.证明中矩形公式的peano核误差公式为:

其中 并由此导出误差形式。

解已知中矩形公式对于一次多项式精确成立,由taylor展开:

又:5. 求系数,使求积公式。

对于次数的一切多项式都是精确成立的。

解:求积公式。

是一个插值型求积公式,令得:

解得: ,第五章线性代数方程组的解法 (习题)

6. 假定已知的三角分解:,试设计一个算法来计算的元素。

解: 记: ,其中:为的第列元素。

由于:, 故: ,矩阵的元素即为

因为有 ,记: 有

解系数矩阵为三角阵的方程组(*)求得,然后解系数矩阵为上三角阵的方程。

组(**即可求得,因而得到的结果。

7.试证对维向量有。

证明: 又故。

18.设有迭代格式。

其中。试证该迭代格式收敛。并取,计算。

解: 可先验证: .

因迭代格式收敛,由于,故:

即 ,记:

于是。即 ,即为精确解。

19.给定方程组。

证明jacobi迭代方法收敛而g-s迭代方法发散。

解:方程组:

jacobi方法:迭代矩阵:

特征方程:

或: jacobi方法收敛。

gauss-seidel迭代方法:

迭代矩阵:特征方程)

或的特征化为下面方程的根:

即: 重根)

故: ,gauss-seidel 迭代方法发散。 □

23.求证矩阵。

当时正定,当时jacobi迭代法解收敛。

证明: 矩阵。

正定的充分必要条件为:

(1)、(2)均满足,由(1)有:

由(2) 有: 记。

可知: 于单调上升,于单调下降,于单调上升, 为的极值点,为的零点。

故: ,当时,合并、,可知当时,矩阵正定。

又若要jacobi迭代法收敛,则要求取,使得和同时正定,此时。

为正定矩阵的充分必要条件为:

、同时成立,由得:

由得: 记:

可知: 于单调下降,于单调上升,又: ,为的零点 ,故: 当 ,综合起来,当时矩阵和同时为正定对称矩阵,故

当jacobi方法解方程是收敛的。 □

第七章非线性方程数值解法 (习题)

2. 为求方程在附近的一个根,设将方程改写为下列等价形式,并建立相应的迭代公式:

1),迭代公式。

2),迭代公式,3),迭代公式,试分析每一种迭代公式的收敛性,并问哪一种迭代收敛得快?

解:取的邻域来考察。

(1), 故迭代公式(1)收敛。

故迭代公式(2)也收敛。

故迭代公式(3)发散。

由于越小,越快地收敛于根,故(2)式收敛最快。□

16.讨论以下序列当时的收敛阶:

对(1)、(2)再用aitken方法计算,至。

解 (1),有: ,线性收敛速度。

(2) ,有: ,线性收敛速度。

(3) ,有。

平方收敛速度。

第八章常微分方程数值解 (习题)

1.用euler方法解初值问题。

并证明其截断误差。

解。将enler方法应用于值问题。

得差分方程初值问题。

这里, 由此得。

迭代得到。而此问题真解为。

于是其截断误差。

2.证明:改进的euler方法是稳定的。

证明:设所考虑的初值问题。

其中关于满足lipschitz条件:

考虑 令,则有。

于是,当时,有。

从而。注意,这里。

取,则有。即改进的euler方法是稳定的。

4.构造形如下面形式的三阶格式。

解将,,,在处泰勒展开,并考虑线性组合。得。令

解得 将此参数代入原格式中便得到三阶格式族。

5.求具有最高阶的三步方法的系数。

解一般三步方法为。

由阶条件。

数值分析报告

数值分析 计算报告。第一次作业。定义矩阵,由矩阵性质可知。根据计算实习说明书要求,所有矩阵的零元素都不存储,所以将矩阵压缩存储为矩阵。在矩阵c中检索矩阵的带内元素的方法为 首先,用幂法求出矩阵按模最大的特征值。然后,根据原点平移法,同样根据幂法求出矩阵按模最大的特征值。当时,当时,当时,或,若,则,...

数值分析考点

第一章。最重要的是数字,有效数字的计算。1.2.2例题的第一题,定理1.2.1有效数字和相对误差限的关系。例1.2.2,基本每年都考。第二章。高斯消去法知道就行了。lu分解知道lu什么意思就行了。迭代法雅格比和高斯s迭代法很重要,基本每年大题目。迭代格式和迭代矩阵。sor不用管。迭代法收敛条件。求1...

数值分析作业

12.求在 0,1 上的一次最佳平方逼近多项式与二次最佳平方逼近多项式。函数 function s zjpfbj n,a,b 创建一个函数,里面填入次数,和区间范围。base inline x j 1 x j 定义多项式。quan inline 1 x 权函数。a zeros n 1 y zeros...