matlab方程求解.docVIP

开课学院、实验室： 数统学院 实验时间 ： 2012年 3月 3日 课程 名称 数学建模与数学实验 实验项目 名 称 飞机如何定价 ——方程求解 实验项目类型 验证 演示 综合 设计 其他 指导 教师 成 绩 实验目的 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 二、实验过程 运行Matlab，利用命令solve()和fsolve()对方程组求解。 运算结果： x1 = 0.56714329040978387299996866221036 x2 = 0.56714329040978387299996866221036 对于方程组（2），键入Matlab命令： [x1,x2,x3]=solve(x1^2-5*x2^2+7*x3^2=-12,3*x1*x2+x1*x3-11*x1=0,2*x2*x3+40*x1=0,x1,x2,x3)；x1=vpa(x1,4)，x2=vpa(x2,4)，x3=vpa(x3,4) 运算结果： x1 = 1.0 0 0 0 0 -0.3145 - 32.7*i - 387.0 32.7*i - 387.0 x2 = 5.0 1.549 -1.549 0 0 2.958 50.81*i - 0.3123 - 50.81*i - 0.3123 x3 = -4.0 0 0 1.309*i -1.309*i 2.126 11.94 - 152.4*i 152.4*i + 11.94 应用实验 一、实验内容 汽车厂生产三种类型的汽车，已知各类型每辆车对钢材、劳动时间的需求，利润及工厂每月的现有量如表1所示，试制订月生产计划，使工厂的利润最大，建立模型并求解。 表1 汽车厂相关参数 小型 中型 大型 现有量 钢材（吨） 1.5 3 5 600 劳动时间（小时） 280 250 400 60000 利润（万元） 2 3 4 二、问题分析 根据题意建立求利润的函数关系式。本题中目标函数便是利润，利用所给的变量的条件，进行线性规划，得出最优解。 三、数学模型的建立与求解(一般应包括模型、求解步骤或思路，程序放在后面的附录中) （1）定义变量：x——小型车月产量；y——中型车月产量；z——大型车月产量； P——总利润。 （2）建立模型：p=2x+3y+4z； ……1) 1.5x+3y+5z=600； ……2） 280x+250y+400z=60000； ……3） （3）Matlab命令： c=[-2 -3 -4]; A=[1.5 3 5;280 250 400]; b=[600;60000]; Aeq=[];beq=[]; vlb=[0;0;0];vub=[]; [x,fval,EXITFLAG]=linprog(c,A,b,Aeq,beq,vlb,vub) （4）运行程序，得解。 四、实验结果及分析 运行后输出的结果： x = 64.5161 167.7419 0.0000 fval = -632.2581 EXITFLAG = 1 该结果显示：造小型车64.5161辆，中型车167.7419辆，不造大型车，获利最高，为：632.2581万。 但实际生活中，车辆数均为整数。分别取x=64或65与y=167或168组合，并结合约束条件分析得： 应加工小型车64辆，中型车168辆，最大利益为632万元每月。 附录 c=[-2 -3 -4]; A=[1.5 3 5;280