数学模型小作业(Malthus_Logistic模型)

最新推荐文章于 2024-05-13 08:33:52 发布
最新推荐文章于 2024-05-13 08:33:52 发布
阅读量453 收藏 7
点赞数 12
文章标签： 算法 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/gisshow135540019/article/details/135122950
版权 
//--------------------------------------------------------------------------------------------------------------------
//数学模型小作业(Malthus_Logistic模型)
//题目描述:
//01.Malthus模型称为指数增长模型,它假设:
//02.	a.增长率r是常数.
//03.	b.单位时间增长量与当时数量成正比.
//04.Malthus模型的表达式为:x(t)=x0*[e^(rt)].
//05.Logistic模型称为阻滞增长模型,它的假设如下:
//06.	a.随着人口增长,资源环境开始对人口起阻滞作用,增长率r是当时数量x的线性递减函数r(x),假设r(x)=r-a*x.
//07.	b.xm为最大人口容量(假设).
//08.	c.r表示固有增长率,即人口很少时的增长率(对应时间单位).
//09.	d.线性递减函数r(x)的表达式:x->0,r(x)=r;x->xm,r(x)->0,即0=r-a*xm,解出:a=r/xm,r(x)=r(1-x/xm).
//10.Logistic模型的表达式为:x(t)=xm/{1+(xm/x0-1)*[e^(-r*t)]}.
//11.最大容量xm是根据当时资源环境推测的,随着科技进步,它也是可以改变的.
//12.预测时间越远,保持同样误差的概率越小(即可靠度越低),故每隔若干年要对模型中的参数作修改.
//程序编制要求:
//1.初始数量x0,固定增长率r,最大容量xm交互读取.
//2.计算两种模型的数量值随时间的变化规律.
//3.绘制时间从(0~10)单位变化时,数量值随时间的变化曲线.
//--------------------------------------------------------------------------------------------------------------------
#include <graphics.h>//包含Easyx模拟TC的BGI绘图库头文件
#include <math.h>//包含数学运算头文件
#define DD 40//宏定义标题边框距离
#define DW 10//宏定义宽度边框距离
#define DH 10//宏定义高度边框距离
#define DC 200//宏定义控件区域宽度
#define nx 10//宏定义x轴刻度数量
#define ny 10//宏定义y轴刻度数量
#define PI 3.1415926//宏定义圆周率
#define TA 0//起始时间单位
#define TB 20//累计时间单位
#define c1 RGB(38,47,86)//宏定义辅助线颜色1
#define c2 RGB(38,47,86)//宏定义辅助线颜色2
#define c3 RGB(255,255,255)//宏定义文字及边框颜色
#define cc1 RGB(255,0,0)//宏定义曲线颜色1(红色)
#define cc2 RGB(0,0,255)//宏定义曲线颜色2(蓝色)
#define cc3 RGB(255,255,0)//宏定义曲线颜色3(黄色)
//定义结构体-------------------------------------------------------------------------------------------------------
typedef struct points
{
	double x,y;//定义双精度变量
}DPOINT;//结束结构体
//定义全局变量-----------------------------------------------------------------------------------------------------
int W = GetSystemMetrics(SM_CXSCREEN);//获取整个屏幕右下角X坐标
int H = GetSystemMetrics(SM_CYSCREEN);//屏幕Y坐标 	
char is1[20]="3.9";//预定义:1.数量初值x0
char is2[20]="0.31";//2.固定增长率r
char is3[20]="200";//3.(Log)最大容量xm
char is4[20]="5";//4.输入时间t
char is5[20],is6[20],is7[20];//定义字符数组(存储字符串)
double d1,d2,d3,d4,d5,d6,d7,d8,d9,d10,d11;//定义双精度变量(存储计算值)
//函数声明---------------------------------------------------------------------------------------------------------
void sub_frame();//绘制框架子函数
void sub_calculate();//计算数据子函数
//主函数-----------------------------------------------------------------------------------------------------------
void main()
{
	sub_frame();//绘图初始化
	MOUSEMSG m;//定义鼠标消息 
	while(true)//循环
	{
		m=GetMouseMsg();//获取一条鼠标消息 
		switch(m.uMsg)//根据获得的消息选择分支
		{
		case WM_LBUTTONDOWN://鼠标左键单击时判断数据输入
			{
				//判断鼠标位置是否在指定的控件区域
				//A.获取原始数据(横坐标在输入区域时输入数据)-------------------------------
				if(m.x>=W-DW-DC && m.x<=W-DW)
				{
					//1.数量初值x0:-----------------------------------------------------
					if(m.y>=H-DH-16*36+1 && m.y<=H-DH-15*36-1)//
					{
						InputBox(is1,9,"请输入数据");//读入数据
						setfillcolor(BLACK);
						fillrectangle(W-DW-DC+DC/2+10,H-DH-16*36+1,W-DW,H-DH-15*36);
						outtextxy(W-DW-DC+DC/2+15,H-DH-16*36+15,is1);//写控件字符串
					}//结束if
					//2.固定增长率r:-----------------------------------------------------
					if(m.y>=H-DH-15*36+1 && m.y<=H-DH-14*36-1)//
					{
						InputBox(is2,9,"请输入数据");//读入数据
						setfillcolor(BLACK);
						fillrectangle(W-DW-DC+DC/2+10,H-DH-15*36+1,W-DW,H-DH-14*36);
						outtextxy(W-DW-DC+DC/2+15,H-DH-15*36+15,is2);//写控件字符串
					}//结束if
					//3.(Log)最大容量xm:-----------------------------------------------
					if(m.y>=H-DH-14*36+1 && m.y<=H-DH-13*36-1)//
					{
						InputBox(is3,9,"请输入数据");//读入数据
						setfillcolor(BLACK);
						fillrectangle(W-DW-DC+DC/2+10,H-DH-14*36+1,W-DW,H-DH-13*36);
						outtextxy(W-DW-DC+DC/2+15,H-DH-14*36+15,is3);//写控件字符串
					}//结束if
					//4.输入时间t:-----------------------------------------------
					if(m.y>=H-DH-13*36+1 && m.y<=H-DH-12*36-1)//
					{
						InputBox(is4,9,"请输入数据");//读入数据
						setfillcolor(BLACK);
						fillrectangle(W-DW-DC+DC/2+10,H-DH-13*36+1,W-DW,H-DH-12*36);
						outtextxy(W-DW-DC+DC/2+15,H-DH-13*36+15,is4);//写控件字符串
					}//结束if					
				}//结束if(横坐标在输入区域时输入数据)
				sub_frame();//重绘框架
				//B.计算并绘图(横坐标在绘图区域时计算绘图)---------------------------------
				if(m.x>=DW && m.x<=W-DW-DC)
				{
					sub_calculate();//调用计算绘图函数
				}//结束if
				break;//
			}//结束case(结束鼠标左键单击事件消息处理)
		case WM_RBUTTONDOWN://鼠标移动的时候画个空心的圆
			return;//返回while
		}//结束switch(结束鼠标消息)
	}//结束while
}//结束主函数
//绘制框架子函数
void sub_frame()
{
	//1.全屏效果---------------------------------------------------------------------------
	initgraph(W,H);//初始化绘图窗口
	HWND hWnd = GetHWnd();//获取窗口句柄
	LONG style = GetWindowLong(hWnd,-16);//获得窗口风格
	style = style & ~WS_CAPTION & ~WS_SIZEBOX; //窗口全屏显示且不可改变大小
	SetWindowLong(hWnd,-16,style);//设置窗口风格
	SetWindowPos(hWnd, NULL,0,0,W,H,SWP_NOZORDER);//改变窗口位置,尺寸和Z序
	//ShowCursor(FALSE);//显示时隐藏鼠标
	//------------------------------------------------------------------------------------
	//2.绘制工作区边框
	setcolor(c3);
	rectangle(DW,DD,W-DW,H-DH);//外边框
	rectangle(W-DW-DC,DD,W-DW,H-DH);//控件区域边框
	//3.输出标题
	setfont(16,0,"黑体");
	outtextxy(5,10,"数学模型小作业(Malthus_Logistic模型)");
	//4.输出简单说明
	RECT r={W-DW-DC+10,DD+10,W-DW-10,H-DH-10};//
	setfont(12,0,"宋体");//
	drawtext("数学模型小作业:\n　\n \
				Malthus_Logistic模型\
				\
				\n\n\
				操作说明:\n\n\
				鼠标左键  : 数据输入\n\
				鼠标右键  : 退出程序\n\
			\n\
			1~4为输入数据   5~6为计算数据\n\
			\n\n\
			\n\n\
			\
			\n", &r, DT_WORDBREAK);//
	//5.调整说明文字下面的矩形控件
	int c_h=36;//定义控件高度
	int c_x1=W-DW-DC,c_x2=W-DW;//控件左右x坐标
	int c_y1=H-DH;//第1个控件下y坐标
	for(int i=1;i<=16;i++)
	{
		rectangle(c_x1,c_y1,c_x2,c_y1-c_h);//
		c_y1=c_y1-c_h;
	}//结束for
	//6.绘制控件说明文字
	outtextxy(c_x1+10,H-DH-16*c_h+15,"1.数量初值x0:");//输出字符串模拟控件文字
	outtextxy(c_x1+10,H-DH-15*c_h+15,"2.固定增长率r:");//输出字符串模拟控件文字
	outtextxy(c_x1+10,H-DH-14*c_h+15,"3.(Log)最大容量xm:");//输出字符串模拟控件文字
	outtextxy(c_x1+10,H-DH-13*c_h+15,"4.输入时间t:");//输出字符串模拟控件文字
	outtextxy(c_x1+10,H-DH-12*c_h+15,"5.(Mal)模型数值:");//输出字符串模拟控件文字
	outtextxy(c_x1+10,H-DH-11*c_h+15,"6.(Log)模型数值:");//输出字符串模拟控件文字
	//7.初始化原始数据
	outtextxy(W-DW-DC+DC/2+60,H-DH-16*36+15,is1);//写控件字符串
	outtextxy(W-DW-DC+DC/2+60,H-DH-15*36+15,is2);//写控件字符串
	outtextxy(W-DW-DC+DC/2+60,H-DH-14*36+15,is3);//写控件字符串
	outtextxy(W-DW-DC+DC/2+60,H-DH-13*36+15,is4);//写控件字符串
	//8.说明曲线颜色
	//------------------------------------------------------------
    setcolor(cc1);//说明曲线颜色
	line(W-DW-DC+10,H-DH-10*c_h+15,W-DW-DC+0.35*DC,H-DH-10*c_h+15);
	setcolor(cc3);//说明曲线颜色
	line(W-DW-DC+10,H-DH-9*c_h+15,W-DW-DC+0.35*DC,H-DH-9*c_h+15);
	//------------------------------------------------------------
	setcolor(c3);//说明文字颜色
	outtextxy(W-DW-DC+0.4*DC,H-DH-10*c_h+10,"1.Malthus模型数值");//写控件字符串
    outtextxy(W-DW-DC+0.4*DC,H-DH-9*c_h+10,"2.Logistic模型数值");//写控件字符串
	//------------------------------------------------------------
}//结束子程序
//计算数据子函数
void sub_calculate()
{
	//1.原始数据转化为双精度数
	d1=atof(is1);//1.数量初值x0
	d2=atof(is2);//2.固定增长率r
	d3=atof(is3);//3.(Log)最大容量xm
	d4=atof(is4);//4.输入时间t
	//2.根据数据计算要求项目---------------------------------------------------------------------------
	d5=d1*exp(d2*d4);//5.(Mal)模型数值
	d6=(1.0*d3)/(1+(d3*1.0/d1-1.0)*exp(0.0-d2*d4));//6.(Log)模型数值
	//3.将计算结果转化为字符串--------------------------------------------------------------------------
	gcvt(d5,6,is5);//
	gcvt(d6,6,is6);//
	//4.绘制工作区坐标线
	//画x轴刻度标线及文字---------------------------
	double value;//定义双精度变量
	char stt[10];//定义字符数组
	double xstep=(W-2*DW-DC)/(1.0*nx);//x方向步长(nx个刻度)(对应的绘图单位)
	double ystep=(H-DH-DD)/(1.0*ny);//y方向步长(ny个刻度)(对应的绘图单位)
	for(int i=0;i<=nx;i++)
	{
		setcolor(c1);//竖直辅助线颜色
		line(DW+xstep*i,H-DH,DW+xstep*i,DD);//画竖直小刻度线
		setcolor(c3);//文字颜色
		value=TA+xstep*i*(TB-TA)/(W-2*DW-DC);//所在的刻度数值
		gcvt(value,5,stt);//将浮点数value转换为字符串stt
		outtextxy(DW+xstep*i-10,H-DH-20,stt);//在规定的位置放置字符串
	}//结束for
	//画y轴刻度标线及文字---------------------------
	for(i=0;i<=ny;i++)
	{
		setcolor(c2);//水平辅助线颜色
		line(DW,DD+ystep*i,W-DW-DC,DD+ystep*i);//画水平小刻度线
		setcolor(c3);//文字颜色
		value=ystep*(ny-i)*(2.0*d3)/(H-DH-DD);//所在的刻度数值
		gcvt(value,3,stt);//将浮点数value转换为字符串stt
		outtextxy(DW+15,DD+ystep*i+5,stt);//在规定的位置放置字符串
	}//结束for
	//5.重新绘制工作区边框
	setcolor(c3);
	rectangle(DW,DD,W-DW,H-DH);//外边框
	rectangle(W-DW-DC,DD,W-DW,H-DH);//控件区域边框
	//6.将字符串写入控件中
	setfillcolor(BLACK);
	//------------------------------------------------------------
	fillrectangle(W-DW-DC+DC/2+10,H-DH-12*36+1,W-DW,H-DH-11*36);
	outtextxy(W-DW-DC+DC/2+15,H-DH-12*36+15,is5);//写控件字符串
	//------------------------------------------------------------
	fillrectangle(W-DW-DC+DC/2+10,H-DH-11*36+1,W-DW,H-DH-10*36);
	outtextxy(W-DW-DC+DC/2+15,H-DH-11*36+15,is6);//写控件字符串
	//------------------------------------------------------------
	//7.获取时间从(TA~TB)变化时的变化曲线------------------------------------------------
	double x_scale=(W-2*DW-DC)*0.1/(TB-TA);//每个数学单位x对应的屏幕单位
	double y_scale=(H-DH-DD)*1.0/(2.0*d3);//每个数学单位y对应的屏幕单位
	double n1,n2;//
	DPOINT *p1;//定义指向点表的指针
	DPOINT *p2;//定义指向点表的指针
	p1=new DPOINT[1000];//定义N1点表
	p2=new DPOINT[1000];//定义N1点表
	for(int j=0;j<=(TB-TA)*10.0;j++)//时间从TA到TB变化
	{
		n1=d1*exp(d2*(0.1*j+TA));//5.(Mal)模型数值
		n2=(1.0*d3)/(1+(d3*1.0/d1-1.0)*exp(0.0-d2*(0.1*j+TA)));//6.(Log)模型数值
		p1[j].x=DW+j*x_scale;
		p1[j].y=H-DH-n1*y_scale;
		p2[j].x=DW+j*x_scale;
		p2[j].y=H-DH-n2*y_scale;
	}//结束for
	//8.绘制时间从(TA~TB)变化时数值的变化曲线---------------------------------------------------
	setcolor(cc1);//曲线颜色1
	for(j=0;j<(TB-TA)*10.0;j++)//5.(Mal)模型数值
	{
		line(p1[j].x,p1[j].y,p1[j+1].x,p1[j+1].y);
	}//结束for
	setcolor(cc3);//曲线颜色3
	for(j=0;j<(TB-TA)*10.0;j++)//6.(Log)模型数值
	{
		line(p2[j].x,p2[j].y,p2[j+1].x,p2[j+1].y);
	}//结束for
}//结束子函数*/
//--------------------------------------------------------------------------------------------------------------------

gisshow135540019
关注
  • 12
    点赞
  • 踩
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
人口指数Malthus增长模型和Logistic模型（附matlab代码）
05-26
人口指Malthus数增长模型和Logistic模型，美国人口做例子方便理解，还附带代码
分别通过BP神经网络,线性多项式模型,Logistic模型以及Malthus模型进行人口增长模型的预测拟合+代码操作视频
05-30
1.领域：matlab，人口增长模型 2.内容：分别通过BP神经网络,线性多项式模型,Logistic模型以及Malthus模型进行人口增长模型的预测拟合+代码操作视频 3.用处：用于BP神经网络,线性多项式模型,Logistic模型以及Malthus模型学习 4.指向人群：本硕博等教研学习使用 5.运行注意事项： 使用matlab2021a或者更高版本测试，运行里面的Runme_.m文件，不要直接运行子函数文件。运行时注意matlab左侧的当前文件夹窗口必须是当前工程所在路径。 具体可观看提供的操作录像视频跟着操作。
人口指数Malthus增长模型和Logistic模型，附带matlab代码
08-21
人口指Malthus数增长模型和Logistic模型，美国人口做例子方便理解，还附带代码
Malthus和Logistic模型及其医学应用.pdf
06-15
Malthus和Logistic模型及其医学应用.pdf
数据拟合Malthus人口指数增长模型.rar_Malthus_matlab_streamxng_拟合、模型_模型拟合
07-15
数据拟合Malthus人口指数增长模型中的参数,并做出拟合图象
人口模型（Malthus+Logistic）
学习笔记
01-19 7207
Malthus模型 6.15很明显可以知道当x(0) = x0，对x(t)对t求导之后导数为r * x0 * e^rt = r * x(t)，符合6.14的条件。 但是很明显，这样的模型肯定是不符合实际规律的，因为人口不可能随时间无限增长的，需要对模型进行修正。 阻滞增长型（Logistics模型） 即将原来人口增长的常数变为随基数而递减的线性关系，表明当人口达到一定数量时，人口便不会再继续增长 下面展示一道例题 直接建立对应模型即可。 clc, clear % 非线性最小二乘估计 data
人口模型（Malthus+Logistic）附Matlab代码
热门推荐
taoxing
04-12 9万+
目录一.Malthus模型（指数模型）（1）提出以及假设（2）影响人口增长的因素（3）建立模型（4）结论二.Logistic模型（阻滞增长模型）（1）背景（2）建立 r 的关系式（3）模型建立（4）结论 一.Malthus模型（指数模型） （1）提出以及假设 指数增长模型，由马尔萨斯在1798年提出 基本假设：人口（相对）增长率r是常数（r很小） 相对增长率 = 出生率 － 死亡率 （2）影响人口增长的因素 人口的基数 出生率和死亡率 年龄结构 性别比例 工农业生产水平 医疗水平 政府出台的政策
数学建模-人口模型Logistic模型与 Malthus模型
华数云图小分队的博客
11-05 1万+
根据题目中的数据只是一个粗略预报人口总数的模型，没有对人口的年龄、性别、地域结构进行预报，故采用Malthus模型或Logistic模型
数模(5)：Malthus和Logistic模型
weixin_43441742的博客
08-15 8208
Malthus拟合和Logistic拟合是两种最典型、应用最广泛的拟合模型。它们分别被用于预测人口增长和物种增长，因此被广泛应用于社会学与生物学领域。 Malthus拟合 通常情况下，每年的人口死亡率与出生率是相同的。因此可以设每年的人口增长率为常数，我们记时刻t的人口数目为x(t)，记初始时刻的人口数目为x0x_0x0​，则有： $ \begin{cases} \frac{dx}{dt}=rx(...
数学建模-Logistic模型
ChadWick的博客
01-25 1万+
文章目录Malthus模型模型假设建模与求解模型评价Logistic模型模型假设建模与求解模型检验 为了更好地理解Logistic模型，我们先看看Malthus模型 Malthus模型 这是英国神父Malthus通过对一百多年人口统计资料的分析之后提出的人口模型假设 模型假设 设x(t)表示t时刻的人口数，且x(t)连续可微. 人口的增长率r是常数(增长数=出生率-死亡率) 人口的数量的变化是封闭的，即人口数量的增加与减少只取决于人口中个体的生育和死亡，且每一个体都具有同样的生育能力与死亡率.
一篇文章带你搞定数学建模中的 Logistic 模型（含代码）
南淮北安的博客
07-19 6637
文章目录一、中国人口情况二、阻滞增长模型三、美国人口预测四、人口模型的概述 一、中国人口情况 二、阻滞增长模型 适用：前期增长迅速，后期增长迟缓 三、美国人口预测 已经学习过：一篇文章带你搞定数学建模中的 Malthus人口模型 知道了 Malthus 模型不适用于长期的预判，所以这里我们使用 Logistic 模型重新拟合美国人口预测 所以这里需要使用 非线性拟合估计模型参数 r 和 Xm % 时间周期 t，这里空出来 2000，用来和预测结果比较 t = [1790 1800 1810 182
人口logistic模型公式_人口指数增长模型和Logistic模型
weixin_39968852的博客
02-22 3569
根据美国人口从1790年到1990年间的人口数据(如下表)，确定人口指数增长模型和Logistic 模型中的待定参数，估计出美国2010年的人口，同时画出拟合效果的图形。1860 1870 1880 1890 1900 1910 1930 1940 1950 1960 1970 1980 指数增长模型：rte x t x 0)(=Logistic 模型：()011mrtm x x t x e x ...
数学建模之微分方程模型详解
左手の明天的博客
03-11 1万+
微分方程知识简介 要掌握常微分方程的一些基础知识，对一些可以求解的微分方程及其方程组，要求掌握其解法，并了解一些方程的近似解法。 微分方程的体系 (1)初等积分法（一阶方程及几类可降阶为一阶的方程） (2)一阶线性微分方程组（常系数线性微分方程组的解法） (3)高阶线性微分方程（高阶线性常系数微分方程解法）。其中还包括了常微分方程的基本定理。 0．常数变易法 常数变易法在上面的（1）（2）（3）三部分中都出现过，它是由线性齐次方程（一阶或高阶）或方程组的解经常数变易后求相应的非齐次方程或方程组的
logistic人口模型python代码_人口模型(马尔萨斯--vs--logistic).ppt
weixin_39626586的博客
12-10 792
人口模型;微分方程模型;求出方程的解 ——求出未知函数的解析表达式 ——利用各种数值解法、数值软件(如Matlab)求近似解不必求出方程的解 ——根据微分方程的理论研究某些性质，或它的变化趋势; 为了保持自然资料的合理开发与利用，人类必须保持并控制生态平衡，甚至必须控制人类自身的增长。本节将建立几个简单的单种群增长模型，以简略分析一下这方面的问题。一般生态系统的...
Java 7大排序
m0_74270127的博客
05-10 810
插入排序、希尔排序、选择排序、堆排序、冒泡排序、快速排序、归并排序
1080：余数相同问题
如果没有特别说明，c++都是基于GCC-4.9.2运行的
05-11 374
【代码】：1080：余数相同问题
找第二大的数（分治法）
2201_75525567的博客
05-11 298
这段代码通过递归的方式，将数组分为两半，并分别找出左半部分和右半部分的两个最大数，然后比较这四个数的大小，得到第二大的数。这段代码是用来找出数组中第二大的数的。
华为OD机试 - 求幸存数之和（Java 2024 C卷 100分）
最新发布
学Java，找哪吒
05-13 754
每一题都有详细的答题思路、详细的代码注释、样例测试，发现新题目，随时更新，全天CSDN在线答疑。
logistic模型人口预测
07-29
Logistic模型是一种用于人口预测的数学模型，它考虑到了人口增长会受到一些限制因素的影响。根据Logistic模型，人口的自然增长率是一个减函数，随着人口数量的增加逐渐下降。当人口数量达到一个上限值K时，自然增长...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值