非线性ode及在工程中的应用第五次作业更正版
第一题:
(1)
题式的平衡方程是:
u
˙
=
p
−
u
2
=
0
\dot u=p-u^2=0
u˙=p−u2=0
解得
u
=
±
p
(
p
>
0
)
u=\pm \sqrt p (p\gt0)
u=±p(p>0)以及
u
在
p
≤
0
时
无
解
u在p\le0时无解
u在p≤0时无解
令 f u ( u , p ) = − 2 u = 0 f_u(u,p)=-2u=0 fu(u,p)=−2u=0,和原平衡方程联立求解可得,奇点为 ( 0 , 0 ) (0,0) (0,0)
将这个点带入 f p ( u , p ) = 1 f_p(u,p)=1 fp(u,p)=1中可以得到, f p f_p fp始终大于零, ( 0 , 0 ) (0,0) (0,0)为转折点(也称为鞍节点)。
应该紧接着有个关于分叉之后稳定性的讨论吧,这个不是很熟,如果有朋友能教教我最好了,实在没有我慢慢思考一下想通了再来补充这部分内容,看这题题干暂时不要求,先把作业写完再说吧。
(2)
题式的平衡方程是:
u
˙
=
p
u
−
u
2
+
u
3
=
0
\dot u=pu-u^2+u^3=0
u˙=pu−u2+u3=0
解得
u
=
0
u=0
u=0或者
u
=
1
/
2
+
(
1
/
2
)
1
−
4
p
u=1/2+(1/2)\sqrt{1-4p}
u=1/2+(1/2)1−4p或者
u
=
1
/
2
−
(
1
/
2
)
1
−
4
p
u=1/2-(1/2)\sqrt{1-4p}
u=1/2−(1/2)1−4p当然对于后面两个解
p
≤
1
/
4
p\le1/4
p≤1/4
令 f u ( u , p ) = p − 2 u + 3 u 2 = 0 f_u(u,p)=p-2u+3u^2=0 fu(u,p)=p−2u+3u2=0,和原平衡方程联立求解可得,奇点为 ( 0 , 0 ) 和 ( 1 / 2 , 1 / 4 ) (0,0)和(1/2,1/4) (0,0)和(1/2,1/4)
将这两个点分别代入 f p ( u , p ) = u f_p(u,p)=u fp(u,p)=u可以得到, ( 0 , 0 ) (0,0) (0,0)为分叉点,而 ( 1 / 2 , 1 / 4 ) (1/2,1/4) (1/2,1/4)为鞍节点
(3)
题式的平衡方程是:
u
˙
=
p
2
u
−
u
3
=
0
\dot u=p^2u-u^3=0
u˙=p2u−u3=0
解得
u
=
0
u=0
u=0或者
u
=
p
u=p
u=p或者
u
=
−
p
u=-p
u=−p
令 f u ( u , p ) = p − 2 u + 3 u 2 = 0 f_u(u,p)=p-2u+3u^2=0 fu(u,p)=p−2u+3u2=0,和原平衡方程联立求解可得,奇点为 ( 0 , 0 ) 和 ( 1 / 3 , 1 / 3 ) (0,0)和(1/3,1/3) (0,0)和(1/3,1/3)
将这两个点分别代入 f p ( u , p ) = 2 p u f_p(u,p)=2pu fp(u,p)=2pu可以得到, ( 0 , 0 ) (0,0) (0,0)为分叉点,而 ( 1 / 3 , 1 / 3 ) (1/3,1/3) (1/3,1/3)为鞍节点
稳定性的部分,实在是不能理解,麻烦读到这里而且理解的朋友和我简单聊聊,可以吗~
第二题:
首先补充一点这个定理相关的知识:
至
于
里
面
的
这
么
多
参
数
至于里面的这么多参数
至于里面的这么多参数
R
e
是
一
个
虚
数
的
实
部
,
而
I
m
为
一
个
虚
数
的
虚
部
Re是一个虚数的实部,而Im为一个虚数的虚部
Re是一个虚数的实部,而Im为一个虚数的虚部
1)在 p = 0 p=0 p=0附近的相图:
p=0处相图
(只制作这一个视频,主要是因为上传比较繁琐,其他视频上传到这篇博文的附属资源里了)
从-0.1到0,以0.05为步长制作相图:
(这一部分的代码以及视频打包上传到这篇博文的资源区了)
定理的部分还差一个参数 a a a实在是不好求,如果有同学会求麻烦联系我一下把这部分补充到这里,不然只能是等着老师讲解了(话说最后会上一节习题课么?至于这个到底是hopf分叉还是极限环,我真不知道了,麻烦懂的同学讲解一下儿~之后的三题都先把相图作好画在这里,代码同样上传至资源区。
2)相图:p(-1:0.2,1)
3)相图:p(-1:0.2,1)
4)相图:p(-1:0.2,1)