最新版ch4_2例题.（完整版）.docx

PAGE PAGE # 4.1 矩阵和代数方程 矩阵运算和特征参数 矩阵运算 表4.2-1矩阵运算含义及相应符号 术语 数学含义 MATL AB表 达 矩阵加 减 A m>n 亠 Bm>n — [ aij 亠 hj ]mxn A+B, A-B 标量与 矩阵加 减 a — Bm>n = [a]m：*n 亠[hj] mxn = [a 亠 hj ] m<n a+B, a-B 矩阵乘 积 Cmxn _[ cij ] _ bij 'mX 'hj〔翅 一 Q aik bkj _ A mxj BI xn -kT mxn A*B 标量与 矩阵相 乘 aBm>h — [a 'hj ]mxh a*B 矩阵的 转置 H B = A = (bq )r + i，(bkj)i =(ajk)r — i，(ajk) i A' 〖说明〗 在MATLAB中，“标量与矩阵加减运算”是一种“行业性”说法， 在教科书中没有这种表述。本表格数学含义栏对 标量与矩阵加 减运算”给出了完整的定义。 MATLAB的矩阵运算是在复数域定义的。在本表所列操作中，矩 阵转置可全称为 矩阵的共轭转置(Hermition transpose ，在 数学表述上米用上标H表示。 【例4.2-1】矩阵乘法算例。 本例演示：矩阵乘法的含义；如何减少和避免循环；如何判断 两个双精度矩阵是否相等。 1)实现矩阵乘的“标量程式” clear ran d('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 min( m ,n ) 瓦 a ii ，即矩阵主对角元素之 i =1 和。 trace(A) 行列式 Determinan t n A俯 Y ei)H1a1j A!j， j m 式中Aij是元素aij对应的子行列式 det(A) 【例4.2-3】矩阵标量特征参数计算示例。 本例演示：rank, det, trace的使用；子行 列式的计算。 【例4.2-4】矩阵标量特征参数的性质。 本例演示：两相乘矩阵位置交换迹不变；两 同阶矩阵相乘位置交换行列式不变。 format short g ran d('state',0) A=ra nd(3,3); B=ra nd(3,3); C=ra nd(3,4); D=ra nd(4,3); %采用默认显示格式 %设置均匀分布随机发生器的初始状态 %生成(3*3)随机阵 %生成另一个(3*3)随机阵  %任何符合矩阵乘法规则的两个矩阵的乘积的“迹”不变 %两个同阶矩阵乘积%两个“内维”相等矩阵的乘积tAB二 %两个同阶矩阵乘积 %两个“内维”相等矩阵的乘积 3.6697 tBA = 3.6697 tCD = 2.1544 tDC = 2.1544 %同阶矩阵乘积行列式等于各矩阵行列式之乘积 d_A_B=de