线性代数方程组的数值解法.docVIP

《数学实验》实验报告 实验序号：　03　　　　　　　　　　　　　　　　日期：　2010　　年　6　月19　日 班　　级 08信息与计算科学班 姓　名 学号 200820020103 实验名称 线性代数方程组的数值解法 问题背景描述： 3、已知方程组Ax=b，其中A∈R，定义为： A= 试通过迭代法求解此方程组，认识迭代法的含义以及迭代初值和方程组系数矩阵性质对收敛速度的影响，实验要求： 选取不同的初始向量x和不同的方程组右端项向量b，给定迭代误差要求，用雅可比迭代法和高斯-赛德尔迭代法计算，观测得到的迭代向量序列是否均收敛？若收敛，记录迭代次数，分析计算结果并给出你的结论。 取定右端向量b和初始向量,将A 的主对角线元素成倍增长若干次，非主对角线元素不变，每次用雅可比迭代法计算，要求迭代误差满足||||10，比较收敛速度，分析现象并得出你的结论。 8、种群的繁殖与稳定收获：种群的数量因繁殖而增加，因自然死亡而减少，对于人工饲养的种群（比如家畜）而言，为了保证稳定的收获，各个年龄的种群数量应维持不变。种群因雌性个体的繁殖而改变，为了方便起见一下种群数量均指其中的雌性。 种群年龄记作k=1,2，…，n，当年年龄k的种群数量记作x,繁殖率记作b(每个雌性个体1年繁殖的数量)，自然存活率记作（=1-，为1年的死亡率），收获量记为，则来年年龄k的种群数量应为=，（k=1,2，…，n-1）。要求各个年龄的种群数量每年维持不变就是要使= x（k=1，…，n），已知，给定收获量，建立求各年龄的稳定种群数量x的模型（用矩阵、向量表示）。 设n=5,,,如果要求～为500,400,200,100,100，求～。 要使～均为500，如何达到？ 实验目的： 通过建立数学模型掌握求解线性方程组的方法. 实验原理与数学模型： 线性方程组，雅可比迭代法，直接法 实验所用软件及版本： Matlab7.0 主要内容（要点）： 通过直接法、雅可比迭代法和高斯-赛德尔迭代法求解线性方程组。 试验过程纪录（含基本步骤、主要程序清单及异常情况记录等）： 3、（1）说明： p为雅可比迭代法的次数，m为高斯-赛德尔迭代法的次数 主函数： n=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+A3+A2+A3; b=(1:20); x0=zeros(20,1); Jaccobi(a,b,x0); GaussS(a,b,x0); 功能函数： 雅可比迭代法： function y=Jaccobi(a,b,x0) A1=sparse(1:20,1:20,3,20,20); C=inv(A1); fj=C*b; L=-tril(a,-1); U=-triu(a,1); Bj=C*(L+U); y=Bj*x0+fj;p=1; while norm(y-x0)0.001 x0=y; y=Bj*x0+fj;p=p+1; end 高斯-赛德尔迭代法： function z=GaussS(a,b,x0) A1=sparse(1:20,1:20,3,20,20); L=-tril(a,-1); U=-triu(a,1); BGS=inv(A1-L)*U; fGS=inv(A1-L)*b; m=1; z=BGS*x0+fGS; while norm(z-x0)0.001 x0=z; z=BGS*x0+fGS;m=m+1; end ①x0=zeros(20,1); b=(1:20); 时 ②x0=zeros(20,1); b=[1,2,3,4,5,7,9,5,7,8,2,9,5,8,3,6,2,1,2];时 取定右端向量b和初始向量,将A 的主对角线元素成倍增长若干次，非主对角线元素不变，每次用雅可比迭代法计算，要求迭代误差满足||||10 第一组不变；第二组扩大2倍；第三组扩大6倍 8、（1） -= （2）程序： a=[0,0,5,3,0; 0.4,0,0,0,0; 0,0.6,0,0,0; 0,0,0.6,0,0; 0,0,0,0.4,0]; b=inv(a-eye(5)); h=[0;500;400;200;100]; x=b*h; X （3）a=[0,0,5,3,0; 0.4,0,0,0,0; 0,0.6,0,0,0; 0,0,0.6,0,0; 0,0,0,0.4,0]; b=inv(a-eye(5)); h=[0;500;500;500;500]; x=b*h; x 实