数学建模实训报告2.pdfVIP

单个实验报告·正文 【实验报告正文】 一、实验目的和要求 实验目的：掌握使用MATLAB 软件实现优化类模型的方法，为进一步学习数学建模打下基础。 实验要求：请分别给出每道题的上机实验程序、实验结果，并上传到 。 二、实验内容 （1）生猪的销售时机模型求解 （p63 ）：用图解法绘制目标函数Q(t) = ( 8 – gt )( 80 + rt ) – 4t – 640 的图形（0≤t ≤20 ）。其中， g=0.1 ， r=2 。从图形上观察出曲线Q(t) 的最大值。 （2 ）设某种规格的钢筋原材料每根长10m，求解如下优化问题： 现需要该种钢筋长度为4m 的28 根，长度为1.8m 的33 根，问至少需要购买原材料几根？ 如何切割？ 装 请对上述问题进行求解，要求建立数学模型，写出使用 Lingo 软件的求解程序，并给出 订 结果分析。 （3 ）附加题：用粒子群算法求解y=1-cos(3*x)*exp(-x) 的在[0,4]最大值。 线 三、实验原理或实验方法 使用MATLAB 进行编程。 四、实验步骤（含实验数据记录处理）或操作设计过程记录 （1）程序如下： t=[0:1:20]; Q=( 8 - 0.1.*t ).*( 80 + 2.*t )-4.*t-640; plot (t,Q, -b) grid t=10; Q=( 8 - 0.1.*t ).*( 80 + 2.*t )-4.*t-640; disp (Q); 运行结果如下： 经过观察发现在t=10 处取到最大值 验证过程 Q=-0.2t*t+4*t 单个实验报告·正文 Q’=-0.4t+4 说明Q 先递增后递减 则Q ‘=0 时取到最大值 此时t=10,Q=20 （2 ）程序如下： model : min=x1+X2+X3; 2*x1+x2=28; x1+3*x2+5*x3=33; @gin(x1);@gin(x2);@gin(x3); 单个实验报告·正文 运行结果如下： （3 ）程序如下： c1=1.4965; c2=1.4965; w=0.7298; num=100; n=5; for i=1:n x(i)=4*rand; v(i)=rand; end for i=1:n p(i)=x(i); end pg=x(1); for i=2:n if 1-cos(3*x(i))*exp(-x(i))1-cos(3*pg)*exp(-pg) pg=x(i); 单个实验报告·正文 end end for t=1:num for i=1:n v(i)=w*v(i)+c1*rand*(p(i)-x(i))+c2*rand*(pg-x(i)); x(i)=x(i)+v(i); if x(i)4||x(i)0 x(i)=4*rand; end end for i=1:N if 1-cos(3*x(i))*exp(-x(i))1-cos(3*p(i))*exp(-p(i)) p(i)=x(i); end end for i=1:n if 1-cos(3*x(i))*exp(-x(i))1-cos(3*pg)*exp(-pg) pg=x