我是靠谱客的博主 传统大侠,这篇文章主要介绍数值分析思考题(钟尔杰版)参考解答——第八章题目:1. 人口模型中马尔萨斯模型与逻辑斯模型有何区别?2. 牛顿谐振动和小阻尼振动的微分方程之间有何区别和联系?3. 单摆的常微分方程如何求近似解?4. 求一阶常微分方程数值解的欧拉法与平面向量场图形有何联系?5. 二阶龙格-库塔法和欧拉法有何联系?6. 一阶常微分方程组和一阶常微分方程在求数值解时用龙格-库塔方法有何区别?7. 求二阶以上的常微分方程数值解时如何应用龙格-库塔方法?8. 蝴蝶效应的洛伦兹方程有何特点?,现在分享给大家,希望可以做个参考。

数值分析思考题(钟尔杰版)参考解答——第八章

  • 题目:
  • 1. 人口模型中马尔萨斯模型与逻辑斯模型有何区别?
  • 2. 牛顿谐振动和小阻尼振动的微分方程之间有何区别和联系?
  • 3. 单摆的常微分方程如何求近似解?
  • 4. 求一阶常微分方程数值解的欧拉法与平面向量场图形有何联系?
  • 5. 二阶龙格-库塔法和欧拉法有何联系?
  • 6. 一阶常微分方程组和一阶常微分方程在求数值解时用龙格-库塔方法有何区别?
  • 7. 求二阶以上的常微分方程数值解时如何应用龙格-库塔方法?
  • 8. 蝴蝶效应的洛伦兹方程有何特点?

题目:


1. 人口模型中马尔萨斯模型与逻辑斯模型有何区别?


马尔萨斯的人口论指出:在没有生存资源限制的情况下,人口或生物种群的数量成指数增长。人口增长的数学模型的共同特征是人口数量在单位时间内增长的百分比 r 是一定的。写成一个微分方程的形式,设 t = 0 时刻的人口数量为 N 0 N_0 N0​,则 t mathrm{t} t 时刻的总人口 N t N_t Nt​ 满足 1 N t d N t d t = r ⇒ N t = N 0 e r t frac{1}{N_t} frac{d N_t}{d t}=r Rightarrow N_t=N_0 e^{r t} Nt​1​dtdNt​​=r⇒Nt​=N0​ert
在资源环境的限制下, 假设资源环境能承受的人口数量为 K K K, 则可以建立一个逻辑斯谛方程: 1 N t d N t d t = r ( 1 − N t K ) frac{1}{N_t} frac{d N_t}{d t}=rleft(1-frac{N_t}{K}right) quad Nt​1​dtdNt​​=r(1−KNt​​) 这个方程将得出仅在人口 N t ≪ K N_t ll K Nt​≪K 时 ,增长 N t ∼ N 0 e r t N_t sim N_0 e^{r t} Nt​∼N0​ert 才是指数的。当 N t N_t Nt​ 接近 K mathrm{K} K 时, 人口的增长明显受到天花板 K mathrm{K} K 的 压制。虽然 N t > K N_t>K Nt​>K 时确实会出现人口的负增长, 但不过是平稳地趋近于平衡水 平 K mathrm{K} K, 而并不会发生马尔萨斯所担心的灾难性人口锐减的行为。


2. 牛顿谐振动和小阻尼振动的微分方程之间有何区别和联系?

牛顿谐振动的微分方程:
d 2 u d t 2 + ω 2 u ( t ) = 0 ⇒ u ( t ) = C 1 cos ⁡ ω t + C 2 sin ⁡ ω t frac{d^2 u}{d t^2}+omega^2 u(t)=0 Rightarrow u(t)=C_1 cos omega t+C_2 sin omega t dt2d2u​+ω2u(t)=0⇒u(t)=C1​cosωt+C2​sinωt 小阻尼振动的微分方程: d 2 u d t 2 + 2 ε d u d t + ω 2 u ( t ) = 0 frac{d^2 u}{d t^2}+2 varepsilon frac{d u}{d t}+omega^2 u(t)=0 dt2d2u​+2εdtdu​+ω2u(t)=0, 设解函数有如下格式 u ( t ) = e ( − ε t ) sin ⁡ t u(t)=e^{(-varepsilon t)} sin t u(t)=e(−εt)sint 代入方程化简得:
d 2 v d t 2 + ( ω 2 − ε 2 ) v = 0 ⇒ frac{d^2 v}{d t^2}+left(omega^2-varepsilon^2right) v=0 Rightarrow dt2d2v​+(ω2−ε2)v=0⇒
v ( t ) = C 1 cos ⁡ k t + C 2 sin ⁡ k t ( ω 2 − ε 2 = k 2 ) mathrm{v}(mathrm{t})=C_1 cos k t+C_2 sin k t quadleft(omega^2-varepsilon^2=k^2right) v(t)=C1​coskt+C2​sinkt(ω2−ε2=k2) 可见小阻尼振动的微分方程和牛顿谐振动的微分方程的形式又可以统一起来。


3. 单摆的常微分方程如何求近似解?

θ ′ ′ = − a sin ⁡ θ  令  y 1 = θ , y 2 = θ ′ { y 2 = y 1 ′ y 2 ′ = − a sin ⁡ y 1 ⇒ { y 2 ( 0 ) = 1 y 2 ′ = − asin ⁡ y 1 ⇒ y n + 1 = y n + h ( − a sin ⁡ y 1 n ) n = 0 , 1 , 2 … . . N − 1 begin{aligned} &theta^{prime prime}=-a sin theta quad text { 令 } y_1=theta, y_2=theta^{prime} \ &left{begin{array} { c } { y _ { 2 } = y _ { 1 } ^ { prime } } \ { y _ { 2 } ^ { prime } = - a operatorname { s i n } y _ { 1 } } end{array} Rightarrow left{begin{array}{c} y_2(0)=1 \ y_2^{prime}=-operatorname{asin} y_1 end{array} Rightarrow y_{n+1}=y_n+hleft(-a sin y_{1_n}right)right.right. \ &n=0,1,2 ldots . . N-1 end{aligned} ​θ′′=−asinθ 令 y1​=θ,y2​=θ′{y2​=y1′​y2′​=−asiny1​​⇒{y2​(0)=1y2′​=−asiny1​​⇒yn+1​=yn​+h(−asiny1n​​)n=0,1,2…..N−1​


4. 求一阶常微分方程数值解的欧拉法与平面向量场图形有何联系?

一阶常微分方程解的几何意义是一族曲线。而满足初始条件 y ( x 0 ) = y 0 yleft(x_0right)=y_0 y(x0​)=y0​ 的 解 y ( x ) y(x) y(x) 对应于其中一条曲线。函数 f ( x , y ) f(x, y) f(x,y) 值, 恰好是曲线 y ( x ) y(x) y(x) 在 x mathrm{x} x 处切线的斜率。曲线族中每一条曲线的切线斜率在平面上确定的方向构成平面区域内的向量场。欧拉法的几何意义是用折线段 P 0 P 1 P 2 … P_0 P_1 P_2 ldots P0​P1​P2​… 来近似代替方程:

{ d y d x = f ( x , y ) , a ≤ x ≤ b y ( a ) = y 0 left{begin{array}{c} frac{d y}{d x}=f(x, y), a leq x leq b \ y(a)=y_0 end{array}right. {dxdy​=f(x,y),a≤x≤by(a)=y0​​
的解曲线???? = ????(????),在每一小段上
y n + 1 − y n x n + 1 − x n = f ( x n , y n ) frac{y_{n+1}-y_n}{x_{n+1}-x_n}=fleft(x_n, y_nright) xn+1​−xn​yn+1​−yn​​=f(xn​,yn​)

即: y n + 1 = y n + h f ( x n , y n ) y_{n+1}=y_n+h fleft(x_n, y_nright) yn+1​=yn​+hf(xn​,yn​)


5. 二阶龙格-库塔法和欧拉法有何联系?

对于 r = 2 r=2 r=2 的 R-K 方法, 可得如下计算公式:
{ y n + 1 = y n + h ( c 1 K 1 + c 2 K 2 ) K 1 = f ( x n , y n ) K 2 = f ( x n + λ 2 h , y n + μ 21 h K 1 ) left{begin{array}{c} y_{n+1}=y_n+hleft(c_1 K_1+c_2 K_2right) \ K_1=fleft(x_n, y_nright) \ K_2=fleft(x_n+lambda_2 h, y_n+mu_{21} h K_1right) end{array}right. ⎩ ⎨ ⎧​yn+1​=yn​+h(c1​K1​+c2​K2​)K1​=f(xn​,yn​)K2​=f(xn​+λ2​h,yn​+μ21​hK1​)​
取 c 1 = c 2 = 1 2 , λ 2 = μ 21 = 1 c_1=c_2=frac{1}{2}, lambda_2=mu_{21}=1 c1​=c2​=21​,λ2​=μ21​=1, 得出改进形式的欧拉法
y n + 1 = y n + h 2 [ f ( x n , y n ) + f ( x n + h , y n + h f ( x n , y n ) ) ] y_{n+1}=y_n+frac{h}{2}left[fleft(x_n, y_nright)+fleft(x_n+h, y_n+h fleft(x_n, y_nright)right)right] yn+1​=yn​+2h​[f(xn​,yn​)+f(xn​+h,yn​+hf(xn​,yn​))]


6. 一阶常微分方程组和一阶常微分方程在求数值解时用龙格-库塔方法有何区别?

  • 求解一阶常微分方程组时,只要把方程组的 y 与 f 理解为向量,就可以代入一阶方程组的龙格-库塔法进行求解。

7. 求二阶以上的常微分方程数值解时如何应用龙格-库塔方法?

求解高阶微分方程的初值问题, 可以将其化为一阶方程组进行求解, 对 m m m 阶微分方程 y ( m ) = f ( x , y , y ′ , … , y ( m − 1 ) ) y^{(m)}=fleft(x, y, y^{prime}, ldots, y^{(m-1)}right) y(m)=f(x,y,y′,…,y(m−1)), 只要引进新的变量 y 1 = y , y 2 = y_1=y, y_2= y1​=y,y2​= y ′ , … , y m = y ( m − 1 ) y^{prime}, ldots, y_m=y^{(m-1)} y′,…,ym​=y(m−1), 即可将 m m m 阶微分方程化为一阶微分方程组
{ y 1 ′ = y 2 y 2 ′ = y 3 ⋮ y m − 1 ′ = y m y m ′ = f ( x , y 1 , y 2 , … , y m ) left{begin{array}{c} y_1^{prime}=y_2 \ y_2^{prime}=y_3 \ vdots \ y_{m-1}^{prime}=y_m \ y_m^{prime}=fleft(x, y_1, y_2, ldots, y_mright) end{array}right. ⎩ ⎨ ⎧​y1′​=y2​y2′​=y3​⋮ym−1′​=ym​ym′​=f(x,y1​,y2​,…,ym​)​
然后利用求解一阶常微分方程组的龙格-库塔方法进行求解。


8. 蝴蝶效应的洛伦兹方程有何特点?

  • 蝴蝶效应的洛伦兹方程,每个方程都是一阶常微分方程, 每个分量初始条件下微小的变化能带动整个系统的长期的巨大的连锁反应,使结果很不稳定。

最后

以上就是传统大侠最近收集整理的关于数值分析思考题(钟尔杰版)参考解答——第八章题目:1. 人口模型中马尔萨斯模型与逻辑斯模型有何区别?2. 牛顿谐振动和小阻尼振动的微分方程之间有何区别和联系?3. 单摆的常微分方程如何求近似解?4. 求一阶常微分方程数值解的欧拉法与平面向量场图形有何联系?5. 二阶龙格-库塔法和欧拉法有何联系?6. 一阶常微分方程组和一阶常微分方程在求数值解时用龙格-库塔方法有何区别?7. 求二阶以上的常微分方程数值解时如何应用龙格-库塔方法?8. 蝴蝶效应的洛伦兹方程有何特点?的全部内容,更多相关数值分析思考题(钟尔杰版)参考解答——第八章题目:1.内容请搜索靠谱客的其他文章。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(281)

相关文章

评论列表共有 0 条评论

立即
投稿
返回
顶部