资源描述:
《matlab实验1matlab初步》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、实验1Matlab初步一、问题已知矩阵A、B、b如下:应用Matlab软件进行矩阵输入及各种基本运算。二、实验目的学会使用Matlab软件构作已知矩阵对应的行(列)向量组、子矩阵及扩展矩阵,实施矩阵的初等变换及线性无关向量组的正交规范化,确定线性相关相关向量组的一个极大线性无关向量组,且将其余向量用极大线性无关向量组线性表示,并能编辑M文件来完成所有的实验目的。三、预备知识1、线性代数中的矩阵及其初等变换、向量组的线性相关性等知识。2、Matlab软件的相关命令提示如下;(1)选择A的第i行做一个行向量:a
2、i=A(i,:);(2)选择A的第j行做一个列向量:ai=A(j,:);(3)选择A的某几行、某几列上的交叉元素做A的子矩阵:A([行号],[列号]);(4)n阶单位阵:eye(n);n阶零矩阵:zeros(n);(5)做一个n维以0或1为元素的索引向量L,然后取A(:,L),L中值为1的对应的列将被取到。(6)将非奇异矩阵A正交规范化,orth(A);验证矩阵A是否为正交阵,只需做A*A'看是否得到单位阵E。(7)两个行向量a1和a2的内积:a1*a2'。(8)让A的第i行与第j列互换可用赋值语句:A([
3、i,j],:)=A([j,i],:);(9)让K乘以A的第i行可用赋值语句:A(i,:)=K*A(i,:);(10)让A的第i行加上第j行的K倍可用赋值语句:A(i,:)=A(i,:)+K*A(j,:);(11)求列向量组的A的一个极大线性无关向量组可用命令:rref(A)将A化成阶梯形行的最简形式,其中单位向量对应的列向量即为极大线性无关向量组所含的向量,其它列向量的坐标即为其对应向量用极大线性无关组线性表示的系数。四、实验的内容与要求1、请在MATLAB直接输入下列常数,看它们的值是多少:aibjcep
4、sdinfenanfpigrealmaxhrealmin>>ians=0+1.0000i>>jans=0+1.0000i>>epsans=2.2204e-016>>infans=Inf>>nanans=NaN>>pians=3.1416>>realmaxans=1.7977e+308>>realmin2、MATLAB的sqrt指令可对任一数值进行开平方的运算。用此指令求出下列各数的平方根,并验算之:a)1-ipb)2*ic)-5+12*i其中i是单位虚数。>>sqrt(1-i)>>ans^2>>sqrt(2
5、*i)>>ans^2>>sqrt(-5+12*i)>>ans^21、输入矩阵A、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]>>a1=A(1,:)a2=A(2,:)a3=A(3,:)a4=A(4,:)a5=A(5,:)a6=A(6,:)1、作出B的列向量组:b1,b2,b3,b4,b5,b6。>>B=[1246-
6、32;7916-58-7;81120155;10152813-19;12193625-723;246-305]>>b1=B(:,1)b2=B(:,2)b3=B(:,3)b4=B(:,4)b5=B(:,5)b6=B(:,6)2、由A的一、三、四行和二、三、五列交叉点上的元素作出子矩阵A3。>>A3=A([1,3,4],[2,3,5])3、做一个12阶矩阵A4其分块形式为A4=。>>A4=[A,eye(6);zeros(6),B]4、由索引向量L产生取A的第二、四、五行所成的子矩阵A5。将A的对应的行向量组正交
7、规范为正交向量组A6,并验证所得的结果。>>L=[2,4,5];A5=A(L,:)>>AL=[0,0,0,0,0,0;1,1,1,1,1,1;0,0,0,0,0,0;1,1,1,1,1,1;1,1,1,1,1,1;0,0,0,0,0,0]A5=reshape(A(find(L)),3,6)>>A6=orth(A)>>A6*A6'5、求a1与a2的内积A7。>>A7=a1*a2'完成以下初等变换:将A的第一、四行互换,再将其第三列乘以6,再将其第一行的10倍加至第五行。>>A>>A([1,4],:)=A([4
8、,1],:)>>A(:,3)=6*A(:,3)>>A(5,:)=10*A(1,:)+A(5,:)6、求B的列向量的一个极大线性无关向量组A9,并将其余的向量用极大线性无关向量组线性表示。>>B2=rref(B)>>A9=B2(:,[1,2,3,6])>>b4=B2(:,4)b5=B2(:,5)>>A9b4>>A9b54、对矩阵A进行下列操作(如不清楚指令可用命令lookfor查找)a)找出矩阵的大小(即行维数