2120121243刘振锋最终修改版.docVIP

数值计算方法实践作业 学 院： 材料学院 姓 名： 刘振锋 学 号： 2120121234 2012年12月 第一题：选用一种编程语言实现列主元Gauss消去法，并把该程序用于一般矩阵的求逆。具体算例可以参考书P61.第1题，也可以自己设计算例。 对算法的叙述 在解答低阶稠密矩阵时，如果计算过程中没有舍入误差，通过有限步四则运算方程组精确解，以便减少计算过程中舍入误差对计算解的影响。 列主元Gauss消去法的基本思想是在每次消元前，在要消去未知数的系数中找到绝对值最大的系数作为主元，通过方程对换将其换到对角线上，然后进行消元。 2.列主元Gauss消去法的步骤： （1）先在增广矩阵[A,b]的第1列中选取绝对值最大的元，该元素位于第k行，则将第1行与第k行替换。记第一次互换后的矩阵为[A（1）,b（1）]，然后进行第一次消元，得矩阵[A（2）,b（2）]。 （2）在经过第一次消元后矩阵[A（2）,b（2）]的第二列中选绝对值最大的主元，比如第j行，然后将[A（2）,b（2）]的第2行与第j行互换，再进行第二次消元，得矩阵[A（3）,b（3）]。 （3）在经过第2次消元后的矩阵[A（i）,b（i）]的第i列中选绝对值最大的主元，比如第l列，然后互换，再进行第i次消元。 （4）如此经过n-1步，[A,b]被化成上三角矩阵，最后由回代过程求解。 算法描述： 对于做到（4）。 按列选主元，即确定使 如果，则为奇异矩阵，停止计算。 如果，则交换第行与第行元素。 消元计算： 回代计算： 3.试用MATLAB软件编程实现矩阵的列主元素高斯消去法，并求Pascal矩阵的逆矩阵A-1。 具体算例： 【1 1 1 1 ；2 1 1 1 ；3 2 1 1 ；4 3 2 1 】 计算结果x1=0.999999;x2=1.000000;x3=1.0000002;x4=1.000000 A-1=【-1,1,0,0；1，-2,1,0；0,1，-2,1；1,0,1，-1】 总结 列主元素消去法clear clc A=[1 1 1 1 ;2 1 1 1 ;3 2 1 1 ;4 3 2 1 ]；B=[4;6;8;11]; n=length(B); X=zeros(n,1); c=zeros(1,n); d1=0; for i=1:n-1 max=abs(A(i,i)); m=i; for j=i+1:n if maxabs(A(j,i)) max=abs(A(j,i)); m=j; end end if(m~=i) for k=i:n c(k)=A(i,k); A(i,k)=A(m,k); A(m,k)=c(k); end d1=B(i); B(i)=B(m); B(m)=d1; end for k=i+1:n for j=i+1:n A(k,j)=A(k,j)-A(i,j)*A(k,i)/A(i,i); end B(k)=B(k)-B(i)*A(k,i)/A(i,i); A(k,i)=0; end end X(n)=B(n)/A(n,n);（忘记点标点“；”） for i=n-1:-1:1 sum=0; for j=i+1:n sum=sum+A(i,j)*X(j); end X(i)=(B(i)-sum)/A(i,i); end X 计算结果为：X =（0.999999；1.000000；1.000002；1.000000） 【命令】：inv(A) ans = -1 1 0 0 1 -2 1 0 0 1 -2 1 0 0 1 -2 第二题：选用一种编程语言实现对一般