草庐IT

常微分初值

全部标签

知识点记录:李群李代数,微分流形,微分几何,图论

李群(Liegroup)是具有群结构的实流形或者复流形,并且群中的加法运算和逆元运算是流形中的解析映射。李代数Liealgebra):一类重要的非结合代数。非结合代数是环论的一个分支,与结合代数有着密切联系。结合代数的定义中把乘法结合律删去,就是非结合代数。微分流形(differentiablemanifold),也称为光滑流形(smoothmanifold),是拓扑学和几何学中一类重要的空间,是带有微分结构的拓扑流形。 微分几何Differentialgeometry是运用微积分的理论研究空间的几何性质的数学分支学科。爱因斯坦的广义相对论就以微分几何中的黎曼几何作为其重要的数学基础。拓扑学(

Sympy 微分方程符号解

常微分方程符号解从数学的角度出发只支持对形如 y'=f(x,y)f(x,y)只能是x(与其指数幂)与y(只支持一次)的非线性组合exampleimportsympyassp#定义变量x=sp.symbols("x")y1=sp.Function("y1")y2=sp.Function("y2")y3=sp.Function("y3")#定义微分方程组eq1=sp.Eq(y1(x).diff(x),2*y2(x)+x**3-(1+y1(x)))eq2=sp.Eq(y2(x).diff(x),3*x+y1(x)+2+y2(x))eq3=sp.Eq(y3(x).diff(x),y2(x))eq=(e

ENVI+ERDAS实现Hyperion叶绿素含量反演:经验比值法、一阶微分法

本文介绍基于ENVI与ERDAS软件,依据Hyperion高光谱遥感影像,采用经验比值法、一阶微分法等,对叶绿素含量等地表参数加以反演的具体操作。目录1前期准备与本文理论部分1.1几句闲谈1.2背景知识1.2.1Hyperion数据介绍1.2.2遥感图像分类方法1.2.3大气校正1.2.4反演算法2基于经验比值法、一阶微分法的叶绿素a含量反演2.1数据导入与波段合成2.2辐射定标与波段合成2.3编辑头文件2.4图像格式转换2.5EDRDAS文件导入与裁剪2.6监督分类2.7水体光谱曲线提取2.8特征波段选取与计算3大气校正及经验比值法波段调整3.1转换文件数据格式3.2FLAASH大气校正3.

python - tensorflow 中的矩阵行列式微分

我对使用TensorFlow计算矩阵行列式的导数很感兴趣。我通过实验可以看出,TensorFlow并没有实现通过行列式求微分的方法:LookupError:Nogradientdefinedforoperation'MatrixDeterminant'(optype:MatrixDeterminant)进一步调查表明,实际上可以计算导数;参见例如Jacobi'sformula.我确定,为了实现这种通过行列式进行区分的方法,我需要使用函数装饰器,@tf.RegisterGradient("MatrixDeterminant")def_sub_grad(op,grad):...但是,我对t

python - 是否有 python 模块来求解/集成随机微分方程组?

我有一个我想求解的随机微分方程组。我希望这个问题已经得到解决。我有点担心构建自己的求解器,因为我担心我的求解器会太慢,并且可能存在数值稳定性问题。是否有针对此类问题的python模块?如果没有,是否有解决此类系统的标准方法。 最佳答案 有一个:http://diffusion.cgu.edu.tw/ftp/sde/网站示例:"""addrequiredPythonpackages"""frompysdeimport*fromsympyimport*"""Variablesacclaimed"""x,dx=symbols('xdx')

python - 实现二阶导数的自动微分 : algorithm for traversing the computational graph?

我正在尝试实现automaticdifferentiation对于Python统计包(问题公式类似于优化问题公式)。计算图是使用运算符重载和用于sum()、exp()等操作的工厂函数生成的。我已经使用反向累加实现了梯度的自动微分。但是,我发现实现二阶导数(Hessian)的自动微分要困难得多。我知道如何进行单独的第二次局部梯度计算,但我很难想出一种智能的方法来遍历图形并进行累加。有谁知道为二阶导数提供自动微分算法的好文章或实现相同算法的开源库,我可能会尝试从中学习? 最佳答案 首先,您必须决定是要计算稀疏的Hessian矩阵还是更接

python - 自动微分

我正在做我的毕业设计,特别是关于流体动力学的,我有一个非线性方程组要求解,我选择牛顿法,所以我必须通过矩阵的雅可比矩阵(实际上是12x12矩阵).这个矩阵中的每个元素都是函数在某一点求值的导数,手工编写所有这些并计算每个导数是非常困难的;系统看起来像:f1(x1,x2,x3,...,x12)=0f2(x1,x2,x3,...,x12)=0...f12(x1,x2,x3,...,x12)=0其中x1、x2、x3是变量(温度、压力...等)我可以自动执行此操作吗?如果在Fortran中不行,我可以使用其他脚本语言作为Python(sympy模块)吗? 最佳答案

python - Tensorflow 中的可微分操作列表

是否有可微分(即自动微分)的Tensorflow操作的主列表?另外两种表达方式:未设置ops.NoGradient的操作列表。不会触发LookupError的操作列表。例如,我假设所有控制流操作都是不可微分的(例如,tf.where)。除了通过tf.gradients手动运行它们以查看它们是否抛出LookupError之外,我该如何找到它。“常识”不是有效答案。谢谢。编辑:tf.where是可微分的,所以我的直觉是错误的。也许这里正确的问题是Tensorflow中的哪些操作不可可微分。谢谢。 最佳答案 我已经使用Python代码设计

初识运算放大器(跟随器,比较器,同向/反向放大器,差分放大器,微分器,积分器)

(本人为电子学生小白,以下是个人学习过程中的归纳总结,如有错误,欢迎指正)虚短与虚断的理解虚断:输入电阻很大虚短:开环线性区,深度负反馈跟随器定义:跟随器是一种电子线路,其输出信号基本等同于输入信号,但提高了带负载能力,广泛存在于各类电子线路中。(来自百度)如图1所示,根据串联电阻分压可得同向端的电位V+=12*(2/(1+2))=8V,由虚断得反向端电位为8V,所以此时万用表显示8V图1跟随器 比较器 比较器是将一个模拟电压信号与一个基准电压相比较的电路。如图2所示,设VDC1为基准电压,VDC3为模拟电压,输出连接一个上拉电阻,此时同向端的模拟电压大于反向端的基准电压,则输出为12V。如图

非线性微分跟踪器

微分器描述离散形式的非线性微分跟踪为其中,h为采样周期;v(k)为第k时刻的输入信号,为决定跟踪快慢的系数;fst()为最速控制综合函数。描述如下仿真分析微分器测试输入信号v(t)=sin2pit,采样周期h=0.001,δ=150;扰动信号为0.05的随机信号;上图为带有噪声的正弦信号;下图为理想信号和微分跟踪器的输出信号;上图为理想微分信号与差分求得的微分信号;下图为理想微分信号和经过微分器输出的微分信号;代码如下所示:clearall;closeall;h=0.001;%Samplingtimedelta=150;r1_1=0;r2_1=0;vn_1=0;fork=1:1:10000ti