- 在线时间
- 1630 小时
- 最后登录
- 2024-1-29
- 注册时间
- 2017-5-16
- 听众数
- 82
- 收听数
- 1
- 能力
- 120 分
- 体力
- 563322 点
- 威望
- 12 点
- 阅读权限
- 255
- 积分
- 174219
- 相册
- 1
- 日志
- 0
- 记录
- 0
- 帖子
- 5313
- 主题
- 5273
- 精华
- 3
- 分享
- 0
- 好友
- 163
TA的每日心情 | 开心 2021-8-11 17:59 |
|---|
签到天数: 17 天 [LV.4]偶尔看看III 网络挑战赛参赛者 网络挑战赛参赛者 - 自我介绍
- 本人女,毕业于内蒙古科技大学,担任文职专业,毕业专业英语。
 群组: 2018美赛大象算法课程 群组: 2018美赛护航培训课程 群组: 2019年 数学中国站长建 群组: 2019年数据分析师课程 群组: 2018年大象老师国赛优 |
中国大学生数学建模竞赛备赛(十三)4 ^9 V- P. r" t8 l o5 y
微分方程问题
$ ]4 Z- g/ P0 [/ t/ b/ j+ `+ n- n% i( c微分方程建模
6 U& Z3 ]% h7 p8 [* [ y h& z, c1、根据实际要求确定要研究的量。. P7 r# o) q- m+ K7 f% a9 d
2、找出这些量所满足的基本规律。
6 |# e' G$ b/ x6 ~! o$ A3、运用这些规律列出方程和定解条件。, u0 L2 q0 \/ J% _
常见的列方程方法:+ h `- E. l M; ^
(1)依照规律直接列;(2)微元分析法与任意区域上取积分的方法;(3)模拟分析法。* l; D) h6 k) j. ]3 u( q
. f# i. k; }8 q% ^ ^
7 u7 J, ?& T: {1 G+ Z
几类微分方程的应用实例 B. P4 Q+ [& s6 S8 n# ^" Q
1、发射火箭使用三级推进器。(P103-P107). T& S+ H+ j- N5 @: ~2 f
2、人口模型:(1)Malthus模型;(2)阻滞增长模型(Logistic模型)。(P107-P110)
, H7 E1 S9 T' S+ U y9 p' M0 r例题:5 v" f$ y2 s; t! i3 W
( I; w1 ` u# [ Q& E. U/ e
- n: |3 X6 i% a( a第一种方法:非线性最小二乘估计,也可以称之为微分方程反问题的求解
+ q9 W$ A2 J$ N' j
1 z: n {9 T* L5 F0 ~; F* n/ Q0 W5 G; Y- {
clc,clear' A4 H" {6 q, a; D8 W
a=textread('data4.txt');//把原始的数据保存在纯文本文件data4.txt中
0 Y$ B; v9 o4 z3 o7 J4 lx=a([2:2:6], ';//提出需要的人口数据
5 e' p5 g/ s; b2 Z8 d& bx=nonzeros(x);//去掉后面的0,变成列向量$ u9 J/ x1 H* A# V5 F
t=[1790:10:2000]';
p; W% d+ X, G' Ft0=t(1);x0=x(1);
G9 F. `; Q2 X5 G1 S: ]; Sfun=@(cs,td)cs(1)./(1+(cs(1)/x0-1)*exp(-cs(2)*(td-t0)));//cs(1)=xm,cs(2)=r
# T N/ k; v$ c! h& Q7 h; tcs=lsqcurvefit(fun,rand(2,1),t(2:end),x(2:end),zeros(2,1));//后向差分
' T( `4 U7 @7 M! T& `xhat=fun(cs,[t,2010]);
0 h( j7 E4 _4 E- {* s% }1
( } n, w9 v" I0 w! r2 N$ r2& \" L# ?2 z; \) Q: F) w3 ?4 Y; ^% A
3 }/ g! I ^- I$ s/ |! c
4
+ L5 L; C( ^- ~7 C5: V1 t# w+ `- I- ]! Z
6
# g( Q( X$ }) M70 v0 M) l) {% s+ X4 H
8* t8 C" I; r5 c- ^1 g( u
9
! |6 Y1 v3 L; y+ L" i第二种方法:线性最小二乘法(参考之前写的最小二乘拟合原理)' L8 ]8 G7 ^/ U! O
- B& H1 o, B/ `1 s. `
& k- r( _5 O8 |5 H; T& S6 e; T' h) Y
1 S$ R9 h# t9 Y& o
clc,clear
# r' n0 H9 o( L; X" D3 xa=textread('data4.txt');9 q' z1 U7 n% g7 {
x=a([2:2:6], ';x=nonzeros(x);
0 f" ?' c2 ~- Z- ?t=[1790:10:2000]';
. p5 m& p) f3 K$ n5 D$ F+ @7 fa=[ones(21,1),-x(2:end)];
2 l' r" ]2 G4 Q$ a$ m9 W7 a6 Wb=diff(x)./x(2:end)/10;
7 Q. u$ E3 b9 s5 G; ucs=a\b;! ]- V s7 |9 w8 L
r=cs(1);xm=r/cs(2);
5 G9 ^$ C: B/ F2 {/ C. s+ O8 [8 \1
$ N( }3 g8 c) W2 u2
?( ^: G9 U2 j/ ~ ~3
% G. W9 D: G' }" L8 P44 x" m7 B6 T v
53 q) |1 J% K- m# J2 C6 o/ c
6* M: \! F, `! x) S
7* e! N! e- j# g# o
82 e. _4 o3 O0 L7 b
' U2 _( f6 R9 [% j
8 h9 B1 W6 C1 l$ H- r% `5 s0 o
; W( \! b' b8 m* s1 S
: S. r2 O5 k* t) S/ z) pclc,clear
5 c" M2 _9 o0 O- p; o! [" aa=textread('data4.txt');# f: R* P% k, i0 Z7 F2 @) M1 b
x=a([2:2:6], ';x=nonzeros(x);9 B4 W0 B2 k* y0 K8 }
t=[1790:10:2000]';
9 _" U3 {) L7 h0 m1 h$ pa=[ones(21,1),-x(1:end-1)];
G3 {5 x: Y4 |: W- S; [( j. |b=diff(x)./x(1:end-1)/10;9 \/ e/ @7 c; j( m7 Q
cs=a\b;( u' q& m% x0 h- t4 [8 ^7 X
r=cs(1);xm=r/cs(2);
) I1 f$ m; G M2 r$ d$ m1
|* T# f; L8 g! W2. c) J2 l( I4 L4 I. _
30 X# g; x8 T3 z ]( g; r
4
; @4 { d) c" V! n5
; Y% y( h0 @. k& ]4 ?6
) F; V6 O8 Y# I4 J, u( p3 _7 y. `. A7& t4 x* F8 J! ^: A: q( \) e" W
8
9 ~9 R7 l& I& @9 K, H: T参考文献; E' C& v" n* ^
司守奎,孙玺菁. 数学建模算法与应用. 北京:国防工业出版社,2011./ l6 h8 h# m9 F/ H. ~
————————————————
3 ~" }8 y" M' X版权声明:本文为CSDN博主「小白成长之旅」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。# [" a9 v, z8 E' k
原文链接:https://blog.csdn.net/weixin_45813658/article/details/107780619/ U+ r" u( D# Q
6 F# u3 V6 t+ a/ p0 L! }" g* u$ i: f$ h8 c8 G
|
zan
|