人口阻滞增长模型的参数估计与程序实现.pdfVIP

● 专 题 研 究 。 静 赫 _-I _． 糖 久 阻滞增长模 参数镶讲 程序杂瑰 ◎郑小洋 罗其莉 姚 蕾 宋江敏 (重庆理工大学数学与统计学院 400054) 【摘要】本文直接从建立阻滞增长模型出发，培养学生 分方程可 以通过 MATLAB软件求数值解 ，这些 内容可 以指 应用MATLAB解决实际问题的计算和估计参数 的能力． 导学生课外去学习、实践． 【关键词】阻滞增长模型；参数估计；最小二乘法；数值 3．估计模型参数 计算 人 口阻滞增长模型的解 (5)里含有两个参数 ，这两个参 引言 在大学生数学建模课程教学 中有人 口预测模型 数可以通过某一 国家或者某一地 区的实 际人 口统计数据作 这部分 内容，最开始导出的是人 口指数增长模型，但是 ，对 参数估计 ，本文采用最小二乘法估计参数．如果直接用教材 一 般情况 的人 口预测效果不怎么好 ，后又建立 了人 口阻滞 的线性二乘法估计参数 ，则需要对解作处理，可 以参 阅教 增长模型．这一过程 的推导和教学能够让学生体会到怎么 材．这里调用 MATLAB的函数 lsqcurvefit就可以完成参数估 样从一个粗略的模型提炼 出一个能描述、反映和解决实际 计 ，该函数的算法原理，可 以指导感兴趣 的学生课外学习． 问题 的数学模 型，从而可 以培养学生从实 际问题抽象 出适 本文以 《数学建模》教材 的美 国每隔 10年 (1790-- 1960)的 当的数学 问题 的建模能力．本文直接从建立阻滞增长模型 人 口统计数据 (×l0)估计该模 型的两个参数 和 rn．求 出发，重在培养学生 的建模思想、建模能力、以MATLAB软 解程序为 ： 件求解模型和估计模型参数的能力． 函数文件 logistic：functionY=logistic(a，x) 1．模型的建立 xO=3．9；Y=a(1)．／(1+(a(1)／xO一1)． exp(一a 人 口阻滞增长模型的原理 ：阻滞增长模型是考虑到 自 (2) x)) 然资源、环境条件等因素对人 口增长 的阻滞作用 ，对指数增 主函数 ： 长模型的基本假设进行修改后得到的．阻滞作用体现在对 X：[01234567891011121314151617]； 人 口增长率r的影响上，使得 r随着人 口数量 的增加而下 Y=[3．9，5．3，7．2，9．6，12．9，17．1，23．2，31．4，38．6， 降．若将 r表示为 的函数 r()，则它应是减函数，可以 50．2，62．9，76．0，92．0，106．5，123．2，131．7，150．7，179．3]； 得到 ： ％ 人 口数量 a0=[5700．03]； ％ 参数初值 =r() ，(0)= 。， (1) 【l [a，resnorm]=lsqcurvefit(@logistic，a