实验3 矩阵的分块求逆及解线性方程组一、 问题化已知矩阵为上三角矩阵，构作范德蒙矩阵，高阶非奇异矩阵的分块求逆，求非齐次线性方程组的通解。

二、 实验目的学会用Matlab 语言编程，实施矩阵的初等变换将已知矩阵化为上三角矩阵；掌握用循环语句由已知向量构造范德蒙矩阵；了解高阶非奇异矩阵用不同分块法求逆矩阵的误差分析；能根据由软件求得的非齐次线性方程组增广矩阵的阶梯型的最简形式写出线性方程组的通解。

三、 预备知识1． 线性代数知识：（1） 向量},,,{21n x x x X =作出的 n 阶范德蒙矩阵为⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---112112222121111n n n n n n x x x x x x x x x（2）分块矩阵⎪⎪⎭⎫⎝⎛=22211211A A A A A ，其中11A 为方的可逆子块，求逆矩阵有如下公式： 设⎪⎪⎭⎫ ⎝⎛=-222112111B B B B A，则221211112112111212222,)(B A A B A A A A B ----=-=， 111211211111111212221,----=-=A A B A B A A B B（3）常用的矩阵范数为Frobenius 范数；21112||||||⎪⎪⎭⎫ ⎝⎛=∑∑==n i n j ij F a A2． 本实验所用Matlab 命令提示：（1）输入语句：input('输入提示')；（2）循环语句：for 循环变量=初始值 ：步长 ：终值 循环语句组 end（3）条件语句： if(条件式1)条件块语句组1elseif(条件式2)条件块语句组2 else条件块语句组3 end（4）矩阵和向量的范数：norm(A)； （5）求矩阵A 的秩：rank （A ）；（6）求矩阵A 的阶梯型的行最简形式：rref(A)。

四、 实验内容及要求1． 在建立的sy31.m 文件中编程将任意给定的n 阶方阵B1，化为上三角矩阵B1；调用时输入：B1=A ，n=6；其A 为实验1[矩阵的基本运算]中的矩阵A 矩阵；2． 在建立的sy32.m 文件中编程用1~6单位增量的行向量产生一个范德蒙矩阵B2； 3． 在建立的sy33.m 文件中编程对任意输入的高阶分块可逆矩阵B3实现分块法求逆；（1）调用sy33.m 文件时输入：B3=A^2 ，输入n1=2求出B3的逆C2 ； （2）调用sy33.m 文件时输入同上的B3，输入n1=4求出B3的逆C4 ； （3）调用sy33.m 文件时输入同上的B3，输入n1=6求出B3的逆C6 ；（4）用norm()函数对上面三种方法所求的逆作误差分析[即作(B3*Ci －E)的范数]； 4． 建立sy34.m 文件，求下列非齐次线性方程组的通解。

⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=-++-+-=+++--=+++-=-+++=-++-+=+-+++246635554254843219526335544242364265432165432642165321654321654321x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x五、思考与练习1． 求下列齐次线性方程组的基础解系⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=-++-+=+++--=+++=-+++=-++-+=+-+++046635504254043209526305544202364265432165432642165321654321654321x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x2．用任意输入的8维行向量构造一个8阶范德蒙矩阵。

六、操作提示1．计算过程（1）sy31.m 文件B1=input('B1=') N=input('N=') for j=1:Nfor i=j+1:NB1(i,:)=-B1(i,j)/B1(j,j)*B1(j,:)+B1(i,:) end end B1调用时输入：Load sy1 A B1=A ，n=6（2）sy32.m文件t=input('t=');tn=input('n=');for i=1:nB2(i,:)=t.^(i-1);endB2调用时输入：T=1:6 ，n=6（3）sy33.m文件B3=input('B3=');[n,m]=size(B3);n1=input('n1=');if(n1>n)disp('n1>n')elseif(n1==n)C1=inv(B3)elseb11=B3(1:n1,1:n1);b12=B3(1:n1,n1+1:n);b21=B3(n1+1:n,1:n1);b22=B3(n1+1:n,n1+1:n);ib11=inv(b11);c22=inv(b22-b21*ib11*b12);c12=-ib11*b12*c22;c21=-c22*b21*ib11;c11=ib11-c12*b21*ib11;C=[c11 c12;c21 c22]End（1）调用时输入：B3=A^2 ，n1=4；转换C4=C；（2）调用时输入：同上的B3 ，n1=2；转换C2=C；（3）调用时输入：同上的B3 ，n1=6；转换C6=C；（4）E=eye(size(B3));norm(B3*C2-E);norm(B3*C4-E);norm(B3*C6-E);（4）sy34.m文件A34=[1 2 4 6 -3 2 4;2 4 -4 5 1 -5 3;3 6 2 0 5 -9 -1;...2 3 0 4 0 1 8;0 -4 -5 2 1 4 -5;5 5 -3 6 6 -4 2]rank(A34)A35=rref(A34)2．计算结果>> sy31B1=AB1 =3 4 -1 1 -9 106 5 07 4 -161 -4 7 -1 6 -82 -4 5 -6 12 -8-3 6 -7 8 -1 18 -4 9 1 3 0N=6N =6B1 =3 4 -1 1 -9 100 -3 2 5 22 -361 -4 7 -1 6 -82 -4 5 -6 12 -8-3 6 -7 8 -1 18 -4 9 1 3 0B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 -5.3333 7.3333 -1.3333 9.0000 -11.33332.0000 -4.0000 5.0000 -6.0000 12.0000 -8.0000-3.0000 6.0000 -7.0000 8.0000 -1.0000 1.00008.0000 -4.0000 9.0000 1.0000 3.0000 0 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 -5.3333 7.3333 -1.3333 9.0000 -11.33330 -6.6667 5.6667 -6.6667 18.0000 -14.6667-3.0000 6.0000 -7.0000 8.0000 -1.0000 1.00008.0000 -4.0000 9.0000 1.0000 3.0000 0 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 -5.3333 7.3333 -1.3333 9.0000 -11.33330 -6.6667 5.6667 -6.6667 18.0000 -14.66670 10.0000 -8.0000 9.0000 -10.0000 11.00008.0000 -4.0000 9.0000 1.0000 3.0000 0 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 -5.3333 7.3333 -1.3333 9.0000 -11.33330 -6.6667 5.6667 -6.6667 18.0000 -14.66670 10.0000 -8.0000 9.0000 -10.0000 11.00000 -14.6667 11.6667 -1.6667 27.0000 -26.6667 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 -6.6667 5.6667 -6.6667 18.0000 -14.66670 10.0000 -8.0000 9.0000 -10.0000 11.00000 -14.6667 11.6667 -1.6667 27.0000 -26.6667 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 1.2222 -17.7778 -30.8889 65.33330 10.0000 -8.0000 9.0000 -10.0000 11.00000 -14.6667 11.6667 -1.6667 27.0000 -26.6667 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 1.2222 -17.7778 -30.8889 65.33330 0 -1.3333 25.6667 63.3333 -109.00000 -14.6667 11.6667 -1.6667 27.0000 -26.6667 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 1.2222 -17.7778 -30.8889 65.33330 0 -1.3333 25.6667 63.3333 -109.00000 -0.0000 1.8889 -26.1111 -80.5556 149.3333 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 0 -14.4706 -21.1471 48.29410 0 -1.3333 25.6667 63.3333 -109.00000 -0.0000 1.8889 -26.1111 -80.5556 149.3333 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 0 -14.4706 -21.1471 48.29410 0 0 22.0588 52.7059 -90.41180 -0.0000 1.8889 -26.1111 -80.5556 149.3333B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 0 -14.4706 -21.1471 48.29410 0 0 22.0588 52.7059 -90.41180 -0.0000 0 -21.0000 -65.5000 123.0000 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 0 -14.4706 -21.1471 48.29410 0 0 0 20.4695 -16.79270 -0.0000 0 -21.0000 -65.5000 123.0000 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 0 -14.4706 -21.1471 48.29410 0 0 0 20.4695 -16.79270 -0.0000 0 0 -34.8110 52.9146 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 0 -14.4706 -21.1471 48.29410 0 0 0 20.4695 -16.79270 -0.0000 0 0 0 24.3566 B1 =3.00004.0000 -1.0000 1.0000 -9.0000 10.00000 -3.0000 2.0000 5.0000 22.0000 -36.00000 0 3.7778 -10.2222 -30.1111 52.66670 0 0 -14.4706 -21.1471 48.29410 0 0 0 20.4695 -16.79270 -0.0000 0 0 0 24.3566t=1:6t =1 2 3 4 5 6n=6B2 =1 1 1 1 1 11 2 3 4 5 61 4 9 16 25 361 8 27 64 125 2161 16 81 256 625 12961 32 243 1024 3125 7776 >> sy33B3=A^2n1=2C =-0.0215 0.0020 -0.0408 -0.0039 -0.0064 0.00100.0115 0.0108 0.0057 0.0065 0.0111 0.02870.0339 0.0003 0.0415 0.0100 0.0115 0.00990.0148 -0.0116 0.0159 0.0167 0.0142 -0.0080-0.0215 -0.0054 -0.0172 -0.0011 0.0158 0.0221-0.0230 -0.0025 -0.0259 -0.0024 0.0034 0.0203>> C2=C;>> sy33B3=A^2n1=4C =-0.0215 0.0020 -0.0408 -0.0039 -0.0064 0.00100.0115 0.0108 0.0057 0.0065 0.0111 0.02870.0339 0.0003 0.0415 0.0100 0.0115 0.00990.0148 -0.0116 0.0159 0.0167 0.0142 -0.0080-0.0215 -0.0054 -0.0172 -0.0011 0.0158 0.0221-0.0230 -0.0025 -0.0259 -0.0024 0.0034 0.0203>> C4=CC4 =-0.0215 0.0020 -0.0408 -0.0039 -0.0064 0.00100.0115 0.0108 0.0057 0.0065 0.0111 0.02870.0339 0.0003 0.0415 0.0100 0.0115 0.00990.0148 -0.0116 0.0159 0.0167 0.0142 -0.0080-0.0215 -0.0054 -0.0172 -0.0011 0.0158 0.0221-0.0230 -0.0025 -0.0259 -0.0024 0.0034 0.0203>> sy33n1=6C1 =-0.0215 0.0020 -0.0408 -0.0039 -0.0064 0.00100.0115 0.0108 0.0057 0.0065 0.0111 0.02870.0339 0.0003 0.0415 0.0100 0.0115 0.00990.0148 -0.0116 0.0159 0.0167 0.0142 -0.0080-0.0215 -0.0054 -0.0172 -0.0011 0.0158 0.0221-0.0230 -0.0025 -0.0259 -0.0024 0.0034 0.0203>> C6=C;>> E=eye(size(B3));>> norm(B3*C2-E)ans =5.1671e-015>> norm(B3*C4-E)ans =4.9230e-015>> norm(B3*C6-E)ans =4.9230e-015>> sy34A34 =1 2 4 6 -3 2 42 4 -4 5 1 -5 33 6 2 0 5 -9 -12 3 0 4 0 1 80 -4 -5 2 1 4 -55 5 -36 6 -4 2ans =5A35 =1.0000 0 0 0 0 6.8068 11.09720 1.0000 0 0 0 -2.2520 -0.41920 0 1.0000 0 0 -0.2041 -1.13840 0 0 1.0000 0 -1.4644 -3.23420 0 0 0 1.0000 -3.1000 -5.90000 0 0 0 0 0 0y=k(－6.8068 2.2520 0.2401 1.4644 3.1000 1)’+(11.0972 －0.4192 －1.1384 －3.2342 －5.9000 0)’。