二阶差分方程预测模型matlab程序,差分方程预测模型.doc
差分方程预测模型
差分方程人口预测模型
一、名词和符号说明
名词解释:
(1)拟合: 对于某个变化过程中的多个相互依赖的变量,可建立适当的数学模型,用于分析预报决策或控制该过程.对于两个变量可通过用一个一元函数去模拟这两个变量的取值.用不同的方法可得到不同的模拟函数.下面使用图表介用Mathematica做曲线拟合。
(2)差分方程:含有自变量,未知函数以及未知函数差分的函数方程,称为差分方程。
(3)迭代法:是牛顿在17世纪提出的一种求解方程f(x)=0.多数方程不存在求根公式,从而求精确根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要。 设r是f(x)=0的根,选取x0作为r初始近似值,过点(,f())做曲线y=f(x)的切线L,L的方程为,求出L与x轴交点的横坐标 ,称为r的一次近似值,过点(,f())做曲线y=f(x)的切线,并求该切线与x轴的横坐标称为r的二次近似值,重复以上过程,得r的近似值序列{Xn},其中,称为r的n+1次近似值。上式称为牛顿迭代公式。 第 k年i岁的女性总人数
女性人口的(按年龄)分布向量
第k年i岁的女性生育率
第k年i岁的女性死亡率
第 k年i岁的女性存活率
i岁女性的生育模式 k年总和生育率(控制人口数量的主要参数)
A 存活率矩阵
B 生育模式矩阵
二、模型假设
针对本题中出现的数据的代表意义和建立模型时能够使问题理想化、简单化,我们应用已知数据,将时间离散化, ,因此本模型考虑女性人口的发展变化假设女性最大年龄为岁,,1年为1个时段,不考虑同一时间间隔内人口数量的变化不考虑生存空间等自然资源的制约,不考虑意外灾难等因素对人口变化的影响生育率仅与年龄有关,存活率也仅与年龄有关我们来建立一个离散的人口增长模型, 女性人口的发展变化Leslie人口模型,用差分方程:第k年i岁的女性生育率; : k年总和生育率,或生育胎次;
:第k年i岁的女性死亡率; :第 k年i岁的女性存活率
: i岁女性的生育模式
,
用表示女性人口的(按年龄)分布向量,记A=
B=则模型应表示为:
=A+B
利用matlab软件编程求解,程序如下:
c=zeros(91);
d1=[ … … ];
for i=1:91
for j=1:91
if i==j
c(i+1,j)=d1(i)
end
end
end
A=c1
a1=[ … … ];
b=zeros(91);
for i=1:35
b(1,i+15)=a1(i)
end
B=b1;
=[ … …] %2001对应初始值
y=zeros(91,n)%n表示要预测年数
y(:,1)= ;
for k=1:19
y(:,k+1)=A*y(:,k)+(k)*B*y(:,k)
end
(一)用此模型预测中短期女性人口变化趋势
考虑到男女性别比例波动不大,所以女性人口数量的发展趋势可以预测全国总人口的发展趋势。
对所给数据进行处理,发现近期(k)变化很小,这里我们取=/5即:市:=1;镇:=1.254;乡:=1.649,代入模型方程,得:
x(k)=……………………………………………………………(3.3-1)
x(k)=………………………………………………………(3.3-2)
x(k)=………………………………………………………(3.3-3)
分别代入k=20,即可算出市、乡、镇从2001年到未来20年的预测数值。
分别取2002、2004年的数据拟合,情况如下:
图3-1 2002、2004拟合趋势图
由上图可看出,拟合情况较好,此模型可用于短期预测,预测趋势图如下:
图3-2
预测数据表为:
表3-1
年份20012005201020152020女性人口数5.97E+086.37E+086.83E+087.34E+087.74E+08通过上面的预测数据和图像,可看出2020年之前女性人口呈增长趋势,全国人口总数也呈增长趋势。
(二)长期预测
进行长期预测时,考虑到国家计划生育一对夫妇只生一个孩子的政策,取=1,则模型可化简为 = 其中为2001年女性人口分布向量。
图3-3模型检验拟合图
利用数据来检验我们建立的差分方程模型,发现数据基本吻合,说明模型是很准确的,可以用此模型进行长期预测。
利用方程预测的女性总人口数据如下:
表3-2
年份20012005201020152020202520302035204020452050女性总人口
二阶差分方程预测模型matlab程序,差分方程预测模型.doc相关推荐
- 威布尔 风速 matlab,风速威布尔分布和ARMA预测模型matlab程序.doc
风速威布尔分布和ARMA预测模型matlab程序 clc clear %% 1.计算风速weibull分布 % 数据处理 load data; mu=mean(speed);%原始数据的统计参数 si ...
- matlab三角形单元,平面三角形单元常应变单元matlab程序的编制.doc
平面三角形单元常应变单元matlab程序的编制.doc 1三角形常应变单元程序的编制与使用有限元法是求解微分方程边值问题的一种通用数值方法,该方法是一种基于变分法(或变分里兹法)而发展起来的求解微分方 ...
- matlab pls rmsecv,偏最小二乘法PLS回归NIPALS算法及Matlab程序及例子.doc
偏最小二乘法PLS回归NIPALS算法及Matlab程序及例子 偏最小二乘法PLS回归NIPALS算法的Matlab程序及例子 function [T,P,W,Wstar,U,b,C,B_pls,.. ...
- 灰色预测方法预测温度matlab,灰色预测模型matlab程序
%by allen @ 红嘴海鸥 %灰色模型预测是在数据不呈现一定规律下可以采取的一种建模和预测方法,其预测数据与原始数据存在一定的规律相似性 %下面程序是灰色模型GM(1,1)程序二次拟合和等维新陈 ...
- svm预测模型 matlab,SVM回归预测模型
我把14章,SVM用于回归建模的例子改动后,结果出现了问题: 首先,我描述一下,我想做的东西: 我要预测输出的压力,当前的输出的压力和前3个时刻的5个变量有关,我也想减少变量,可是预测的效果就不好了) ...
- Caputo 分数阶微分方程-慢扩散方程初边值问题基于L1 逼近的空间二阶方法及其Matlab程序实现
2.3.3 Caputo 分数阶一维问题基于 L1 逼近的空间二阶方法 考虑如下时间分数阶慢扩散方程初边值问题 { 0 C D t α u ( x , t ) =
- Caputo 分数阶微分方程-慢扩散方程初边值问题基于快速 L2-1σ 逼近的空间二阶方法及其Matlab程序实现
2.3.6 Caputo 分数阶一维问题基于 L2-1σ逼近的快速差分方法 考虑如下时间分数阶慢扩散方程初边值问题 { 0 C D t α u ( x , t ) =
- 标号法matlab程序,matlab程序大全
上例 Matlab 程序如下: %使用表8.4 给出的小波基 ld=[0 0.0... matlab经典代码大全_计算机软件及应用_IT/计算机_专业资料.哈哈哈 MATLAB 显示正炫余炫图:plo ...
- 组合预测模型 | PSO-ELM、ELM极限学习机数据回归预测对比(Matlab程序)
组合预测模型 | PSO-ELM.ELM极限学习机数据回归预测对比(Matlab程序) 目录 组合预测模型 | PSO-ELM.ELM极限学习机数据回归预测对比(Matlab程序) 预测结果 基本介绍 ...
- gm11模型 MATLAB源代码,灰色预测模型 GM(1,1)- 级比检测不通过 - 平移变换常数c的确定(内含代码)...
灰色预测模型 GM(1,1) - 级比检测不通过 - 平移变换常数c的确定(内含代码) 谢谢! (哇-------没注意就5600多访问量了,激动!!!!我的初衷只是想把blog作为我的学习笔记,回首 ...
最新文章
- 好用的Span(富文本)api
- Java的新项目学成在线笔记-day10(三)
- Python Django 多表设计OneToOneField/ManyToManyField/ForeignKey
- Java黑皮书课后题第5章:*5.31(金融应用:计算CD价值)假设你用10 000美元投资一张CD,年利率为5.75%。编写程序,提示由用户输入一个金额数、年获利率、月份数,然后显示一个表格
- jquery in action 学习笔记
- 洛谷 P2040 打开所有的灯-dfs
- html:(16):a标签
- Qt工作笔记-QXmlStreamReader中的字符编码的坑
- mysql 命令 utf8_Mysql 统一设置utf8字符
- HttpInvoker-----客户端实现
- loss函数之TripletMarginLoss与TripletMarginWithDistanceLoss
- 拓端tecdat|r语言有限正态混合模型EM算法的分层聚类、分类和密度估计及可视化
- 博途调试g120_【免费资料】西门子变频器调试软件汇总
- 微分中的dx和delta x
- 自由技艺 (Liberal arts)
- 获取百度首页的源代码
- 三进制 四进制计算机原理,基因编码为何选择了“四进制”,而不是计算机系统的“二进制”?...
- 数据传输加密——非对称加密算法RSA+对称算法AES(适用于java,android和Web)
- 少儿编程----颜值打分系统
- 广告业务系统 之 核心通道 —— “日志中心-s2s监测上报”
热门文章
- 5行python代码帮你搞定百度文库复制问题
- 《移动通信原理与系统》——第二章《移动通信电波传播与传播预测模型》——笔记
- java面试宝典2019
- java velocity 语法_Velocity的基本语法
- 最佳的75个安全检测工具
- Python实现本地翻译API
- 微信小程序位置定位php,【微信开发】微信小程序通过经纬度计算两地距离php代码实现...
- Landsat卫星MSS/TM/ETM数据(转自ESRI社区)(二)
- 设置SQL Server 2008 以允许远程连接
- 怎么用Wireshark抓包QQ的到对方的IP地址