模型参数灵敏度分析.docxVIP

。 返回主页 模型参数灵敏度分析 建立数学模型的主要目的之一是增进我们对系统的了解，而模型参数的灵敏度分析是对数 学模型的参数动态变化过程 , 即瞬时变化过程进行分析。 因此，通过模型参数的灵敏度分析可以明确哪些参数对系统的总体输出和动态影响较大。 方法简介 下面考虑一个经验模型，其模型的输出 y 可以是一种作物的产量，也可以是一头奶牛的总 泌乳量等， 假定对该模型已经圆满地进行了检验与评价， 包括对试验数据进行了适度的拟合。 模 型中有些参数是生理指标，而有些是环境指标，另有一些参数， 如对作物投入肥料 xi ，按其对产 量的相对效应进行排序， 则无论目标函数怎样， 都可以得到较为客观的度量。 目标函数 y 对参数 x i 的灵敏度 ( , x i ) 的定义为： S y 式中 为边际函数 ( 偏导数 ) ， 为平均的投入产出效应。 S( y, xi ) 表示目标 函数 y 对输入参数 xi 的灵敏度。 灵敏度分析的计算机处理 在分析模型参数变化速率之前 , 先给出模型方程， 并将模型及其参数按规定格式定义成公 式块，按系统要求将待分析的公式放在公式块的第 1 行，这些待分析的公式系由变量和参数组合 起来的表达式。在公式中用 x1 , x2, ? , xp 分别代表 p 个变量 ( 必须从１开始按顺序给出 ) 。公式 中可使用本系统的全部标准函数，最后的公式在形式上必须是合法的数学表达式，定义格式为 : 方程表达式 变量 1 的起始值 , 终止值 , 间隔值或变量 1 的取值水平。 变量 2 的起始值 , 终止值 , 间隔值或变量 2 的取值水平。 ? ? ? 变量 p 的起始值 , 终止值 , 间隔值或变量 p 的取值水平 例如 : 21.5+4.29 x1-3.77 x2-0.059 x1^2-0.015 x2^2+0.0044 x1x2 12, 18, 0.1 50 上述公式块中的第 2 行表示变量 x1 从 12 开始到 18，每隔 0.1 分析一次，第 3 行表 示将 x2 固定为 50。 例如研究产量随某种肥料用量变化的规律，求出的肥料反应经验方程为 ，根据该方程进行模型的灵敏度分析， 先按图 22-7 的方式编辑定义公式： 3x1+2x1x1-0.1x1^3 1 16 0.5 22-7 产量函数灵敏度分析公式定义图 然后进入菜单操作，选择“模型参数变化速率分析”菜单，执行后系统即输出分析结果， 包括灵敏度、导数 ( 偏导数 ) 、平均效应及目标函数的计算结果，详见如下输出单“ y x1 灵敏度 导数 ( 边际值 ) 平均效应 y/ x 目标函数 0.0000 ...3.0000... 0.0000 0.5000 1.2390 4.9250 3.9750 1.9875 1.0000 1.3673 6.7000 4.9000 4.9000 。 1 。 1.5000 1.4416 8.3250 5.7750 8.6625 2.0000 1.4848 9.8000 6.6000 13.2000 2.5000 1.5085 11.1250 7.3750 18.4375 3.0000 1.5185 12.3000 8.1000 24.3000 3.5000 1.5185 13.3250 8.7750 30.7125 4.0000 1.5106 14.2000 9.4000 37.6000 4.5000 1.4962 14.9250 9.9750 44.8875 5.0000 1.4762 15.5000 10.5000 52.5000 5.5000 1.4510 15.9250 10.9750 60.3625 6.0000 1.4211 16.2000 11.4000 68.4000 p class=MsoNormal align=center style=text-align: center; line-height: 10.0pt; mso-line-height-rule: exactly; mso-border-alt: solid windowtext .75pt; mso-padding-alt: 1.0pt 1.0pt 1.0pt 1.0pt; word-spacing: 0; border-style: none; border-width: medium; margin: 0; padding 。 2 。 欢迎您的下载， 资料仅供参考！ 致力为企业和个人提供合同协议， 策划案计划书， 学习资料等等 打造全网一站式需求 。 3