实验5-线性代数方程组的数值解法.docxVIP

实验5 线性代数方程组的数值解法实验目的学会用MATLAB数值求解线性代数方程组，对迭代法的收敛性和稳定性做初步分析；通过实例学习通过用线性代数方程组解决简化实际问题。实验内容3.已知方程组（略）选取不同的初始向量和不同的方程组右端项ｂ，给定迭代误差要求，用雅克比迭代法和高斯-赛德尔迭代法计算，观察得到的迭代向量序列是否均收敛？若收敛，记录迭代次数，分析计算结果并得出结论；雅克比迭代法的函数：jacobi.mfunction [x,k]=jacobi(x0,a,b,ep)ang=diag(a);n=length(ang);for i=1:nfor j=1:n a(i,j)=a(i,j)/ang(i);end b(i)=b(i)/ang(i);endcoefficient=a-eye(n);x=[x0];cha=100;n=1;while chaep^2 storage=-coefficient*(x(n,:))+b; x=[x;storage]; var=(x(n,:)-x(n+1,:)).^2; cha=max(var); n=n+1;endk=n;高斯-赛德尔迭代法的函数：gauss.mfunction [x,k]=gauss(x0,a,b,ep)ang=diag(a);n=length(ang);for i=1:nfor j=1:n a(i,j)=a(i,j)/ang(i);end b(i)=b(i)/ang(i);endcoef=a-eye(n);x=[x0];cha=100;n=1;while chaep^2 xn=x(n,:);for i=1:length(ang) storage=-coef(i,:)*(xn)+b(i); xn(i)=storage;end x=[x;xn]; var=(x(n,:)-x(n+1,:)).^2; cha=max(var); n=n+1;endk=n;实现任务的脚本：Iteration.mn=20;A1=sparse(1:n,1:n,3,n,n);A2=sparse(1:n-1,2:n,-1/2,n,n);A3=sparse(1:n-2,3:n,-1/4,n,n);A=A1+A2+A2+A3+A3;AA=full(A);x0=zeros(1,20);b=ones(1,20);ep=0.001;[x1,k1]=jacobi(x0,AA,b,ep);x1k1[x2,k2]=gauss(x0,AA,b,ep);x2k2输出结果如下：x1 = 1 至 11 列 0 0 0 0 0 0 0 0 0 0 0 0.3333 0.3333 0.3333 0.3333 0.3333 0.3333 0.3333 0.3333 0.3333 0.3333 0.3333 0.4167 0.4722 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.4537 0.5278 0.5718 0.5810 0.5833 0.5833 0.5833 0.5833 0.5833 0.5833 0.5833 0.4689 0.5527 0.6046 0.6184 0.6236 0.6248 0.6250 0.6250 0.6250 0.6250 0.6250 0.4758 0.5638 0.6196 0.6362 0.6430 0.6451 0.6457 0.6458 0.6458 0.6458 0.6458 0.4789 0.5689 0.6266 0.6445 0.6523 0.6549 0.6559 0.6562 0.6562 0.6562 0.6562 0.4804 0.5713 0.6298 0.6485 0.6568 0.6598 0.6609 0.6613 0.6614 0.6614 0.6614 0.4810 0.5724 0.631