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