养猪策略.docx

养猪策略摘要本文根据某个养猪场的生猪从第一天到第24天间的生猪的体重数据建立数学模型对该养猪场未来猪的体重作出相应的预测，首先我们利用plot将题目所给的生猪的体重与时间的关系利用图形将其表现出来；接着通过综合各方面因素，我们决定建立Logistic模型来解决问题。通过matlab数据拟合并求解出第25天的生猪的体重预测为*，通过数据检验和误差分析，发现相对误差在合理的范围内。关键字养猪策略，阻滞增长模型（Logistic模型），合理范围 一、问题重述3. 养猪策略：一个养猪场每天投入5元的饲料喂养每头猪，其中某头猪的体重增长如下表： 时间(天)12345678重量(千克)10.111.41314.816.919.121.524.1时间(周)910111213141516重量(天)26.93033.336.740.143.947.651.3时间(天)1718192021222324重量(千克)55.158.762.365.769.172.175.177.8（1）试建立数学模型描述猪的体重与饲养时间的关系，并预测猪在25天时的体重。（2）已知在第24天时的市场价格为每千克10元，但是预测每天会降低 0.1元，问何时把猪卖出去，才能获得最大利润？我们先从图像上来观察该养猪场时间和体重的大致关系，利用matlab的绘图功能把上表绘制如下：图一从上图我们可以看出猪的体重从第1天到第24天随时间的推移而呈现上升趋势，但其增长速率渐渐的有所下降。其原因是-由于受到自然环境、养猪场自身环境等因素的影响，通过各方面的分析，我们得出猪的体重的增长率不是一个常数，所以我们建立了阻滞增长模型（Logistic模型）。阻滞增长模型（Logistic模型）二、模型假设、模型变量和函数定义1.猪的体重增长率随着时间的增加而线性下降，即是r(g)=r-sg 其中2.Gm为一般猪的体重所能达到的最大值，显然r(Gm)=03. 为固有人口增长率，即：三、符号说明1 人口固有增长率2Gm 为一般猪的体重所能达到的最大值3 人口增长系数4g 当前猪的体重5r(g) 当前增长率四、模型分析与模型建立猪的体重的增长受到自然资源、养猪场环境条件、等因素影响，当猪的体重达到一定时候，人口增长将出现阻滞作用。阻滞作用体现在对猪体重增长率的影响上。依照上面的定义和假设，有： r(g)=r-sg r(Gm)=0 由上面的关系式推出： r(g)=r(1-g/Gm) 将r(g)的表达式代入指数增长模型中的微分方程中得：dx/dg=r(1-g/Gm)gg(0)=g0 把g(1)=10.1代入上述方程二式解之得：g(t)=-(Gm*(tan(-(Gm*(r + 2*atan(101*i)/(5*Gm) - i)*i)/Gm - (r*t)/Gm)*i)/2) + i)*i)/2五、模型求解及检验通过拟合求解，我们得到Gm=100.58kg，r=-0.1488 。通过这个模型可以预测x(25)=80.34六、误差检验为了检验这个模型的预测效果，我们求出预测值与真实值的相对误差。拟合后真实值与预测值以及相对误差表如下：(表一)：时间1234567891011真实值10111.413014.816.919.121.524.126.93033.3预测值10.1011.5313.1414.9416.9319.1321.5424.1727.0130.0533.27相对误差00.130.140.140.030.030.040.070.110.050.03时间12131415161718192021真实值36.740.143.947.651.355.158.762.365.769.1预测值36.6640.2043.8447.5451.2855.0158.6962.2765.7369.04相对误差0.040.10.060.060.020.090.010.030.030.06 时间222324真实值72.175.177.8预测值72.1775.1077.83相对误差0.0100.03平均相对误差 0.0546 原有数据点与拟合曲线如下图：图二由上表可以看出预测值与真实值的相对误差较小。七、结果分析我们先通过matlab利用plot将题目所给的数据制作成图形，从而对生猪的体重与时间关系有了较为直观的了解，通过综合各方面的因素（如自然环境因素等）分析，我们决定用阻滞增长模型（Logistic模型）对该养猪场的猪的重量进行预测，发现整体的相对误差比较小情况，而且后面几年的数据也比较吻合。 八、模型评价模型跟根据生猪体重增长的规律，做出了合理的假设，建立了符合实际的阻滞增长模型（Logistic模型），该模型无论是对已知数据的拟合还是对未知数据的预测都相当合理。这样的模型符合事物发展规律，预测效果不错。问题二：简单优化模型：一、符号说明1p 出售单价（元/千克）2C 投入的资金（元）3Q 纯利润4R 出售的收入二、模型分析与模型建立投入资金可使生猪的体重随时间的增长，但销售单价随时间减少，应该存在一个最佳出售时机，使得获得利润最大，根据分析建立以下模型： p=10-0.1*(t-24)g=100.5807./(1+(100.5807/10.1-1)*exp(-0.1488*(t-1)R=p.*gC=5tQ= R-c*t得到目标函数的（纯利润）为：Q= 10-0.1*(t-24) 100.5807./(1+(100.5807/10.1-1)*exp(-0.1488*(t-1)- 5t三、模型求解及检验通过MATLAB得到如下图形：由图形分析可得当t=30天 时所得利润最大，最大利润为Q= 694.3705元 所以最适合出售的时间是第30天。 四、 结果分析敏感性分析：由于模型假设中的参数（生猪每天的价格降低g）是估计和预测的，所以应该研究它有所变化时对模型结果的影响。由模型一，我们可清楚的看到g与t的关系，接下来我们用相对改变量衡量结果对参数的敏感程度。T对g的敏感度记作s(t,g);定义为：s(t,g)=dt/dg*(g/t)=9.9/(g-50)当g=0.1.时s(t,g)=-0.2即生