肥料效应方程求解施肥量和作物产量的方法.doc

1、利用肥料效应方程求解施肥量和作物产量的方法 1概述 利用逐步回归这一数学方法，将试验中测得的产量数据与肥料施用量拟合成回归方程，一般叫做肥料效应方程， 或叫肥效反应方程，或叫效应方程，或直接叫做产量曲线。这一方法在当前的土壤肥料研究中得到越来越广泛的应用。 以下简要介绍如何利用肥效回归方程来求出最高产量施肥量、理论最高产量；最佳经济施肥量、最佳经济产量；经济 合理施肥量、经济合理产量等计算关键性施肥量和产量的方法。 将施肥量选在合适的范围内布置田间试验，所得作物产量可以与施肥量模拟出抛物型二次回归方程式。在施肥量 较低时，作物产量随施肥量增加而提高，直到达到最高产量。此后，作物产量随施肥量增加

2、而降低。能够取得最高产 量的施肥量叫做最高产量施肥量。 在商品生产中，需要根据作物产品的产值和肥料投入的价值比较来计算经济效益，这时首先要考虑最后增加的 那一个单位的肥料投入量所能增加的作物产量，这就叫做边际产量。利用高等数学中微分求导方法，很容易得到边际 产量。根据已知回归方程，即可求出肥料用量对作物产量的一阶导数，就是最后一个单位肥料投入量的作物产量增加 量，这就是边际产量方程。 一般情况下，随着施肥量增加，边际产量递减，当达到最高产量时，边际产量为零，就是说这时再增加施肥，作 物产量不再增加，这是重要的一点。根据这一原理，令边际产量方程为 0 时，就可求出最高产量施肥量，将其代回肥 效方

3、程，可求出理论最高产量。 因为肥料是有价格的，所以达到最高产量时，往往在经济上已经不合算，所以要找到最佳经济产量和最佳经济产 量施肥量。作物产量乘以作物产品价格等于作物产值，肥料投入量乘以肥料价格等于肥料投入值。当最后一个单位肥 料投入增加值与作物增产值相等，即两者比例为 1 时，就是最佳经济产量，此时的施肥量即为最佳经济施肥量。在简 单的情况下，可以利用求出的边际产量方程乘以肥料作物价格比得到边际产值方程，当边际产值方程等于 1 时，最后 一个单位肥料投入增加值等于作物产值增加值，这时的产量为最佳经济产量，施肥量为最佳经济施肥量。这时的边际 产值，通俗地说即施肥投入的最后一元钱，作物增产值只

4、能收回一元钱，利润为零。 为了提高利润率，我们可以引入边际利润率的概念，设边际利润率为 R，当边际产量方程等于 1+R 时，就会使最 后一个单位肥料投入值产生一定的利润，其值为 R。 2 介绍几种不同肥料和作物之间的效应方程 2.1 单一肥料单一作物的肥料效应方程 下面分析的是施用一种肥料，种植一种作物的情况，此时的肥料效应二次回归方程为 （式中 y 为作物产量，x 为肥料施用量） 肥料 x 对作物产量 y 的一阶导数，即边际产量方程为： dy/dx = b + 2cx 当边际产量方程为 0 时，解出的 x 值就是最高产量施肥量 xmax b + 2cx max = 0 x max = -b/

5、2c 将 xmax代回肥料效应方程，即可求得理论最高产量 ymax y max = a ? (b2/4c) 在这种简单情况下，令作物产品价格为 Py，肥料价格为 P x，边际产值等于边际投入时就是经济上所允许的最大 产量，用数学公式表示为： dyPy = dxPx 因此边际产值方程为边际产量方程与作物肥料价格比的乘积。当边际产值方程等于 1 时，可求出最佳经济施肥量 xe，则 整理得： b + 2cx e = Px/Py xe =( Px/Py) ? b/2c 将 xe代回肥料效应方程，可求得最佳经济施肥量 ye 以下推出的回归方程系数形式特定产量表达式的书写形式都比较复杂，本文不再给出。 下

6、面讨论经济合理施肥量和经济合理产量。首先，我们回顾一下边际利润率 R。就是最后一个单位肥料投入增加 值所能得到的利润，当 R 大于 0 时，所投入的最后一个单位肥料价值，又叫边际肥料价值能赚钱；当 R 小于 0 时就赔 钱。当 R 为 0.20.5 时；相应的施肥量 xR为 将 xR代回肥效反应方程，可求得经济合理产量 yR。此时肥料投入的最后一元钱仍可赚回 1.21.5 元，既 0.20.5 元的利润，或 20%50%的利润率。 2.2 两种肥料一种作物的肥料效应方程 我们再来分析在一种作物上同时施用两种肥料的情况，此时肥料效应二次方程为： 式中 y 为作物产量，x 1为第一种肥料施用量，x

7、 2为第二种肥料施用量，分别用 x1和 x2对产量 y 求一阶偏导数， 得出它们的边际产量方程： 当两个边际产量方程都等于 0 时，解方程组可得到最高产量施肥量： 将(x 1)max和(x 2)max 代入肥效方程中可得到作物的理论最高产量 ymax如下。 令作物产品价格为 Py，肥料价格分别为 Px1和 Px2。用边际产值方程来求最佳经济施肥量(x 1)e和(x 2)e。 或 b 1 +2c1x1 +c12x2 = Px1/Py b 2 +2c2x2 +c12x1= Px2/Py 求得： 将(x 1)e和(x 2)e代回肥效反应方程，可求出最佳经济产量 ye。 求解经济合理施肥量时先设第一种

8、肥料投入期望得到边际利润率 R1，第二种肥料边际利润率 R2（当然 R1可以等 于 R2）。代入边际产值方程可得到 或 b 1 +2c1x1 +c12x2 = Px1/Py（1 + R 1） b 2 +2c2x2 +c12x1= Px2/Py（1 + R 2） 当 R1和 R2为 0.20.5 时可解得经济合理施肥量 2.3 交互作用、肥料增产趋势比较和肥料主效应方程 这里介绍从肥效反应方程衍生出来的其它几种定性观察的方法。 2.3.1 交互作用 从肥效方程中两种肥料交互作用项 c12x1x2的系数 c12的正负可以看出两种肥料在这种作物大田试验当时当地特定 条件下是表现为协同的正交互作用，还

9、是拮抗的负交互作用。另外从 c12的绝对值大小可看出交互作用的强弱程度。 2.3.2 两种肥料增产趋势比较 用两种肥料的边际产量方程可比较它们的增产趋势大小（同样也可用边际产值方程，因为它们只差系数 Px/Py， 所以此处视同等价）。 比较两个边际产量方程中常数项 b1和 b2（即肥效方程中一次项系数），可以知道两种肥料最初施用时的起始增 产幅度。数值大的起始增产幅度高。 比较两个边际产量方程中一次项系数 c1和 c2（即肥效方程中二次项系数），可以知道两种肥料随施肥量增加衰 减的快慢程度。绝对值大者衰减得快，施用时要注意不要过量。 这种比较类似于几个不同地块上一种作物施用单一肥料，分别得出几

10、个肥效回归方程时，比较方程间对应的各项 系数可知肥料用在哪块地上更划算。 2.3.3 肥料主效应方程 当一种作物上施用两种肥料时，为了了解其中任一种肥料对产量的增产作用，将肥效方程改写后减去截距 a（a 为除肥料外其它因素的产量效应），可以用来考虑在固定另一种肥料用量时该种肥料的主效应方程。其数学表达式为 Ex 1 = (b2x2 +c2x22) + (b1 + c12x2)x1 + c1x12 其中第二种肥料用量 x2是已给定值，因此式中 b2x2+c2x22是常数项，b 1+c12x2是一次项系数。给定不同的 x2，可 得到不同的主效应方程式。其中较有意义的 x2值有(x 2)max、(x

11、 2)e、(x 2)R、(x 2)0等。(x 2)0=0，即第二种肥料用量为 0。同理，可得到 Ex 2 = (b1x1 +c1x12) + (b2 + c12x1)x2 + c2x22 肥料主效应方程还有一种重要作用。当两种肥料中其中任一种肥料投入量受到限制时，可用另一种肥料的主效应 方程来求出另一种肥料的最高产量施肥量、最佳经济施肥量和经济合理施肥量等重要数据。 2.4 单一作物有限制地施用单一肥料的肥料效应方程 由于土壤肥力限制因素不同、水分条件不同、光热条件不同、土壤物理限制因素不同，同一种作物在不同地块和 /或不同季节种植时，虽然施用同一种肥料，也会得到不同的产量反应方程。例如在 A

12、、B 两块地上种植一种作物时的 两条产量曲线分别为 在肥料使用不受数量限制时，可按第二节所介绍的方法分别计算出 A、B 两块地的(x A)max、(x B)max、(x A)e、(x B) e、(x A)R、(x B)R等。 但是如果没有充裕的资金购买足量肥料达到每块地的理想施肥量时，如何将有限的肥料分配在 A、B 两块地上使 用呢？ 先从 A、B 两块地上所能得到的产量总和考虑，即 y = y A+yB。另外还知道两块地上使用的肥料总量 Q = xA+xB。 y = a A + bAxA + cAxA2 + aB + bBxB + cBxB2 代入 xB = Q - xA，并求 xA对 y

13、的一阶导数，得到边际产量方程，令其等于 0，可计算出 A 地块最高产量施肥 量 (b A ? bB ? 2cBQ) + 2(cA + cB)(xA)max = 0 同理，如果代入 xB = Q - xA，得到 B 地块最高产量施肥量 当(x A)max或(x B)max小于 0 时，令其等于 0 后重新计算另一个的最高产量施肥量值（应为 Q）。 将最高产量施肥量(x A)max和(x B)max代回各自的肥效反应方程可得到理论最高产量(y A)max和(y B)max以及 y max 。 在得知产品价格 Py 和肥料价格 Px 后，利用产值总和公式 yP y = (yA+yB)Py求出边际产值

14、公式，令其等于 1，可 计算出最佳经济施肥量。即 (b A ? bB ? 2cBQ) + 2(cA + cB)xA = Px/Py 解出 当(x A)e或(x B)e小于 0 时令其为 0。将求出的最佳经济施肥量代入各自化肥效益反应方程可求出最佳经济产量 (yA)e、(y B)e和 ye。 当边际利润率=0.20.5 时用边际产值方程或可求出经济合理施肥量。 (b A ? bB ? 2cBQ) + 2(cA + cB)xA = Px/Py(1 + R) 解出 当(x A)R或(x B)R小于 0 时，令其为 0，代入各自肥效方程可求得经济合理施肥量 yR、(y A)R、(y B)R。 2.5

15、两种作物施用一种肥料的肥效方程 当两种作物连作或在不同地块同时种植时，不同作物会对施用的肥料有不同的反应。得到不同的肥效反应方程。 比如种 A、B 两种作物，它们的肥效反应方程为 在肥料用量可以同时满足两种作物的最高需求时，最高产量、最佳经济产量和经济合理产量及其相应的施肥量 可根据 A、B 两种作物各自的肥效反应方程分别求出。当肥料总用量受到限制时，则需要按第 2.4 节描述的方法重新 考虑这些问题，通过求 A、B 两种作物总产量的一阶导数，得到边际产量方程，可求得 A 作物和 B 作物的最高产量施 肥量分别为： 当(x A)max或(x B)max小于 0 时，令其为 0 ，将它们分别代入

16、各自的肥效反应方程中可以求得最高产量 ymax、(y A) max、(y B)max。 在得知 A 作物产品价格为 PA、B 作物产量价格为 PB、肥料价格为 PX后利用产值总和公式， yP y = yAPA+ yBPB A、B 两种作物的肥效反应方程代入此式，两种作物所用肥料各用 xA和 xB表示 yP y = (aA + bAxA + cAxA2)PA + (aB + bBxB + cBxB2)PB 两种作物用肥总量为 Q= xA+xB， 将 xB=Q?xA代入上式 yP y = aAPA + (aB + bBQ + cBQ2)PB+ bAPA - (bB + 2cBQ)PB xA+ (c

17、APA + cBPB) xA2 求导得到边际产值公式，令其等于 1，得到最佳经济施肥量。 得 同理可得 B 作物的最佳经济施肥量 当(x A)e或(x B)e小于 0 时，令其为 0。将(x A)e和(x B)e分别代入各自的肥效方程，得到最佳经济产量 ye、(y A)e 和(y B)e。 令边际利润为 0.20.5，代入边际产值方程后得到两种作物的经济合理施肥量。 当(x A)R或(x B)R小于 0 时，令其为 0 。将它们代回各自的肥效方程，得到经济合理产量 yR、(y A)R和(y B)R。 我们可以进一步设想，如果两块地上种植同一种作物，但由于作物产品质量不同或收获季节不同而产生差价

18、时， 也可利用这里介绍的方法求得最佳经济产量和经济合理产量。 2.6 两种作物有限制地施用两种肥料 当两种作物连作或在不同地块同时种植时，不同作物会对施用的两种肥料有不同的反应。以 A、B 两种作物为例， 它们的肥效反应方程为 在两种肥料供应都充足的情况下可分别按计算出的最高产量或最佳经济或经济合理施用量，但如果受到资金 限制，只能得到有限的肥料供应，则要重新安排两种作物的用肥量，假设第一种肥料总量 Q1 = x1A+x1B，式中 x1A和 x1B分别为第一种肥料在 A、B 两种作物上的施用量。第二种肥料总用量 Q2 = x2A+x2B，式中 x2A和 x2B分别为第二种肥 料在 A、B 两种

19、作物上的施用量。先求出两种作物产量总和 y = y A+yB，然后将 x1B = Q1?x1A和 x2B = Q2?x2A代入 y = a A + bA1x1A + bA2x2A + cA1x1A2 + cA2x2A2 + cA12x1Ax2A + aB + bB1 (Q1 - x1A) + bB2 (Q2 ? x2A) + cB1(Q1 - x1A)2 + cB2(Q2 ? x2A)2 + cB12(Q1 - x1A)(Q2 ? x2A) = a A + aB + bB1Q1 + bB2Q2 + cB1Q12 + cB2Q22 + cB12Q1Q2 + (bA1- bB1-2cB1Q1 -

20、cB12Q2) x1A + (bA2- bB2-2cB2Q2 - cB12Q1) x2A + (cA1 + cB1)x1A2 + (cA2 + cB2)x2A2 + (cA12 + cB12)x1A x2A 两种肥料用量 x1A和 x2A分别对 y 求导，得到边际产量方程组，令其等于 0，得 如果求出的数值小于 0，令其等于 0。依此方法，或用 x1B = Q1-x1A； x 2B = Q2-x2A，可进一步求出两种 肥料在 B 种作物上的最高产量施肥量，代回肥效反应方程，可得到理论最高产量(y A)max、(y B)max和 y max。 求两种作物产值总和 yP y = yAPA+yBPB

21、，式中 PA为 A 种作物产品价格，P B为 B 种作物产品价格。另外 PX为小麦 价格，将第一种肥料总用量 Q1 = x1A+x1B，和第二种肥料总用量 Q2 = x2A+x2B代入 yP y然后求出 x1A和 x2A对 y 的一 阶导数。y 的一阶导数等于 1 时，可得到最佳经济施肥量；y 的一阶导数等于 1+R 时，可得到经济合理施肥量等， 如果求出的数值小于 0，令其等于 0。代回肥效反应方程，可求得最佳经济产量和经济合理产量。因以系数表示的 xmax 表示式书写十分复杂，故略去。 2.7 三种肥料一种作物的肥料效应方程 在一种作物上同时施用三种肥料与施用两种肥料的情况相似，此时肥料效

22、应二次方程为： 式中 y 为作物产量，x 1为第一种肥料施用量，x 2为第二种肥料施用量，x 3为第种肥料施用量，分别用 x1、x 2和 x3对产量 y 求一阶偏导数，得出它们的边际产量方程： 当三个边际产量方程都等于 0 时，解方程组可得到最高产量施肥量（x 1） max、（x 2） max 和（x 3） max。 如果(x 1)max、(x 2)max 和(x 3)max小于 0 时令其等于 0。将(x 1)max、(x 2)max 和(x 3)max代入肥效方程中可得到作 物的理论最高产量 ymax。 令作物产品价格为 Py,肥料价格分别为 Px1、Px 2和 Px3。用边际产值方程来求最佳经济施肥量(x 1)e、(x 2)e和 (x3)e。进一步可求出经济合理施肥量(x 1)R、(x 2)R和(x 3)R。因系数形式最佳经济施肥量和经济合理施肥量书写形式很 复杂，故此处不予给出。