准备工作

题目

人口增长模型的确定
1790-1980年间美国每隔10年的人口记录如下表所示。
1.试用以上数据建立马尔萨斯(Malthus)人口指数增长模型,并对接下来的每隔十年预测五次人口数量,并查阅实际数据进行比对分析。
2.如果数据不相符,再对以上模型进行改进,寻找更为合适的模型进行预测,并对两次预测结果进行对比分析。
3.查阅资料找出中国人口与表1同时期的人口数量,用以上建立的两个模型进行人口预测与分析。

理解算法

http://www.doc88.com/p-773826839689.html
http://www.docin.com/p-938750966.html
https://wenku.baidu.com/view/df497b3410661ed9ad51f36f.html

参考代码

http://www.doc88.com/p-414778587173.html

难点、问题

  1. http://www.docin.com/p-938750966.html 中公式4到公式5的推导 — 未解决
  2. 已知公式5,和表格数据,怎么求出r与Xm。— 已知公式+数据,求参数(即拟合问题)
  3. 预测是指数增长,真实数据不是指数增长。

实现算法

##马尔萨斯

clear;
clc;%% ---------------------------------------------------
% 马尔萨斯1:求模型
% ------------------------------------------------------------
t = (0:19);                 % 1790-1980 共20年
x = [3.9, 5.3,  7.2,  9.6,   12.9,  17.1,  23.2,  31.4,  38.6,  50.2, 62.9, 76.0, 92.0, 106.5, 123.2, 131.7, 150.7, 179.3, 204.0, 226.5];
subplot(1, 2, 1);
plot(t, x, ' o');
y = log(x);
p = polyfit(t, y, 1)        % 求出p, y = p[1] * t + p[2]
% 得到,p[1] = 0.2142, p[2] = 1.7213
% 所以,x0 = e^p[2] = e^1.7213 = 5.5917
x0 = exp(p(2))
% 所以,x(t) = 5.5917 * e^0.2142t%% ---------------------------------------------------
% 马尔萨斯2:预测
% ------------------------------------------------------------
tPredict = (20:24);
xPredict = 5.5917 * exp(0.2142 * tPredict);
subplot(1, 2, 2);
plot([t,tPredict], [x , xPredict], ' *');
% xPredict = [405, 502, 622, 771, 955]

##logistic 阻滞增长
已知公式+数据,求参数(即拟合问题)

clear
clc% 准备数据
t = (0:19);                 % 1790-1980 共20年
x = [3.9, 5.3,  7.2,  9.6,   12.9,  17.1,  23.2,  31.4,  38.6,  50.2, 62.9, 76.0, 92.0, 106.5, 123.2, 131.7, 150.7, 179.3, 204.0, 226.5];
plot(t, x, ' o');
hold on;% 已知x(0) = 3.9
myfunc = @(a, t)[a(1) / ((1 + a(1) / 3.9) - 1) * exp(-a(2) * (t - 0))];
a0 = [500, 1];
a = lsqcurvefit(myfunc, a0, t, x);
disp(['a= ' num2str(a)]);atPredict = (20:24);
xPredict = myfunc(a, tPredict);
plot([t, tPredict], [x, xPredict], ' *r');

题目其他要求

对比真实数据

真实数据1
真实数据2
对比以后发现并不理想,下午继续看。

数学建模2-美国人口增长模型的确定相关推荐

  1. 【Python数学建模】SEIR传染病模型模型延伸-SEIDR模型(一),加入疫苗接种、政府管控、病毒变异等因素的影响

    目录 一. SEIR传染病模型 二. SEIR模型的延伸--SEIDR模型 三. 模型延伸--影响因素1:疫苗接种 四. 模型延伸--影响因素2:政府管控 五. 模型延伸--影响因素3:病毒变异 写在 ...

  2. 数学建模学习记录——数学规划模型

    数学建模学习记录--数学规划模型 一.线性规划问题 MatLab中线性规划的标准型 MatLab中求解线性规划的命令 二.整数线性规划问题 三.非线性规划问题 MatLab中非线性规划的标准型 Mat ...

  3. Python小白的数学建模课-09.微分方程模型

    小白往往听到微分方程就觉得害怕,其实数学建模中的微分方程模型不仅没那么复杂,而且很容易写出高水平的数模论文. 本文介绍微分方程模型的建模与求解,通过常微分方程.常微分方程组.高阶常微分方程 3个案例手 ...

  4. dna序列分类数学建模matlab,数学建模DNA序列分类模型(终稿).doc

    您所在位置:网站首页 > 海量文档 &nbsp>&nbsp高等教育&nbsp>&nbsp生物学 数学建模DNA序列分类模型(终稿).doc32页 本文 ...

  5. 【数学建模】常用微分方程模型 + 详细手写公式推导 + Matlab代码实现

    文章目录 一.学习内容 二.学习时间 三.学习产出 3.1 微分方程基本概念 3.2 微分方程在数学建模中的应用 3.3 微分方程常用模型 3.3.1 人口增长模型 3.3.1.1 指数增长模型(马尔 ...

  6. 【数学建模】种群竞争模型(最优化)

    文章目录 一.算法介绍 二.适用问题 三.算法总结 1.竞争模型 四.应用场景举例 五.MATLAB代码 六.实际案例 七.论文案例片段(待完善) 种群竞争模型主要针对数学建模问题中的一些小的子问题进 ...

  7. 2023电工杯数学建模B题思路模型

    文章目录 0 赛题思路 1 竞赛信息 2 竞赛时间 3 组织机构 4 建模常见问题类型 4.1 分类问题 4.2 优化问题 4.3 预测问题 4.4 评价问题 0 赛题思路 (赛题出来以后第一时间在C ...

  8. 2023电工杯数学建模A题思路模型代码

    占个位置吧,开始在本帖实时更新电工杯数学建模赛题思路代码,文章末尾获取! A题思路分析 问题一:1典型住户电采暖负荷用电行为分析 (1)在满足温控区间约束条件下,分析典型房间温变过程微分方程稳态解的性 ...

  9. 2023第十五届电工杯数学建模AB题思路模型

    文章目录 0 赛题思路 1 竞赛信息 2 竞赛时间 3 组织机构 4 建模常见问题类型 4.1 分类问题 4.2 优化问题 4.3 预测问题 4.4 评价问题 0 赛题思路 (赛题出来以后第一时间在C ...

  10. 数学建模--时间序列分析、模型预测

    学习自:b站 清风数学建模 第1部分_时间序列分析的概念与时间序列分解模型 时间序列 时间序列的基本概念 区分时期和时点序列 时期序列适用于灰色预测模型 时间序列分解 时间要小于1年,才能进行时间序列 ...

最新文章

  1. centos7 安装 wireshark
  2. 手持机设备公司(WINCE/ANDROID/LINUX)
  3. 不同操作系统之间的网络配置
  4. Effective Java之请不要在新代码中使用原生态类型(二十三)
  5. AcWing 703. 数独检查
  6. 【Python】Python3编码规范
  7. 通用采集器Modbus协议应用
  8. 制作一个简单HTML静态网页(HTML+CSS)
  9. 推荐育儿书《正面管教》
  10. 电源输出的Overshoot和Undershoot 测试
  11. R语言中的缺失值处理
  12. 将git项目提交到其他的git仓库
  13. 郭全中:对大数据的认识该正本清源
  14. 关系抽取(二)远程监督方法总结
  15. 如果结束进程拒绝访问,可以尝试以下
  16. MCU-51:单片机蜂鸣器播放音乐和提示音
  17. 一位金融工程小硕的华丽逆袭人生!超真实Quant菜鸟的修行路
  18. python查看包的依赖关系
  19. 1.1 网络安全的前景
  20. 国际原子时、协调世界时、世界时和GPS时间系统之间的关系

热门文章

  1. CAPL学习之路-诊断函数
  2. 画手机原型图工具推荐
  3. 百度编辑器(ueditor)魔改:5、自定义一键排版(按照需求保留标签和样式)
  4. linux命令行启动向日葵无法编辑,向日葵linux版无法启动
  5. 如何划分机器学习的训练集和测试集
  6. Turbo码 - 接近完美的编码
  7. 非线性调频信号matlab仿真,线性调频信号回波仿真
  8. 计算机无法访问权限,局域网共享无法访问权限不足的解决方法
  9. android 存储盘 dcim,Android上的DCIM目录路径 – 返回值
  10. 使用Android Studio向SVN上传新项目