线性规划问题求解----数学建模实验报告.docxVIP

084实验报告实验目的：（1)学会用matlab软件解决线性规划问题的最优值求解问题。（2）学会将实际问题归结为线性规划问题用MATLAB软件建立恰当的数学模型来求解。（3）学会用最小二乘法进行数据拟合。（4）学会用MATLAB提供的拟合方法解决实际问题。实验要求：（1）按照正确格式用MATLAB软件解决课本第9页1.1、1.3，第100页5.1、5.3这几个问题，完成实验内容。（2）写出相应的MATLAB程序。（3）给出实验结果。（4）对实验结果进行分析讨论。（5）写出相应的实验报告。实验步骤：（1)、对于习题1.1:a.将该线性规划问题首先化成MATLAB标准型b．用MATLAB软件编写正确求解程序：程序如下:c=[3,-1,-1];a=[4,-1,-2;1,-2,1]; b=[-3;11]aeq=[-2,0,1]; beq=1;[x,y]=linprog(-c,a,b,aeq,beq,zeros(3,1))x,y=-y（2)、对于习题1.3:a.建立适当的线性规划模型：对产品I来说，设以A1，A2完成A工序的产品分别为，件，转入B工序时，以B1，B2，B3完成B工序的产品分别为，，件；对产品II来说，设以A1，A2完成A工序的产品分别为，件，转入B工序时，以B1完成B工序的产品为件；对产品III来说，设以A2完成A工序的产品为件，则以B2完成B工序的产品也为件。由上述条件可得，.由题目所给的数据可建立如下的线性规划模型：Min z(1.250.25)()(20.35)(2.80.5)-(510)(7912)(68) (411)-7s.t.b.运用MATLAB软件编写程序求解：程序如下：c=[0.75,1-(321*7*0.0001),-16*6,(-783*4)/7000,-7/20,-0.5,-321*9*0.0001,1.15,2.3-(321*12*0.0001-(783*11)/7000)];a=[-5,0,0,0,0,-10,0,0,0;0,-7,0,0,0,0,-9,0,-12;0,0,-6,0,0,0,0,-8,0;0,0,0,-4,0,0,0,0,-11;0,0,0,0,-7,0,0,0,0];b=[-6000;-10000;-4000;-7000;-4000];aeq=[1,1,-1,-1,-1,0,0,0,0;0,0,0,0,0,1,1,-1,0];beq=[0;0];[x,y]=linprog(c,a,b,aeq,beq,zeros(3,1))（3)、对于习题5.1:用MATLAB中的三次函数，二次函数，四次函数进行数据拟合，然后与原来结果进行比较。程序如下：x=-5:0.3:5;L=length(x);a=[1,-6,5,-3];y=polyval(a,x);plot(x,y,.-)no=randn(1,L);hold onplot(x,y+no,*);b1=ployfit(x,y+no,3);y1=ployval(b1,x);plot(x,y1,-)b2=ployfit(x,y+no,2);y2=ployval(b2,x);plot(x,y2,-)b3=ployfit(x,y+no,4);y3=ployval(b3,x);plot(x,y3,rp-)（4)、对于习题5.3:用MATLAB中最小二乘法求拟合表中的数据。程序如下：x=[1:8];y=[15.3,20.5,27.4,36.6,49.1,65.6,87.87,117.6];xishu=[ones(8,1),x];%构造系数矩阵cs=xishu\log(y);%线性最小二乘法拟合参数cs(1)=exp(cs(1));%把lna变换成a实验结果与分析:（1)、对于习题1.1:最优解为x = 4.0000 1.000 9.0000对应的最优值为 y = 2.0000（2)、对于习题1.3:x11200，x2230.0493，x30，x4858.6207，x5571.4286，x60，x7500，x8500，x9324.1379， 最优值为z1146.567元。对应整数规划的最优解为x11200，x2230，x30，x4859，x5571，x60，x7500，x8500，x9324， 最优值为z1146.414元。（3)、对于习题5.1:相比较而言，三次拟合和四次拟合结果都较好，二次拟合结果相对较差。（4)、对于习题5.3：最后拟合的函数为