第四章 线性方程组的应用和矩阵的应用.docVIP

第四章：线性方程组的应用和矩阵的应用 一、线性方程组的应用 MATLAB命令： rref(A) 将矩阵A化为最简行阶梯形矩阵 例１、闭合经济问题 一个木工，一个电工，一个油漆工，三个人相互同意彼此装修他们自己的房子。在装修之前，他们达成如下协议：（1）每人总工作十天（包括给自己家干活在内）；（2）每人的日工资根据一般的市价在60~80元之间；（3）每人的日工资数应使得每人的总收入与总支出相等。下表示他们协商后的制定出的工作天数分配方案： 木工 电工 油漆工 在木工家的工作天数 2 1 6 在电工家的工作天数 4 5 1 在油漆工家的工作天数 4 4 3 问题分析与数学模型 根据协议中每人总支出与总收入相等的原则，分别考虑木工、电工及油漆工的总收入和总支出。设木工的每日工资为x1，电工的日工资为x2，油漆工的日工资为x3。则木工的10个工作日总收入应该是10 x1，而木工、电工、油漆工三人在木工家工作的天数分别为：2天，1天，6天，按日工资累计木工的总支出2 x1+ x2+ 6x3。于是木工的收支平衡可描述为等式 2 x1+ x2+ 6x3=10x1 同理，可建立描述电工、油漆工各自的收支平衡关系的另外两个等式，将三个等式联立，可描述实际问题的方程组。 整理后为： 这是一个齐次线性方程组问题。 算法与数学模型求解 写出齐次线性方程组的系数矩阵如下 为了求出齐次线性方程组的基础解系，将方程组的系数矩阵化为最简行 阶梯形，在MATLAB环境下输入系数矩阵A，然后用命令rref将其化简，键入命令 A=[-8 1 6;4 -5 1;4 4 -7] format rat rref(A) 可得： A = -8 1 6 4 -5 1 4 4 -7 ans = 1 0 -31/36 0 1 -8/9 　　 0 0 0 由此得等价的齐次方程组 根据齐次方程组基础解系的理论，齐次方程组的同解可以表示为 其中，k为任意实数。最后，为了确定满足条件 的方程组的解。即选择恰当的k以确定木工、电工及油漆工每人的日工资：60~80元。取 k=72满足题意得 x1=62,x2=64,x3=72 问题解答 尽管这一问题是在方程组的无穷多组解中寻求解答，但是由于题目条件限制，对于参数k,没有更多的选择余地。为了使日工资为整数值，可确定k=72,使得 木工日工资：62元／日 电工日工资：64元／日 油漆工日工资：72元／日 例２、世界人口预测问题 据统计，20世纪60年代世界人口数据如下表所示，有人根据表中数据，预测2000年世界人口会超过60亿。这一结论在20世纪60年代末令人难以置信，但现在已成为事实。作出这一预测结果所用的方法就是数据拟和方法，根据数学模型，构造出能逼近下表中数据的拟和函数。正是拟和函数反映了人口增长的趋势。根据下表中数据构造拟和函数，预测2000年时的世界人口 世界人口统计表 年 1960 1961 1962 1963 1964 1965 1966 1967 1968 人口 29.72 30.61 31.51 32.13 32.34 32.85 33.56 34.02 34.83 问题分析与数学模型 据人口增长的统计资料和人口理论数学模型知，当人口总数N不是很大时，在不太长的时期内，人口增长接近于指数增长。因此，采用指数函数 N=ea+bt 对数据进行拟和。为了计算方便，将上式两边同取对数，得lnN=a+bt,令 y=lnN或N=ey 变换后的拟和函数为 y(t)=a+bt 由人口数据对数(y=lnN)计算，得表如下 t 1960 1961 1962 1963 1964 1965 1966 1967 1968 y 3.3918 3.4213 3.4503 3.4698 3.4763 3.4920 3.5133 3.5322 3.5505 算法与数学模型求解 根据表中数据及等式a+btk=yk(k=1,2,……,9)可列出关于两个未知数a、b的９个方程的超定方程组（方程数多于未知数个数的方程组） 输入下面命令 format short; t=1960:1968 y=[3.3918 3.4213 3.4503 3.4698 3.4763 3.4920 3.5133 3.5322 3.5505] A=[ones(9,1) t] A\y 直接求解方程组