数学实验第三章A.pptVIP

* * 第三章 矩阵代数 解 设初始分布X(1)=(0.9 0.1 0)’, 第n代分布为X(n)= A = B= 则 X(n) = An-1X(1) X(n) = Bn-1X(1) 分别是 两种情况下第n代的基因型分布 * * 第三章 矩阵代数 结论 计算方法：取n足够大n=20 方案一：An-1X(1) ?[0.95, 0, 0.05]’ 存在5%的显性患者 方案二： Bn-1X(1) ?[1, 0, 0]’ 显性患者和隐性患者都消失 本例说明了杂交的优势 * * 第三章 矩阵代数 习题 P59 ex2, ex3, ex4, ex5, ex6, ex10 第三章 矩阵代数 第三章 矩阵代数 MATLAB数学实验 第三章 矩阵代数 * * 第三章 矩阵代数 第三章 矩阵代数 3.1 预备知识：线性代数 3.2 矩阵代数的MATLAB指令 3.3 计算实验：线性方程组求解 3.4 建模实验：投入产出分析和基因遗传 * * 第三章 矩阵代数 3.1 预备知识：线性代数 线性方程组 记为 A x = b * * 第三章 矩阵代数 3.1 预备知识：线性代数 线性方程组 若秩(A) ? 秩(A,b)，则无解； 若秩(A) = 秩(A,b) = n, 存在唯一解； 若秩(A) = 秩(A,b) n, 存在无穷多解； 通解是齐次线性方程组 Ax=0 的基础解系与 Ax=b 的一个特解之和。 * * 第三章 矩阵代数 3.1 预备知识：线性代数 逆矩阵 方阵A称为可逆的，如果存在方阵B，使A B = B A = E，记 B = A-1 方阵A可逆的充分必要条件:?A??0 A-1 =A*/|A| 这里A*为A的伴随矩阵 （A E） 行变换 （E A-1） * * 第三章 矩阵代数 3.1 预备知识：线性代数 特征值与特征向量 对于方阵A，若存在数?和非零向量x 使 A x = ? x，则称?为A的一个特征值，x 为A 的一个对应于特征值?的特征向量。 特征值计算归结为特征多项式的求根。 特征向量计算：齐次线性方程组 (A - ?E) x = 0 的所有一组线性无关解。 * * 第三章 矩阵代数 3.2 矩阵代数的MATLAB指令 运算符 A’ (共轭)转置, A.’ 转置 A+B与A-B 加与减 k+A与k-A 数与矩阵加减 k*A或A*k 数乘矩阵 A*B 矩阵乘法 A^k 矩阵乘方 左除A\B 为AX=B的解 右除B/A 为XA=B的解 与数组运算不同 * * 第三章 矩阵代数 3.2 矩阵代数的MATLAB指令 矩阵运算与数组运算的区别 数组运算按元素定义，矩阵运算按线性代数定义 矩阵的加、减、数乘等运算与数组运算是一致的 矩阵的乘法、乘方和除法与数组乘法、乘方和除法不同 数与矩阵加减、矩阵除法在数学上是没有意义的。但在MATLAB中有定义。 例子 P45-46 * * 第三章 矩阵代数 3.2 矩阵代数的MATLAB指令 特殊矩阵生成 zeros(m,n) m行n列的零矩阵; ones(m,n) m行n列的元素全为1的阵; eye(n) n阶单位矩阵; rand(m,n) m行n列[0,1]上均匀分布随机数矩阵 randn ：产生均值为0，方差为1的标准正态分布随机矩阵。 * * 第三章 矩阵代数 zeros生成全部元素为0的零矩阵 A=zeros(n) 生成n×n零矩阵 A=zeros(m,n) 或 zeros([m,n]) 生成m×n零矩阵 A=zeros(m,n,p,… ) 生成m×n×p的零矩阵B=zeros(size(A)) 生成和矩阵A大小相等的全零矩阵。 举例: * * 第三章 矩阵代数 例2-3 分别建立3×3、3×2和与矩阵A同样大小的零矩阵。(1) 建立一个3×3零矩阵。zeros(3) (2) 建立一个3×2零矩阵。zeros(3,2) (3) 设A为2×3矩阵，则可以用zeros(size(A))建立一个与矩阵A同样大小零矩阵。A=[1 2 3;4 5 6]; %产生一个2×3阶矩阵Azeros(size(A)) %产生一个与矩阵A同样大小的零矩阵 * * 第三章 矩阵代数 3.2 矩阵代数的MATLAB指令 矩阵处理 trace(A) 迹(对角线元素的和) diag(A) A对角线元素构成的向量; diag(x) 向量x的元素构成的对角矩阵. tril(A) A的下三角部分 triu(A) A的上三角部分 flipud(A) 矩阵上下翻转 fliplr(A)