实验1矩阵方法实验.doc

PAGE  PAGE 36 实验1 矩阵方法实验 实验内容： 矩阵运算和基于矩阵的代数运算 二、实验目的： 熟悉有关矩阵运算的各种命令。 能熟练地进行代数运算，包括计算矩阵的加、减、乘、逆和方阵的行列式等。 能求解线性方程组。 能运用矩阵方法求解代数问题。 三、基本知识 MATLAB的操作对象是矩阵,标量为1×1的矩阵,向量为1×n矩阵,多项式也可用矩阵表示. 1．矩阵输入 矩阵输入有两种方式： (1)用中括号表示,每行元素间用逗号或空格分开,行与行之间用分号隔开； (2)a=初始值:步长:终值,可输入行矩阵. 例如输入 a=[1 2 3; 4 5 6; 7 8 9] 或 a=[1,2,3;4,5,6;7,8,9] 会得到同样的结果: a= 1 2 3 4 5 6 7 8 9 输入 a=2:2:10 得到从2到10的以2为公差的数组: a = 2 4 6 8 10 2．多项式的表示法和运算 MATLAB用行向量来表示一个多项式.例如 可用矩阵表示为: p=[1 0 –1 6] MATLAB提供的多项式运算函数名见(表3-1). 表3-1 函数名含义 root求多项式的根 poly已知方程的根,还原方程 polyval计算多项式的值  solve求方程的解polyder多项式的微分polyfit数据的多项式曲线拟合factor因式分解例1 求解方程 . 解 输入如下命令 syms x s=solve(x^3-2*x^2=x-1) 或 syms x s= x^3-2*x^2-x+1; solve(s) 求多项式的微分. 解 输入 p=[1 -2 0 -1]; q=polyder(p) 得 q = 3 -4 0 3．矩阵运算 1.MATLAB提供的一些特殊矩阵,见(表3-2). 表3-2 函数含义 Eye单位矩阵 Company伴随矩阵 Rand产生均匀分布随机数矩阵 zeros零矩阵 ones元素均为1的矩阵 randn产生正态分布随机矩阵例如rand(3,9)产生一3行9列的矩阵,其元素数值范围为(0,1). 2.运算符,见(表3-3). 表3-3 操作符含义+加-减*矩阵乘\矩阵左除/矩阵右除^矩阵的冥.*数组乘.\数组左除./数组右除.^数组乘方3.关系运算符,见(表3-4). 表3-4 操作符含义= =等于~=不等于大于小于=小于等于=大于等于逻辑和~逻辑非|逻辑或4.矩阵运算的命令符,见(表3-5). 表3-5 语法含义AˊA的转置det(A)A的行列式rank(A)A的秩A(i,:)=k*A(i,:)k乘以A的第i行A([i,j],:) =A([j.i],:)A的第i行和第j行互换aj=A(:,j)选择A的第j列做一列向量inv(A)A的逆eig(A)A的特征值poly(A)A的特征多项式Ai=A(i,:)选择A的第i行做一个行向量A(i,:)=A(i,:)+k*A(j,:)A的第i行加上第j行的k倍a1*a2ˊ两行向量的内积若方程组为,且存在,则.因此可用或计算线性方程组的解. 例3 已知矩阵A、B、b如下 (1)输入矩阵A、B、b; (2)求X1=A的转置,X2=A+B,X3=AB; (3)求X4=A的行列式,X5=A的秩; (4)求X6=A的逆; (5)求方程组AX=b的解向量. 解 完成上述求解的命令为： (1)A=[3 4 –1；6 5 0;1 –4 7] B=[1 3 4;7 9 16;8 11 20] b=1:2:5 (2)x1=A x2=A+B x3=A*B 得 x1 = 3 6 1 4 5 -4 -1 0 7 x2 = 4 7 3 13 14 16 9 7 27 x3 = 23 34 56 41 63 104 29 44 80 (3)x4=det(A) x5=rank(A) 得 x4 = -34 x5 = 3 (4)x6