- 在线时间
- 1630 小时
- 最后登录
- 2024-1-29
- 注册时间
- 2017-5-16
- 听众数
- 82
- 收听数
- 1
- 能力
- 120 分
- 体力
- 565610 点
- 威望
- 12 点
- 阅读权限
- 255
- 积分
- 174906
- 相册
- 1
- 日志
- 0
- 记录
- 0
- 帖子
- 5313
- 主题
- 5273
- 精华
- 3
- 分享
- 0
- 好友
- 163
TA的每日心情 | 开心 2021-8-11 17:59 |
|---|
签到天数: 17 天 [LV.4]偶尔看看III 网络挑战赛参赛者 网络挑战赛参赛者 - 自我介绍
- 本人女,毕业于内蒙古科技大学,担任文职专业,毕业专业英语。
 群组: 2018美赛大象算法课程 群组: 2018美赛护航培训课程 群组: 2019年 数学中国站长建 群组: 2019年数据分析师课程 群组: 2018年大象老师国赛优 |
中国大学生数学建模竞赛备赛(十三)
- V: `3 j$ N6 r# O8 j2 A+ o+ D微分方程问题 o. t) f: @" s; d
微分方程建模. B1 B3 ]" t( b: V) C7 p* K
1、根据实际要求确定要研究的量。
- o8 ~( }2 j5 Z1 q2、找出这些量所满足的基本规律。
5 C/ @& V; ^9 q2 Y3、运用这些规律列出方程和定解条件。
- v; P& o9 n R S常见的列方程方法:- d. t$ f, M# ^: b$ I
(1)依照规律直接列;(2)微元分析法与任意区域上取积分的方法;(3)模拟分析法。/ O/ U. S! ]& n: O+ J( N$ g
h: A- `, [! g- c; l
" z* E5 J$ t7 c, |几类微分方程的应用实例
, m# g$ i5 {, X) H/ o1、发射火箭使用三级推进器。(P103-P107)8 T* Z4 s7 E& h, o) c
2、人口模型:(1)Malthus模型;(2)阻滞增长模型(Logistic模型)。(P107-P110); }2 X" R' m9 r( C& X$ C M; H
例题:! z' g, ~" \: @2 j+ y6 |
) u- L; ^& |6 _$ e
1 R- w/ s: \: M! W
第一种方法:非线性最小二乘估计,也可以称之为微分方程反问题的求解: f4 S7 R9 x9 O) {$ x
" z& L) L, h+ H; V/ v8 Z6 D z0 N8 z6 W+ n) o; H) H5 V
clc,clear( p# \4 q: B- ~5 |- N
a=textread('data4.txt');//把原始的数据保存在纯文本文件data4.txt中
, g; b# o( h# j$ [x=a([2:2:6], ';//提出需要的人口数据/ C& ]; x( A! V: h8 A, c2 q: }
x=nonzeros(x);//去掉后面的0,变成列向量
3 b3 Q* E2 Z$ T( P' xt=[1790:10:2000]';/ d' E4 J) X, D2 |; J
t0=t(1);x0=x(1);7 t3 b, n+ F1 S6 X5 b, X# Q6 Z
fun=@(cs,td)cs(1)./(1+(cs(1)/x0-1)*exp(-cs(2)*(td-t0)));//cs(1)=xm,cs(2)=r2 m0 W# [2 v8 u4 @5 @3 ^! n
cs=lsqcurvefit(fun,rand(2,1),t(2:end),x(2:end),zeros(2,1));//后向差分4 _1 K7 t( p) Q7 v7 q. [/ E# |
xhat=fun(cs,[t,2010]);
+ b4 @2 e! ?0 R5 d$ @+ ?7 o" Z1! N, n- B) H. u
20 W" a* ~( ^0 d9 }! u
32 ?* a* z3 v. ]1 Q& L
4
& o4 n Y! o" [; P5 i3 c8 W8 I5
7 H h# _. P2 _1 g9 n8 v% P% j- e6
5 C. B- O \+ ~7
9 {! M: [1 C n) l! P! T% b: ^' V8$ w* @1 P+ `% e
90 m& p3 b; a3 |9 Y$ Z( Z3 f( [
第二种方法:线性最小二乘法(参考之前写的最小二乘拟合原理)
; T7 t: o* W2 l6 a' v2 I' j) F& m4 U; B( n5 O
7 k! C. E& E# _. G8 H
3 o1 R( r- R2 } _8 x4 n) k
- s7 E8 M' {0 a, g d$ zclc,clear
# p9 @3 h% U% h0 T7 ]4 U1 fa=textread('data4.txt');! e4 \* c. s6 a) I
x=a([2:2:6], ';x=nonzeros(x);
E* B& X* }* k7 N2 |2 B$ n- i3 zt=[1790:10:2000]';
$ k% V% p$ p5 {% D7 y- da=[ones(21,1),-x(2:end)];5 d5 _- D0 G O' U9 @9 o, f
b=diff(x)./x(2:end)/10;4 E/ D& p: D3 Y0 \' l, k
cs=a\b;5 b& \* t. E1 \3 B1 N: H
r=cs(1);xm=r/cs(2);
$ O0 e+ j) a! ?5 u/ |. z* C1. c/ F( J# ^& F p3 ^
2
: J/ D# @% v! }: \5 q/ l+ ~3
" a7 l, t5 F: }* A4 y: t4
9 X9 V; f# _) ?" ~% ^9 o) m4 d5
! M9 ]2 l0 G! S' r6
# m1 `! J" m% F1 b! L, E7
& x8 @6 @* E9 a9 f' Y4 D9 T8
( l" Q; H8 R. v
* F/ T5 w1 J$ {. x; Y# p" H$ z: P
! m! l M- p0 Q- R
- S. v7 V' h4 e" J3 m: N' C9 D, _5 L* g/ l6 p. c& \* Q% C
clc,clear2 M3 D$ Z6 o$ a6 h
a=textread('data4.txt');8 o$ E* i8 k m5 e! H! w# P
x=a([2:2:6], ';x=nonzeros(x);
# M+ o* n. ~, E& L8 H2 g9 ~t=[1790:10:2000]';
' E$ Y, G+ m* s5 ?' V7 O+ o% Z$ Z# r9 Qa=[ones(21,1),-x(1:end-1)];$ E# D7 f, z. {$ G' ]7 Q3 b
b=diff(x)./x(1:end-1)/10;7 ~) i( |, X. {" K4 k; y# F2 W
cs=a\b;+ t; ^$ l/ [4 O6 l7 t
r=cs(1);xm=r/cs(2);
% ~* m) |- f9 Y+ v1 F- W" z1
8 h+ X+ i3 O& Q' f! c; k22 f! u' x( r% D
34 Y! ~9 q' B( p \0 R6 C5 y7 @
4
: x" D5 s" ^$ L4 `& B. I8 v5
2 ^8 ~, j$ j% Y/ m6
$ S/ T0 v" N: p. M3 k7
, i9 `. o! |# H+ N% @8; @* N p; X3 h% Q8 O1 V5 p
参考文献7 k' l2 Z# W* E/ Z; \
司守奎,孙玺菁. 数学建模算法与应用. 北京:国防工业出版社,2011.& x1 Z! u! J1 T. u+ t; j2 r- |
————————————————% D5 m% }/ k H e. p+ Q
版权声明:本文为CSDN博主「小白成长之旅」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
* l% v: |( c* R) X原文链接:https://blog.csdn.net/weixin_45813658/article/details/107780619, ~+ V4 m9 t6 r$ A2 ?6 ~% N
1 `( a% i! n9 B! ` x* Z- \& N6 w8 e" {7 R" D1 m
|
zan
|