小访客的博客

原创 1,2 多维运动公式

原创 1,1 一维运动公式

已知条件：加速度a，时间t，初始速度V0，初始位置X0

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

一，预备知识非线性自治微分方程组：{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)​图像是一个速度场：F⃗=fi^+gj^\vec{F}=f\widehat{i}+g\widehat...

原创 第三十二讲 极限环

一，预备知识：非线性一阶微分自治方程组的一般形式：{x′=f(x,y)y′=g(x,y)\left\{\begin{matrix}{x}'=f(x,y)\\ {y}'=g(x,y)\end{matrix}\right.{x′=f(x,y)y′=g(x,y)​等式右边不显含变量t等式右边是非线性函数（如三角函数、二次项）速度场：F⃗=fi^+gj^=x′...

原创 第三十一讲 非线性微分自治方程组及图解

一，非线性一阶微分自治方程组的一般形式：{x′=f(x,y)y′=g(x,y)\left\{\begin{matrix}{x}'=f(x,y)\\ {y}'=g(x,y)\end{matrix}\right.{x′=f(x,y)y′=g(x,y)​等式右边不显含变量t等式右边是非线性函数（如三角函数、二次项）二，例题（含阻尼的非线性摆）：如图，一根木...

原创 第二十三讲 解一阶微分方程组

例题{u1′=−u1+2u2u1′=u1−2u2\left\{\begin{matrix}{u_{1}}'=-u_{1}+2u_{2}\\ {u_{1}}'=u_{1}-2u_{2}\end{matrix}\right.{u1​′=−u1​+2u2​u1​′=u1​−2u2​​一，将微分方程组化为矩阵形式[u1′u2′]=A[u1u2]=[−121−2][...

原创 第三十讲 解耦

一，线性的变量代换：耦合的方程组：{u=ax+byv=cx+dy\left\{\begin{matrix}u=ax+by\\ v=cx+dy\end{matrix}\right.{u=ax+byv=cx+dy​，求u，v解耦后的两个一阶方程：{u′=k1uv′=k2v\left\{\begin{matrix}{u}'=k_{1}u \\ {v}'=k_{2}...

原创 第二十九讲 求方程组通解和特解的公式（矩阵指数）

一，齐次方程组的通解：通解形式：x⃗=c1x1⃗+c2x2⃗\vec{x}=c_{1}\vec{x_{1}}+c_{2}\vec{x_{2}}x=c1​x1​​+c2​x2​​用基本矩阵简化为：x⃗=[x1⃗x2⃗][c1c2]=Xc⃗\vec{x}=\begin{bmatrix}\vec{x_{1}} & \vec{x_{2}}\end{bmatrix}\begin{bmat...

原创 第二十八讲 解非齐次线性方程组

一，关于二阶方程组x⃗′=Ax⃗{\vec{x}}'=A\vec{x}x′=Ax的理论（对n阶方程也成立）：（假设A是常数矩阵）定理A：x⃗′=Ax⃗{\vec{x}}'=A\vec{x}x′=Ax的通解是x⃗=c1x1⃗+c2x2⃗\vec{x}=c_{1}\vec{x_{1}}+c_{2}\vec{x_{2}}x=c1​x1​​+c2​x2​​（x1⃗...

原创 第二十七讲 微分方程组解的图像

一，竞争模型（含参数的微分方程组）：{x′=−x+byy′=cx−3y\left\{\begin{matrix}{x}'=-x+by \\ {y}'=cx-3y\end{matrix}\right.{x′=−x+byy′=cx−3y​这个模型是上海和北京之间的旅游竞争，两个地方都忙于做广告，想要吸引游客x′{x}'x′代表上海，y′{y}...

原创 第二十六讲 有特殊特征值的微分方程组

一，温度扩散问题：假设有三个房间，房间初始温度分别为：高、中、低，但经过一段时间后，三个房间互相交换热量，最终稳定在同一个温度上。假设房间温度分别为：X1,X2,X3X_{1},X_{2},X_{3}X1​,X2​,X3​，如图：温度Xi(t)X_{i}(t)Xi​(t)是时间的函数二，建立微分方程组：根据温度变化率，建立数学模型：X1′=a(X3−X1)+a(X2−X1){X_{1}}...

原创 第二十二讲 对角化分解和幂公式

一，对角化分解A=SΛS−1A=S\Lambda S^{-1}A=SΛS−1A表示具有n个线性无关的x（特征向量）的矩阵S表示由x组成的可逆方阵，称作特征向量矩阵Λ\LambdaΛ表示由A的λ（特征值）作为对角元素的对角矩阵比较：A=LU（消元化分解），A=QR（正交化分解）二，矩阵的幂...

原创 第二十一讲 特征值和特征向量

我个人认为麻省理工线性代数这门课，到二十一讲才真正进入有用的部分，因此从这一讲开始做笔记。一，概念满足条件：Ax=λx解释：当向量x经过矩阵A变换后，效果等于向量x乘上任意常数λ则：x是矩阵A的特征向量，λ是矩阵A的特征值二，性质性质1：如果A是奇异矩阵，且Ax=0，则x是0空间的非0向量，λ=0注：奇异=不可逆=线性相关，非奇异=可逆=线性无关性质2：λ的和=A的迹trace解...

原创 第二十五讲 用线性代数解微分方程组

一，上一讲的例题，如图：设x=T1x=T_1x=T1​，y=T2y=T_2y=T2​方程组为·：{x′=−2x+2yy′=2x−5y\left\{\begin{matrix}{x}'=-2x+2y\\ {y}'=2x-5y\end{matrix}\right.{x′=−2x+2yy′=2x−5y​用消元法求出的通解为：{x=c1e−t+c2e−6ty=1...

原创 第二十四讲 一阶常微分方程组

一，常微分方程组：必须满足：自变量只有1个，因变量有多个 一阶形式：，，x和y是因变量，t是自变量二，什么样的常微分方程组是线性方程组？必须满足：x和y以线性形式出现 例如：，，a、b、c和d可以是t的函数三，什么样的线性方程组是常系数方程组?必须满足：a、b、c和d都是常数 例如：，四，什么样的常系数方程组是齐次的？必须满足：， 例如：，五，初始条件：需要初始...

原创 第二十三讲 狄拉克函数（冲激函数）（补充）

四，例题： 如图，假设没有阻尼，弹簧常数，在的瞬间，输入一个冲量A，单位冲量为 建立数学模型： 初始条件：（初始位置），（初始速度） 第一步：两边拉式变换左边： 代入初始条件：， 左边： 右边：（利用延迟定理1） 左边等于右边：第二步：解出第三步：拉式逆变换（查表） （查表） 因为中含有指数函数，因此要用到逆变换的唯一性 利用延迟定理1： 因此 当时， 当...

原创 第二十三讲 狄拉克函数（冲激函数）

一，脉冲及建立模型在一个时间区间上的作用等于，如果是恒定的，则 如图： 在小车上作用一个脉冲（一个力使小车左移或右移一段时间），假设时间从0到h，冲量是1（曲线下的面积是1）如图： 假设没有阻尼，弹簧常数，建立数学模型（参照第九讲）： 因为，是单位方框函数（参考第二十二讲） 所以二，用拉普拉斯变换解这个微分方程（参考第二十讲）：初始条件：， 第一步：两边拉氏变换左边...

原创 第二十二讲 延迟定理（补充）

八，计算：已知 因为当时，，因此可以用延迟定理 根据延迟定理1： 同理： 因此九，计算：因为当时，，因此可以用延迟定理 根据延迟定理2： 因此 图像见视频35:00~36:40十，计算： 当表达式中存在指数时，如，就要用到逆变换的唯一性： 利用延迟定理1： 因此 当时，， 当时，， 结果...

原创 第二十二讲 延迟定理

一，单位阶跃函数：如图：，它有三个定义： 设为单位阶跃函数，那么表示单位阶跃函数向右平移a个单位，如图： 二，单位方框函数： 如图： 意义：去掉了在区间以外的部分。三，单位阶跃函数的拉普拉斯变换： 因为当时，，所以： ， （查表）四，逆变换的唯一性：如图： 因为拉普拉斯变换只关注这段区间，所以在这段区间内相等的函数，变换后的结果相等（无法区别）...

原创 第二十一讲 卷积公式

一，卷积公式：已知：， 设： 求： 因为拉氏变换是由幂级数变过来的，所以上面的问题可以转换为下面的问题方便计算：已知：， 设： 求：，（求解过程省略）解得卷积公式： 文字解读：两个函数的乘积，等于分别将它们变换后的乘积，再逆变换的结果，由于被变换卷在了一起，因此称为卷积。 满足交换律：二，例1：求： 代入卷积公式： 验证：因为：，， 所以：三，例2：求...

原创 第二十讲 拉普拉斯变换求解线性ODE

一，拉氏变换公式：，二，保证的拉氏变换存在的条件：只有在收敛的情况下，拉氏变换才存在 增长条件：保证拉氏变换存在的前提是的增长速度不能超过的收敛速度， 指数型表达法：，，， 正例： 反例：，因为 ，因为，，三，导数的拉氏变换：解微分方程： 设初始条件：， 一阶导数的拉氏变换：（分部积分） 要使，必须满足增长条件：，，，，或者说 代入初始条件：...

原创 第十九讲 拉普拉斯变换引入

一，从傅里叶变换到拉普拉斯变换：傅里叶变换： 拉普拉斯变换：，其中为衰减因子， 为什么要在傅里叶变换中乘上衰减因子？ 因为当非周期函数随时间单调递增或单调递减，趋于无穷大（直男）时，无法使用傅里叶变换，如图：  因此给乘上衰减因子，使其在远处逐渐衰减下来（被掰弯），就可以用傅里叶变换了，如图（红色）： 令，s表示拉普拉斯算子 拉普拉斯变换标准形式： 收敛域，如图（红色...

原创 可降阶的高阶微分方程

一，n阶微分方程的一般形式或 通解：二，n阶微分方程思路：这种n阶微分方程可以通过n次积分求出通解 三，特殊的二阶微分方程第一种：，不显含y 思路：令，则原方程化为，解出p的通解，再还原成，得y的通解 第二种：，不显含x 思路：令，，原方程化为，解出p的通解，再还原成，得y的通解 ...

原创 第十讲 幂级数展开式的应用

一，对误差函数的近似计算 二，对正弦积分函数的近似计算 三，复数项级数的性质定义：复数项级数收敛，当且仅当实部级数和虚部级数都收敛，如图： 绝对收敛定理：若复数项级数各项模的级数收敛，则也收敛。证明如图： 四，欧拉公式 ...

原创 第九讲 函数间接展开成幂级数

一，间接展开法（利用已知的展开式，展开其他函数）已知的展开式：， ， ， ，二，例题1，将展开成x的幂级数，三，例题2，将展开成x的幂级数，四，例题3，将展开成x的幂级数思路：先求导，再积分 ， 如果，那么求就得后面加个，如图： 五，例题4，将展开成x的幂级数，六，例题5，将展开成x-1的幂级数，七，二项展开式 八，利用函数的展开式求高阶导...

原创 第八讲 函数直接展开成幂级数

一，泰勒级数定义：若在点处有各阶导数，，则可求出在点处的泰勒级数： 这是个的幂级数，收敛中心在，收敛区间为 如果设，则得到麦克劳林级数：二，函数的泰勒展开式是唯一的：若，则三，的充要条件：泰勒公式的余项，在x处的极限为0： 若，则 泰勒公式如图： 四，将函数展开成麦克劳林级数的步骤：直接展开法求出在处的各阶导数：,，如果处某阶导数不存在，则不能展开 写出幂级...

原创 第七讲 幂级数的运算与和函数

一，幂级数的四则运算设幂级数的收敛半径为，幂级数的收敛半径为 ，收敛域为 二，幂级数的连续性、导数和积分设幂级数，收敛区间为（-R，R）则幂级数在收敛区间（-R，R）内连续 则幂级数在收敛区间（-R，R）内有各阶导数，并且可以逐项求导：，收敛区间不变但收敛域可能改变 则幂级数在收敛区间（-R，R）内可积，并且可以逐项积分：，收敛区间不变但收敛域可能改变三，求幂级数的和函数...

原创 第六讲 幂级数的收敛半径和收敛域

一，函数项级数定义：， 部分和： 收敛点：使函数项级数收敛的点 发散点：使函数项级数发散的点 收敛域：D={  收敛}，即所有收敛点的集合 和函数：， 余项：二，幂级数及其收敛性幂级数 有一个明显的收敛点：x=0 幂级数是否收敛，跟的形式相关 幂级数的收敛有三种可能：处处绝对收敛，如 仅在x=0处收敛，如 在以x=0为中心，两边对称的区间内收敛，如，收敛域（-1,...

原创 第五讲 交错级数、绝对收敛和条件收敛

一，负项级数可以通过添加负号转化为正项级数，因此不用讨论二，交错级数定义：，（）三，交错级数的审敛法（又称Leibniz判别法）若，（），即逐项递减 且，即通项趋于0 则收敛，其和 余项 截断误差：余项四，数列的极限理论一个数列，如果它的偶次项子数列收敛于s，奇次项子数列也收敛于s，那么整个数列将收敛于s。如图： 五，变号级数的审敛法任意变号级数：非正项级数，非...

原创 第四讲 比值、根值和积分审敛法

一，正项级数的比值审敛法设是正项级数，且， 若，则收敛 若，则发散 若或不存在，则无法判断敛散性 例题，如图： 二，正项级数的根值审敛法设是正项级数，且， 若，则收敛 若，则发散 若或不存在，则无法判断敛散性 ， 例题，如图： 三，正项级数的积分审敛法设函数在区间连续，非负，单调递减 则级数与广义积分具有相同的敛散性 例题，如图： 四，正项级数...

原创 第三讲 正项级数的比较审敛法

一，用比较审敛法的原因对于大多数级数，很难得到部分和的表达式，因此很难用定义来研究其敛散性 比较审敛法不能用于非正项级数二，正项级数定义：，其中 部分和为单调递增数列 收敛的充要条件：有上界 推论：若无界，则发散，且 正项级数结合律：发散+发散=发散三，正项级数的比较审敛法设和是正项级数，且，（） 若收敛，则收敛 若发散，则发散四，p-级数的敛散性定义：，（）...

原创 第二讲 无穷级数的性质

一，收敛级数的分配率若收敛，则也收敛，且 推论：若，则和具有相同的敛散性二，收敛级数的结合律若和收敛，则也收敛，且 推论1：若收敛，发散，则必发散 推论2：若发散，发散，则未必发散三，收敛级数的线性性质四，前有限项与敛散性无关改变（减少或增加）级数的有限项，不会改变该级数的敛散性 推论：决定级数敛散性的不是前有限项，而是级数的余项五，收敛级数内部的结合律收敛级数...

原创 第一讲 无穷级数的概念

一，无穷级数的定义 通项： 部分和： 余项： 离散的和： 连续的和（反常积分）：二，级数的收敛和发散收敛：，s为任意数余项 ，其中截断误差为发散：，或者不存在三，等比级数，（）的敛散性展开： 通项：，公比： 假设：等比数列求和公式 当时，，，级数收敛 当时，，，级数发散假设：，级数发散 假设：展开： 当n为奇数时，， 当n为偶数时，， ...

原创 第十八讲 傅里叶变换

一，傅里叶级数的复数形式：一般形式：，周期， ，n为任意整数 ，n为任意整数利用逆向欧拉公式（第十讲第二节）展开： 原式： 第一项： 第三项：当时， 代入原式： 当时： 当时： 当时： 总结：，二，傅里叶变换：用表示周期为T的， 对于任意，其中和都是不变的 唯一改变的是 频谱如图（图中的和意思相同）： 两个不同角速度之间的间距，如图（图中的和意思相...

原创 第十七讲 利用傅里叶级数求特解

一，几何变换法，求傅里叶级数：假设是周期，的函数，求它的傅里叶级数。如图1： 第一步，求周期的函数的傅里叶级数。如图2： 因为是奇函数， 当，；当， 因此，前提： 第二步，将压缩成，周期改变，变，如图3： 的周期，，振幅是的一半： 基频率， 将替换：，前提： 用上一节拓展1的公式验证答案：因为是奇函数， 当，；当， 因此，前提：第三步，将图3上移1...

原创 第十六讲 傅里叶级数拓展

一，推论：       一个函数不能有两个不同的傅里叶级数，因为傅里叶系数公式表明对应唯一的和。二，奇偶性（缩减计算）：如果是偶函数，即，y轴对称，那么，，证明： 如果是偶函数，那么是偶函数， 因此，，因为是偶函数如果是奇函数，即，原点对称，那么，，证明：同上 如果是奇函数，那么也是偶函数（奇X奇=偶）， 因此，，因为是奇函数三，傅里叶级数与泰勒级数的不同之...

原创 第十五讲 傅里叶级数引入

注：在学习这章之前建议看看《漫画傅里叶解析》和可汗学院的傅里叶级数课程，那里对傅里叶级数有很基础的讲解。一，二阶非齐次常系数线性ODE的标准形式：       二，将输入项变成傅里叶级数：             任何以为周期的周期函数f(t)都可以用傅里叶级数表示三，输入项和输出项（通解）的关系：输入（驱动）项f(t) 输出（响应）项y(t) ...

原创 第十四讲 共振

一，共振的数学模型： 弹簧—质量—阻尼系统 如图 标准形式： 设阻尼常数，弹性常数，表示弹簧振荡的角速度， 原方程化为： 假设阻尼系数，方程右侧输入的驱动项，驱动角速度 （注：视频中的，这里用表示，以便区别第三节的阻尼角速度） 方程为： 方程左边换成线性算子式： 方程右边复数化：是的实部，将代替 方程化简为：，表示复数解， 因为 指数输入定理 取出实部： 幅度...

原创 第十三讲 二阶非齐次常系数线性ODE的特解

一，二阶非齐次常系数线性ODE的标准形式：二，通解： 三，将方程化成特殊形式：设方程右边的输入项为“纯振荡”：，，表示复数 方程左边换成线性算子式： 四，代换法则： 证明： D表示对函数求导：， 五，指数输入定理（当时）：的特解为，是常数 证明：将特解代入原方程，利用代换法则六，求特解例题1： 将复数化：是的虚部，将代替， 方程左边换成线性算子式：...

原创 第十二讲 二阶非齐次线性ODE解的结构

一，二阶非齐次线性ODE的标准形式： 表示输入、驱动  表示输出、响应二，通解：  叫补充解，是齐次方程的通解 （作为的相伴方程） 是非齐次方程的特解三，应用：弹簧—质量—阻尼系统（第九讲） 如图 在系统中，对小车施加额外的一个力 标准形式： 被动系统：当时，系统没有外力干涉，只是被动地对初始条件进行响应 受迫系统：当时，系统一直有外力干涉，强迫它运动 电...