日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

第三十三讲 非线性方程组化为一阶方程

發(fā)布時(shí)間:2025/4/16 编程问答 15 豆豆
生活随笔 收集整理的這篇文章主要介紹了 第三十三讲 非线性方程组化为一阶方程 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

一,預(yù)備知識(shí)
非線性自治微分方程組:{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
圖像是一個(gè)速度場(chǎng):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:
只需將方程組中的兩個(gè)方程相除:dydx=g(x,y)f(x,y)\frac{dy}{dx}=\frac{g(x,y)}{f(x,y)}dxdy?=f(x,y)g(x,y)?
此時(shí)方程組變?yōu)榱艘浑A常微分方程
圖像中,消去了自變量t,也就沒有了速度(向量的長度和方向),只剩下該點(diǎn)的斜率,速度場(chǎng)變成了斜率場(chǎng)。
解不再是一對(duì)參數(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表示沒有小魚,鯊魚會(huì)消亡(假設(shè)鯊魚只吃小魚)
bxy表示鯊魚吃了小魚
cy表示沒有鯊魚,小魚會(huì)增多(這里不是邏輯斯蒂增長方程)
-dxy表示小魚被鯊魚吃了
為了簡(jiǎn)化方程組,假設(shè)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?
第一步,找臨界點(diǎn)
計(jì)算{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}ozvdkddzhkzd=1\\ y=\frac{a}{b}=1\end{matrix}\right.{x=dc?=1y=ba?=1?
第二步,對(duì)每個(gè)臨界點(diǎn)附近,線性化方程組,并畫出軌跡
當(dāng){x=0y=0\left\{\begin{matrix}x=0\\ y=0\end{matrix}\right.{x=0y=0?時(shí):xy是兩個(gè)無窮小的乘積,可以忽略
方程組變?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的特征向量朝向原點(diǎn),λ2=1\lambda _{2}=1λ2?=1的特征向量朝向無窮遠(yuǎn)
圖像:

當(dāng){x=1y=1\left\{\begin{matrix}x=1\\ y=1\end{matrix}\right.{x=1y=1?時(shí):沒有無窮小了
計(jì)算雅克比矩陣: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ū)域”
圖上的點(diǎn)(0,1)對(duì)應(yīng)跡=0,行列式=1的方程組,如果矩陣的系數(shù)(方程組的系數(shù))有微小的變化,點(diǎn)(跡,行列式)就不會(huì)在中心線上,原來的同心圓圖像,就會(huì)變成螺旋源圖像或螺旋匯聚圖像。如圖:

因?yàn)橹行氖欠蔷€性方程組的近似(實(shí)際的系數(shù)肯定有微小的變化),所以無法判斷這個(gè)臨界點(diǎn)是螺旋匯聚還是螺旋源,或者仍是中心。
這個(gè)問題叫沃爾泰拉問題。
解決辦法,消去tdydx=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ù)圖像:

當(dāng)u→0u\rightarrow 0u0時(shí),ue?u→uue^{-u}\rightarrow uue?uu
當(dāng)u→∞u\rightarrow \inftyu時(shí),ue?u→0ue^{-u}\rightarrow 0ue?u0
對(duì)ue?uue^{-u}ue?u求導(dǎo),當(dāng)u=1時(shí),導(dǎo)數(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)這個(gè)方程的最大值在點(diǎn)(x=1,y=1)處,當(dāng)x=0或y=0時(shí)h(x,y)=0h(x,y)=0h(x,y)=0,在最大值和0值之間環(huán)繞著等值線,如圖:

通常每一條水平線和一條等值線只有兩個(gè)交點(diǎn),但在最大值時(shí)只有一個(gè)交點(diǎn),因此等值線不可能是螺旋。
因?yàn)闆]有小魚,鯊魚會(huì)消亡;沒有鯊魚,小魚會(huì)增多;小魚增多,鯊魚開始增多……因此可以判斷等值線的方向是順時(shí)針。

以恒定速率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?
原來的臨界點(diǎn)是{x=cd=1y=ab=1\left\{\begin{matrix}x=\frac{c}ozvdkddzhkzd=1\\ y=\frac{a}{b}=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}ozvdkddzhkzd\\ y=\frac{a+k}{b}\end{matrix}\right.{x=dc?k?y=ba+k??

捕魚的結(jié)果:降低了鯊魚的數(shù)量,增加了小魚的數(shù)量。
這個(gè)現(xiàn)象叫沃爾泰拉法則。

完結(jié)撒花!

總結(jié)

以上是生活随笔為你收集整理的第三十三讲 非线性方程组化为一阶方程的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。

主站蜘蛛池模板: 国产超碰91 | av手机天堂网 | 六月综合激情 | 亚洲av区无码字幕中文色 | 一区二区三区欧美在线 | 极品美女高潮出白浆 | 9i看片成人免费看片 | 欧美一级黄色大片 | 亚洲一区二区在线免费 | 精品欧美在线 | 一二三不卡视频 | av网站免费观看 | 孕妇疯狂做爰xxxⅹ 国产精品乱码久久久久久 99久久久成人国产精品 | 日本一二三区视频 | 不卡av在线 | 日韩中字在线 | 亚洲精品在线中文字幕 | 国产精品毛片一区二区在线看舒淇 | 久久尤物视频 | 欧美10p | 国产大片一区 | 国产经典自拍 | 国产精品免费电影 | 日本不卡一区在线观看 | 一区二区三区资源 | 日韩1页| 亚洲一区二区三区国产 | 视频在线a | juliaann第一次和老师 | 免费看黄色一级片 | 性做久久久久久 | 精品视频在线播放 | 极品美女被c | 爆操白虎| 亚洲男女激情 | 国产精品腿扒开做爽爽爽挤奶网站 | 中文字幕视频一区二区 | xxxx69视频| 在线观看成人免费 | 美女试爆场恐怖电影在线观看 | 看片久久 | 欧美永久精品 | 超碰av免费 | 夜色视频网 | 中文字幕精品一二三四五六七八 | 免费啊v在线观看 | 99re在线精品视频 | 亚洲天堂美女 | 激情综合丁香五月 | 人妻少妇偷人精品久久性色 | 精品国产免费看 | 对白刺激theporn| 久久久久无码国产精品一区李宗瑞 | 调教少妇视频 | 亚洲av无码精品色午夜果冻不卡 | 日日夜夜网 | 在线国产黄色 | 日日日插插插 | 日本人性爱视频 | 香蕉久久久久久 | 黄色裸体片 | 香蕉成视频人app下载安装 | 亚洲精品视频一区 | 丰满少妇在线观看资源站 | 国产超碰人人爽人人做人人爱 | 灌满闺乖女h高h调教尿h | 激情小说激情视频 | 久青草视频 | 涩涩涩av | 精品久久久久久久久久久国产字幕 | 国产精品熟妇一区二区三区四区 | 欧美变态另类刺激 | 不卡av中文字幕 | 日韩精品视频网 | 怡红院av亚洲一区二区三区h | 久久电影一区二区 | 国产成人无码精品亚洲 | 亲嘴扒胸摸屁股免费视频日本网站 | 久草操| 五月婷婷综 | 国产麻豆午夜三级精品 | 天天综合色网 | 99爱在线观看 | 国产精品一区二区三区线羞羞网站 | 极品少妇在线 | 欧美日韩免费一区二区三区 | 欧美极品一区二区三区 | 人妻妺妺窝人体色www聚色窝 | 国产啊v在线 | 啪视频在线观看 | 成人欧美一区二区三区 | 黄色片视频免费看 | 特级av| 欧美精品网址 | 亚洲一级片网站 | 久久久久久久久久久久久久免费看 | 五月天婷婷影院 | 毛片免费播放 | 日韩孕交|