matlab矩阵代数.pptVIP

第三章 矩阵代数 MATLAB数学实验 A 矩阵代数 对于方程Ax＝b，A 为n×m矩阵，求解时分三种情况： ? 当n=m时，此方程成为“恰定”方程 ? 当nm时，此方程成为“超定”方程 ? 当nm时，此方程成为“欠定”方程 线性方程组的解 若秩(A) ? 秩(A,b)，则无解； 若秩(A) = 秩(A,b) = n, 存在唯一解； 若秩(A) = 秩(A,b) n, 存在无穷多解； 通解是齐次线性方程组 Ax=0 的基础解系与 Ax=b 的一个特解之和。 方程Ax=b A=[1 2;2 3];b=[8;13]; ?x=inv(a)*b ? x=A\b x = x = 2.00 2.00 3.00 3.00 2.超定方程组的解 方程 Ax=b ,mn时此时不存在唯一解。 方程解 (AA)x=A b x=(A‘ A)-1A’ b —— 求逆法给出一个基本解 x=A\b —— matlab用最小二乘法找一个基本解。 例: x1 +2x2=1 2x1+ 3x2=2 3x1+ 4x2=3 a=[1 2;2 3;3 4];b=[1;2;3]; 解1 x=A\b 解2 x=inv(A?A) ? A ? b x = 1.00 0 3.欠定方程组的解 当方程数少于未知量个数时,即不定情况,有无穷多个解存在。 matlab可求出两个解： 具有最多零元素的解(除法):x=A\b 具有最小长度或范数的解，这个解是基于伪逆pinv求得的: x=pinv(A)*b x1+2x2+3x3=1 2x1+3x2+4x3=2 A=[1 2 3;2 3 4];b=[1;2]; x=A\b x=pinv(A)?b x = x = 1.00 0.83 0 0.33 0 -0.17 例2 线性方程组通解 用rref化为行最简形以后求解 用除法求出一个特解，再用null求得一个齐次组的基础解系 相似对角化 如果n阶方阵A有n个线性无关的特征向量，则必存在可逆矩阵P, 使得 P-1AP= ?, 其中?是A的特征值构成的对角矩阵。P的列向量是对应的n个正交特征向量。（[P, ?]=eig(A)） 使用MATLAB函数eig求得的每个特征向量都是单位向量，并且属于同一特征值的线性无关特征向量已正交化。 例3 用相似变换矩阵P将A相似对角化，并求 ，这里 三 建模举例 设有n个经济部门，xi为部门i的总产出，cij为部门j单位产品对部门i产品的消耗，di为外部对部门i的需求，fj为部门j新创造的价值。 分配平衡方程组 消耗平衡方程组 i,j =1,2,…,n 令 C =（cij），X = (x1, …, xn),D = (d1, …, dn)’，F= (f1, …, fn)’, 则 X=CX+D 令 A = E－C，E为单位矩阵，则 AX = F C称为直接消耗矩阵 A称为列昂杰夫（Leontief）矩阵。 例 某地有三个产业，一个煤矿，一个发电厂和一条铁路，开采一元钱的煤，煤矿要支付0.25元的电费及0.25元的运输费; 生产一元钱的电力，发电厂要支付0.65元的煤费，0.05元的电费及0.05元的运输费; 创收一元钱的运输费,铁路要支付0.55元的煤费和0.10元的电费。在某一周内煤矿接到外地金额50000元定货，发电厂接到外地金额25000元定货，外界对地方铁路没有需求。问三个企业间一周内总产值多少才能满足自身及外界需求？三个企业间相互支付多少金额？三个企业各创造多少新价值？ 例 x1+2x2=8 2x1+3x2=13