第三十三讲 非线性方程组化为一阶方程
一,預備知識
非線性自治微分方程組:{dxdt=f(x,y)dydt=g(x,y)\left\{\begin{matrix}\frac{dx}{dt}=f(x,y)\\ \frac{dy}{dt}=g(x,y)\end{matrix}\right.{dtdx?=f(x,y)dtdy?=g(x,y)?
等式右邊不顯含變量t
圖像是一個速度場:F?=fi^+gj^\vec{F}=f\widehat{i}+g\widehat{j}F=fi+gj?
方程組的解為:{x=x(t)y=y(t)\left\{\begin{matrix}x=x(t)\\ y=y(t)\end{matrix}\right.{x=x(t)y=y(t)?
圖像是一條軌跡:
二,從方程組中消除t:
只需將方程組中的兩個方程相除:dydx=g(x,y)f(x,y)\frac{dy}{dx}=\frac{g(x,y)}{f(x,y)}dxdy?=f(x,y)g(x,y)?
此時方程組變?yōu)榱艘浑A常微分方程
圖像中,消去了自變量t,也就沒有了速度(向量的長度和方向),只剩下該點的斜率,速度場變成了斜率場。
解不再是一對參數(shù)方程,而是y′=y(x){y}'=y(x)y′=y(x)(顯函數(shù),隱函數(shù)都有可能),不再是軌跡,而是曲線。
消除t的好處:有可能使原來無法解的方程組變得可解。
三,線性方程組例題:
線性方程組:{x′=yy′=?x\left\{\begin{matrix}{x}'=y\\ {y}'=-x\end{matrix}\right.{x′=yy′=?x?
通解:[xy]=c1[cos(t)?sin(t)]+c2[sin(t)cos(t)]\begin{bmatrix}x\\ y\end{bmatrix}=c_{1}\begin{bmatrix}cos(t)\\ -sin(t)\end{bmatrix}+c_{2}\begin{bmatrix}sin(t)\\ cos(t)\end{bmatrix}[xy?]=c1?[cos(t)?sin(t)?]+c2?[sin(t)cos(t)?]
圖像:
消去t:dydx=?xy\frac{dy}{dx}=\frac{-x}{y}dxdy?=y?x?
分離變量再積分,得通解:x2+y2=cx^{2}+y^{2}=cx2+y2=c
圖像也是圓。
四,非線性方程組例題:
鯊魚(x)-小魚(y)方程組:{x′=?ax+bxyy′=cy?dxy\left\{\begin{matrix}{x}'=-ax+bxy\\ {y}'=cy-dxy\end{matrix}\right.{x′=?ax+bxyy′=cy?dxy?,(a,b,c,d>0a,b,c,d>0a,b,c,d>0)
-ax表示沒有小魚,鯊魚會消亡(假設鯊魚只吃小魚)
bxy表示鯊魚吃了小魚
cy表示沒有鯊魚,小魚會增多(這里不是邏輯斯蒂增長方程)
-dxy表示小魚被鯊魚吃了
為了簡化方程組,假設a,b,c,d=1a,b,c,d=1a,b,c,d=1
方程組變?yōu)?#xff1a;{x′=?x+xyy′=y?xy\left\{\begin{matrix}{x}'=-x+xy\\ {y}'=y-xy\end{matrix}\right.{x′=?x+xyy′=y?xy?
第一步,找臨界點:
計算{x(?1+1y)=0y(1?1x)=0\left\{\begin{matrix}x(-1+1y)=0\\ y(1-1x)=0\end{matrix}\right.{x(?1+1y)=0y(1?1x)=0?
解得:{x=0y=0\left\{\begin{matrix}x=0\\ y=0\end{matrix}\right.{x=0y=0?,{x=cd=1y=ab=1\left\{\begin{matrix}x=\frac{c}ze8trgl8bvbq=1\\ y=\frac{a}=1\end{matrix}\right.{x=dc?=1y=ba?=1?
第二步,對每個臨界點附近,線性化方程組,并畫出軌跡:
當{x=0y=0\left\{\begin{matrix}x=0\\ y=0\end{matrix}\right.{x=0y=0?時:xy是兩個無窮小的乘積,可以忽略
方程組變?yōu)?#xff1a;{x′=?xy′=y\left\{\begin{matrix}{x}'=-x\\ {y}'=y\end{matrix}\right.{x′=?xy′=y?
矩陣A:[?1001]\begin{bmatrix}-1 & 0\\ 0 & 1\end{bmatrix}[?10?01?]
特征值:λ1=?1\lambda _{1}=-1λ1?=?1,λ2=1\lambda _{2}=1λ2?=1
依此可以判斷,圖像是鞍形(不穩(wěn)定),λ1=?1\lambda _{1}=-1λ1?=?1的特征向量朝向原點,λ2=1\lambda _{2}=1λ2?=1的特征向量朝向無窮遠
圖像:
當{x=1y=1\left\{\begin{matrix}x=1\\ y=1\end{matrix}\right.{x=1y=1?時:沒有無窮小了
計算雅克比矩陣:J0=[?1+yx?y1?x]=[01?10]J_{0}=\begin{bmatrix}-1+y & x\\ -y & 1-x\end{bmatrix}=\begin{bmatrix}0 & 1\\ -1 & 0\end{bmatrix}J0?=[?1+y?y?x1?x?]=[0?1?10?]
即方程組:{x′=yy′=?x\left\{\begin{matrix}{x}'=y\\ {y}'=-x\end{matrix}\right.{x′=yy′=?x?
圖像:
里面和外面有很多同心圓(沒畫在圖上),這種叫中心。
現(xiàn)在,陷入了邊界線情形的困境:
邊界線情形:
中心隔開了“螺旋源區(qū)域”和“螺旋匯聚區(qū)域”
圖上的點(0,1)對應跡=0,行列式=1的方程組,如果矩陣的系數(shù)(方程組的系數(shù))有微小的變化,點(跡,行列式)就不會在中心線上,原來的同心圓圖像,就會變成螺旋源圖像或螺旋匯聚圖像。如圖:
因為中心是非線性方程組的近似(實際的系數(shù)肯定有微小的變化),所以無法判斷這個臨界點是螺旋匯聚還是螺旋源,或者仍是中心。
這個問題叫沃爾泰拉問題。
解決辦法,消去t:dydx=y(1?x)x(?1+y)\frac{dy}{dx}=\frac{y(1-x)}{x(-1+y)}dxdy?=x(?1+y)y(1?x)?
分離變量再積分,(過程見視頻),得通解:xe?x?y?e?y=cxe^{-x}\cdot y\cdot e^{-y}=cxe?x?y?e?y=c
通解是方程xe?x?y?e?y=h(x,y)xe^{-x}\cdot y\cdot e^{-y}=h(x,y)xe?x?y?e?y=h(x,y)的一條等值線
ue?uue^{-u}ue?u的函數(shù)圖像:
當u→0u\rightarrow 0u→0時,ue?u→uue^{-u}\rightarrow uue?u→u
當u→∞u\rightarrow \inftyu→∞時,ue?u→0ue^{-u}\rightarrow 0ue?u→0
對ue?uue^{-u}ue?u求導,當u=1時,導數(shù)=0,得最大值e?1e^{-1}e?1
因此,xe?x?y?e?y=h(x,y)xe^{-x}\cdot y\cdot e^{-y}=h(x,y)xe?x?y?e?y=h(x,y)這個方程的最大值在點(x=1,y=1)處,當x=0或y=0時h(x,y)=0h(x,y)=0h(x,y)=0,在最大值和0值之間環(huán)繞著等值線,如圖:
通常每一條水平線和一條等值線只有兩個交點,但在最大值時只有一個交點,因此等值線不可能是螺旋。
因為沒有小魚,鯊魚會消亡;沒有鯊魚,小魚會增多;小魚增多,鯊魚開始增多……因此可以判斷等值線的方向是順時針。
以恒定速率k捕魚的影響:
方程組變?yōu)?#xff1a;{x′=?ax+bxy?kxy′=cy?dxy?ky\left\{\begin{matrix}{x}'=-ax+bxy-kx\\ {y}'=cy-dxy-ky\end{matrix}\right.{x′=?ax+bxy?kxy′=cy?dxy?ky?
整理:{x′=?(a+k)x+bxyy′=(c?k)y?dxy\left\{\begin{matrix}{x}'=-(a+k)x+bxy\\ {y}'=(c-k)y-dxy\end{matrix}\right.{x′=?(a+k)x+bxyy′=(c?k)y?dxy?
原來的臨界點是{x=cd=1y=ab=1\left\{\begin{matrix}x=\frac{c}ze8trgl8bvbq=1\\ y=\frac{a}=1\end{matrix}\right.{x=dc?=1y=ba?=1?,現(xiàn)在變?yōu)?#xff1a;{x=c?kdy=a+kb\left\{\begin{matrix}x=\frac{c-k}ze8trgl8bvbq\\ y=\frac{a+k}\end{matrix}\right.{x=dc?k?y=ba+k??
捕魚的結果:降低了鯊魚的數(shù)量,增加了小魚的數(shù)量。
這個現(xiàn)象叫沃爾泰拉法則。
完結撒花!
總結
以上是生活随笔為你收集整理的第三十三讲 非线性方程组化为一阶方程的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 第三十一讲 非线性微分自治方程组及图解
- 下一篇: 1,1 一维运动公式