Matlab实验1 Matlab初步.doc

问题已知矩阵A、B、b如下:应用Matlab软件进行矩阵输入及各种基本运算。实验目的学会使用Matlab软件构作已知矩阵对应的行(列)向量组、子矩阵及扩展矩阵,实施矩阵的初等变换及线性无关向量组的正交规范化,确定线性相关相关向量组的一个极大线性无关向量组,且将其余向量用极大线性无关向量组线性表示,并能编辑M文件来完成所有的实验目的。预备知识线性代数中的矩阵及其初等变换、向量组的线性相关性等知识。Matlab软件的相关命令提示如下;选择A的第i行做一个行向量:ai=A(i,:);选择A的第j行做一个列向量:ai=A(j,:);选择A的某几行、某几列上的交叉元素做A的子矩阵:A([行号],[列号]);n阶单位阵:eye(n);n阶零矩阵:zeros(n);做一个n维以0或1为元素的索引向量L,然后取A(:,L),L中值为1的对应的列将被取到。将非奇异矩阵A正交规范化,orth(A);验证矩阵A是否为正交阵,只需做A*A'看是否得到单位阵E。两个行向量a1和a2的内积:a1*a2'。让A的第i行与第j列互换可用赋值语句:A([i,j],:)=A([j,i],:);让K乘以A的第i行可用赋值语句:A(i,:)=K*A(i,:);让A的第i行加上第j行的K倍可用赋值语句:A(i,:)=A(i,:)+K*A(j,:);求列向量组的A的一个极大线性无关向量组可用命令:rref(A)将A化成阶梯形行的最简形式,其中单位向量对应的列向量即为极大线性无关向量组所含的向量,其它列向量的坐标即为其对应向量用极大线性无关组线性表示的系数。四、实验的内容与要求1、请在MATLAB直接输入下列常数,看它们的值是多少:aibjcepsdinfenanfpigrealmaxhrealmin>>ians=0+>>jans=0+>>epsans=-016>>infans=Inf>>nanans=NaN>>pians=>>realmaxans=+308>>realminans=-3082、MATLAB的sqrt指令可对任一数值进行开平方的运算。用此指令求出下列各数的平方根,并验算之:a)1-ipb)2*i-5+12*i其中i是单位虚数。a)>>sqrt(1-i)ans=->>ans^2ans=-)>>sqrt(2*i)ans=+>>ans^2ans=0+)>>sqrt(-5+12*i)ans=+>>ans^2ans=-+、B作出A的行向量组:a1,a2,a3,a4,a5,a6。>>A=[3,4,-1,1,-9,10;6,5,0,7,4,-16;1,-4,7,-1,6,-8;2,-4,5,-6,12,-8;-3,6,-7,8,-1,1;8,-4,9,1,3,0]A=34-11-91065074-161-47-16-82-45-612-8-36-78-118-49130>>a1=A(1,:)a2=A(2,:)a3=A(3,:)a4=A(4,:)a5=A(5,:)a6=A(6,:)a1=34-11-910a2=65074-16a3=1-47-16-8a4=2-45-612-8a5=-36-78-11a6=8-49130作出B的列向量组:b1,b2,b3,b4,b5,b6。>>B=[1246-32;7916-58-7;81120155;10152813-19;12193625-723;246-305]B=1246-327916-58-781120**********-1912193625-723246-305>>b1=B(:,1)b2=B(:,2)b3=B(:,3)b4=B(:,4)b5=B(:,5)b6=B(:,6)b1=17810122b2=291115194b3=4162028366b4=6-511325-3b5=-385-1-70b6=2-759235由A的一、三、四行和二、三、五列交叉点上的元素作出子矩阵A3。>>A3=A([1,3,4],[2,3,5])A3=4-1-9-476-4512做一个12阶矩阵A4其分块形式为A4=。>>A4=[A,eye(6);zeros(6),B]A4=34-11-910**********-160100001-47-16-80010002-45-612-8000100-36-78-110000108-4913000000**********-320000007916-58-70000008112015500000010152813-1900000012193625-72300000024