查看: 3164|回复: 0

中国大学生数学建模竞赛备赛（十三）

字体大小: 正常放大

杨利霞

5273 主题	82 听众	17万积分

TA的每日心情

	开心 2021-8-11 17:59

签到天数: 17 天

[LV.4]偶尔看看III

网络挑战赛参赛者

自我介绍: 本人女，毕业于内蒙古科技大学，担任文职专业，毕业专业英语。

群组: 2018美赛大象算法课程

群组: 2018美赛护航培训课程

群组: 2019年数学中国站长建

群组: 2019年数据分析师课程

群组: 2018年大象老师国赛优

电梯直达

1^#

发表于 2021-5-24 15:59 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

中国大学生数学建模竞赛备赛（十三）
微分方程问题
微分方程建模
1、根据实际要求确定要研究的量。
2、找出这些量所满足的基本规律。
3、运用这些规律列出方程和定解条件。
常见的列方程方法：
（1）依照规律直接列；（2）微元分析法与任意区域上取积分的方法；（3）模拟分析法。

几类微分方程的应用实例
1、发射火箭使用三级推进器。（P103-P107）
2、人口模型：（1）Malthus模型；（2）阻滞增长模型（Logistic模型）。（P107-P110）
例题：

第一种方法：非线性最小二乘估计，也可以称之为微分方程反问题的求解

clc,clear
a=textread('data4.txt');//把原始的数据保存在纯文本文件data4.txt中
x=a([2:2:6],

';//提出需要的人口数据
x=nonzeros(x);//去掉后面的0，变成列向量
t=[1790:10:2000]';
t0=t(1);x0=x(1);
fun=@(cs,td)cs(1)./(1+(cs(1)/x0-1)*exp(-cs(2)*(td-t0)));//cs(1)=xm,cs(2)=r
cs=lsqcurvefit(fun,rand(2,1),t(2:end),x(2:end),zeros(2,1));//后向差分
xhat=fun(cs,[t,2010]);
1
2
3
4
5
6
7
8
9
第二种方法：线性最小二乘法（参考之前写的最小二乘拟合原理）

clc,clear
a=textread('data4.txt');
x=a([2:2:6],

';x=nonzeros(x);
t=[1790:10:2000]';
a=[ones(21,1),-x(2:end)];
b=diff(x)./x(2:end)/10;
cs=a\b;
r=cs(1);xm=r/cs(2);
1
2
3
4
5
6
7
8

clc,clear
a=textread('data4.txt');
x=a([2:2:6],

';x=nonzeros(x);
t=[1790:10:2000]';
a=[ones(21,1),-x(1:end-1)];
b=diff(x)./x(1:end-1)/10;
cs=a\b;
r=cs(1);xm=r/cs(2);
1
2
3
4
5
6
7
8
参考文献
司守奎，孙玺菁. 数学建模算法与应用. 北京：国防工业出版社，2011.
————————————————
版权声明：本文为CSDN博主「小白成长之旅」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/weixin_45813658/article/details/107780619

风控大赛.jpg (209.44 KB, 下载次数: 0)