LSEC小陆的博客

原创 MATLAB绘图如何使用原点居中的标准坐标系（坐标系为十字带箭头）

目录初衷目标实现效果不足初衷在进行文档写作的时候，我们往往会利用 matlab 绘制几张图片插入到文章中。所有 matlab 画出来的图千篇一律，以至于看文章的人一眼就看出来，你这张图使用 matlab 画的。matlab 所画的图有个特色就是坐标刻度和标签在左边和下边，如下图所示的样子：十分刻板。我们现在想要把 matlab 做的图坐标轴移到图形的中间部分去(与数学的做图习惯一致)，也...

原创 Gronwall 不等式

Gronwall 不等式文章目录Gronwall 不等式定理一定理二定理三参考文献维基百科和百度百科上具有对 Gronwall 不等式的描述，但是有点凌乱，为了方便查询，我整理了一份干净的。Gronwall 引理是对满足一定微分不等式的一个实变量的（负）函数的基本估计。引理广泛用于研究演化问题（例如偏微分方程和常微分方程，连续动力系统）以限制依赖于时间的有界量的数学领域。不等式的最基本形式如下。定理一令ϕ:[0,T]→R\phi:[0, T] \rightarrow \mathbb{R}ϕ:

原创 [有限元方法基础理论] 质量集中有限元（谱元）

[有限元方法基础理论] 质量集中有限元（谱元）集中质量有限元，也称谱元，要求具有对角质量矩阵，且对角线上的值等于行和。可以通过用数值积分公式以及其他方法来构造谱元。不妨考虑一个初边值问题，ut−Δu=f in Ω,t>0u=0 on ∂Ω,t>0, with u(⋅,0)=v in Ω\begin{aligned}&u_{t}-\Delta u=f \quad \text { in } \Ome

原创 科研工具链汇总

科研工具链汇总（持续补充，欢迎收藏）SCI 期刊查询LetPub 网站查询：http://www.letpub.com.cn/index.php?page=journalappWeb Of Science 查询：http://apps.webofknowledge.com/Journal Citation Reports 官方平台上面查询：https://jcr.incites.thomsonreuters.com/JCRLandingPageAction.action兄弟单位

原创 [有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序（Helmholtz 方程）

[有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序（Helmholtz 方程）文章目录[有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序（Helmholtz 方程）前言一些链接动机问题描述关于特征值问题亥姆霍兹方程程序与结果问题一问题二写在后面的话前言一些链接之前写过三篇基础的有限元基础入门级别的程序介绍，引起比较大的反响，很多人在 CSDN 私信和添加我微信讨论了相关问题，有限元入门见如下链接。有限元方法入门：有限元方法简单的一维算例有限元方法入门：有限元方法简单的二维算例（三角形

原创 [MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合

[MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合文章目录[MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合前言引述MATLAB 曲线拟合函数简述一二维数据非线性拟合一维数据拟合例子二维数据拟合例子高维数据非线性拟合lsqcurvefitnlinfit 函数数据拟合工具箱的使用写在后面的话前言引述所谓曲线拟合，就是给定数据和数据可能满足的带参函数表达，希望能找到一组最好的参数，是的这个函数能够“最大限度”地刻画给定的数据。在一维（自变量维数）的情况下，换句话说，

原创 数据科学和机器学习中的优化理论与算法（下）

数据科学和机器学习中的优化理论与算法（下）数据科学和机器学习当前越来越热，其中涉及的优化知识颇多。很多人在做机器学习或者数据科学时，对其中和优化相关的数学基础，包括随机梯度下降、ADMM、KKT 条件，拉格朗日乘数法、对偶问题等，不太了解，一知半解地用，用着用着就出错了。本文希望从基础知识的角度，尽可能全地对最简单的优化理论和算法做一个小结。内容涵盖以下几个方面：优化简介、无约束优化、线搜索方法、信赖域方法、共轭梯度方法、拟牛顿方法、最小二乘问题、非线性方程、约束优化理论、非线性约束优化算法、二次规

原创 数据科学和机器学习中的优化理论与算法（上）

数据科学和机器学习中的优化理论与算法（上）数据科学和机器学习当前越来越热，其中涉及的优化知识颇多。很多人在做机器学习或者数据科学时，对其中和优化相关的数学基础，包括随机梯度下降、ADMM、KKT 条件，拉格朗日乘数法、对偶问题等，不太了解，一知半解地用，用着用着就出错了。本文希望从基础知识的角度，尽可能全地对最简单的优化理论和算法做一个小结。内容涵盖以下几个方面：优化简介、无约束优化、线搜索方法、信赖域方法、共轭梯度方法、拟牛顿方法、最小二乘问题、非线性方程、约束优化理论、非线性约束优化算法、二次规划

原创 数分习题的若干问题说明（持续扩展补充）

数分习题的若干问题说明文章目录数分习题的若干问题说明原函数可导，导函数未必连续参量方程问题关于向径和切线的夹角原函数可导，导函数未必连续 证 明 函数 \begin{array}{ll}\text { 证 明 函数 }\end{array} 证 明 函数 ​f(x)={e−1x2,x≠00,x=0f(x)=\left\{\begin{array}{ll}\mathrm{e}^{-\frac{1}{x^{2}}}, &

原创 数学分析中的点荟萃（持续补充修缮）

数学分析中的点荟萃（持续补充）文章目录数学分析中的点荟萃（持续补充）导数篇常用求导公式导数技巧微分中值定理基本认知积分篇积分技巧常用积分表处理原则积分表格林公式和高斯公式格林公式高斯公式积分等式积分中值定理积分第一中值定理积分第二中值定理广义积分柯西判别法无穷限积分瑕积分极限和连续篇极限技巧基本认知二阶曲面方程函数篇泰勒展开定义常见函数泰勒级数一些重要的特定函数符号函数黎曼函数狄利克雷函数分段三角函数分段指数函数其他函数双曲函数二元函数Γ\GammaΓ和β\betaβ函数常用等价无穷小一般stirling

原创 关于两个周期函数的和的周期性的讨论

关于两个周期函数的和的周期性的讨论因为排版和敲数学公式的局限性，很多地方写得并不是非常严格，或者有些跳跃，望海涵。初衷想这个问题的初衷是在给同学们习题课的时候（华东师大版的数学分析），里面有一道题，如下：求下列函数的周期:cos⁡x2+2sin⁡x3\cos \frac{x}{2}+2 \sin \frac{x}{3}cos2x​+2sin3x​。这道题本身比较简单，显然 12π12\pi12π 是它的一个周期，如果这里的周期理解为基本周期（最小正周期）的话（有同学发问了），我们还得 chec

原创 分布式优化和去中心化优化概述

分布式优化和去中心化优化概述陆嵩中科院 数学与系统科学研究院 科学与工程计算国家重点实验室文章目录分布式优化和去中心化优化概述简介预备知识优化基础理论梯度下降方法（GD）Gradient Descent步长选取收敛率约束优化问题的拉格朗日乘数法次梯度算法图谱理论基础平均一致性问题模型问题梯度算法收敛性结果对于去中心化的梯度下降方法凸问题的 DGD 算法一些结论和讨论凸问题的较新方法ADMM（原始-对偶）对偶上升对偶分解增广拉格朗日方法ADMM 及其收敛性分布式 ADMM 算法EXTRA和其他方法非

原创 最优控制与深度学习

最优控制与深度学习这是 2020 北京大学“应用数学专题讲习班”课程的一个极简的总结，省略了很多符号描述、定理条件描述等等，最后一部分，关于本征正交分解和动态模态分解，以及 Koopman 分析的一些东西，不太想写了，权当做个记录吧。文章目录最优控制与深度学习最优控制方程PMPHJB平均场最优控制公式最优性性条件在平均场中的迁移最优控制诱导的网络学习算法神经网络稳定性稳定性不同控制方程离散诱导的神经网络结构PDE 神经网络架构动力系统逼近的 NN 数学结果最优控制方程考虑最优控制的 Bolza

原创 泛函和变分：从最速降线谈起

泛函和变分：从最速降线谈起处理最速降线等类似的问题，大家都喜欢直接对推导出来的极小化问题，直接代入欧拉-拉格朗日方程求解，我觉得这不易于对“变分”的理解，对于大家理论分析功底的增近，并无裨益。通过思考，我觉得直接从 Frechet 导数出发，进行变分推导得到更具化的欧拉-拉格朗日方程，是更为广泛的一种变分方法。本质上是一样的，只不过这个过程更容易推广。文章目录泛函和变分：从最速降线谈起简介变分Frechet 导数最速降线问题变分求解最速降线最速降线的物理意义简介所谓的泛函，就是“函数的函数”，自

原创 CTeX 套装替代方案：Texlive+Texmaker

cTeX 套装替代方案：Texlive+Texmaker文章目录cTeX 套装替代方案：Texlive+Texmaker容易混淆的概念说明LaTeX 编译系统LaTeX 文本编辑器编译引擎CTeX 套装和 ctex 宏包CTeX or TeXlive？？？对于 CTeX 套装或者 WinEdt 有执念一些下载链接容易混淆的概念说明LaTeX 编译系统编译系统用于提供编译环境，主流的编译系统有，如 TEXLive（MacOS 下叫 MacTeX）和 MiKTEX，包括他们的一些七七八八的衍生品。La

原创 通过例子理解头文件库文件的编译链接

原力计划

通过例子理解头文件库文件的编译链接文章目录通过例子理解头文件库文件的编译链接多 cpp 和 h 文件编译链接例子关于方法实现位置的一个说明打包成静态链接库和动态链接库进行链接小小的总结网传头文件的包含和静动态链接库的链接无非就是文件的赋值黏贴而已。是耶非耶？是的话，是如何复制，又是如何黏贴，今天通过这篇文章，通过简单的例子，带你们了解一下多个 cpp 文件和 h 头文件之间是如何组织在一块的。多 cpp 和 h 文件编译链接例子为了说明这个问题，首先我写了以下几个文件，命名分别为 1main.c

原创 MPI 通讯开销测算

一个简单的pingpong程序测试mpi消息通讯的开销及并行计算通讯启动时间测算一个简单的pingpong程序测试mpi消息通讯的开销随着科技的进步，集群单节点计算能力的提高，似乎通讯开销成了并行计算中dominant，再提高计算能力对于并行的增益似乎效果不明显，限制性能的瓶颈从处理器计算能力上转移到通讯开销上。显然，此时设法降低MPI消息通讯带来的时间消耗，成为了当务之急。因此，写了一个极...

原创 网页弹窗命令的使用

下载下图到工作目录，并保存为名为pc.jpg复制以下代码到记事本，并保存为html到同图片目录下。<!DOCTYPE html><html><head><meta http-equiv="Content-Type" content="text/html; charset=GBK"><meta itemprop="name" content="你被耍了额"><meta itemprop="ima

原创 C++ 高性能计算之多线程简单基础入门教程

C/C++ 高性能计算之多线程简单基础入门教程比起别人的盲目罗列函数接口，鹦鹉学舌式的解释每一个输入参数和输出参数，一味求全而无重点，我的文章更侧重于入门知识的讲解，宁缺毋滥，只有一些最简单的入门用法，先亮代码，让你以最快地速度感受到这个东西原来是这个样子。更深入的知识，可以查看文档以及其他参考。不管什么东西，都是尽可能地用新不用旧，所以，文章涉及的库和工具，都是基于 C++11 标准的。前面有相关的博文已经介绍过 MPI，OpenMP 等相关的并行技术，以及 hadoop 生态圈相关的一些知识。

原创 汉诺塔游戏玩法介绍（攻略和编程实现）

汉诺塔游戏玩法介绍（攻略和编程实现）汉诺塔游戏介绍汉诺塔是根据一个传说形成的数学问题。描述如下：有三根杆子 A、B、C 。 A 杆上 N 个 穿孔圆盘，盘的尺寸由下到上依次变小。要求按下列规则将所有圆盘移至 C 杆：每次只能移动一个圆盘大盘不能叠在小盘上面可将圆盘临时置于 B 杆，也可将从 A 杆移出的圆盘重新移回 A 杆，但都必须遵循上述两条规则。问：如何移？最少要移动多少次？最...

原创 C++ 基础之公有、私有和保护（一句话说清楚）

C++ 基础之公有、私有和保护公有继承、私有继承和保护继承，什么七七八八的，一句话说清楚。数据成员和成员函数的公有、私有和保护类的定义我们知道，定义一个类，类里面的成员包括数据成员（属性）和成员函数（方法）。数据成员是描述类的属性的一些数据，成员函数是描述类的行为的一些数据。类的成员根据它们的访问权限，在类内被分成了三类：公有（public）、私有（private）和保护（prote...

原创 C++ 基础之虚函数和虚基类

C++ 基础之虚函数和虚基类C++ 有别于 C ,比较重要的就是其抽象、封装、继承和多态等特性。多态又具体地体现在重载、虚函数和模板等等之上。今天就和大家聊一聊虚函数和虚基类。虚函数我们来思考一个问题。假如子（派生）类 B 继承自父（基）类 A，现在有一个 A 类型的指针指向了 B，这时候指针只能调用 B 中继承自 A 的那一部分东西。那么，现在假设 B 中有一个和 A 中同名、同传入参数...

原创 C++ 基础之函数模板和类模板

C++ 基础之函数模板和类模板C++ 有别于 C ,比较重要的就是其抽象、封装、继承和多态等特性。多态又具体地体现在重载、虚函数和模板等等之上。今天就和大家聊一聊函数模板和类模板。函数模板我们来看函数重载：观察发现，这几个函数代码基本相同，只有参数类型不同，功能相同。问题就来了：写多个相同的函数体，是重复劳动，而且还可能因为代码的冗余造成不一致性。就不能找个东西去通用地表示和替代这些类...

原创 C++ 基础之函数重载和内联函数

C++ 基础之函数重载和内联函数C++ 有别于 C ,比较重要的就是其抽象、封装、继承和多态等特性。多态又具体地体现在重载、虚函数和模板等之上。今天就和大家聊一聊函数重载和内联函数。函数重载所谓的函数重载，就是说我们可以定义不同功能的同名函数，它们传入参数类型是不一样的，当你调用这个函数名的时候，编译器会根据你给的参数类型，去匹配和调用相对于的函数。举例如下，这里交换两个数的值，使用了...

原创 交替方向乘子法（ADMM）的数学基础

交替方向乘子法（ADMM）网上的一些资料根本就没有把ADMM的来龙去脉说清楚，发现只是一个地方简单写了一下流程，别的地方就各种抄，共轭函数，对偶梯度上升什么的，都没讲清楚，给跪了。下面我来讲讲在机器学习中用得很多的ADMM方法到底是何方神圣。共轭函数给定函数f:Rn→Rf: \mathbb{R}^{n} \rightarrow \mathbb{R}f:Rn→R，那么函数f∗(y)=max⁡...

原创 三道题入门C语言

三道题入门C语言判断输入天数为当年的第几天输入一个日期，格式为xxxx-xx-xx，判断这一天为当年的第几天？输入：2019-1-2输出：2#include<stdio.h>#include<string.h>#include<stdlib.h>#include<time.h>int main()/*用譬如2019-12-...

原创 从一个简单的示例入门Cmake

从一个简单的示例入门Cmake文件树形目录结构 ├─meshes│ └─poisson├─param│ ├─geom│ │ └─unspecified│ ├─levelset│ │ ├─brick_transp│ │ ├─film│ │ ├─prJump│ │ ├─twophasedrops│ │ └─unspecified│ ├─os...

原创 Ubuntu(linux)知识琐碎汇总（持续更新）

Ubuntu 18.xx 安装和使用systemback：Ubuntu经常玩崩溃，时常备份以防不测。另外，有些软件安装起来极其复杂，安装好后，制作成iso镜像文件，分享给别人，也是一桩好事。systemback可用于系统的备份还原以及iso镜像制作。终端安装，代码如下：sudo add-apt-repository --remove ppa:nemh/systembacksudo add-...

原创 压缩感知基础入门知识

压缩感知压缩感知是2006年才开始兴起的研究方向，它主要是借助自然信号的规律性，从而可大大减少观测次数。这在很多领域都有很好的应用前景。对于自然信号的规律性，用数学语言可以做很多种描述，比较流行的一种就是自然信号在一组基底表示下是稀疏的。压缩感知说的是对于方程组Ax=bAx=bAx=b，AAA是m×Nm\times Nm×N，这里m<Nm<Nm<N的矩阵，什么条件下，能保证x...

原创 函数型数据主成分分析（FPCA）

本文主要介绍了以下几个方面的内容：简单介绍了经典的主成分分析方法，包括其数学推导，算法步骤，和几个实际算例；简单介绍了其它的数据降维方法，譬如局部线性嵌入以及它的简单算例；更近一步，我们介绍了函数型主成分分析方法（FPCA），包括其基本思想、数学推导、算法描述等，最为重要的是，我们将该方法和本领域进行结合，有了一些新的思考，感谢"数据科学与矩阵优化"课程给带来的灵感。文章目录前言主成分分析（...

原创 基于nips改进的LaTex中文模板

之前写了一份中文模板。但是标题摘要那一块趋于扁平化，没有比较酷的style。另外，如果用latex编译，是不支持除了eps之外的图片的插入。故而，做了这一份基于nips模板改进的模板。首先，下载nips官方提供的模板，使用nips14submit_e.sty文件。Ctex中文套装，编辑器是WinEdt，采用pdflatex编译。源文件如下。\documentclass[12pt,a4pap...

原创 关于评论、私信和邮件的一个声明

原创 曲面偏微分方程：参数化有限元方法

曲面偏微分方程：参数化有限元方法前面介绍的P\mathbf{P}P和Pd\mathbf{Pd}Pd，及其对应的同参等值延伸及投影自然延伸，对于曲面上有限元方法的理论分析起着至关重要的作用。利普希茨参数曲面上的FEM利普希茨参数化曲面因为同参映射P\mathbf{P}P的双利普希茨性，存在一个常数LLL，使得 (64) L−1∣x1−x2∣≤∣x~1−x~2∣≤L∣x1...

原创 曲面扰动理论之C2曲面的H2扩张

曲面扰动理论之C2曲面的H2扩张我们在后续提到的窄带（narrow band）方法和迹方法（trace）的分析中，对解要求：u~∈H2(γ)\widetilde{u} \in H^{2}(\gamma)u∈H2(γ)∥u∥H2(N(δ))≲δ12∥u~∥H2(γ)\|u\|_{H^{2}(\mathcal{N}(\delta))} \lesssim \delta^{\frac{1}{2...

原创 流体力学引论

流体力学引论简介流体力学以流体为研究对象，流体就是在剪切外力作用下会发生流动的物体。流体力学的研究方法包括实验观测和理论分析。物质结构包括气体、固体、液体等等。所谓的流体力学，是建立在连续介质假设基础上的流体力学。所谓的连续介质假设，是说我们把研究的对象当成一个整体，考虑大量分子的平均行为，不考虑物质的分子结构和单个分子的运动细节。应力张量流体的作用力有体积力和面积力，体积力一般指的是外...

原创 曲面扰动理论之C2曲面的扰动

曲面扰动理论之C2曲面的扰动依然假设邻域：(48)N={x∈Rn+1:∣d(x)∣<12K∞}(48) \quad \mathcal{N}=\left\{\mathrm{x} \in \mathbb{R}^{n+1}:|d(\mathrm{x})|<\frac{1}{2 K_{\infty}}\right\}(48)N={x∈Rn+1:∣d(x)∣<2K∞​1​}假设有...

原创 曲面扰动理论之C1曲面的扰动

曲面扰动理论之C1曲面的扰动定义曲面χΓ\chi_\GammaχΓ​，定义uΓ∈H#1(Γ)u_{\Gamma} \in H_{\#}^{1}(\Gamma)uΓ​∈H#1​(Γ)，它是如下弱形式的解：(33)∫Γ∇ΓuΓ⋅∇Γv=∫ΓfΓv∀v∈H#1(Γ)(33) \quad \int_{\Gamma} \nabla_{\Gamma} u_{\Gamma} \cdot \nabla_{\...

原创 曲面上的微积分（二）

曲面上的微积分（二）从上一次讨论结束，回去之后，我花了一点时间走了一遍证明。心中依然遗留了两个问题。对于Laplace-Beltrami方程的弱形式，∫γφ~Δγv~=−∫γ∇γφ~⋅∇γv~\int_{\gamma} \tilde{\varphi} \Delta_{\gamma} \tilde{v}=-\int_{\gamma} \nabla_{\gamma} \tilde{\varph...

原创 曲面上的微积分（一）

曲面上的微积分（一）参数化方法、迹方法和窄带方法是三种主要方法去求解Laplace-Beltrami问题。首先，我会介绍一些关于曲面上微积分的知识。参数化曲面定义曲面的参数化表示，χi:Vi→Ui∩γ⊂Rn+1\chi_{i} : \mathcal{V}_{i} \rightarrow \mathcal{U}_{i} \cap \gamma \subset \mathbb{R}^{n+1}...

原创 linux下.bashrc的使用和个人用户家目录下的软件安装

linux下.bashrc的使用和个人用户家目录下的软件安装.bashrc文件简介:环境配置.bashrc文件主要保存个人的一些个性化设置，如命令别名替代、设置系统的环境变量(PATH)等。在公用的集群环境下，个人用户登录到该账户下的家目录，此只对这个用户的个性化设置起效。它们是一些隐藏文件，ls -a可查看。可以看到有四个相关的文件：.bashrc：当你登入shell时以及每次打开新的...