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