matlab在数学中的应用

《matlab在数学中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第六讲MATLAB在数学建模中的应用一、在线性规划中的应用1、数学原理：线性规划是处理线性目标函数和线性约束的一种较为成熟的方法，目前已经广泛应用于军事、经济、工业、教育、商业和社会科学等许多方面。线性规划的标准形式是线性规划的标准形式要求目标函数最小化，约束条件取等式，变量非负，不符合这几个条件的线性规划要首先转化为标准形式。线性规划的求解方法主要是单纯形法（simpleMethod),此法由Dantzig于1947年提出，以后经过多次改进，2、线性规划的MATLAB求解：linprog函数数学模型：其中：f，x，b，beq，lb，ub为向量，A，Aeq为矩阵。使用形式：x=

2、linprog(f,A,b)x=linprog(f,A,b,Aeq,beq,lb,ub)x=linprog(f,A,b,Aeq,beq,lb,ub,x0)x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options)[x,fval]=linprog(…)[x,fval,exitflag]=linprog(…)[x,fval,exitflag,output]=linprog(…)[x,fval,exitflag,output,lambda]=linprog(…)注：以上几种形式在使用时根据具体的模型适当选用。每一种形式都有特定的涵义。可查有关书籍3、实例例1.

3、投资问题：某单位有一批资金用于4个工程项目的投资，用于各个工程项目得到的净收益（投入资金的百分比）如表所示：工程项目ABCD收益（%）1510812由于某种原因用项目A的总投资不大于其他各各项目的和，用用项目B和C的投资要大于项目D的投资。试确定使该单位收益最大的投资分配方案。分析问题建立模型：用x1,x2,x3,x4分别代表用于项目A、B、C、D的投资百分数，由于各项目的投资百分数之和为100%，所以x1+x2+x3+x4=1根据题意可建立下面的数学模型：把它转化为标准形式为首先输入下列系数：f=[-0.15;-0.1;-0.08;-0.12]；A=[1-1-1-1;0-1-

4、11]；b=[0;0]；Aeq=[1111]；beq=[1]；lb=zeros(4,1)；调用linprog函数[x,fval,exitflag,output,lambda]=linprog(f,A,b,Aeq,beq,lb)结论：4个项目的投资百分数分别为0.50、0.25、0.00和0.25时可获得最大的收益，最大收益为13%二、其它应用动物的繁殖的规律问题：某农场饲养的某种动物能达到的最大年龄为15岁，将其分为三个年龄组：第一组0~5岁，第二组6~10岁，第三组11~15岁。动物从第二个年龄组开始繁殖后代第二个年龄组的动物在其中年龄段平均繁殖4个后代，第三年龄组的动物在其

5、年龄段平均繁殖3个后代。第一年龄组和第二年龄组的动物能顺利进入下一年龄组的成活率分别为0.5和0.25.假设农场现有的三个年龄组的动物各1000头，根据有关生物学研究结果，对于足够大的时间值k,有其中是莱斯利矩阵L的唯一的正特值请检验这一结果是否正确，如果正确给出适当的k值计算5年后、10年后、15年后各年龄段的动物数量。20年后农场三个年龄段的动物情况会怎样？以五年为一个年龄段，则某一时刻三个年龄段的动物数量可以用一个向量X(k)=[x1(k)x2(K)x3(k)]T为第k个时间段动物数分布向量。问题分析：由题设，在初始时刻0~5岁、6~10岁、11~15岁的三个年龄段动物数

6、量分别为x1(0)=1000，x2(0)=1000，x3(0)=1000如果每五年平均向市场供应动物数c=[sss]T，在20年后农场动物不至于灭绝的前提下，c为多少为好？当k=0,1,2,3时，X(k)表示现在、五年后、十年后、十五年后的动物数分布向量。根据第二年龄段和第三年龄段的繁殖能力，在第k个时间段，第二年龄组动物在其年龄段平均繁殖4个后代，第三年龄组动物在其年龄段平均繁殖3个后代。由此得第一个年龄组在第k+1个时间段的数量如下：x1(k+1)=4x2(k)+3x3(k)同理，根据第一年龄组和第二年龄组的存活率，可得等式X2(k+1)=0.5x1(k)X3(k+1)=0

7、.25x2(k)K=0,1,2,3可得数学模型如下：或写成矩阵形式：X（k+1)=LX(k)其中是莱斯利矩阵由此得X（k+1)=Lk+1X(０)程序和计算X0=[1000;1000;1000]A=[043;0.500;00.250]X1=A*X0X2=A*X1X3=A*X2X4=A*X3为了计算L的特征值，输入下面的命令eig(A)得到特征值为ans=1.5000-1.3090-0.1910这说明只有一个正特征值1.5为了验证运行下面程序；x=[1000;1000;1000];d1=1.5;