四种线性代数模型
- 格式:doc
- 大小:217.00 KB
- 文档页数:5
常见的线性代数求解方法
1.列主元消去法
列主元消去法是一种经典的求解线性方程组的方法。
它通过将
方程组转化为上三角矩阵的形式来求解。
这个方法的关键在于选取
主元的策略。
一种常见的选取主元的策略是选择当前列中绝对值最
大的元素作为主元,然后进行消去操作,直到将矩阵转化为上三角
矩阵。
2.高斯-约当消去法
高斯-约当消去法是另一种常见的线性方程组求解方法。
它通
过消去矩阵的下三角部分来将线性方程组转化为上三角矩阵的形式。
这个方法也需要选择主元,常见的选择策略是选取当前行中绝对值
最大的元素作为主元,然后进行消去操作。
3.LU分解法
LU分解法是将矩阵分解为一对矩阵的乘积的方法。
这个方法的思想是先将矩阵分解为一个下三角矩阵和一个上三角矩阵,然后通过求解上三角矩阵和下三角矩阵的两个方程组来求解原始的线性方程组。
4.Jacobi迭代法
Jacobi迭代法是一种迭代求解线性方程组的方法。
它通过将原始的线性方程组转化为一个对角矩阵和另一个矩阵的乘积的形式,然后通过迭代求解这个对角矩阵和另一个矩阵的方程组来逼近线性方程组的解。
5.Gauss-Seidel迭代法
Gauss-Seidel迭代法是另一种迭代求解线性方程组的方法。
它与Jacobi迭代法类似,但是在每一次迭代中,它使用前一次迭代得到的部分解来更新当前的解。
这个方法通常比Jacobi迭代法收敛得更快。
以上是一些常见的线性代数求解方法。
每种方法都有其特点和适用范围,我们可以根据具体情况选择合适的方法来求解线性方程组的问题。
线性分类模型(一):线性判别模型分析前言前几篇文章介绍了线性回归算法,线性分类模型分为判别式模型和生成式模型,本文首先简单复习了与算法相关的数学基础知识,然后分析各线性判别式分类算法,如最小平方法,Fisher线性判别法和感知器法,最后总结全文。
目录1、相关的数学知识回顾2、判别式模型和生成性模型3、最小平方法4、Fisher线性判别函数5、感知器算法6、总结相关数学知识回顾1、直线方程和平面方程拓展到分类思想:直线l为分类决策方程,坐标点落在直线l上方时,则分类为C1;坐标点落在直线l下方时,则分类为C2(如下图)。
平面方程类似,在这里不展开。
2、点到直线和点到平面的距离点到直线的距离:点到平面的距离拓展到分类思想:平面方程为决策方程,正确分类的情况下,当点P到决策方程的距离越大,则分类模型越好;错误分类的情况下,点P到决策方程的距离作为损失函数,损失函数最小化过程即是模型参数最优化过程。
3、向量内积的数学意义几何意义:向量A与向量B的内积等于向量A在向量B的投影与向量B的乘积,当向量B是单位向量时,则等于向量A在单位向量方向的投影,单位向量类似于基函数或者可以理解成坐标轴,即向量A在向量B的投影可理解成向量A在向量B方向的坐标,如下图,B'是B 在OA坐标轴方向的投影。
拓展到分类思想:C1与C2属于不同的类,给定一条决策性直线l,当C1与C2在直线L2的投影间距越大,则分类效果越好。
增加不同类间的距离可以作为模型参数优化的方向。
如下图,C1和C2的在直线L2的投影距离|C1'C2'|大于|C1'C2'|,因此决策方程直线L2优于直线L1.4、梯度下降法函数f(x0,x1,...,xn)在梯度方向是函数值变化(增加或减少)最快的方向(本文只给出结论,后续文章会有详细的说明)。
拓展到分类思想:损失函数最小化过程即是模型参数最优化过程,损失函数最小化可通过梯度下降法来实现,当迭代到一定程度,损失函数收敛,则迭代结束,参数w即是要求的最优参数。
线性代数是高等学校理工科和经济类学科相关专业的一门重要基础课,它不仅是其他数学课程的基础,也是物理、力学、电路等专业课程的基础。
作为处理离散问题工具的线性代数,也是从事科学研究和工程设计的科研人员必备的数学工具之一。
实验一 生物遗传模型1.工程背景设一农业研究所植物园中某植物的基因型为AA 、Aa 和aa 。
常染色体遗传的规律是:后代是从每个亲体的基因对中个继承一个基因,形成自己的基因对。
如果考虑的遗传特征是由两个基因A 、a 控制的,那末就有三种基因对,记为AA 、Aa 和aa 。
研究所计划采用Aa(AA)型的植物与每一种基因型植物相结合的方案培育植物后代。
问经过若干年后,这种植物的任意一代的三种基因型分布如何?2.问题分析分析双亲体结合形成后代的基因型概率,如表6-4所示。
表6-4基因型概率矩阵 后代 基因对 父体—母体的基因对AA —AAAA —Aa AA —aa Aa —Aa Aa —aa aa —aa AA 1 1/2 0 1/4 0 0 Aa 0 1/2 1 1/2 1/2 0 aa1/41/213.模型建立与求解设,,n n n a b c 分别表示第n 代植物中基因型AA 、Aa 、aa 型的植物占植物总数的百分率。
则第n 代植物的基因型分布为()n n n n a x b c ⎛⎫ ⎪= ⎪ ⎪⎝⎭,0(0)00a x b c ⎛⎫ ⎪= ⎪ ⎪⎝⎭表示植物型的初始分布。
依据上述基因型概率矩阵,有1112n n n a a b --=+,1112n n n b b c --=+,0n c =,1n n n a b c ++=,表示为矩阵形式11111/2001/21000n n n n n n a a b b c c ---⎛⎫⎛⎫⎛⎫⎪ ⎪⎪= ⎪ ⎪⎪ ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭记11/2001/21000M ⎛⎫ ⎪= ⎪ ⎪⎝⎭,则()(1)2(2)3(3)(0)n n n n n x MxM x M x M x ---=====。
数学建模四大模型总结1优化模型1.1 数学规划模型线性规划、整数线性规划、非线性规划、多目标规划、动态规划。
1.2 微分方程组模型阻滞增长模型、SARS传播模型。
1.3 图论与网络优化问题最短路径问题、网络最大流问题、最小费用最大流问题、最小生成树问题(MST)、旅行商问题(TSP)、图的着色问题。
1.4 概率模型决策模型、随机存储模型、随机人口模型、报童问题、Markov链模型。
1.5 组合优化经典问题l 多维背包问题(MKP)背包问题:个物品,对物品,体积为,背包容量为。
如何将尽可能多的物品装入背包。
多维背包问题:个物品,对物品,价值为,体积为,背包容量为。
如何选取物品装入背包,是背包中物品的总价值最大。
多维背包问题在实际中的应用有:资源分配、货物装载和存储分配等问题。
该问题属于难问题。
l 二维指派问题(QAP)工作指派问题:个工作可以由个工人分别完成。
工人完成工作的时间为。
如何安排使总工作时间最小。
二维指派问题(常以机器布局问题为例):台机器要布置在个地方,机器与之间的物流量为,位置与之间的距离为,如何布置使费用最小。
二维指派问题在实际中的应用有:校园建筑物的布局、医院科室的安排、成组技术中加工中心的组成问题等。
l 旅行商问题(TSP)旅行商问题:有个城市,城市与之间的距离为,找一条经过个城市的巡回(每个城市经过且只经过一次,最后回到出发点),使得总路程最小。
l 车辆路径问题(VRP)车辆路径问题(也称车辆计划):已知个客户的位置坐标和货物需求,在可供使用车辆数量及运载能力条件的约束下,每辆车都从起点出发,完成若干客户点的运送任务后再回到起点,要求以最少的车辆数、最小的车辆总行程完成货物的派送任务。
TSP问题是VRP问题的特例。
l 车间作业调度问题(JSP)车间调度问题:存在个工作和台机器,每个工作由一系列操作组成,操作的执行次序遵循严格的串行顺序,在特定的时间每个操作需要一台特定的机器完成,每台机器在同一时刻不能同时完成不同的工作,同一时刻同一工作的各个操作不能并发执行。
线性代数知识点总结(第4章)(一)方程组的表达形与解向量1、解的形式:(1)一般形式(2)矩阵形式:Ax=b;(3)向量形式:A=(α1,α2,…,αn)2、解的定义:若η=(c1,c2,…,c n)T满足方程组Ax=b,即Aη=b,称η是Ax=b的一个解(向量)(二)解的判定与性质3、齐次方程组:(1)只有零解←→r(A)=n(n为A的列数或是未知数x的个数)(2)有非零解←→r(A)<n4、非齐次方程组:(1)无解←→r(A)<r(A|b)←→r(A)=r(A)-1(2)唯一解←→r(A)=r(A|b)=n(3)无穷多解←→r(A)=r(A|b)<n5、解的性质:(1)若ξ1,ξ2是Ax=0的解,则k1ξ1+k2ξ2是Ax=0的解(2)若ξ是Ax=0的解,η是Ax=b的解,则ξ+η是Ax=b的解(3)若η1,η2是Ax=b的解,则η1-η2是Ax=0的解【推广】(1)设η1,η2,…,ηs是Ax=b的解,则k1η1+k2η2+…+k sηs为Ax=b的解(当Σk i=1)Ax=0的解(当Σk i=0)(2)设η1,η2,…,ηs是Ax=b的s个线性无关的解,则η2-η1,η3-η1,…,ηs-η1为Ax=0的s-1个线性无关的解。
变式:①η1-η2,η3-η2,…,ηs-η2②η2-η1,η3-η2,…,ηs-ηs-1(三)基础解系6、基础解系定义:(1)ξ1,ξ2,…,ξs是Ax=0的解(2)ξ1,ξ2,…,ξs线性无关(3)Ax=0的所有解均可由其线性表示→基础解系即所有解的极大无关组注:基础解系不唯一。
任意n-r(A)个线性无关的解均可作为基础解系。
★7、重要结论:(证明也很重要)设A施m×n阶矩阵,B是n×s阶矩阵,AB=O(1)B的列向量均为方程Ax=0的解(2)r(A)+r(B)≤n(第2章,秩)8、总结:基础解系的求法(1)A为抽象的:由定义或性质凑n-r(A)个线性无关的解(2)A为数字的:A→初等行变换→阶梯型自由未知量分别取1,0,0;0,1,0;0,0,1;代入解得非自由未知量得到基础解系(四)解的结构(通解)9、齐次线性方程组的通解(所有解)设r(A)=r,ξ1,ξ2,…,ξn-r为Ax=0的基础解系,则Ax=0的通解为k1η1+k2η2+…+k n-rηn-r (其中k1,k2,…,k n-r为任意常数)10、非齐次线性方程组的通解设r(A)=r,ξ1,ξ2,…,ξn-r为Ax=0的基础解系,η为Ax=b的特解,则Ax=b的通解为η+ k1η1+k2η2+…+k n-rηn-r (其中k1,k2,…,k n-r为任意常数)(五)公共解与同解11、公共解定义:如果α既是方程组Ax=0的解,又是方程组Bx=0的解,则称α为其公共解12、非零公共解的充要条件:方程组Ax=0与Bx=0有非零公共解←→有非零解←→13、重要结论(需要掌握证明)(1)设A是m×n阶矩阵,则齐次方程A T Ax=0与Ax=0同解,r(A T A)=r(A)(2)设A是m×n阶矩阵,r(A)=n,B是n×s阶矩阵,则齐次方程ABx=0与Bx=0同解,r(AB)=r(B)。
解线性方程组的迭代法Haha送给需要的学弟学妹摘要:因为理论的分析表明,求解病态的线性方程组是困难的,但是实际情况是否如此,需要我们来具体检验。
系数矩阵H 为Hilbert 矩阵,是著名的病态问题。
因而决定求解Hx b =此线性方程组来验证上述问题。
详细过程是通过用Gauss 消去法、J 迭代法、GS 迭代法和SOR 迭代法四种方法求解Hx b =线性方程组。
关键词:病态方程组、Gauss 消去法、J 迭代法、GS 迭代法、SOR 迭代法目录:一、问题背景介绍二、建立正确额数学模型 三、求解模型的数学原理1、Gauss 消去法求解原理2、Jacobi 迭代法求解原理3、G-S 迭代法求解原理4、SOR 迭代法求解原理5、Jacobi 和G-S 两种迭代法收敛的充要条件 四、计算过程(一)Hilbert 矩阵维数n=6时1、Gauss 消去法求解2、Jacobi 迭代法求解3、G-S 迭代法求解4、SOR 迭代法求解(二)Hilbert 矩阵维数n=20、50和100时1、G-S 迭代法求解图形2、SOR 迭代法求解图形 五、编写计算程序 六、解释计算结果1、Gauss 消去法误差分析2、G-S 迭代法误差分析3、SOR 迭代法误差分析G-S 迭代法与SOR 迭代法的误差比较 七、心得体会正文:一、问题背景介绍。
理论的分析表明,求解病态的线性方程组是困难的。
实际情况是否如此,会出现怎样的现象呢?二、建立正确的数学模型。
考虑方程组Hx b =的求解,其中系数矩阵H 为Hilbert 矩阵,,,1(), , ,1,2,,1i j n n i j H h h i j n i j ⨯===+-这是一个著名的病态问题。
通过首先给定解(为方便计算,笔者取x 的各个分量等于1),再计算出右端,b Hx =这样Hx b =的解就明确了,再用Gauss 消去法、J 迭代法、GS 迭代法和SOR 迭代法四种方法分别求解,Hx b =将求解结果与给定解比较,而后求出上述四种方法的误差,得出哪种方法比较好。
线性代数是高等学校理工科和经济类学科相关专业的一门重要基础课,它不仅是其他数学课程的基础,也是物理、力学、电路等专业课程的基础。
作为处理离散问题工具的线性代数,也是从事科学研究和工程设计的科研人员必备的数学工具之一。
实验一 生物遗传模型1.工程背景设一农业研究所植物园中某植物的基因型为AA 、Aa 和aa 。
常染色体遗传的规律是:后代是从每个亲体的基因对中个继承一个基因,形成自己的基因对。
如果考虑的遗传特征是由两个基因A 、a 控制的,那末就有三种基因对,记为AA 、Aa 和aa 。
研究所计划采用Aa(AA)型的植物与每一种基因型植物相结合的方案培育植物后代。
问经过若干年后,这种植物的任意一代的三种基因型分布如何2.问题分析分析双亲体结合形成后代的基因型概率,如表6-4所示。
表6-4基因型概率矩阵 后代 基因对 父体—母体的基因对AA —AA AA —Aa AA —aa Aa —Aa Aa —aa aa —aa AA 1 1/2 0 1/4 0 0 Aa 0 1/2 1 1/2 1/2 0 aa1/41/213.模型建立与求解设,,n n n a b c 分别表示第n 代植物中基因型AA 、Aa 、aa 型的植物占植物总数的百分率。
则第n 代植物的基因型分布为()n n n n a x b c ⎛⎫ ⎪= ⎪ ⎪⎝⎭,0(0)00a x b c ⎛⎫ ⎪= ⎪ ⎪⎝⎭表示植物型的初始分布。
依据上述基因型概率矩阵,有1112n n n a a b --=+,1112n n n b b c --=+,0n c =,1n n n a b c ++=,表示为矩阵形式11111/2001/21000n n n n n n a a b b c c ---⎛⎫⎛⎫⎛⎫⎪ ⎪⎪= ⎪ ⎪⎪ ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭记11/2001/21000M ⎛⎫ ⎪= ⎪ ⎪⎝⎭,则()(1)2(2)3(3)(0)n n n n n x MxM x M x M x ---=====。
线性模型的标准形式
线性模型是统计学中常见的一种模型,它在各个领域都有着广泛的应用。
线性模型的标准形式是指模型的数学表达式,通常包括自变量、因变量和参数。
在本文中,我们将详细介绍线性模型的标准形式,包括线性回归模型和线性分类模型。
首先,我们来介绍线性回归模型的标准形式。
线性回归模型用于建立自变量和因变量之间的线性关系。
其标准形式可以表示为:
Y = β0 + β1X1 + β2X2 + ... + βpXp + ε。
其中,Y表示因变量,X1、X2、...、Xp表示自变量,β0、β1、β2、...、βp表示模型的参数,ε表示误差。
在这个模型中,我们希望通过最小化误差来估计参数的取值,从而建立自变量和因变量之间的线性关系。
接下来,我们将介绍线性分类模型的标准形式。
线性分类模型用于将样本划分到不同的类别中,其标准形式可以表示为:
Y = β0 + β1X1 + β2X2 + ... + βpXp。
与线性回归模型不同的是,线性分类模型不需要考虑误差项,因为它的目标是将样本正确地分类,而不是建立自变量和因变量之间的精确关系。
在实际应用中,线性模型的标准形式可以通过最小二乘法、梯度下降法等方法来估计参数的取值,从而得到模型的最优解。
此外,线性模型还可以通过添加交互项、多项式项等方式进行扩展,以适应更复杂的数据模式。
总之,线性模型的标准形式是建立模型的基础,它描述了自变量和因变量之间的线性关系。
通过对模型参数的估计,我们可以得到一个能够描述数据特征的线性模型,从而进行预测、分类等任务。
希望本文对您理解线性模型的标准形式有所帮助。
`线性代数是高等学校理工科和经济类学科有关专业的一门重要基础课, 它不单是其余数学课程的基础,也是物理、力学、电路等专业课程的基础。
作为办理失散问题工具的线性代数,也是从事科学研究和工程设计的科研人员必备的数学工具之一。
实验一 生物遗传模型1.工程背景AA 、Aa 和aa 。
常染色体遗传的规律是:设一农业研究所植物园中某植物的基因型为后辈是从每个亲体的基因对中个继承一个基因,形成自己的基因对。
假如考虑的遗传特点是由两个基因A 、a 控制的,那么就有三种基因对, 记为AA 、Aa 和aa 。
研究所计划采纳Aa(AA)型的植物与每一种基因型植物相联合的方案培养植物后辈。
问经过若干年后,这栽种物的任意一代的三种基因型散布怎样 ?问题剖析剖析双亲体联合形成后辈的基因型概率,如表6-4所示。
表6-4基因型概率矩阵 后辈 父体—母体的基因对基因对AA —AA AA —Aa AA —aa Aa —Aa Aa —aaaa —aaAA 1 1/2 0 1/4 0 0 Aa 0 1/2 1 1/2 1/2 0 aa 0 0 0 1/41/21模型成立与求解设a n ,b n ,c n 分别表示第n 代植物中基因型AA 、Aa 、aa 型的植物占植物总数的百分率。
a na 0则第n 代植物的基因型散布为x(n)b n ,x(0)b 0 表示植物型的初始散布。
依照上述基c nc 0因型概率矩阵,有a n an1 1b n1,b n1b n1c n1,c n0,a n b n c n1,表示为22矩阵形式a n 1 1/2 0 an1b n0 1/2 1 bn1c n0 00 cn11 1/2 0记M0 1/21 ,则x (n) Mx(n1)M2x (n2)M3x (n3)LM n x (0)。
0 0于是问题归纳为怎样计算 M n,可将M 对角化。
易于计算 M 的特点值为1、1/2、0,其相应的特点向量为(1,0,0)T,(0,1,0)T,(1,2,1)T。
数学建模模型常用的四大模型及对应算法原理总结四大模型对应算法原理及案例使用教程:一、优化模型线性规划线性回归是利用数理统计中回归分析,来确定两种或两种以上变量间相互依赖的定量关系的一种统计分析方法,在线性回归分析中,只包括一个自变量和一个因变量,且二者的关系可用一条直线近似表示,这种回归分析称为一元线性回归分析。
如果回归分析中包括两个或两个以上的自变量,且因变量和自变量之间是线性关系,则称为多元线性回归分析。
案例实操非线性规划如果目标函数或者约束条件中至少有一个是非线性函数时的最优化问题叫非线性规划问题,是求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。
建立非线性规划模型首先要选定适当的目标变量和决策变量,并建立起目标变量与决策变量之间的函数关系,即目标函数。
然后将各种限制条件加以抽象,得出决策变量应满足的一些等式或不等式,即约束条件。
整数规划整数规划分为两类:一类为纯整数规划,记为PIP,它要求问题中的全部变量都取整数;另一类是混合整数规划,记之为MIP,它的某些变量只能取整数,而其他变量则为连续变量。
整数规划的特殊情况是0-1规划,其变量只取0或者1。
多目标规划求解多目标规划的方法大体上有以下几种:一种是化多为少的方法,即把多目标化为比较容易求解的单目标,如主要目标法、线性加权法、理想点法等;另一种叫分层序列法,即把目标按其重要性给出一个序列,每次都在前一目标最优解集内求下一个目标最优解,直到求出共同的最优解。
目标规划目标规划是一种用来进行含有单目标和多目标的决策分析的数学规划方法,是线性规划的特殊类型。
目标规划的一般模型如下:设xj是目标规划的决策变量,共有m个约束条件是刚性约束,可能是等式约束,也可能是不等式约束。
设有l个柔性目标约束条件,其目标规划约束的偏差为d+, d-。
设有q个优先级别,分别为P1, P2, …, Pq。
在同一个优先级Pk中,有不同的权重,分别记为[插图], [插图](j=1,2, …, l)。
数学建模_四大模型总结四类基本模型1 优化模型1.1 数学规划模型线性规划、整数线性规划、非线性规划、多目标规划、动态规划。
1.2 微分方程组模型阻滞增长模型、SARS 传播模型。
1.3 图论与网络优化问题最短路径问题、网络最大流问题、最小费用最大流问题、最小生成树问题(MST)、旅行商问题(TSP)、图的着色问题。
1.4 概率模型决策模型、随机存储模型、随机人口模型、报童问题、Markov 链模型。
1.5 组合优化经典问题● 多维背包问题(MKP)背包问题:n 个物品,对物品i ,体积为i w ,背包容量为W 。
如何将尽可能多的物品装入背包。
多维背包问题:n 个物品,对物品i ,价值为i p ,体积为i w ,背包容量为W 。
如何选取物品装入背包,是背包中物品的总价值最大。
多维背包问题在实际中的应用有:资源分配、货物装载和存储分配等问题。
该问题属于NP 难问题。
● 二维指派问题(QAP)工作指派问题:n 个工作可以由n 个工人分别完成。
工人i 完成工作j 的时间为ij d 。
如何安排使总工作时间最小。
二维指派问题(常以机器布局问题为例):n 台机器要布置在n 个地方,机器i 与k 之间的物流量为ik f ,位置j 与l 之间的距离为jl d ,如何布置使费用最小。
二维指派问题在实际中的应用有:校园建筑物的布局、医院科室的安排、成组技术中加工中心的组成问题等。
● 旅行商问题(TSP)旅行商问题:有n 个城市,城市i 与j 之间的距离为ij d ,找一条经过n 个城市的巡回(每个城市经过且只经过一次,最后回到出发点),使得总路程最小。
● 车辆路径问题(VRP)车辆路径问题(也称车辆计划):已知n 个客户的位置坐标和货物需求,在可供使用车辆数量及运载能力条件的约束下,每辆车都从起点出发,完成若干客户点的运送任务后再回到起点,要求以最少的车辆数、最小的车辆总行程完成货物的派送任务。
TSP 问题是VRP 问题的特例。
线性代数是高等学校理工科和经济类学科相关专业的一门重要基础课,它不仅是其他数学课程的基础,也是物理、力学、电路等专业课程的基础。
作为处理离散问题工具的线性代数,也是从事科学研究和工程设计的科研人员必备的数学工具之一。
实验一 生物遗传模型1.工程背景设一农业研究所植物园中某植物的基因型为AA 、Aa 和aa 。
常染色体遗传的规律是:后代是从每个亲体的基因对中个继承一个基因,形成自己的基因对。
如果考虑的遗传特征是由两个基因A 、a 控制的,那末就有三种基因对,记为AA 、Aa 和aa 。
研究所计划采用Aa(AA)型的植物与每一种基因型植物相结合的方案培育植物后代。
问经过若干年后,这种植物的任意一代的三种基因型分布如何?2.问题分析分析双亲体结合形成后代的基因型概率,如表6-4所示。
表6-4基因型概率矩阵 后代 基因对 父体—母体的基因对AA —AAAA —Aa AA —aa Aa —Aa Aa —aa aa —aa AA 1 1/2 0 1/4 0 0 Aa 0 1/2 1 1/2 1/2 0 aa1/41/213.模型建立与求解设,,n n n a b c 分别表示第n 代植物中基因型AA 、Aa 、aa 型的植物占植物总数的百分率。
则第n 代植物的基因型分布为()n n n n a x b c ⎛⎫ ⎪= ⎪ ⎪⎝⎭,0(0)00a x b c ⎛⎫ ⎪= ⎪ ⎪⎝⎭表示植物型的初始分布。
依据上述基因型概率矩阵,有1112n n n a a b --=+,1112n n n b b c --=+,0n c =,1n n n a b c ++=,表示为矩阵形式记11/2001/21000M ⎛⎫⎪= ⎪ ⎪⎝⎭,则()(1)2(2)3(3)(0)n n n n n x MxM x M x M x ---=====。
于是问题归结为如何计算nM ,可将M 对角化。
易于计算M 的特征值为1、1/2、0,其相应的特征向量为(1,0,0)T ,(0,1,0)T -,(1,2,1)T-。
令101012001P ⎛⎫ ⎪=-- ⎪ ⎪⎝⎭,则111/2001/21000M P P -⎛⎫ ⎪= ⎪ ⎪⎝⎭。
于是()(0)1(0)11/2001/21000nn n xM x P P x -⎛⎫⎪== ⎪ ⎪⎝⎭1001001(1/2)(1/2)(1/2)(1/2)0n n n n b c b c --⎛⎫-- ⎪=+ ⎪ ⎪⎝⎭。
当n →∞,1,0n n a b →→,因此,可以认为经过若干年后,培育出的植物基本上呈现AA 型。
实验二 员工培训问题 1.工程背景某试验性生产线每年一月份进行熟练工与非熟练工的人数统计,然后将1/6熟练工支援其他生产部门,其缺额由招收新的非熟练工补齐。
新、老非熟练工经培训及实践至年终考核有2/5成为熟练工。
若记第n 年一月份统计的熟练工与非熟练工所占比例分别为n n x y ⎛⎫⎪⎝⎭。
2.问题问题1:第n+1年熟练工与非熟练工所占比例11n n x y ++⎛⎫⎪⎝⎭与第n 年熟练工与非熟练工所占比例n n x y ⎛⎫⎪⎝⎭的关系。
问题2:若第1年熟练工与非熟练工所占比例为111212x y ⎛⎫⎪⎛⎫= ⎪ ⎪ ⎪⎝⎭ ⎪⎝⎭,求11n n x y ++⎛⎫ ⎪⎝⎭。
3.模型建立与求解 依据题意,有1521()656n n n n x x x y +=++,131()56n n n y x y +=+。
整理化简得119210513105n n nn n n x x y y x y ++⎧=+⎪⎪⎨⎪=+⎪⎩,即119210513105n n n n x x y y ++⎛⎫ ⎪⎛⎫⎛⎫=⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭ ⎪⎝⎭,记9210513105A ⎛⎫⎪= ⎪ ⎪ ⎪⎝⎭,亦有11n n n n x x A y y ++⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭。
由问题1结果,有112111212n n n n n n n x x x A A A y y y +-+-⎛⎫ ⎪⎛⎫⎛⎫⎛⎫==== ⎪⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎪⎝⎭。
问题归结为求n A ,可将A 对角化。
易于计算1、1/2是矩阵A 的两个特征值,且相应的特征向量为()()4,1,1,1TT-。
记4111P -⎛⎫= ⎪⎝⎭,则1921010511302105P P -⎛⎫⎛⎫⎪ ⎪=⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭11104()44()411111221111411550()1()14()222n n n n n n A ⎛⎫+-⎛⎫ ⎪-⎛⎫⎛⎫ ⎪==⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭ ⎪-+ ⎪⎝⎭⎝⎭。
因此111183()122111023()22n n n n n x A y ++⎛⎫⎛⎫- ⎪ ⎪⎛⎫== ⎪⎪ ⎪ ⎪ ⎪⎝⎭+ ⎪ ⎪⎝⎭⎝⎭。
实验三 多金属分选流程计算1. 工程背景设,j γγ—原矿产率及第j 种产品产率,%,100%γ=;i α—原矿中第i 种金属品位,%;ij β—第j 种产品中第i 种金属品位,%;ijβε—第j 种产品中第i 种金属的理论回收率,%;按照金属平衡和产率平衡进行计算。
为了计算方便,尾矿视为产品。
金属平衡, 产品平衡,其中,尾矿产率及金属品位为,n i in θγγθβ== 解次多元线性方程组求出产品产率。
各产品任一金属回收率1100%ijj ijnjijj βγβγβε=⨯=∑。
2. 问题某铅锌矿选矿厂生产的产品为铅、锌、硫精矿和尾矿,已化验知各产品的金属品位(见下表),试计算各产品产率和回收率。
表6-5各产品的化验品位原矿3.143.6315.41铅 71.04 3.71 15.70 锌 1.20 51.50 30.80 硫 0.38 0.35 42.38 尾矿0.340.101.403. 模型建立与求解设铅、锌、硫和尾矿的产率为123,,x x x 和4x ,按照金属平衡与产率平衡,可建立以下线性方程组:MATLAB 源代码:A=[71.04 1.20 0.38 0.34;3.71 51.50 0.35 0.10;15.70 30.80 42.38 1.40;1 1 1 1] %创建系数矩阵b=[314 363 1541 100]’; %常数列矩阵 x=A\b %利用x=inv(A)*b x =3.8659 6.4590 28.2046 61.4706又 x0=repmat(x,[1,4]); %创建多维数组B0=repmat(b,[1 4])’;s=x0.*A ’./B0 %计算各产品的理论回收率,最后一列为产率 s=87.4623 3.9511 3.9386 3.8659 2.4684 91.6361 12.9096 6.4590 3.4133 2.7194 77.5671 28.2046 6.6560 1.6934 5.5846 61.4706 将计算结果填入下表实验四 交通流量模型1. 问题图6-8.200300 500 400 500 x9 x10600700图6-8假设:(1)全部流入网络的流量等于全部流出网络的流量.(2)全部流入一个节点的流量等于全部流出此节点的流量. 2. 模型的建立与求解由假设可知,所给问题满足如下线性方程组:3. Matlab 程序实现A=[0,1,-1,1,0,0,0,0,0,0;0,0,0,1,1,0,0,0,0,0;0,0,0,0,0,-1,1,0,0,0;1,1,0,0,0,0,0,0,0,0;1,0,0,0,1,0,0,0,0,0;0,0,0,0,0,0,1,1,0,0;0,0,0,0,0,0,0,0,1,0;0,0,0,0,0,0,0,0,-1,1;0,0,0,0,0,0,0,0,0,1;0,0,1,0,0,1,0,1,0,0] % 矩阵A b=[300;500,200;800;800;100;400;200;600,100]B=[A,b] % 增广矩阵B Rank(A) % 计算矩阵A 的秩Rank(B) % 计算增广矩阵B 的秩,若秩相等,则有解 rref(B) % 将增广矩阵B 化为最简型4.结果分析 增广矩阵系数矩阵的秩Rank(A)=8增广矩阵的秩Rank(B)=8<10,说明该非齐次线性方程组有无穷多个解. 增广矩阵的最简型为: 其对应的齐次同解方程组为:以85,x x 做为自由变量,将最简形方程转化为⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧==+-=+-=+-==+=+-=600400100800500200080010987865435251x x x x x x x x x x x x x 求得其通解为⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛+⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--+⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛6004000100080005002000800001110000000000110112110987654321C C x x x x x x x x x x。