数值分析试卷(2020-1)解答.docxVIP

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

PAGE3

数值分析(开卷,09/01/2020)

(2019级研究生用)

学院:_______姓名:__________学号:____________

(15分)求函数在区间上关于权函数wx=(1-x2)-1/2的最佳Chebyshev多项式

这里

T

T

T

T

T

?

分别计算当时的逼近函数,请将所得到的逼近函数和原函数绘在一张图上,比较逼近的效果。

解:由于基函数在区间上是正交的,但不是归一化的。事实上,

T

所以有

a

当n=2时

a

a

=

=2

a

=0

a

=

=-

=

s

s

(10分)选定一种矩阵范数,计算矩阵11+ε1

解例如选取一范数。记

A=

A

A

因此,条件数为A1A-11=2+ε2ε

3.(15分)下表给出了函数在不同点的函数值。分别利用二次、三次Lagrange插值多项式逼近函数值。你能估计这两种插值逼近的大致精度吗?为什么?

1.1

0.7652

1.3

0.6201

1.5

0.4554

1.8

0.2818

2.1

0.1104

解估计1.4点值,如果用二次Lagrange插值,可用的三个点x0=1.1,x1

L

其中

l

L

利用三次插值时的基点为x0=1.1,x1=1.3

L

其中

l

所以

L

这时,二次估计的误差大约为0.5371

4.(15分)用Newton迭代法求出方程的至少一个解,选择你认为合适的初始点,计算方程的根,使得近似解的相对误差不超过。请从理论上估计达到精度所需的迭代次数。

解记fx=2x-cosx=0.不难发现,fx在[0,1]区间上异号,且

为了简化相应的Newton迭代公式的迭代函数,我们对原问题进行简化。事实上,原方程等价于

4x=

8x=1+

4

g

这时,

g

=

不难证明

L

取初始值t0=0.5,N

t

t

t

t

收敛。

迭代次数估计公式

L

解出k=

5.(15分)用SOR(松弛因子取ω=1.1)迭代法解方程组

对于你所给定的初始值,估计精度达到需要的迭代次数,并实际计算之。计算该迭代的渐进收敛速度。

解为了给出SOR迭代公式,我们首先给出G-S迭代的分量形式,即

x

B

其中

B=

x

而SOR迭代

x

其中x(k)表示由G-S

x

所以SOR迭代公式为

x

=

这里的迭代矩阵

T

ω

取x0=0

x

x

x

x

而迭代次数的估计公式如下

T

这里我们取∞范数Ts∞=0.309,

6.(15分)利用Broyden方法解非线性方程组

选择初始值x10x20=00计算

解如果该题目不限制方法,也可以将非线性方程组化简为非线性方程求解。事实上,第二个方程可以转化为

x

再将上述表达式代入方程组的第一个方程,得

4

化简得,

4

不难发现这是一个关于变量x2的6

x

7.(15分)设fx=x3-2x,x∈[1,3]。已知数据点xj

解由于函数fx的定义域[1,3],与基函数{

t=(x-2)π

x=

这时原函数变为

g

变换后的数据点:tj

考虑到正交性,我们有

a

b

解出

a

b

g

f

=

文档评论(0)

oyj + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档