实验3解线性方程组实验.docVIP

实验3 解线性方程 姓名 报告日期 4.5 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：进一步熟练掌握高斯消去法解线性方程解线性方程 实验原理 高斯消去法 上三角线性方程组的回代解法 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 实验总结 实验报告：根据实验情况和结果撰写并递交实验报告。 参考程序 1.书上的参考程序如下： ①建立backsub.m函数文件，内容如下： function X backsub A,B n length B ; X zeros n,1 ; X n B n /A n,n ; for k n-1:-1:1 X k B k -A k,k+1:n *X k+1:n /A k,k ; end ②建立uptrbk.m函数文件，内容如下： function X uptrbk A,B [N N] size A ; X zeros N,1 ; C zeros 1,N+1 ; Aug [A B]; for p 1:N-1 [Y,j] max abs Aug p:N,p ; C Aug p,: ; Aug p,: Aug j+p-1,: ; Aug j+p-1,: C; if Aug p,p 0 A was singular. No unique solution break end for k p+1:N m Aug k,p /Aug p,p ; Aug k,p:N+1 Aug k,p:N+1 -m*Aug p,p:N+1 ; end end X backsub Aug 1:N,1:N ,Aug 1:N,N+1 ; ③在command window里按如下操作，得输出结果： A [1 2 1 4;2 0 4 3;4 2 2 1;-3 1 3 2] A 1 2 1 4 2 0 4 3 4 2 2 1 -3 1 3 2 B [13;28;20;6] B 13 28 20 6 X uptrbk A,B X 3 -1 4 2 2.自己编写改进的程序： ①用高斯消去法解线性方程Gauss.m函数文件，内容如下： function x Gauss A,B Aug [A B]; n length B ; [row,col] size A ; if n~ row disp I am sorry!The length of B must be equal to rows of A! ; end if det A 0 disp I am sorry!A must be unsignlar! ; end for i 1:n-1 for k i+1:n m -Aug k,i /Aug i,i ; Aug k,: m*Aug i,: +Aug k,: ; end end x n Aug n,n+1 /Aug n,n ; for k n-1:-1:1 s 0; for j k+1:n s s+Aug k,j *x j ; end x k Aug k,n+1 -s /Aug k,k ; end ②在command window里按如下操作，得输出结果： A [1 2 1 4;2 0 4 3;4 2 2 1;-3 1 3 2] A 1 2 1 4 2 0 4 3 4 2 2 1 -3 1 3 2 B [13;28;20;6] B 13 28 20 6 x Gauss A,B x 3 -1 4 2