《矩阵与数值分析》数值实验报告.docVIP

大连理工大学 2013级工科硕士研宄生 《矩阵与数值分析》数值实验报告 教学班号：04 任课教师：金光日 所在院系： 学 号： 姓 名： 完成日期：2013年12月20日 1 - 1 1 - 1 - ..-2- ..-2- ??-2- ..-3- ..-5- ..-5- ..-5- ..-5- ..-7- ..-8- ..-8- ..-8- ..-8- ..-9- -10- -10- -12- -12- -12- -12- -12- -13- -13- -13- -13- -13- -15- -19- -19- -19- -21 - 目录 一、 求数列的前N项和 题目一 1、 问题 2、 程序 3、 结果与分析 二、 解线性方程组 题目一 1、 问题 2、 程序 3、 结果与分析 题目二 1、 问题 2、 程序 3、 结果与分析 三、 非线性方程的迭代解法 题目一 1、 问题 2、 程序 3、 结果与分析 题目二 1、问题 3、结果与分析 四、 数值积分 题目一、 1、 问题 2、 程序 3、 结果与分析 五、 插值与逼近 题目一、 K 问题 2s程序 3、结果与分析 题目二 1、 问题 2、 程序 3、 结果与分析 一、求数列的前N项和 题目一、 1、问题 N 1Q6 S s,v = —，分别编制从小到大和从大到小的顺序程序分别计算兄画，\胃）0 y=2 j -1 并指出两种方法计算结果的有效位数。 2、程序 (1)从小到大(从2到N) clear: format long; N= input C 请输入N= ’)： s=0; for j=2:N s=s+i(r(6)/(j*j-i): end%求近似值 for k=0:100 if s=l(T(k) break end end 3=0.5*1*10*(6)*(1. 5-l/N-l/(N+l)); e=abs(s-a) ;%求误差界 for i=0:100 if e=l/2*10*(-i) break end end dispf结果是 s if e=0 dispC有效教字的位数是:?’)； k+i-1 end (2)从大到小(从N到2) clear: format long; N= input C 请输入N= ’)： s=0; for j=2:N s=s+l(T (6)/[(N-j+2厂2-1]: end%求近似值 for k=0:100 if s=l(T(k) break end end a=0.5*l0A (6)*(1. 5-l/N-l/(N+l)); e=abs(s-a) :%求误差界 for i=0:100 if e=l/2*10(-i) break end end dispf结果是 s if 广=0 dispC’有效数学的位数是: k+i-1 end if e==0 disp C误差为零): end 3、结果与分析 (1)从小到大 (2)从大到小 请输入N=1OOOO 请输入N=1OOOO 结果是： 结果是： s = s = 7,499000049994994e+05 7.499000049995000e+05 有效教字的位教是： 误差为零 请输入N=1OOOOOO ans = 结果是： 14 s = 请输入N=1OOOOOO 结果是： 7.499990000005000e+05 s- 误差为零 7.4999900000052006+05 有效教字的位教是： ans = 13 I 分析：使用format long，显示15位双精度。当按照从2到N方法叠加时，a??an+i,会出 现“大数吃小数“现象，导致有效数字位数损失，有效数字位数少于使用=1000,或N=1000000 逆序开始求和。 二、解线性方程组 题0—、 1、问题 分别利用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组Ax = /)，其屮常向S为if 维随机生成的列向量，系数矩阵A具有如下形式 ■ ? ? ?- 为《-1阶矩阵，为n-1阶单位矩阵，迭代法计算停止的条件 ...-1 -1 2 j 为：k+i 一 给出” =10,20,30吋的不同迭代步数. 2、程序 (1) Jacobi迭代法 %jacobi迭代法 clear: n= input (’ input n二?) a=sparse ([1 :n-b 1 :n-2, 2:n-1, 1:n-5, 5:n-1]，[1 :n-l, 2:n-1, 1:n-2, 5:n-l，l:n-5]，. [4*ones (b n-1), -l*ones (1，n-2)-l*ones (1，n-2), -l*ones(1, n-5)-l*ones (b n-5) ]: n-1