响应曲面法

反应曲面法（Response surface methodology，简写RSM）为结合数学与统计而延生出的方法，为最适实验设计或作业条件的有利工具，于1951年，Box 和 Wilson 共同进行数学模式的建立与推导，而后普遍应用于电子、机械、农业、化学工业、生物科技、材料科学、食品科学及工业制程改善等各项研究领域中。

**说明**
反应曲面法在协助研究人员对科学系统或工业制程中最佳产品设计、制程改善、系统最佳化等问题提供一套分析、求解程序，大部分应用时机均属工业性研究，尤其是当系统特性受大量变数影响状况下最为适当。

**历史**

**概念**
反应曲面法之研究问题，一般假设问题为限制性之最佳化问题，目标函数的确切型式是未知的 ， 为误差，反应曲面法一般在此前提的假设与应用系统的限制下，可有效地求得最佳实验或作业变数值。
一般来说，执行反应曲面法大致分为两阶段:
*第一阶段称为反应曲面设计（response surface design）
*第二阶段称为反应曲面最佳化（response surface optimization）

**反应曲面设计**
为探讨独立变数与反应变数之间的数学模式关系，因此欲对于反应和独立变数之间找出一个适当的近似函数。通常利用独立变数在一些范围里的低阶多项式近似，即为一阶回归模型 （first-order model），
如果系统中有曲率，则必须利用较高阶的多项式，如二阶模型（second-order model）。
获得最适化实务模型便是本阶段最重要的议题。收集资料后以最小平方法 （least squares estimation, LSE） 配适，以寻找出一个适当近似的函数，采用回归分析的显著性检定 （general linear test approach） 来了解独立变数与反应变数间的关系强弱，并检定配适的模式是否恰当 （statistical adequacy）。当实验区域接近最佳反应值附近时，真实反应曲面的曲率 （curvature） 会增加，则考虑二阶模型，同样的，我们需要检定二阶模式的适当性。当这个二阶回归模式配适良好时，便可以利用这二阶模式求得最适操作点及特征化反应曲面。

**反应曲面最佳化**

**优点**
#经济性原则：反应曲面法可以使用部分因子设计或特殊反应曲面设计（如混种设计等 （hybrid design）），以较少的实验成本及时间获得不错且有效的资讯。
#深入探讨因子间交互作用影响：反应曲面法可以经由分析与配适模式来研究因子间的交互作用，并且进而讨论多因子对反应变数影响的程度。
#获得最适化的条件：根据数学理论求得最适的实验情况，同时利用配适反应方程式绘出模式三度空间曲面图与等高线图，观察并分析出最适的操作条件。
#减少模拟时间：可获得模拟独立变数与响应变数关系之数学模型，借此将实验次数及时间降低。

**限制**
在应用上主要存在下列二项限制：
#只适用于连续性的系统，是假设所有反应值与独立变数的量测刻度是连续性的。
#影响系统之独立变数（可控制和不可控制变数）是属于计量性。

**应用**

**背景**

**问题与方法**
而这里主要探讨的是行动式随意网路里 IEEE 802.11 中的 媒体存取控制 （Medium access control, MAC） 协定和随意随选距离向量 （Ad hoc On demand Distance Vector, AODV） 路由协议间的通讯。
由于目前没有一个通用的方法论来确认和最适化协定之间的通讯，这里利用反应曲面法来解决这些问题。
考虑三个可能影响的因子：
A：Active route timeout
B：Max route request wait timeout
C：Max retransmissions
将三个因子分别设定成二个水准（利用 ns-2 模拟器中 AODV 软体分配的预设值）
而网路性能上考虑以下二个变数做为反应变数：
Average throughput：目的地节点每秒收到的数量（bytes/sec）
Average packet delay：一个封包从来源送到目的节点所花的平均时间（sec）
因此我们将网路分析的问题转变为二水准三因子实验设计（response surface design）
接下来利用反应曲面法寻求最佳化 （response surface optimization） ，使得 average throughput 最大，average packet delay 最小。
并且在某个限制式下，寻求最佳解，例如，延迟时间在 400ms 之内。
可利用 ns-2 模拟器来模拟出 20 次不同的流动情境，拿预设值和用 RSM 找出的最佳解做比较。

**因子动态**
在不同的网路设定下，有可能导致不同因子对于反应变数的影响程度。因此对于因子在不同网路环境下之表现行为的了解，对 RSM 是很重要的。
不仅可以决定超初的搜寻范围，亦可决定最陡上升（下降）法的步伐大小。
某些应用中，平均网路速度及节点密度是已知的条件，在自然世界的动态中可能会不合适，因此在分析时亦需将当下的网路条件随时做改变。
关于这个问题，可将因子配上一个权重表示，设定所有因子的权重起始值为1，而权重将随着网路设定的不同而改变；0表示对反应变数没有影响。

**相关工作**
此例子为在 IEEE 802.11 中的媒体存取控制协定和随意随选距离向量路由协定中的通讯，其他的协定下的补充研究可参考
*AODV：E.M. Royer and C.E. Perkins
*DSR：Maltz, J. Broch, and D.B. Johnson
*ODMRP：S.H. Bae, S.-J. Lee, and M. Gerla
*ABR：C.K. Toh and M. Delawar

**参考文献**
K.K. Vadde, V.R. Syrotiuk, D.C. Montgomery, “Optimizing Protocol Interaction Using Response Surface Methodology,”IEEE Transations On Mobile Computing, Vol. 5, NO. 6, June 2006.
*E.M. Royer and C.E. Perkins, “An Implementation Study of the AODV Routing Protocol,” Proc. IEEE Wireless Comm. and Networking Conf. （WCNC）, vol. 3, pp. 23-28, Sept. 2000.
*D.A. Maltz, J. Broch, and D.B. Johnson, “Experiences Designing and Building a Multi-Hop Wireless Ad-Hoc Network Testbed,” Technical Report CMU-CS-99-11, Carnegie Mellon Univ., Mar. 1999.
*D.A. Maltz, J. Broch, and D.B. Johnson, “Lessons from a Full-Scale Multihop Wireless Ad Hoc Network Testbed,” IEEE Personal Comm. Magazine, vol. 8, no. 1, pp. 8-15, Feb. 2001.
*S.H. Bae, S.-J. Lee, and M. Gerla, “Unicast Performance Analysis of the ODMRP in a Mobile Ad-Hoc Network Testbed,” Proc. IEEE Conf. Computer Comm. and Networks （ICCCN ’00）, pp. 148-153, Oct. 2000.
*C.K. Toh and M. Delawar, “Implementation and Evaluation of an Adaptive Routing Protocol for Infra-Structureless Mobile Networks,” Proc. IEEE Int’l Conf. Computer Comm. and Networks, pp. 16- 18, Oct. 2000.

**外部连结**
  Response surface designsdate=20201203031554 -  Matlab SUrrogate MOdeling Toolbox - SUMO Toolbox - Matlab code for Response Surface Modeling