参数拟合题目.docVIP

习题1 酶促反应是使用酶作催化剂的化学反应，反应物又称作底物，描述反应速度和底物浓度关系的数学模型是Michaelis-Menten模型：，其中是底物浓度，是反应速度，和 是待定参数. 请根据表1.15的数据，完成下列数据拟合问题： 表1.15 酶促反应试验中的反应速度和底物浓度数据 底物浓度 0.02 0.06 0.11 0.22 0.56 1.10 反应速度 76 47 97 107 123 139 159 152 191 201 207 200 请用MATLAB统计工具箱的函数nlinfit计算待定参数和 的最小二乘估计值，要求写出程序，给出拟合参数和误差平方和的计算结果，并展示拟合效果图； 通过变量替换，可以将属于非线性模型的Michaelis-Menten模型转化成线性模型，并用MATLAB函数polyfit进行计算，请说明转化成线性模型的详细过程，然后写出程序，给出拟合参数和误差平方和的计算结果，并展示拟合效果图； 请分析Michaelis-Menten模型非线性拟合和线性拟合的结果有何区别？原因是什么？ 土豆生长所需的主要营养素是氮(N).磷(P).钾(K). 某作物研究所在某地对土豆做了一定数量的实验，取得的实验数据如表1.16所示,其中ha表示公顷，t表示吨，kg表示公斤. 当一个营养素的施肥量变化时，总将另外两个营养素的施肥量保持在第七个水平上，如对土豆产量关于N的施肥量做实验时，P与K的施肥量分别取为196kg/ha与372kg/ha. 请分析土豆的施肥量与产量之间的关系，要说明选择什么函数模型，为什么选择这些函数模型；要给出拟合参数、误差平方的计算结果，并展示拟合效果图. 表1.16 土豆的施肥量和产量实验数据 N施肥量 /(kg/ha) 产量 /(t/ha) P施肥量 /(kg/ha) 产量 /(t/ha) K施肥量 /(kg/ha) 产量 /(t/ha) 0 15.18 0 33.46 0 18.98 34 21.36 24 32.47 47 27.35 67 25.72 49 36.06 93 34.86 101 32.29 73 37.96 140 38.52 135 34.03 98 41.04 186 38.44 202 39.45 147 10.09 279 37.73 259 43.15 196 41.26 372 38.43 336 43.46 245 42.17 465 43.87 404 40.83 294 40.36 558 42.77 471 30.75 342 42.73 651 46.22 1928年，美国经济学家C.Cobb和P.Douglas在他们关于1899年至1922年美国经济增长的研究报告中提出了生产函数模型. 他们使用美国政府发表的经济数据（见表1.18），以1899年为基准，1899年的Q（产值），K（投资），L（劳动力）都设为100，其他年份的数据表示成1899年数据的百分数，用最小二乘法拟合出生产函数模型中的待定参数c和α. 表1.18 C.Cobb和P.Douglas使用的政府数据 年 份 1899 1900 1901 1902 1903 1904 1905 1906 Q 100 101 112 122 124 122 143 152 K 100 107 114 122 131 138 149 163 L 100 105 110 117 122 121 125 134 年 份 1907 1908 1909 1910 1911 1912 1913 1914 Q 151 126 155 159 153 177 184 169 K 176 185 198 216 216 226 236 244 L 140 123 143 147 148 155 156 152 年 份 1915 1916 1917 1918 1918 1920 1921 1922 Q 189 225 227 223 218 231 179 240 K 266 298 335 366 387 407 417 431 L 156 183 198 201 196 194 146 161 请用MATLAB统计工具箱的函数nlinfit计算生产函数的数据拟合问题，要求写出陈旭，给出拟合参数和误差平方和的计算结果，并展示拟合效果图. 通过变量替换，可以将属于非线性模型的生产函数转化成线性模型，并用MATLAB函数polyfit进行计算，请说明转化成线性模型的详细过程，然后写出程序，给出拟合参数和误差平方和的计算结果，并展示拟合效果图.