线性方程组的数值解法.docxVIP

《大学数学实验》作业 线性方程组的数值解法 班级 姓名 学号 日期 目录 TOC \o 1-5 \h \z \o Current Document 目录 2 \o Current Document 【实验目的】 3 \o Current Document 【实验内容】 3 \o Current Document 【题目1］（课本习题第五章第3题） 3 \o Current Document 【第⑴问求解】 3 \o Current Document 【第（2）问求解】 11 \o Current Document 【本题小结】 12 \o Current Document 【题目21 （课本习题第五章第5题） 13 \o Current Document 【第（1）问求解】 14 \o Current Document 【第（2）问求解】 15 \o Current Document 【本题小结】 16 \o Current Document 【题目31 （课本习题第五章第7题） 16 【第（1）问求解】 16 【第（2）问求解】 18 \o Current Document 【本题小结】 19 \o Current Document 【实验心得、体会】 19 注： 本实验作业matlab脚本文件均以ex5_3形式命名，其中ex代表作业，5_3表示第五章第三 小题 【实验目的】 学会用MATLAB软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定 性作初步分析； 通过实例学习用线性代数方程组解决简化的实际问题。 【实验内容】 【题目11《课本习题第五章第3题〉 已知方程组Ax二b,其屮Ae /?20x2°,定义为 -1/4 -1/4 -1/2 3 -1/2 -1/4 3 -1/2 -1/2 3 -1/4 -1/2 1/4 一 1/4 -1/2 -1/4 3 ??. TOC \o 1-5 \h \z ?? ?. 3 -1/2 -1/4 ???-1/2 3 -1/2 _ -1/4 -1/2 3 _ 试通过迭代法求解此方程组，认识迭代法收敛的含义以及迭代初值和方程组系数 矩阵性质对收敛速度的影响。实验要求： ⑴选取不同初始向量兀⑹和不同的方程组右端向量b,给定迭代误差要求，用雅 克比和高斯一一赛徳尔迭代法计算，观测得到的迭代向量序列是否均收敛？若收 敛，记录迭代次数，分析计算结果并得出你的结论； (2)取定右端向量b和初始向量x(0),将A的主对角元素成倍增长若干次，非主 对角元素不变，每次用雅克比迭代法，要求迭代误差满足 比较收敛速度，分析现象并得出你的结论。 【第⑴问求解】 【模型分析】 对于矩阵A,将其分解为： A=D-L-U, 其中:_3003? ? ?00? ? ?00D= 其中: _3 0 0 3 ? ? ? 0 0 ? ? ? 0 0 D= ■ ■ ? 0 ■ ■ ■ 0 0 0 ■ ■ ■ 0 3 0 _0 0 ? ? ? 0 3 \o Current Document % % 0 0 % ? ? ? ? ? ? ? ? 0 0 0 … 0 0 0 %% 0 严）=?1（厶 高斯一赛德尔迭代公式为： 兀（Z二矿】（厶兀伙+1）+伙））+ ￡）, k = 0,1,2...选取初始向量*二［0,0,0???0］丁，b二［1,1,1…1］丁，迭代误差为：卜曲）一严L 10= 用雅可比迭代公法和高斯一赛徳尔迭代法对方程组Ax二b进行迭代计算，在 k = 0,1,2... ［matlab程序】 % 作业题5_3脚本M文件源程序 clear all;clc; n 二 20; Al=sparse(l: n, 1: n, 3, n, n); A2=sparse (1:n~l, 2: n, -0. 5, n, n); A3二sparse (1: n-2, 3: n, -0. 25, n, n); A二A1+A2+A3+A2 +A3; D=diag(diag(A)); L=-tril (A, -1); U=-triu(A, 1); b= ones (20, 1); x(:, 1)二zeros(n, 1); for i=2:1000 x(:, i)=D\ (L+U)*x(:, i-l)+D\b; if no rm (x (:, i) - x(:, i-1), inf )le-5 break; end encl Bg_s= (D-L)\U; %输入A的对角元素%输入A的（上）次对角元素%输入（上）次次对角元素 %输入A的对角元素 %输入A的（上）次对角元素 %输入（上）次次对角元素 %得到A矩阵 %提取D矩阵 %提取L矩阵 %提取U矩阵 %设定b初值（进行变化） %设