高等数学实验第3章.pptVIP

第三章 矩阵代数;3.1 预备知识：线性代数;1、线性方程组;若秩(A) ? 秩(A,b)，则无解； 若秩(A) = 秩(A,b) = n, 存在唯一解； 若秩(A) = 秩(A,b) n, 存在无穷多解； 通解是齐次线性方程组 Ax=0 的基础解 系与 Ax=b 的一个特解之和。;2、逆矩阵;对于方阵A，若存在数?和非零向量x使Ax=?x，则称?为A的一个特征值，x为A 的一个对应于特征值?的特征向量。;3.2 矩阵代数的MATLAB函数;1、矩阵运算符; a=[1+2i 2 3;4 5 6];a. ans = 1.0000 + 2.0000i 4.0000 2.0000 5.0000 3.0000 6.0000 ;2、特殊矩阵的生成;3、矩阵处理;4、行列式、逆矩阵和特征值;解： a=[4 1 -1;3 2 -6;1 -5 3];det(a),inv(a),[v,d]=eig(a);3.3 计算实验：线性方程组求解;1、矩阵除法;例1、应用矩阵除法解下列方程组;? A=[1 2;3 -2:1 -1]; ? B=[1;4;2];x=A \B 求得一最小二乘近似解;2、线性方程组通???;?a=[1 -1 1 -1;-1 1 1 -1;2 -2 -1 1]; b=[1;1;-1]; ? [rank(a),rank([a,b])] ans = 2 2 ?rref([a,b]) ％方法一（行最简形） ans = 1 -1 0 0 0 0 0 1 -1 1 0 0 0 0 0;?x0=a\b, x=null(a) ％方法二(基础解系) x0 = 0 0 1 0 x = -0.7071 0 -0.7071 0 -0.0000 0.7071 -0.0000 0.7071;3、方阵的相似对角化及应用;例3、设;1、投入产出分析;令 C =（cij）， X = (x1, …, xn) , D = (d1, …, dn)’， F= (f1, …, fn)’ 则 X=CX+D；令 A = E－C，E为单位矩阵，则 AX = D;3.4 建模实验：投入产出分析和基因遗传;例4 某地有三个产业，一个煤矿，一个发电厂和一条铁路，开采一元钱的煤，煤矿要支付0.25元的电费及0.25元的运输费; 生产一元钱的电力，发电厂要支付0.65元的煤费，0.05元的电费及0.05元的运输费; 创收一元钱的运输费,铁路要支付0.55元的煤费和0.10元的电费，在某一周内煤矿接到外地金额50000元定货，发电厂接到外地金额25000元定货，外界对地方铁路没有需求。;;则原方程为 X=CX+D ;2、基因遗传;例5、设金鱼某种遗传病染色体的正常基因为A，不正常基因为a, 那么AA，Aa，aa分别表示正常金鱼，隐性患者，显性患者。设初始分布为90%正常金鱼，10%的隐性患者，无显性患者。考虑下列两种配种方案对后代该遗传病基因型分布的影响;解 设初始分布X(1)=(0.9 0.1 0)’,;3.5 扩展实验：病态方程组和矩阵分解;3.6 习题