非线性ode及在工程中的应用第五次作业更正版

第一题:

第五次作业第一题

(1)

题式的平衡方程是: u ˙ = p − u 2 = 0 \dot u=p-u^2=0 u˙=pu2=0
解得 u = ± p ( p > 0 ) u=\pm \sqrt p (p\gt0) u=±p (p>0)以及 u 在 p ≤ 0 时 无 解 u在p\le0时无解 up0

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˙=puu2+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)14p 或者 u = 1 / 2 − ( 1 / 2 ) 1 − 4 p u=1/2-(1/2)\sqrt{1-4p} u=1/2(1/2)14p 当然对于后面两个解 p ≤ 1 / 4 p\le1/4 p1/4

f u ( u , p ) = p − 2 u + 3 u 2 = 0 f_u(u,p)=p-2u+3u^2=0 fu(u,p)=p2u+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˙=p2uu3=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)=p2u+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)为鞍节点

稳定性的部分,实在是不能理解,麻烦读到这里而且理解的朋友和我简单聊聊,可以吗~

第二题:

第二题图片,在课本的P168
首先补充一点这个定理相关的知识:
PAhopf定理
至 于 里 面 的 这 么 多 参 数 至于里面的这么多参数
R e 是 一 个 虚 数 的 实 部 , 而 I m 为 一 个 虚 数 的 虚 部 Re是一个虚数的实部,而Im为一个虚数的虚部 ReIm

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)

在这里插入图片描述

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值