大连理工矩阵与数值分析上机作业.doc
文本预览下载声明
矩阵与数值分析
（上机作业）
Matrix and numerical analysis
学 院（系）： ~~~~~~~~~~
学 生 姓 名： ~~
指 导 老 师： ~~~
学 号： ~~~~~~~~~
完 成 日 期： 2012.11.24
大连理工大学
Dalian University of Technology
1.给定n阶方程组，其中
，
则方程组有解。对和，分别用Gauss消去法和列主元消去法解方程组，并比较计算结果。
1）程序代码：
输入A，b
function [X,A,b,ep] = gauss(A,b)
n=length(A);
%选主元
for k=1:n-1
[~,p]=max(abs(A(k:n,k)));
p=p+k-1;
if pk
y=A(k,:);A(k,:)=A(p,:);A(p,:)=y;
y2=b(k);b(k)=b(p);b(p)=y2;
end
if abs(A(k,k))1e-10
disp(高斯消去法求解失效)
break
end
for i=k+1:n
m=A(i,k)/A(k,k);
A(i,k:n)=A(i,k:n)-m.*A(k,k:n);
b(i)=b(i)-m*b(k);
end
end
%求解X
X(n)=b(n)/A(n,n);
for k=n-1:-1:1
X(k)=(b(k)-sum(A(k,k+1:n).*X(k+1:n)))/A(k,k);
end
end
输入数组，A，b
调用函数 [X,a,c,ep]=gauss(A,b);
2）结果分析：
n=10，两种方法运行结果均近似为x=[1,1,…,1];
n=84，选主元的gauss消去法运算后的右下角元素变得很小，约为10^-24,矩阵变得开始病态。最大误差为2.80*10^-6
不选主元，求解的最大误差为5.37*10^8,可见求解结果误差极大。
2. 设有方程组，其中，
(a) 选取不同的初始向量和不同的右端向量，给定迭代误差要求，用Jacobi和Gauss-Seidel迭代法计算，观测得出的迭代向量序列是否收敛。若收敛，记录迭代次数，分析计算结果并得出你的结论。
(b) 选定初始向量初始向量和不同的右端向量，如取。将的主对角线元素成倍增长若干次，非主对角元素不变，每次用Jacobi法计算，要求迭代误差满足，比较收敛速度，分析现象并得出你的结论。
1）程序代码：
function [X,k,t] = e2(A,b,ep,X0)
U=-triu(A,1);
L=-tril(A,-1);
D=diag(diag(A,0));
B=D^-1*(L+U);
f=D^-1*b;
%B=(D-L)^-1*U;
%f=(D-L)^-1*b;
k=0;
t=1;
while tep
X=B*X0+f;
t=max(abs(X-X0));
X0=X;
k=k+1;
if k1000
disp(不收敛)
break
end
end
end
2）结果分析：
情况一：对于b=（1,1,…,1），X0=(0,0,…，0);Jacobi迭代收敛，迭代次数20；G-S收敛，迭代次数13；X0=（1，2，3，…，10），Jacobi迭代收敛，迭代次数23；G-S收敛，迭代次数16
情况二：对于b=（1，2，…，10），X0=(0,0,…，0) ；Jacobi迭代收敛，迭代次数23；G-S收敛，迭代次数16；X0=（1，2，3，…，10），Jacobi迭代收敛，迭代次数22；G-S收敛，迭代次数15
在这两种情况下，两种方法均收敛，且G-S收敛速度较快
3. 用迭代法求方程的全部根，要求误差限为。
1）程序代码：
二分法：
clear
t=1;
ep=0.5*10^-8;
%求根区间
a=-10000;b=10000;
m=0;p=3;c(10000)=0;cc(10000)=0;c(1)=a;c(2)=b;c(3)=(a+b)/2;d(6)=0;
while m3
m=0;
for i=1:p-1
cc(2*i-1)=c(i);
cc(2*i)=(c(i)+c(i+1))/2;
end
cc(2*i+1)=c(p);
c=cc;
p=2*
显示全部