农作物杀虫剂的合理使用.pptVIP

农药的合理使用 规划模型的应用 问题的描述 农作物生长期间遭到虫害，农作物减产； 使用农药可减少虫害，但会给农作物的安全带来隐患； 使用农药必须考虑到农药残留量的限制； 问题的实例 设一地块为5亩，种植了某农作物，其生长期为100天，每次农药的使用量每亩不得超过1千克。假定在农作物生长初期，已发现农作物的害虫密度为每平方米1.5只，平均每只害虫能吞食农作物0.001千克，该地块害虫的横向增加率约为200只/天（即害虫由其它区域向该地块迁入与迁出的数量之差），害虫每天的增长率为2%，每次喷洒农药需固定设备费50元，人工费10元，农药每千克10元，农作物收购价每千克1元，问应如何使用农药，能降低损害，又使花费最小？ 提出假设 农作物的数量足以满足害虫的消耗，害虫的自然增长率符合Mathus模型； 单位时间内害虫的固有增长率为r； 由其它地块向该地块的净迁入率为?； 害虫的自然死亡率可忽略不计； t时刻害虫的数量； t时刻使用了剂量为u的农药后，害虫的剩余量。即 时刻害虫的残存量。 提出假设 农作物的被损毁的速率与害虫数成正比； 每个害虫在单位时间内对农作物的吞食量是一个常数； 喷洒农药的每次费用是关于农药剂量的一个线性函数； 单位面积内使用每单位农药，其能杀死害虫的比例为一个常数?。 引入变量符号 农作物的生长期长度; 期初该地块的害虫数量; 单位时间内害虫的自然增长率; 每千克农药对害虫的杀死率; 横向增长率； 每只害虫在单位时间内对农作物的破坏量； 农作物的单价； 生长期内喷洒农药的固定费用； 农药的单价； 每次喷洒农药的人工费用； 单位面积内每次农药使用量的最大值。 模型的建立 害虫数量的变化规律，可由微分方程描述： 初始条件: 使用农药后,农作物上残留的害虫数量为: 可得微分方程的解： 模型的建立 考虑农药的使用对害虫数量的影响： 单位数量的农药对害虫的杀死率为一个常数?。 解此方程，得： 使用了农药后，害虫的数量为： 第k次农药的使用量。 以上确定了害虫数量的变化规律。 模型的建立 被害虫破坏后农作物的损失费用： 在整个生长期内被害虫破坏的农作物的损失费用： 喷洒农药的费用： 农药使用的模型 整个问题可归结为如下的规划问题： 是一个多元动态规划模型，求解困难。 模型的简化 喷洒农药策略的简化； 每次喷洒农药使用相同的剂量； 喷洒农药的时间间隔相同； 每次喷洒农药后，害虫的数量都降为z。 初始时刻，考虑使用农药数量： 可确定 模型的简化 经过一个周期（h时间），则害虫的数量为： 此时使用农药数量u，使害虫数量降到： 农作物的损失费用： 模型的简化 整个生长期内的总费用： 约束条件为： 简化后的模型 将某些约束条件代入，整个规划问题为： 代入求解。 模型计算 应用数学软件求解，可得最优解。 模型计算 function y = Nongyaoobj( x ) u = x(1); h = x(2); arfa = 0.5; beta = 100; x0 = 5000; r = 0.02; b = 0.001; c0=50; c1 = 10; c2 = 10; p = 1; T = 100; y = c0+c1/arfa*log(r*x0*(exp(arfa*u)-exp(r*h))/beta/(exp(r*h)-1))... +c1*(T/h-1)*u+c2*T/h+p*b*T*beta/h*(exp(arfa*u)-1)*(exp(r*h)-1)/(r^2*... (exp(arfa*u)-exp(r*h)))-p*b*T*beta/r; function [ C,Ceq ] = Nongyaocon( x ) u = x(1); h = x(2); arfa = 0.5; beta = 100; x0 = 5000; r = 0.02; C = [ beta*(exp(r*h)-1)/r/(exp(arfa*u)-exp(r*h))-x0; r*h-arfa*u]; Ceq = []; xt=[10,20];lb=[0 0];ub=[5,100]; [x,f]=fmincon(@Nongyaoobj,xt,[],[],[],[],lb,ub,@Nongyaocon)