《微分方程与动力系统（Differential Equations and Dynamical Systems）》听课笔记3

为了深入理解fMRI分析的原理，还是要学习一些微分方程与动力系统的东东，这是我的学习笔记，内容还没有很好的整理。

视频在这里1-Differential Equations and Dynamical Systems Overview_哔哩哔哩_bilibili

课程的原始地址在这里：ME 564 - Mechanical Engineering Analysis (washington.edu)

17. 2x2 阶ODEs源与汇

17-2x2 Systems of ODEs Sources and Sinks_哔哩哔哩_bilibili

此处的源（source）我理解是系统的固定点（fixed point）当系统处于其他状态时，由于系统是不稳定的，因此会向其他方向(可以理解为向相平面上的其他状态进行转移)发散。
如下图所示，圆球是可以立于另一个球的顶端的，但是当球处于其他位置时，会迅速滚走。最顶点的平衡状态，即可认为是系统的source。

image.png

此处的汇（sink）我理解是系统的稳定点，当系统处于其他状态时，会迅速的向稳定点状态运动。

image.png

$\underline{\dot x} = A \underline{x}$

若
$A= \begin{bmatrix} 3 & -1 \\ -1 & 3 \\ \end{bmatrix}$
则
$T = \begin{bmatrix} 1 & 1 \\ 1 & -1 \\ \end{bmatrix}$
$D = \begin{bmatrix} 2 & 0 \\ 0 & 4\\ \end{bmatrix}$
$T^{-1} = \begin{bmatrix} 0.5 & 0.5 \\ 0.5 & -0.5\\ \end{bmatrix}$

image.png

如图所示，图中的黑色坐标为 $\underline{x}$ 的坐标，蓝色的表示两个特征向量的方向，在这两个方向上，得到的 $\underline{z}$ 是彼此解耦的，在 $(1,1)$ 方向上的运动速度是在 $(1,-1)$ 方向上速度的 $\frac{1}{2}$ （特征值是2倍的关系）

18. 2x2系统的ode鞍点及其不稳定性

18-2x2 Systems of ODEs Saddle Points and Instability_哔哩哔哩_bilibili

此处建议看视频，了解一下鞍点的含义

19. 2x2 ode的虚特征值和中心不动点

19-2x2 Systems of ODEs Imaginary Eigenvalues and Center Fixed Points_哔哩哔哩_bilibili

特征值的实数部分给出了增加或减少的趋势，特征值的虚部则代表了周期震荡

20. 特征值与稳定性

20-Stability and Eigenvalues What does it mean to be a stable eigenvalue_哔哩哔哩_bilibili

任意特征值的实数部分 $Re(\lambda)>0$ 则，系统不稳定；若所有 $Re(\lambda)<0$ ，则系统稳定。

21. 线性微分方程

21-What is a Linear Differential Equation_哔哩哔哩_bilibili

线性的定义为：
$f(\alpha x+ \beta y)=\alpha f(x)+ \beta f(y)$

22. 非线性微分方程在固定点附近的线性化

22-Linearizing Nonlinear Differential Equations Near a Fixed Point_哔哩哔哩_bilibili

非线性系统是一个比较复杂的问题，我们一般采用首先找到其固定点，之后在固定点附近进行线性化处理。

对于微分方程：
$\dot x = f(x)$
对于任意点 $\overline x$ 而言，若 $f(\overline x) = 0$ ，则 $\overline{x}$ 是一个固定点。

对于 $\overline x$ 附近的点 $x=\overline x + \Delta x$ ，则利用泰勒公式展开可得：

$\begin{align} \dot x = f(x) =&f(\overline x+\Delta x) \\ =&f(\overline x) + \frac{Df}{Dx}( \overline x )\Delta x + \frac{D^2f}{Dx^2}( \overline x)\Delta x^2 + \cdots \end{align}$

与此同时，考虑到：
$\begin{align} \because& f(\overline{x}) = 0 \\ \therefore& \dot x \approx \frac{Df}{Dx}( \overline x )\Delta x \\ \therefore& \dot x = \frac{d}{dt}( \overline x+\Delta x) = \frac{d}{dt}\Delta x=\frac{Df}{Dt}(\Delta x) \end{align}$

此处得到了一个关于 $\Delta x$ 的线性常微分方程，已知：

$\frac{Df}{Dx}= \begin{bmatrix} \frac{\partial f_1}{x_1} & \frac{\partial f_1}{x_2} & \cdots & \frac{\partial f_1}{x_n}\\ \frac{\partial f_2}{x_1} & \frac{\partial f_2}{x_2} & \cdots & \frac{\partial f_2}{x_n}\\ \vdots & \vdots & \ddots & \vdots \\ \frac{\partial f_n}{x_1} & \frac{\partial f_n}{x_2} & \cdots & \frac{\partial f_n}{x_n}\\ \end{bmatrix}$

一个例子，已知：
$\underline {f}(\underline x)= \begin{bmatrix} f_1(x_1, x_2) \\ f_2(x_1, x_2) \\ \end{bmatrix} = \begin{bmatrix} x_1-x_1^2 \\ x_1+x_2 \\ \end{bmatrix}$
所以固定点为
$\overline x = \begin{bmatrix} 0 \\ 1 \\ \end{bmatrix}， \overline x = \begin{bmatrix} 1 \\ -1 \\ \end{bmatrix}$

雅可比为：
$\frac{Df}{Dx}= \begin{bmatrix} 1-2x_1 & 0 \\ 1 & 1 \end{bmatrix}$

对于固定点1,
$\frac{Df}{Dx}(\overline x)= \begin{bmatrix} 1 & 0 \\ 1 & 1 \\ \end{bmatrix}$
特征值为
$\lambda=1,1$

对于固定点2,
$\frac{Df}{Dx}(\overline x)= \begin{bmatrix} -1 & 0 \\ 1 & 1 \\ \end{bmatrix}$
特征值为
$\lambda=\pm1$

最后编辑于：2023.12.14 11:45:04

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 159,835评论 4赞 364
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 67,598评论 1赞 295
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 109,569评论 0赞 244
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 44,159评论 0赞 213
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 52,533评论 3赞 287
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 40,710评论 1赞 222
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 31,923评论 2赞 313
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 30,674评论 0赞 203
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 34,421评论 1赞 246
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 30,622评论 2赞 245
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 32,115评论 1赞 260
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 28,428评论 2赞 254
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 33,114评论 3赞 238
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 26,097评论 0赞 8
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 26,875评论 0赞 197
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 35,753评论 2赞 276
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 35,649评论 2赞 271

《微分方程与动力系统（Differential Equations and Dynamical Systems）》听课笔记3

17. 2x2 阶ODEs源与汇

18. 2x2系统的ode鞍点及其不稳定性

19. 2x2 ode的虚特征值和中心不动点

20. 特征值与稳定性

21. 线性微分方程

22. 非线性微分方程在固定点附近的线性化

推荐阅读更多精彩内容