二阶差分方程预测模型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作为我的学习笔记,回首 ...
最新文章
- 教你在Tableau中绘制蝌蚪图等带有空心圆的图表(多链接)
- 研究生数学建模比赛日程安排
- 重磅!全网最全13010本SCI2020最新影响因子下载!CA+四大神刊!预警期刊最新影响因子一览!
- php中的getdate+函数,PHP中的getdate()函数
- 转载爱哥自定义View系列--Paint详解
- ASP.NET Core 认证与授权[5]:初识授权
- python 人脸关键点检测_opencv+python+dlib人脸关键点检测、实时检测
- linux批量筛选序列变异位点,使用bedtools获取指定坐标上下游的序列
- ucgui移植到rt_thread simulator
- Rust: mod文件、main文件调用
- 用PS将照片背景变成白色
- intent-filter属性介绍
- 《一万小时理论》读书笔记
- Classic界面chatter中的子选项卡配置
- 如何学会做一名优秀的下属
- 如何给Windows计算机加域
- 《汇编语言》第十章 call 和 ret 指令
- uni.showLoading,uni.hideLoading,uni.showToast,uni.showModal
- AVI和WAV格式文件分析
- 英文邮件:greetings
热门文章
- [android]qq输入法手机版 v1.0 build0580,腾讯Android版手写输入法惊艳亮相 手写更畅快...
- Unity3D 脚本模板插件
- [转]PKM-个人知识体系建设
- 计算机音乐apple,Apple Music苹果音乐电脑版
- 如何将Excel工作表另存为PDF
- python arp 网关_python arp欺骗伪造网关代码详解
- vue zxing 实现一维码、二维码扫描可移动端h5使用代码亲测可用
- java学习手册下载_java学习手册
- 《细说PHP》读书笔记
- 几行Java解决图片提取文字功能