数值分析报告上机作业2.doc

word
word
1 / 17
word
数值实验
数值实验综述：线性代数方程组的解法是一切科学计算的根底与核心问题。求解方法大致可分为直接法和迭代法两大类。直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的准确解的方法，因此也称为准确法。当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时，Gauss消去法是目前最根本和常用的方法。如假如系数矩阵具有某种特殊形式，如此为了尽可能地减少计算量与存储量，需采用其他专门的方法来求解。
Gauss消去等同于矩阵的三角分解，但它存在潜在的不稳定性，故需要选主元素。对正定对称矩阵，采用平方根方法无需选主元。方程组的性态与方程组的条件数有关，对于病态的方程组必须采用特殊的方法进展求解。
实验一
实验名称：矩阵的LU分解.
实验目的：用不选主元的LU分解和列主元LU分解求解线性方程组Ax=b, 并比拟这两种方法.
实验内容与要求
〔1〕用所熟悉的计算机语言将不选主元和列主元LU分解编成通用的子程序，然后用编写的程序求解下面的84阶方程组
将计算结果与方程组的准确解进展比拟，并就此谈谈你对Gauss消去法的看法。
〔2〕写出追赶法求解三对角方程组的过程，并编写程序求该实验中的方程组
〔1〕①不选主元高斯消去法求解方程组
function [L,U]=gauss1(A,b)
n=length(A);
for k=1:n-1
A(k+1:n,k)=A(k+1:n,k)/A(k,k);
A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);
end
L=tril(A(:,1:n),-1)+eye(n);
U=triu(A);
主程序为：
Clear；clc；
word
word
2 / 17
word
a=ones([84,1])*6;
b=ones([83,1]);
c=ones([83,1])*8;
A=diag(a)+diag(b,1)+diag(c,-1);
d=ones([82,1])*15;
b=[7;d;14];
[L U]=gauss1(A,b);
n=length(A);
y=ones(n,1);
y=L\b;
x=ones(n,1);
x=U\y
解为：
x=
1





























word
word
3 / 17
word












**********