矩阵与数值分析上机作业.doc_第1页
矩阵与数值分析上机作业.doc_第2页
矩阵与数值分析上机作业.doc_第3页
矩阵与数值分析上机作业.doc_第4页
矩阵与数值分析上机作业.doc_第5页
全文预览已结束

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1矩阵与数值分析上机实习1. 考虑计算给定向量的范数:输入向量x = (x1; x2; ; xn)T ,输出 x 1, x 2, x 1。请编制一个通用程序,并用你编制的程序计算如下向量的范数:x = (1;111)T; ; y = (1; 2; ; n)T :23n对n = 10,100,1000甚至更大的n计算其范数,你会发现什么结果?你能否修改你的程序使得计算结果相对精确呢?2. 考虑y = f(x) = ln(1+x),其中定义f(0) = 1,此时f(x)是连续函数。用此公式计算x当x10 15; 10 15时的函数值,画出图像。另一方面,考虑下面算法:d = 1 + xif d = 1 then y = 1elsey = ln d=(d 1) end if用此算法计算x10 15; 10 15时的函数值,画出图像。比较一下发生了什么?3. 首先编制一个利用秦九韶算法计算一个多项式在给定点的函数值的通用程序,你的程序包括输入多项式的系数以及给定点,输出函数值。利用你编制的程序计算p(x) = (x 2)9 = x9 18x8 + 144x7 672x6 + 2016x5 4032x4 + 5376x3 4608x2 + 2304x 512在x = 2邻域附近的值。画出p(x)在x1:95; 20:5上的图像。4. 编制计算给定矩阵A的LU 分解和P LU 分解的通用程序,然后用你编制的程序完成下面两个计算任务:(1) 考虑10. . 01. .n n;A =.1 . . . . . .1111R.011111 自己取定x Rn,并计算b = Ax。然后用你编制的不选主元和列主元的Gauss消去法求解该方程组,记你计算出的解为x。对n从5到30估计计算解的精度。(2) 对n从5到30计算其逆矩阵。25. 编制计算对称正定阵的Cholesky分解的通用程序,并用你编制的程序计算Ax = b,其中A = (aij )Rn n,aij =1。b可以由你自己取定,对n从10到20验证程序的可靠性。i+j16. (1) 编制程序House(x),其作用是对输入的向量x,输出单位向量u使得(I2uuT )x =x 2e1。(2)编制Householder变换阵H = I 2uuT Rn n乘以A Rn m的程序HA,注意,你的程序并不显式的计算出H。(3)考虑矩阵12341323A =22e;102730275=2用你编制的程序计算H使得HA的第一列为 e1的形式,并将HA的结果显示。7. 用Jacobi和Gauss-Seidel迭代求解下面的方程组,输出迭代每一步的误差 xk x:5x1x2 3x3 = 2x1+ 2x2 + 4x3 = 13x1+ 4x2 + 15x3 = 108. 取不同的初值用Newton迭代以及弦截法求方程x3 + 2x2 + 10x 100 = 0的实根,列表或者画图说明收敛速度。9. 用二分法求方程ex cos x + 2 = 0在区间0; 4 上的所有根。10. 考虑函数f(x) = sin( x),x 0; 1。用等距节点作f(x)的Newton插值,画出插值多项式以及f(x)的图像,观察收敛性。11. 对函数f(x) = 1+1x2 ,x 5; 5,取不同的节点数n,用等距节点作lagrange插值,观察Runge现象。12. 令f(x) = e3x cos( x),考虑积分02 f(x)dx。区间分为50,100,200,500,1000等,分别用复合梯形以及复合Simpson积分公式计算积分值,将数值积分的结果与精确值比较,列表说明误差的收敛性。13. 分别用2点,3点以及5点的Gauss型积分公式计算如下定积分:1x2sin x2(1)dx;(2)0dx:1x1x214. 考虑微分方程初值问题:dx1dt=(ty y2);(t+1)2x(0) = 2:分别用Euler法,改进的Euler法,Runge-Kutta法求解该方程。分别取步长为0:1,0:01,0:001,计算到x(1),画图说明结果。3要求:1. 考试前提交作业(以Word形式提交结果,按照题目要求提交代码以及数值结果),主题写学号+姓名,发送至邮箱:张宏伟老师:zhuke 2015163.com孟兆良老师

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论