山东大学数值计算实验报告2.pdfVIP

山东大学 计算机科学与技术 学院 数值计算 课程实验报告 学号： 姓名： 鸿武 班级：计机14.1 实验题目： 1、 Computer Problems P101 2.2 计算向量b 的1 范数、2 范数、无穷范数。 计算A 矩阵的1 范数、2 范数、无穷范数、条件数、F-范数。 2、 Cholesky 分解 Computer Problems P101 2.6 (1) 写出Cholesky 分解函数 （参考教材上算法） (2) 解方程组Hx=b （取n=2，5，10，分别使用Cholesky 分解及三角形方程组解法， 和调用直接使用软件环境提供的函数解法） (3) 计算H 矩阵的1 范数、2 范数、无穷范数、条件数（取n=2，5，10，可以直接使 用软件环境提供的函数完成）   1 2  2 1     3、 （1）方程组Ax=b，A 1 1 1 ，b 3 ，用简单迭代法（Jacobi 迭代法）和G-         2 2 1 5     S 迭代法分别解方程组。精确到小数点后5 位，最大迭代次数N=100，说明是否收敛   2  1 1  2     （2）方程组Ax=b，A 1 1 1 ，b  3 。         1 1  2 0     实验学时：4 实验日期：2016 年11 月5 日 实验目的： Cholesky 分解，范数，条件数，迭代法解方程组 硬件环境： PC 机 软件环境： MATLAB、JAVA、C++ 实验步骤与内容： 2 4 −2 2 1、A = [ 4 9 −3] b = [ 8 ] −2 −1 7 10 程序代码如下 test1.m clc,clear b = [2 8 10]; A = [2 4 -2 4 9 -3 -2 -1 7]; fprintf(%s\n,向量b的1范数) fanshub1 = abs(b(1)) + abs(b(2)) + abs(b(3)) fprintf(%s\n,向量b的2范数) fanshub2 = sqrt((b(1)).^2 + (b(2)).^2 + (b(3)).^2) fprintf(%s\n,向量b的无穷范数) fanshub3 = max(abs(b)) fprintf(%s\n,矩阵A的1范数) a = zeros(1,3); sum = 0; for i = 1:3 for j = 1:3 sum = sum + abs(A(j,i)); end a(i) = sum; sum = 0; end