实验1 Matlab初步(1).doc

实验1 Matlab初步 问题 已知矩阵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 直接输入下列常数，看它们的值是多少： a i i ans = 0 + 1.0000i b j j ans = 0 + 1.0000i c eps eps ans = 2.2204e-016 d inf inf ans = Inf e nan nan ans = NaN f pi pi ans = 3.1416 g realmax realmax ans = 1.7977e+308 h realmin realmin ans = 2.2251e-308 2、MATLAB 的 sqrt 指令可对任一数值进行开平方的运算。用此指令求出下列各数的平方根，并验算之： A? sqrt(1-i) ans = 1.0987 - 0.4551i b) 2*i sqrt(2*i) ans = 1.0000 + 1.0000i c) -5+12*i sqrt(-5+12*i) ans = 2.0000 + 3.0000i 其中 i 是单位虚数。 输入矩阵A、B 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 = 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 B=[1 2 4 6 -3 2;7 9 16 -5 8 -7;8 11 20 1 5 5;10 15 28 13 -1 9;12 19 36 25 -7 23;2 4 6 -3 0 5;] B = 1 2 4 6 -3 2 7 9 16 -5 8 -7 8 11 20 1 5 5 10 15 28 13 -1 9 12 19 36 25 -7 23 2 4 6 -3 0 5 作出A的行向量组：a1，a2，a3，a4，a5，a6。 a1=A(1,:) a1 = 3 4 -1 1 -9 10 a2=A(2,:) a2 = 6 5 0 7