环境系统分析上机实验报告Matlab优化工具箱的用(I)――无条件极值的求解.docVIP

环境系统分析 上机实验报告 上机主题 Matlab优化工具箱的使用(I)――无条件极值的求解 上机时间 2013年9月22日 班 级 1127201 学 号 1112720113 姓 名 魏晓邢 评 分 　 市政环境学院 环境科学与工程系 指导教师 郑彤 上机实验目的 1、了解数量化理论，学习模型参数的估值方法，初步学习利用matlab最优化方法建立函数模型。 2、学习运用优化工具箱，运用fminunc语句和fminsearch语句求无条件极值。 3、掌握模型参数估计、灵敏度分析和修正的方法 4、进一步熟悉matlab的编程环境。 二．上机实验题目 题目3. 下表为一有机污染物进入水体后，其浓度的时间数据序列 时间（h） 0 1 3 5 7 9 23 27 31 浓度(mg/l) 2.30 2.22 1.92 1.6 1.52 1.07 0.73 0.50 0.45 试确定其浓度～时间关系表达式的形式，并估计其中的参数，给出估参所用的Matlab程序及最后的目标函数值（计算值与实测值离差的平和） 三．上机程序要点 1.首先根据列表画出时间与浓度的图形 程序： 2．经验公式的选择 由于题目给出的数据是污染物（有机物）进入水体后的浓度变化，通常采用公式 公式中含有三个参数，运用无条件极值法估参，确定函数模型。 3．估参 1）建立离差平方和的函数 命名为‘von_f’（function语句） 程序： 2）引用函数f，求极值（fminunc语句），得到三个参数值和函数模型。 程序： 运行结果： 故，估得得参数分别为0.3742、1.9787和-0.0923，中值误差为0.0631 4．灵敏度计算 1）c对ai1的灵敏度 图像： 2）c对ai2的灵敏度 图像： 3）c对ai3的灵敏度 图像： 5.结果分析 由以上可知：c对ai1，ai2，ai3的灵敏度均在[-1 1]范围内，满足误差要求。则目标函数为： 6.拟合结果： 四．上机实验心得 1.关于优化工具箱 在实际解题中因为使用了优化工具箱中的函数（fminunc）而使问题得到了简化。fminunc是一种便捷的极值求法，运用中可以节省时间。 2.关于此类题的解题方法 在解决此类问题时，应当把握好思路，即分析数据—作图—选择经验公式—估参—检验灵敏度—确定函数模型。只有清楚解决方法和思路，才能高效地完成此类问题。另外，此题也可以用网格法，梯度法等方法进行求解。