z=0;Р for q=1:k-1Р z=z+L(i,q)*U(q,k);Р endР L(i,k)=(A(i,k)-z)/U(k,k);Р endР end Р2、输入代码Р(第一大题):Р1(1)、在matlab主窗口中打入代码(选列主元)РA=[3e-16,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]Рb=[59.17,46.78,1,2];Р[x,det,flag]=Gauss (A,b)Р1(2)、在matlab主窗口中打入代码(不选列主元)РA=[3e-16,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]Рb=[59.17,46.78,1,2];Р[x,det,flag]=Gauss_1 (A,b)Р2(1)、在matlab主窗口中打入代码(选列主元)РA=[10,-7,0,1;-3,2.099999,6,2;5,-1,5,-1;0,1,0,2];Рb=[8,5.909901,5,1]';Р[x,det,flag]=Gauss (A,b)Р2(2)、在matlab主窗口中打入代码(不选列主元)РA=[10,-7,0,1;-3,2.099999,6,2;5,-1,5,-1;0,1,0,2];Рb=[8,5.909901,5,1]';Р[x,det,flag]=Gauss_1 (A,b)Р(第一大题):在matlab主窗口中打入代码РA=[4,2,1,5;8,7,2,10;4,8,3,6;12,6,11,20]Р[L,U,flag]=lu(A)Р3、输出结果分别为A B C D EР Р A BР C DР РEР四、实验小结Р对一部分线性方程组求解时,使用选列主元的方法比直接求解得到的数值解更接近真实解,因此在用Gauss消元法解此类方程时,我们应选择选主元的思想。
全文预览
