ch4_2例题.docVIP

矩阵和代数方程 矩阵运算和特征参数 矩阵运算 表4.2-1 矩阵运算含义及相应符号 术 语 数 学 含 义 MATLAB表达 矩阵加减 A+B, A-B 标量与矩阵加减 a+B, a-B 矩阵乘积 A*B 标量与矩阵相乘 a*B 矩阵的转置 A 〖说明〗 在MATLAB中，“标量与矩阵加减运算”是一种“行业性”说法，在教科书中没有这种表述。本表格数学含义栏对“标量与矩阵加减运算”给出了完整的定义。 MATLAB的矩阵运算是在复数域定义的。在本表所列操作中，矩阵转置可全称为“矩阵的共轭转置（Hermition transpose）”，在数学表述上采用上标H表示。 【例 4.2-1】矩阵乘法算例。本例演示：矩阵乘法的含义；如何减少和避免循环；如何判断两个双精度矩阵是否相等。 （1）实现矩阵乘的“标量程式” clear rand(state,12) A=rand(2,4);B=rand(4,3); %---------用三重循环体实现矩阵乘法--------- C1=zeros(size(A,1),size(B,2)); %为C预分配内存，可加快计算 4 for ii=1:size(A,1) for jj=1:size(B,2) for k=1:size(A,2) C1(ii,jj)=C1(ii,jj)+A(ii,k)*B(k,jj); end end end % 13 （2）实现矩阵乘的“saxpy程式”（saxpy是LAPACK术语） %---------用二重循环体实现矩阵乘法--------- C2=zeros(size(A,1),size(B,2)); for jj=1:size(B,2) for k=1:size(B,1) C2(:,jj)=C2(:,jj)+A(:,k)*B(k,jj); end end （3）MATLAB中矩阵乘的最好实现方式 C3=A*B; %直接运用MATLAB的矩阵乘 （3）计算结果比较 C3_C3=norm(C3-C1,fro) %采用计算矩阵F范数的方式比较两个矩阵之间的“差” C3_C2=norm(C3-C2,fro) 〖说明〗 “标量程式”的效率最差，“saxpy”是MATLAB提高计算效率一系列措施中的一个，直接采用MATLAB中的矩阵乘算符“*”实现矩阵乘最好。 由于数值计算是有限精度计算，不同计算程式所得结果的比较一般采用矩阵范数进行。 【例 4.2-2】观察矩阵的转置操作和数组转置操作的差别。本例演示：共轭转置和非共轭转置；有理表示的显示格式 p.4上的思考题 本章习题都要做，下次课堂抽查。 矩阵的标量特征参数 表4.2-2 计算矩阵标量特征参数的指令 术 语 数 学 含 义 MATLAB指令 秩 Rank 可采用以下任何一种表述： 矩阵A中线性无关列（或行）向量组中最大向量数； 矩阵A中最高非零子行列式的阶数； 矩阵A中最高非奇异子矩阵的维数。 rank(A) 迹 Trace ，即矩阵主对角元素之和。 trace(A) 行列式Determinant ， 式中是元素对应的子行列式 det(A) 【例4.2-3】矩阵标量特征参数计算示例。本例演示：rank, det, trace的使用；子行列式的计算。 【例4.2-4】矩阵标量特征参数的性质。本例演示：两相乘矩阵位置交换迹不变；两同阶矩阵相乘位置交换行列式不变。 format short g %采用默认显示格式 rand(state,0) %设置均匀分布随机发生器的初始状态 A=rand(3,3); %生成(3*3)随机阵 B=rand(3,3); %生成另一个(3*3)随机阵 C=rand(3,4); D=rand(4,3); %任何符合矩阵乘法规则的两个矩阵的乘积的“迹”不变 tAB=trace(A*B) %两个同阶矩阵乘积 tBA=trace(B*A) tCD=trace(C*D) %两个“内维”相等矩阵的乘积 tDC=trace(D*C) tAB = 3.6697 tBA = 3.6697 tCD = 2.1544 tDC = 2.1544 %同阶矩阵乘积行列式等于各矩阵行列式之乘积 d_A_B=det(A)*det(B) dAB=det(A*B) dBA=det(B*A) d_A_B = 0.0846 dAB = 0.0846 dBA = 0.0846 %对于非同阶矩阵，两个矩阵的乘积行列式随相乘次序不同而不同