- 博客(407)
- 资源 (17)
- 收藏
- 关注
RSS订阅
原创 MATLAB 提取图片曲线
文章目录我用 MATLAB 提取图片曲线行文动机图像的读入与裁剪颜色拾取颜色转换与色差计算分离曲线二值化，提取数据数据点分类与排序后话我用 MATLAB 提取图片曲线给你一张图片，如何提取里面曲线的数据，从而利用这些数据进行图像重绘、加工处理、测距、拟合得到函数表达式等操作呢？行文动机前段时间，有个朋友问了我一个问题，大概意思就是要给图像的流线测距离，在我的印象里面，MATLAB 是似乎没有这种直接的功能的。那么换个角度来理解一下这个问题，如果给你一张图像，如何提取里面点的数据？其实，有了曲
2021-08-16 23:09:19
9893
19
原创 反赌科普宣传——赌徒的谬误
反赌科普宣传——赌徒的谬误文章目录反赌科普宣传——赌徒的谬误游戏简单介绍不公平的游戏彩民的谬误求神拜佛，是邪？非邪？隐形的陷阱，无底的深渊在赌博中，有一种“荷兰赌式”的“倍赌”策略，即输了加赌注加倍，直到赢了为止。似乎在所有人看来，只要本金足够，似乎此法无往而不利。那么这种做法可不可取？背后潜藏的玄机又在哪呢？我们从一个简单的游戏说起。所用的概率论知识，一个正常上过大学的理工科的人都应该会懂。游戏简单介绍在我们的家乡广泛流行着一个“机会游戏”，叫做“六合彩”，其游戏规则如下：人的年龄（虚岁）对
2021-06-13 11:01:25
2915
11
原创 线性方程组的求解方法的选择和加速
线性方程组的求解方法的选择和加速文章目录线性方程组的求解方法的选择和加速迭代法概述矩阵分裂方法Krylov 子空间方法线性方程组迭代方法的选择共轭梯度类残差类更牛逼的加速手段：多重网格和并行（并发）计算其他加速方法迭代法概述求解线性方程组一般有直接解法（分解）、迭代解法和启发式方法。如果系数矩阵 A 很大并且是稀疏矩阵，分解方法一般情况下将不会有效。大部分情况下，小规模计算喜欢用直接解法，不存在时间问题了，就不需要加速。大规模问题的计算，我们青睐于用迭代解法，也有的求解器做出了适用于大规模问题的直接分
2021-05-29 21:29:16
4422
8
原创 MATLAB 求函数极值的内置函数一览表（实则优化算法函数汇总）
MATLAB 求函数极值的内置函数一览表收集了几乎所有的 MATLAB 内置的优化函数，可收藏，需要时查阅。文章目录MATLAB 求函数极值的内置函数一览表简介概览求函数极值（或最大值最小值）求零点（解方程）最小二乘问题求极值fminbnd：单变量fmincon：约束、非线性、多变量fminunc：无约束、多变量fminsearch：无约束、多变量、无导数linprog：线性规划quadprog：二次规划fminimax：minmax 问题fgoalattain：目标达到fseminf：半无限约束、
2021-05-22 17:38:09
3886
4
原创 [有限元方法基础理论] 质量集中有限元（谱元）
[有限元方法基础理论] 质量集中有限元（谱元）集中质量有限元，也称谱元，要求具有对角质量矩阵，且对角线上的值等于行和。可以通过用数值积分公式以及其他方法来构造谱元。不妨考虑一个初边值问题，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
2021-01-12 22:49:50
1508
原创 [有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序（Helmholtz 方程）
[有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序（Helmholtz 方程）文章目录[有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序（Helmholtz 方程）前言一些链接动机问题描述关于特征值问题亥姆霍兹方程程序与结果问题一问题二写在后面的话前言一些链接之前写过三篇基础的有限元基础入门级别的程序介绍，引起比较大的反响，很多人在 CSDN 私信和添加我微信讨论了相关问题，有限元入门见如下链接。有限元方法入门：有限元方法简单的一维算例有限元方法入门：有限元方法简单的二维算例（三角形
2020-12-26 02:05:39
2965
8
原创 [MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合
[MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合文章目录[MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合前言引述MATLAB 曲线拟合函数简述一二维数据非线性拟合一维数据拟合例子二维数据拟合例子高维数据非线性拟合lsqcurvefitnlinfit 函数数据拟合工具箱的使用写在后面的话前言引述所谓曲线拟合，就是给定数据和数据可能满足的带参函数表达，希望能找到一组最好的参数，是的这个函数能够“最大限度”地刻画给定的数据。在一维（自变量维数）的情况下，换句话说，
2020-12-22 01:07:39
5122
原创 数据科学和机器学习中的优化理论与算法（下）
数据科学和机器学习中的优化理论与算法（下）数据科学和机器学习当前越来越热，其中涉及的优化知识颇多。很多人在做机器学习或者数据科学时，对其中和优化相关的数学基础，包括随机梯度下降、ADMM、KKT 条件，拉格朗日乘数法、对偶问题等，不太了解，一知半解地用，用着用着就出错了。本文希望从基础知识的角度，尽可能全地对最简单的优化理论和算法做一个小结。内容涵盖以下几个方面：优化简介、无约束优化、线搜索方法、信赖域方法、共轭梯度方法、拟牛顿方法、最小二乘问题、非线性方程、约束优化理论、非线性约束优化算法、二次规
2020-12-21 15:37:53
790
原创 数据科学和机器学习中的优化理论与算法（上）
数据科学和机器学习中的优化理论与算法（上）数据科学和机器学习当前越来越热，其中涉及的优化知识颇多。很多人在做机器学习或者数据科学时，对其中和优化相关的数学基础，包括随机梯度下降、ADMM、KKT 条件，拉格朗日乘数法、对偶问题等，不太了解，一知半解地用，用着用着就出错了。本文希望从基础知识的角度，尽可能全地对最简单的优化理论和算法做一个小结。内容涵盖以下几个方面：优化简介、无约束优化、线搜索方法、信赖域方法、共轭梯度方法、拟牛顿方法、最小二乘问题、非线性方程、约束优化理论、非线性约束优化算法、二次规划
2020-12-21 15:28:46
1026
2
原创 函数型数据主成分分析（FPCA）
本文主要介绍了以下几个方面的内容：简单介绍了经典的主成分分析方法，包括其数学推导，算法步骤，和几个实际算例；简单介绍了其它的数据降维方法，譬如局部线性嵌入以及它的简单算例；更近一步，我们介绍了函数型主成分分析方法（FPCA），包括其基本思想、数学推导、算法描述等，最为重要的是，我们将该方法和本领域进行结合，有了一些新的思考，感谢"数据科学与矩阵优化"课程给带来的灵感。文章目录前言主成分分析（...
2019-11-24 19:57:27
14708
18
原创 五线谱中蕴含的数学问题
降调判断成什么，我们一般读成降什么调，比如说降C大调，降D大调，降E大调，降G大调，降A大调，降B大调。调号称呼前面的升降，表示主音在原有的基础上升降半个音，用于正推，即告诉你是调号称呼，让你用“全全半全全全半”的方法推断升降在哪些调上的时候，如果是升（升C，升F），主音首先升半个调，然后来推断，降调亦然。将音名C、D、E、F、G、A、B，分别唱做do、re、mi、fa、sol、la、si，固定不变，不依调性调式的变化而改变唱法，这种唱（谱）法叫做“固定唱名法”（也就是说，各音名的唱名是固定的。
2019-07-25 16:27:08
4993
15
原创 强化学习（强化主干知识，看这篇就够了）
强化学习强化学习到底是什么？若是想略微了解强化学习大概的意思和操作，把握一个主干知识，而不做深入研究，看这篇文章就很够了。文章目录强化学习核心内容元素简介优化目标：状态值函数状态-动作值函数策略迭代和值迭代策略迭代值迭代Q学习、SARSA、时序差分学习Q学习SARSA和时序差分学习深度强化学习一些学习资料核心内容元素简介强化学习考察的是对于一个对象（智能体），从一个状态执行某种动作到...
2019-07-23 11:56:22
934
原创 一维激波管（Lax shock tube）问题的数值求解
一维激波管（Lax shock tube）问题的数值求解文章目录一维激波管（Lax shock tube）问题的数值求解@[toc]问题描述一般的守恒算法格式精确解常用守恒格式（先介绍标量）Lax-Friedrichs格式Roe格式Lax-Wendroff格式ENO格式简介（先介绍标量）5阶WENO格式空间离散Lax-Friedrichs通量分裂原函数方法重构左值右值空间算子和时间三阶Rung...
2019-05-28 11:52:59
12495
47
原创 函数型数据主成分分析
函数型数据主成分分析基本思想函数型主成分分析（FPCA，Functional Principal Components Analysis）是传统的PCA的一种推广。考虑我们已经从数据中得到拟合曲线xi(s),s∈T,i=1,⋯ ,nx_{i}(s), s \in \mathcal{T}, i=1, \cdots, nxi(s),s∈T,i=1,⋯,n，所谓的第一...
2019-04-15 08:59:29
9308
25
原创 快速行进算法（Fast Marching）
快速行进算法（Fast Marching）Fast Marching方法简介快速行进算法 （Fast Marching Method） 是求解程函方程 （Eikonal Equation）F∣∇T∣=1F|\nabla T|=1F∣∇T∣=1 的一种高效数值算法，程函方程属于非线性偏微分方程，可以认为是一种近似波动方程。在界面演化问题中的求解的结果 TTT 的物理含义是曲线以速度 F(x⃗)F...
2019-03-11 21:01:07
21700
43
原创 FFT的并行实现
关键词：快速傅里叶变换 高维FFT 并行计算快速傅里叶变换简介离散傅里叶变换离散傅里叶变换（DFT）一般定义为： Fn≡∑k=0N−1fke−2πink/NF_n \equiv \sum_{k=0}^{N-1}f_ke^{-2\pi ink/N}Fn≡k=0∑N−1fke−2πink/N离散傅里叶逆变换可以定义为： fk=1N∑k=0N−1Fne−2πink/Nf_...
2018-12-21 11:06:23
11539
5
原创 方程组求解的直接法与迭代法实现
方程组求解的直接法与迭代法实现问题描述我们的目的在于求解如下所示的方程组：其中的A11、A12、A21、A22A_{11}、A_{12}、A_{21}、A_{22}A11、A12、A21、A22以及右端项fff以不同的的稀疏存储方式存在文件中。我们需要完成的问题包括：理解并描述A11、A12、A21、A22A_{11}、A_{12}、A_{21}、A_{22}A11、A12...
2018-12-14 20:25:01
5317
9
原创 傅里叶谱方法及其工程实现
傅里叶谱方法求解PDE快速傅里叶变换用一维举例，在(−∞,+∞)(-\infty,+\infty)(−∞,+∞)有定义且绝对可积，它的傅里叶变换及其逆变换为：u^(k)=∫−∞+∞u(x)e−ikxdxu(k)=12π∫−∞+∞u^(x)eikxdx\hat u(k)=\int_{-\infty}^{+\infty}u(x)e^{-ikx}dx\\u(k)=\frac{1}{2\pi}\...
2018-11-19 20:28:11
14320
25
原创 深度学习中的一些基本概念
深度学习中的一些基本概念一维卷积（离散）对于一维卷积的理解，是在一次坐火车的时候。所谓的一维卷积，无非就是两列火车相向而驰，从头碰头，到尾离尾的过程。我们把两个向量比作两列火车，向量的元素比作每个车厢，那么两个向量做卷积得到的值的每个分量就是计算每个时刻火车重叠部分的点乘，所谓的点乘，就是对应位置做乘法最后求和。花了一点时间做了一个gif如下：这里需要注意的是，两列火车的长度不定，更不一定...
2018-11-06 16:35:37
4375
2
原创 多重网格方法
多重网格方法基本思想一般的迭代法是在一种固定的网格上进行迭代，当网格比较细时，计算量十分大。多重网格说的是，在计算细网格上的精确解时，其初值是比它粗一些网格上的精确解构造的，因而迭代次数少。当然，求较粗的网格上的精确解，它的初值就是有更粗一些的网格上的精确解所构造的，如此往复，直到一个相当粗的网格位置。多重网格法的求解过程是一个递归的过程。算法过程预平滑选定一个初值，先采用一般方法，比...
2018-10-05 02:44:41
18474
27
原创 图像与视频处理中的优化方法
问题与初步思考问题重述问题一：视频的前后背景分离给定一段视频（也可以说给定很多帧连续的图像），视频中有一些人或者别的东西在动（前景），也有一些不动的景象（背景）。要求使用背景建模的相关技术（最好使用优化方法），区分出前景和背景，输出为视频的每一帧对应的前景和背景。如图[pic1]所示：图像，左边是输入，右边是输出。问题二：主成分提取输入为一张由若干张小图放在一块...
2018-07-26 23:07:35
10146
原创 有限元方法入门：有限元方法简单的二维算例（三角形剖分）
有限元方法简单的二维算例（三角形剖分）算例描述我们对下述椭圆边值问题 \label{eq1}{−Δu=fu|∂Ω=0{−Δu=fu|∂Ω=0\left \{\begin{aligned} & -\Delta u = f\\& u|_{\partial \Omega} = 0 &\end{aligned}\right. 其中，Ω=(0,1)2Ω=(0,1)2\Omega = (...
2018-07-20 01:09:08
30873
72
原创 有限元方法入门：有限元方法简单的二维算例（矩形剖分）
有限元方法简单的二维算例（矩形剖分）算例描述我们对下述椭圆边值问题 \label{eq1}{−Δu=fu|∂Ω=0{−Δu=fu|∂Ω=0\left \{\begin{aligned} & -\Delta u = f\\& u|_{\partial \Omega} = 0 &\end{aligned}\right. 其中，Ω=(0,1)2Ω=(0,1)2\Omega = (0...
2018-07-20 00:52:28
17032
31
原创 有限元方法入门：有限元方法简单的一维算例
有限元方法简单的一维算例，入门，基础，且介绍清晰，带你走进有限元的世界。算例描述我们对下述边值问题\label{eq1} {u′′(x)+u(x)=(1−π2)sinπxu(0)=u(1)=00≤x≤1{u″(x)+u(x)=(1−π2)sinπx0≤x≤1u(0)=u(1)=0\left \{\begin{aligned} & u''(x)+u(x) = (1-\pi^2)\text{sin}\pi x & 0\leq x ...
2018-07-20 00:43:45
23688
41
原创 机器学习算法概括
前言随着时代的发展，机器学习被用到了工业、学习和生活等各个方面，它在在模式识别、自然语言处理、数据挖掘、计算机视觉等方面发挥着不可替代的作用。机器学习主要分为监督学习、无监督学习以及最近兴起的强化学习。监督学习为带标签的学习，是最常见的一部分，它以回归和分类为主。无监督学习为不带标签的学习，主要分为聚类和降维。强化学习是带奖励机制的一种学习。我们主要来介绍回归分类和聚类降维中最为核心的算法...
2018-07-08 19:13:02
1046
原创 求极小值的线搜索方法应用（SD，Newton）
求极小值的线搜索方法应用（SD，Newton） 最优化问题中，线搜索 是一种寻找目标函数f:Rn→Rf:Rn→R {\displaystyle f:\mathbb {R} ^{n}\to \mathbb {R} } 的局部最小值 x∗x∗{\displaystyle \mathbf {x} ^{*}} 的近似方法。它是最基础的迭代近似方法之一，另一种是信赖域方法。 线搜索近似首先找到...
2018-06-05 00:25:54
3909
9
原创 约束优化问题的罚函数求解方法
约束优化问题的罚函数求解方法（外罚） 惩罚函数法是求解有约束的最优化问题的一种算法。 惩罚函数法的要旨是将一个有约束的最优化问题转化为一系列的无约束问题；这些无约束问题由原问题及罚函数，再加上惩罚因子组成；而且，这些无约束问题的解会收敛于所求问题的解。算法和步骤 程序和示例写一个程序用罚函数方法求解： 计算前两部迭代的结果，也就是x1x1x_1和x2x2...
2018-06-03 20:01:25
22029
18
原创 求解非约束优化问题的拟牛顿方法（BFGS、DFP）
求解非约束优化问题的拟牛顿方法（BFGS、DFP） 拟牛顿法是一种以牛顿法为基础设计的，求解非线性方程组或连续的最优化问题函数的零点或极大、极小值的算法。当牛顿法中所要求计算的雅可比矩阵或Hessian矩阵难以甚至无法计算时，拟牛顿法便可派上用场。数学原理考虑模型问题最小化之，可得拟牛顿步长：那么，新的迭代为：问题是，这里的αkαk\alpha_k 和...
2018-06-03 19:10:57
2522
6
原创 信赖域狗腿（dogleg）方法
信赖域狗腿方法 信赖域方法（Trust-region methods）又称为TR方法，它是一种最优化方法，能够保证最优化方法总体收敛。现今，信赖域算法广泛应用于应用数学、物理、化学、工程学、计算机科学、生物学与医学等学科。相信在不远将来，信赖域方法会在更广泛多样的领域有着更深远的的发展。本文简单介绍信赖域狗腿（dogleg）方法，希望能为其他专业诸如机器学习的同学在使用优化工具时提供一个参...
2018-06-03 18:14:53
19661
8
原创 PCA算法的数学推导及其简单算例
PCA算法的数学推导及其简单算例PCA算法操作简述假设XXX是一个m*n矩阵，表示n个对象的m个特征表示数据，即每一列表示一个对象，每一行表示一个特征。我们希望将特征降为d维，d远小于m。输出结果为YYY，一个d*n的矩阵。 - 记X=[x1,x2...xn]X=[x1,x2...xn]X=[x_1,x_2...x_n]，计算每个对象点的平均值x0=1n∑i=1nxix0=1n∑i=...
2018-05-30 23:33:57
3230
原创 浅谈大数据和深度学习和计算数学的一点关系
浅谈大数据和深度学习和计算数学的一点关系专业介绍计算数学是数学的一个分支，研究的内容包括设计和分析算法以及数学建模等，目的是为了在实际工程中利用快速稳定的算法得到精确值的近似值。在计算机科学高度发展的今天，其基础计算理论的发展使计算数学进入现代化阶段。（维基百科 ）计算数学，我个人的理解是，不管是生活中，还是工程上，只要是涉及到计算，都可以划到计算数学的研究范围内，自然包括人工智能...
2018-05-05 19:45:34
17136
6
原创 国科大（UCAS）课件自动批量下载python3脚本
国科大课程网站于18年4月9日进行了换血升级，李博伟学长原有的课件下载工具不能用了（不知道他后来更新没，可以follow他的github）。特此新写一份脚本，供广大校友使用。初衷最初的想法就是因为上的课比较多，每个老师上传的课件都很多，手动下载很费事；国科大课程网站下载课件只能一个一个下载，如果有文件夹，还要先点进去，下载十分麻烦；老师上传课件不进行通知，想要知道是否有新课件...
2018-04-18 18:44:51
3625
5
原创 数据降维之局部线性嵌入（LLE）
LLE（局部线性嵌入）数据降维原理解释当数据具备某些非线性结构，如流形结构时，我们希望降维后的数据仍然保持这些结构。那么就提出了LLE降维算法。 LLE (Locally linear embedding)：在数据降维后仍然保留原始高维数据的拓扑结构，这种拓扑结构表现为数据点的局部邻接关系。此算法我们首先要寻求每个数据点的k个最近邻，然后将当前数据点用k个最近邻线性表出，那么就有相对...
2018-03-25 03:08:30
3763
2
原创 国科大教务选课系统课程余量监测python3脚本
国科大教务选课系统课程余量监测python脚本初衷:果壳同学都知道，在国科大选课，很多课还没找到就被抢光了。所以，有一个想法，可以写一个python自动选课的脚本，到了时间自动秒选课程。可是，谁能冒这个风险呢？万一出bug了呢？所以，建议选课还是自己手动去选。 每到选课时间中午12:00一过，很多想选的课程几秒钟就满了，导致没选上。在接下来的日子里，很有可能有一些人因为某一些原因退课...
2018-02-22 15:45:12
4357
原创 额，双曲守恒律
双曲守恒律的ENO格式和WENO格式问题描述考虑下述Burgers方程的初值问题 {ut+(u2/2)x=0−1≤x≤1,t>0u(x,0)=u0(x)=sin(πx)periodic B.C. for x" role="presentation">⎧⎩⎨⎪⎪
2018-01-15 23:36:03
2979
23
原创 常用守恒格式求解双曲问题（Lax-Friedrichts、local Lax-F、Roe、Engquist-Osher、Godunov、Lax-Wendroff）
Burgers方程的守恒格式求解（Proj No.1）陆嵩 计算数学问题描述我们对下述Burgers方程的初边值问题 ⎧⎩⎨⎪⎪⎪⎪ut+(u2/2)x=0u(x,0)=u0(x)=sin(πx+π)periodicB.C.forx−1≤x≤1,t>0\label{q}\left \{\begin{aligned} & u_t + (u^2/2)_x = 0 & -1\leq x
2017-12-09 03:15:46
10352
2
原创 Matlab如何制作和保存gif动图
在Mathematica当中，有Export函数可以很方便地将你绘出的多张图片合成为一张gif动图，在Matlab当中有没有类似的函数呢？当然，没有用起来那么方便的。但是Matlab可以使用imwrite函数进行gif的复写叠加形成gif动图，这也是最为普遍的一种方法。动机来源于很多时候，你想要向别人讲解或者展示某个函数的变化过程，或者要说明某个变化的物理过程（比如说流体），利用matlab的类型
2017-11-25 16:07:32
73078
50
原创 共轭梯度（CG）算法
共轭梯度（CG）方法计算数学与科学工程计算研究所 陆嵩简单介绍共轭梯度方法也是一种迭代方法，不同于Jacobi，Gauss-Seidel和SOR方法，理论上只要n步就能找到真解，实际计算中，考虑到舍入误差，一般迭代3n到5n步，每步的运算量相当与矩阵乘向量的运算量，对稀疏矩阵特别有效。 共轭梯度方法对于求解大型稀疏矩阵是很棒的方法，但是这个方法看起来总不是太靠谱。这个方法也不是越迭代精度越高，有
2017-11-16 14:27:38
84690
26
原创 SOR和SSOR迭代
SOR和SSOR迭代计算数学与科学工程计算研究所 陆嵩基本思想对于一般的大型稀疏方程组的求解，我们一般使用迭代方法求解，而古典的迭代方法中，松弛方法最优。对于Jacobi迭代，Gauss-Seidel迭代和（对称）超松弛方法迭代，我们可以简单地用矩阵分裂的思想来考虑这个问题，也可以从方程组的角度来理解这个问题。Jacobi迭代是将第i个方程中的第i个分量单独挪到一边，每个方程单独迭代的过程。Gau
2017-11-16 09:55:05
9123
1
原创 方程组求解的切比雪夫半迭代加速方法
方程组求解的切比雪夫半迭代加速方法计算数学与科学工程计算研究所 陆嵩背景介绍解方程组的迭代算法有Jacobi迭代，SOR方法等，但是对于一般的矩阵，这类算法不一定收敛，即使收敛，也有可能收敛得很慢。所以我们试图找到一个方法，来加速迭代算法的收敛速度。基本思想考虑迭代方法如下，GG为迭代矩阵 xk+1=Gxk+c\label{1} x_{k+1} = Gx_k+c 这里计算xk+1x_{k+1
2017-11-15 00:29:50
3283
5
各种课程大作业、项目等 latex 报告模板，有精美的 logo 封面、标题、摘要等。
2022-03-30
MATLAB 曲线提取分离数据生成傻瓜工具
2021-08-16
LaTeX模板.rar
2021-07-05
fast marching 方法程序.rar
2021-06-07
XMtool2.0.rar
2021-04-30
MATLAB微分方程高效解法.rar
2021-03-10
函数型数据主成分分析数据与代码全部.rar
2021-02-03
song_data.csv
2020-03-13
单位分解有限元方法(PUFEM)_PPT
2019-12-14
Fast Marching.rar
2019-12-14
双曲守恒律的ENO格式和WENO格式.zip
2019-12-14
一维激波管问题的数值求解代码.rar
2019-12-14
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人