数学原理
关注
分享
扫一扫
文章平均质量分 91
陆嵩
中国科学院数学与系统科学研究院 计算数学与科学工程计算研究所 科学与工程计算国家重点实验室
State Key Laboratory of Scientific and Engineering Computing,
Institute of Computational Mathematics and Scientific/Engineering Computing,
Academy of Mathematics and Systems Science,
Chinese Academy of Sciences
Email: lusong@lsec.cc.ac.cn
展开
-
MATLAB 工具箱傻瓜式求解 NS（Navier Stoke）方程
可惜的是，MATLAB 的自带的工具箱无法对它进行很好地求解，主要原因在于非线性项（对流项）无法处理，如果把它和右端载荷项放在一块处理，显然 MATLAB 工具箱是没办法支持这件事情的。打开工具箱，选择圆柱绕流算例，就是这个例子，点击 Run，它会自动地给你演示如何做。这个工具够傻瓜够强大，能求解很多 CFD 的问题，能和 FEniCS 和 Openfoam 联动，求足够完整的文档，以及完备的社区。安装就不用说，一搜就有了。需要指明的是，这份代码的运行需要在工具箱打开，变量载入空间的时候才能 work。原创 2022-10-18 13:05:46 · 5372 阅读 · 3 评论 -
核酸检测多少人为一组混检合适？
但是对于被检测者来说，不管是混检还是单检，做一次核酸真的就是要付出做一次核酸的成本，排一次队，刮一次鼻子，遭一次罪。有人说，可以再采用的时候，把每个人的样本分两份，当检测出患病的时候，再把备份拿出来单独测，这样就不需要喊人来车第二次了。考虑到患病率极低，这种分两份的人工成本是极高的，已经远远超过了把一起混检的人喊来再测一次的成本，更远远超过了不采用最优参数。其实很简单，2个管子，6 个人，每个管子 3 个人，2 个患病的人进入同一个管子的概率不是 1/2，而是。是一件非常麻烦的事情，不是初等数学可以做的。原创 2022-10-14 10:13:32 · 15623 阅读 · 22 评论 -
MATLAB Tensor 和 N-Way 工具箱张量补全功能评测
既然现有的方法这么垃圾，我们就基于比较快的分解方法，做一个改进。一个小改进就是在做 cp_als 最小二乘交替步的时候，直接在已知数值的位置上，每步迭代用已知值去替代它。Nway 工具箱实现的非常 stupid，使得速度非常慢，为什么 Nway 实现得这么 stupid，举个例子，它算一堆矩阵 KR 除某个之外的乘积的时候，一个一个算的，一点都不明智。ALS 方法就是先固定 B 和 C 以求解 A，然后再固定 A 和 C 求解 B，接着固定 A 和 B 以求解 C，并继续重复整个过程，直到收敛标准被满足。原创 2022-10-11 13:16:58 · 1954 阅读 · 2 评论 -
欠定方程组的最小范数解
当然，这些条件都有一些奇奇怪怪的名字，比如“互补性条件”之类的，可以不用管。当然，如果只有等式约束，那么其实只有前两条，很简单，也很好记。对于这个问题，我们知道这里未知量的个数比方程的个数多，也就是解是存在不唯一的。了，如果是低于三维的情况，高中的知识很快就能解出来了。除了上述的分裂法推导外，我们也可以用拉格朗日乘子法得到这个最小范数解。容易证明，它也是满足最小二乘误差相同的最小范数解。几何的提法是，超平面。是可逆的，我们有更加漂亮的写法，从第一块我们可以得到。就是最小范数解，它极小化了如下的优化问题，..原创 2022-08-26 14:14:48 · 4588 阅读 · 7 评论 -
MATLAB 张量工具箱的使用
中的所有列都拿出来，做 Kronecker 内积之后求和，求和的系数来自于张量$ \mathcal{G}$ 的。运行的结果上可以看出，原来缺失的位置补全的效果可能不是非常地好。MATLAB 的 Tensor Toolbox 是一个强大的张量工具，它提供了一些强大的张量函数。是完全的，没被噪声污染的数据，我们称之为精确解，用上述的张量形式存储。给出了补全矩阵和原来的矩阵的一个误差表示，它不仅仅考虑缺失的部分。表示的是 CP 分解求和符号的上标，也是因子的个数。的 分量，并且可以执行许多操作，例如。......原创 2022-08-04 08:12:29 · 7902 阅读 · 2 评论 -
关于向量值函数方程变分的一点注记
一般来说，我们喜欢把梯度定义为列向量，一列一列地来，但是上述的关于向量值函数的梯度，却是对每个分量求梯度，然后按行排的，这是反直觉的，所以导致了这个地方的。有人说，自己做的都是标量解函数的偏微分方程变分，当遇到解函数是向量值函数时，就不知道有限元空间怎么取，变分也不知道怎么做了。再 一个对于矩阵的散度和梯度算子，是支持分块矩阵的操作的。有了以上的准备，我们举个简单的例子来看向量值解的方程的变分，头上加转置变回原来的定义方式，就得到了最开始的上述的表达。的地方，是可以随意进出梯度和散度算子的，即。...原创 2022-08-03 04:51:48 · 1062 阅读 · 2 评论 -
模型降阶方法之张量方法应用举例
TABLE5.3-5.5不同的参数离散、不同的reducedbasis维数、不同的投影方式下，不同方法的相对增益，相对增益指的是相对于POD而言的相对误差，FIG5.3说的是样本外预测的投影误差随着reducedbasis的个数、参数的采样间隔、张量分解精度的变化而变化的情况。同样地，用某种格式，比如CN格式得到snapshots，之后用张量方法得到reducedbasis，再对原方程进行降阶，得到了，用标准的有限元方法，在空间上做离散，就得到这样一个时间连续的方程，...原创 2022-07-29 06:04:32 · 920 阅读 · 0 评论 -
模型降阶方法之张量方法
POD 方法的介绍如下所示：https://blog.csdn.net/lusongno1/article/details/125944587经典的 POD 方法无法处理非线性项，所以发展了很多的方法，比如 DEIM：https://blog.csdn.net/lusongno1/article/details/125955245如果我们要考虑的方程里面有一组参数，不同的参数都能获得一个 Matrix，那么我们应该如何把多组参数下的得到的数据放在一块作为先验，来给新到来的一组参数做模型降阶呢？这就是这篇原创 2022-07-27 13:49:49 · 1200 阅读 · 0 评论 -
模型降阶方法之 DEIM
为什么叫插值，为什么叫greedy，可以看得出来，我们选的就是投影之后，距离最大的一个点，然后把它固定。也就是说，我们对于$\mathbf{F}$选行，我们只要计算里面必要分量作用即可，也就是可以把选行操作放到$\mathbf{F}$里面去。POD方法固有的局限性，使得无法很好地处理非线性项，所以必须有一个合适的方法处理非线性项。这里看起来，似乎还需要把非线性项施加于每个分量算出来，其实并不然，这里的。同理，对于稳态相关的问题，我们需要计算其雅克比，非线性项的雅克比矩阵为，...原创 2022-07-24 02:52:18 · 3667 阅读 · 4 评论 -
模型降阶方法之 POD
我们想要找到解所在的超平面，那么事实上，我们只要找到这个超平面上足够多个“典型”的点，那么事实上，我们就获得了平面的信息。在更高维的空间中，一般的ODEs的解的分布区域也不是充满整个空间的，它是位于一个特定的平面或者流形上的。这揭示一种共性，即虽然是解空间是一个三维的空间，但是他们的解都位于一个二维的超平面里面，也就是说，用两个变量，就可以把三维的向量值解表达出来了。那么问题来了，上述的两个例子很简单，我们用肉眼都能观察到方程之间的关系，从而进行降阶，对于更一般的问题，我们应该怎么处理呢？...原创 2022-07-23 10:39:27 · 6969 阅读 · 0 评论 -
曲面曲率驱动的斑图传播
前几天 PRL 杂志发布了一篇文章，是关于图灵斑图的。大概意思是说，对于在曲面的上斑图演化，受到曲面的形状的影响。这个 Letter 首次报告了在平面上的静态图灵斑图在曲面上是会传播的，而我们以往的研究里，都是假设图灵斑图在曲面上也是静态的。这篇文章研究了反对称曲面上的反应扩散系统。结果表明，曲面的对称性和斑图的样子都会影响斑图的传播。作者研究了曲面的曲率对斑图传播影响的一个机制。一些研究已经表明了，当图灵斑图放到曲面上时，不稳定性条件是如何变化的，并且图案的位置如何受到曲面曲率的不均匀性的影响，这被称为原创 2022-06-12 13:12:38 · 626 阅读 · 0 评论 -
曲面抛物方程的质量集中法
曲面抛物方程的质量集中法该篇内容来自文章 Xiao X, Feng X, Yuan J. The lumped mass finite element method for surface parabolic problems: error estimates and maximum principle[J]. Computers & Mathematics with Applications, 2018, 76(3): 488-507.部分内容来自张林波老师 slides。特此说明。文原创 2022-05-29 10:08:04 · 1100 阅读 · 0 评论 -
共轭梯度法（CG）详解
共轭梯度法（CG）详解文章目录共轭梯度法（CG）详解线性共轭梯度法共轭方向共轭方向法CG 方法预条件非线性共轭梯度法FR 方法其他非线性 CGPR+ 方法重启动之前写过几个关于共轭梯度法的注记，譬如：https://blog.csdn.net/lusongno1/article/details/78550803https://blog.csdn.net/lusongno1/article/details/68942821但事实上很多人反应，看得一头雾水，基于此，本篇文章旨在对于共轭梯度方法从优原创 2022-05-26 17:12:54 · 9502 阅读 · 2 评论 -
特征值问题误差分析：Babuska–Osborn紧算子谱逼近理论
特征值问题误差分析：Babuska–Osborn紧算子谱逼近理论文章目录特征值问题误差分析：Babuska–Osborn紧算子谱逼近理论从一个简单的例子开始一维例子P1 元逼近一些观测和思考一些基本定义变分特征值问题特征值问题特征值刻画误差的刻画Babuska-Osborn 理论收敛性通常的框架收敛阶对称变分特征值问题上午好，各位。非常荣幸站在这里。且，今天我想给大家分享一些关于特征值问题的误差分析。关于我们的讨论班，我提议用说英文。就是，关于我们的讨论班，我提一个建议。我觉得大家可以英文来讲。事实上原创 2022-05-18 22:57:37 · 612 阅读 · 0 评论 -
紧算子的谱逼近理论及其在拉普拉斯特征值问题上的应用
紧算子的谱逼近理论及其在拉普拉斯特征值问题上的应用文章目录紧算子的谱逼近理论及其在拉普拉斯特征值问题上的应用紧算子的谱逼近理论预备知识收敛阶对称变分特征值问题拉普拉斯特征值问题简述T=TVT=T_VT=TV 时的分析T=THT=T_HT=TH 时的分析紧算子的谱逼近理论预备知识V1V_1V1 和 V2V_2V2 是复希尔伯特空间，我们需要找 λ∈C\lambda \in \mathbb{C}λ∈C 和 u∈V1u \in V_{1}u∈V1 且 u≠0u\neq 0u=0，使得，a(原创 2022-05-18 22:53:58 · 887 阅读 · 0 评论 -
密度泛函理论平面波基组展开
密度泛函理论：平面波基组展开文章目录密度泛函理论：平面波基组展开基本理论与概念原胞绝热近似、 HF 方法DFT 与 KS 方程自洽迭代交换关联势赝势DFT 应用赝势平面波 DFT 方法各项能量表达Fourier 展开动量空间总能量及发散问题动量空间 KS 方程截断能平面波 DFT 方法概述参考书基本理论与概念密度泛函理论 (Density Functional Theory, DFT) 是一种研究多电子体系电子结构的量子力学方法。密度泛函理论在物理、 化学、材料上都有广泛的应用, 特别是用来研究分子和原创 2022-05-18 22:44:05 · 2128 阅读 · 0 评论 -
DFT 平面波方法
DFT 平面波方法文章目录DFT 平面波方法变分法DFT 平面波方法DFT 平面波迭代解法abacus 源码阅读变分法考虑薛定谔方程，H^ψ=Eψ\hat{H} \psi=E \psiH^ψ=Eψ这里的波函数我们一般不能精确得到，所以我们一般找一个数学上可以处理的函数去逼近，ϕ≈ψ\phi \approx \psiϕ≈ψ特征值理论告诉我们，极小化能量泛函E~=⟨ϕ∣H^∣ϕ⟩⟨ϕ∣ϕ⟩\tilde{E}=\frac{\langle\phi|\hat{H}| \phi\rangle}原创 2022-05-18 22:29:40 · 970 阅读 · 0 评论 -
小议移动有限元方法（转载）
小议移动网格有限元方法（转载）转载一篇提到移动网格有限元的文章。网格是偏微分方程数值解法的基础，网格体系的好坏直接影响计算结果的精度，甚至影响计算的成败。网格方法的研究经历了从结构化到非结构化，从单一网格到混合网格的过程。经过几十年的发展，这些网格方法已经很好地用于各种问题的计算，并不断出现新的针对不同情况的网格生成技术，而且形成了一些好的网格生成软件。近三十年来，自适应网格方法（主要有移动网格方法和局部细化或粗化的网格方法）一直引起了国际学术界和各类应用部门的高度重视，并且成为网格方法研究的热点问原创 2022-05-18 22:10:48 · 304 阅读 · 0 评论 -
电子结构计算简介
电子结构计算简介文章目录电子结构计算简介问题背景基本概念密度泛函理论方法波函数方法最基础的 Hartree-Fock 算法薛定谔方程Hartree-Fock 方法软件包求解薛定谔方程是定义在 R3N\mathbb{R}^{3 N}R3N 高维线性特征值问题. 除极少数体系外, 无法直接求解. 因此, 人们需要寻找与其等价或简化的可计算模型. 这些等价或简化的可计算模型中最典型的是 Hartree-Fock 方程和 Kohn-Sham 方程. 求解这两种方程，其对应了波函数方法和密度泛函理论方法。前原创 2022-05-04 23:08:28 · 2953 阅读 · 0 评论 -
循环制比赛要赢几场可能（一定）晋级
循环制比赛要赢几场可能（一定）晋级最近，世界杯正在如火如荼地进行。它的晋级规则胜平负分别积3、1、0分。每个小组积分的前两名球队出线进入淘汰赛阶段的1/8决赛，共16支队，即“16强”。积分相同的接下来看总净胜球数，还相同看总进球数，再相同只能抽签。16强后是对分PK剩下一半，一直选出冠军和亚军，前四强中的第三名和第四名要额外比一场，选出季军。今天，我一直在想一个问题。在小组中，要赢了几场才能保证出线？有赢几场才有可能出线？显然，4 进 2 的情况比较简单。下面考虑一般情况。假如有 n 支队伍，两原创 2022-03-28 21:09:46 · 3151 阅读 · 0 评论 -
《剑指 Offer（第 2 版）》系列刷题
剑指offer刷题文章目录剑指offer刷题03 [数组中重复的数字]04 [二维数组中的查找]05 [替换空格]06 [从尾到头打印链表]07 [重建二叉树]0原创 2022-03-04 04:30:00 · 670 阅读 · 0 评论 -
关于（广义）代数特征值问题的一点注记
关于（广义）代数特征值问题的一点注记感谢王同学和王同学提供的一些信息和资料。文章目录关于（广义）代数特征值问题的一点注记常用特征值计算方法概述对称特征值问题非对称特征值问题Krylov 子空间方法对于大规模问题小小的总结常用特征值计算方法概述对称特征值问题Jacobi 迭代 、Rayleigh 商迭代 、对称 QR 迭代方法、 分而治之法 、对分法和反迭代法对称 QR 迭代：对称矩阵的 QR 方法，首先将矩阵约化为 Herssenberg 矩阵，然后在求解该矩阵的特征值。对于对称三角原创 2022-02-26 22:06:00 · 2523 阅读 · 0 评论 -
【深度学习入门基础】一、从线性代数和微积分的角度看神经网络
【深度学习入门基础】从深度学习和微积分的角度看神经网络这是深度学习入门系列文章，我们企图用最简洁的语言、最干净的表达，让读者快速获取到他所想要的。本系列文章持续更新。一些网上能查到的基础知识，诸如激活函数有哪些等，不再介绍。导言大多数介绍深度学习的资料，一开始就从生物神经元、轴突、树突、刺激等等讲起，本来很简单的东西，却被介绍得云里雾里。本文将从矩阵乘法和导数的角度入手，让你一下子 get 到什么是神经网络。矩阵乘法看神经网络问题是这样一个问题：假设有一个列向量 x\mathbf{x}x，希望原创 2022-02-26 21:57:33 · 1132 阅读 · 0 评论 -
过年回家和小朋友玩什么
过年回家和小朋友玩什么文章目录过年回家和小朋友玩什么尼姆游戏的必胜策略尼姆游戏简介初体验玩法必胜法反推二进制如何保持平衡态二者区别必胜法记忆反推法平衡态记忆二进制心算的作法尼姆类似游戏减法游戏多排尼姆游戏循环尼姆Grundy 游戏尼姆游戏的必胜策略三堆数量随机的石头，两个人轮流从其中一堆中拿走任意数量的石头，拿走最后一块的人输，必胜法是什么？尼姆游戏简介初体验三年级的时候，来探亲的姑丈和用扑克牌玩过一个游戏，大概意思的是：三五七张数的三垛牌，每人每次轮流从一垛牌里面拿不少于一张牌，谁拿到最后原创 2022-01-09 22:14:30 · 2890 阅读 · 0 评论 -
水平集和符号距离函数
水平集和符号距离函数零水平集定义： 对于一个函数 ϕ(x⃗)：Rn→R\phi(\vec x)：{\mathbf{R}^n}\rightarrow \mathbf{R}ϕ(x)：Rn→R（其中x⃗∈Rn\vec x \in {\mathbf{R}^n}x∈Rn ，下同），取其值域为零部分对应的定义域：Γ={x⃗∣ϕ(x⃗)=0}(1)\Gamma = \{ \vec x|\phi (\vec x) = 0\} \tag{1}Γ={x∣ϕ(x)=0}(1)这里， Γ∈Rn−1\Gamma \in原创 2021-12-19 15:27:33 · 3658 阅读 · 1 评论 -
图灵与图灵斑图
图灵与图灵斑图文章目录图灵与图灵斑图图灵简介图灵成就与贡献战争中的英雄计算机科学之父人工智能之父图灵斑图同性恋、迫害与平反怪人和运动天赋英镑和电影图灵斑图（Turing Pattern）自然和物理现象形成机理形而上的理解反应扩散系统形成条件图灵斑图的数值方法之参数化有限元方法简介演化表面上反应扩散方程的推导曲面梯度演化表面上的反应扩散系统模型说明曲面有限元方法变分形式演化曲面有限元方法时间离散肿瘤增长模型建立数值方法数值结果网格自适应和重新划分网格图灵斑图的 Trace 方法图灵简介图灵成就与贡献原创 2021-12-05 17:01:18 · 3430 阅读 · 0 评论 -
有限元分析中的常识（持续更新）
有限元分析中的常识（持续更新）介绍一些学过有限元的人都不得不超级熟练掌握的基本常识。文章目录有限元分析中的常识（持续更新）通用符号基本不等式Cauchy-Schwarz 不等式Hölder's 不等式Young's 不等式Poincaré 不等式离散版本的 Gronwall 不等式Sobolev 不等式嵌入定理嵌入定理紧嵌入定理Cea 引理证明技巧Aubin-Nitsche 对偶技巧Lax-Milgram 定理通用符号∥⋅∥m,p\|\cdot\|_{m, p}∥⋅∥m,p 表示 Sobol原创 2021-11-20 19:45:55 · 3008 阅读 · 0 评论 -
梯度流是个什么玩意儿
梯度流是个什么玩意儿举个最最最最简单的例子。考虑一个经典的极小化问题：minxf(x)\min _{\mathbf{x}} f(\mathbf{x})xminf(x)其中 x∈Rn,f:Rn→R\mathrm{x} \in \mathbb{R}^{n}, f: \mathbb{R}^{n} \rightarrow \mathbb{R}x∈Rn,f:Rn→R 是一个二次可微的函数。x∗\mathrm{x}^{*}x∗ 是上述问题的解的必要条件是∇f(x∗)=0\nabla f\left原创 2021-11-15 10:40:26 · 6071 阅读 · 3 评论 -
肿瘤增长数学模型
曲面偏微分方程的一个实际应用——肿瘤增长文章目录曲面偏微分方程的一个实际应用——肿瘤增长简介演化表面上反应扩散方程的推导曲面梯度演化表面上的反应扩散系统模型说明曲面有限元方法变分形式演化曲面有限元方法时间离散肿瘤增长模型建立数值方法数值结果网格自适应和重新划分网格简介基本想法是利用表面有限元数值方法求解演化曲面上的反应扩散方程。不断增长的生物表面上的模式形成，是通过求解曲面上的反应扩散方程得到的，可以参考图灵的经典文章 1。这种方程表现出空间均匀结构的扩散驱动不稳定性，导致空间不均匀的图案。这项工作的原创 2021-11-06 00:00:00 · 3465 阅读 · 3 评论 -
机器学习基础知识点
机器学习基础知识点文章目录机器学习基础知识点监督学习回归线性回归岭回归lasso回归分类k最近邻分类朴素贝叶斯分类logistic回归支持向量机其他随机梯度下降线性判别分析决策树无监督学习聚类k均值分层次聚类谱聚类高斯混合模型降维PCA降维LLE降维MDS和t-SNE独立成分分析等距特征映射谱嵌入深度学习（一些基本概念）一维卷积（离散）二维卷积（离散）padding反卷积（转置卷积）学习率反向传播正则化softmaxdropout历史随机梯度下降（mini-batch）万有逼近定理几种激活函数白化动量方法原创 2021-10-24 14:41:14 · 5067 阅读 · 0 评论 -
三万字长文带你了解电子设计自动化（EDA）
三万字长文带你了解电子设计自动化（EDA）文章目录三万字长文带你了解电子设计自动化（EDA）EDA 的工业视角简介物理设计不断发展的设计环境早期物理设计：交互式支持自动设计路由布局布线中的关键支持子系统经济影响较小的物理设计工具关于物理设计EDA的未来VSLI 设计的仿真和验证事件周期模型编译代码硬件描述语言测试平台语言加速和仿真事件驱动的加速器基于周期的加速器和仿真器大规模并行通用处理器基于 FPGA 的仿真器分析模型检查器等效检查器未来：SoC综合开端逻辑优化和设计编译器RTL 综合顺序映射和优化扩大原创 2021-10-11 22:13:52 · 5351 阅读 · 5 评论 -
零基础使用 MATLAB 求解偏微分方程（建议收藏）
零基础使用 MATLAB 求解偏微分方程（建议收藏）文章目录零基础使用 MATLAB 求解偏微分方程（建议收藏）偏微分开源工具介绍PDE 工具箱函数汇总介绍0 基础：GUI 界面操作示例问题工具箱求解导出为代码形式代码导出相关数据0.1 基础：编程调用 PDE 工具箱偏微分开源工具介绍百分之九十以上的重要的工程和数学科学研究，和偏微分方程都脱不开关系。在所有的偏微分方程中，百分之九十九都是没有解析解的。没有解析解怎么办，我们只能通过有限元或者有限差分等方法，求解偏微分方程数值解。如果您有一些代码基础，原创 2021-10-02 01:30:07 · 45542 阅读 · 4 评论 -
匹配算法告诉你为什么要找女（男）朋友一定要主动？
稳定匹配：该如何给男孩女孩们配对？文章目录稳定匹配：该如何给男孩女孩们配对？定义算法：延迟接受（盖尔沙普利算法）越主动配偶越好代码实现对于恋爱的启示纯文字描述稳定匹配（stable matching）。定义有这样一个问题，大概意思说是给一群男的和一群女的配对，每个男的对所有女的喜欢程度都有一个排序，所有女的对所有男的喜欢程度也有一个排序，男女人数相等。我们给一个不稳定对（blocking pair）的定义。假设小明的女朋友更喜欢小东，比起小东自己的女朋友，小东更喜欢小明的女朋友，这时候，（小明原创 2021-09-27 19:04:57 · 1345 阅读 · 3 评论 -
LeetCode 刷题汇总（第一阶段，随机游走）
从 0 开始 LeetCode刷题面临毕业，为了谋求一份工作，不得已走上了刷题之路。不追求解法的性能和代码的优雅，只追求在最短的时间内写出来。很多人喜欢给各种各样的方法起一个奇奇怪怪的名字，比如说“回溯”、“滑动窗口”等等，不懂这些，也可以做呀。你最后自己想了一个方法做出来了，其实默默地也就用了这些所谓的方法。这些题目我做完了，别人和我聊起什么回溯，我也不懂什么意思。我只知道，暴力出奇迹，干就完了。我们做计算数学的人的长处不在于刷题，奈何这个世界这么无奈。其实，我一直在思考，用这些小算法题去考验一原创 2021-09-15 14:07:28 · 742 阅读 · 0 评论 -
xin片设计的中的数学问题
xin片设计的中的数学问题（为了过审只能用拼音替代）文章目录xin片设计的中的数学问题（为了过审只能用拼音替代）概述EDA 简介EDA 的主要工具和相关数学问题数字全流程工具数学问题定制电路全流程工具数学问题晶圆制造工具数学问题写在后面的话概述EDA 与光刻机均被称为“xin片之母”，一个是xin片设计的必备工具，一个是xin片制造不可或缺的“神器”。恰恰在这两个关键领域，国产 EDA 及光刻机全面落后于国际大厂。目前最先进的 EUV 光刻机被荷兰阿斯麦尔垄断。国内最先进的光刻机制程仅有 90nm原创 2021-09-12 20:40:09 · 748 阅读 · 0 评论 -
MATLAB 提取图片曲线
文章目录我用 MATLAB 提取图片曲线行文动机图像的读入与裁剪颜色拾取颜色转换与色差计算分离曲线二值化，提取数据数据点分类与排序后话我用 MATLAB 提取图片曲线给你一张图片，如何提取里面曲线的数据，从而利用这些数据进行图像重绘、加工处理、测距、拟合得到函数表达式等操作呢？行文动机前段时间，有个朋友问了我一个问题，大概意思就是要给图像的流线测距离，在我的印象里面，MATLAB 是似乎没有这种直接的功能的。那么换个角度来理解一下这个问题，如果给你一张图像，如何提取里面点的数据？其实，有了曲原创 2021-08-16 23:09:19 · 11675 阅读 · 20 评论 -
Pinsker 不等式的简单证明
Pinsker 不等式的简单证明网上有很多很多关于 Pinsker 不等式的证明方法，但是我没有看到一个用数学归纳法证明的，也没有看到一个不加先验定义的自包含的证明。下面我给出一个关于一个极简的证明。任何的引用请注明本出处。Pinsker 不等式请证明如下不等式：∑i=1nailnaibi≥∑i=1n(ai−bi)2\sum_{i=1}^{n}a_i \ln \frac{a_i}{b_i}\geq \sum_{i=1}^{n}(a_i-b_i)^2i=1∑nailnbiai≥i=1原创 2021-08-01 22:26:10 · 1554 阅读 · 2 评论 -
学数学有个球用？
学数学有个球用？我们学数学的，不能有丝毫的骄傲，却要有无比的自信。仅以此文致敬中国科学院数学与系统科学研究院的 2021 届毕业生们。文章目录学数学有个球用？如何调好收音机AFBT 算法收敛性分析梯度下降方法（GD）Gradient Descent步长选取收敛率学好了数学，你能找一份工作，找到了一份不错的工作，你就有钱了，你有钱了，就可以买各种各样的球，足球篮球羽毛球，网球棒球羽毛球，手球水球曲棍球，桌球排球扔链球。结论是，学数学当然有球用，想要什么球都有，想要多少有多少。学数学不仅有球用，还有原创 2021-07-03 20:48:51 · 816 阅读 · 2 评论 -
【科普】有趣“小学”数学题，做出一道即可成名（持续补充）
【科普】有趣“小学”数学题，做出一道即可成名，赶紧让你孩子想想（持续补充）小学数学老师问了我一道数学题，“任何大于2的偶数都可以写成两个质数的和”，“对吗”，“怎么论证”。我吃惊于这种著名的猜想普及率却如此值低，故做此文，普及数学中的一些著名猜想和定理。文章目录【科普】有趣“小学”数学题，做出一道即可成名，赶紧让你孩子想想（持续补充）哥德巴赫猜想费马大定理（亦名费马最后定理）孪生素数猜想3n+1 猜想四色定理黎曼猜想P=NP?P 问题NP 问题NPC 问题NP-Hard 问题NP=P？问题庞加莱猜原创 2021-06-27 19:32:45 · 5553 阅读 · 2 评论 -
反赌科普宣传——赌徒的谬误
反赌科普宣传——赌徒的谬误文章目录反赌科普宣传——赌徒的谬误游戏简单介绍不公平的游戏彩民的谬误求神拜佛，是邪？非邪？隐形的陷阱，无底的深渊在赌博中，有一种“荷兰赌式”的“倍赌”策略，即输了加赌注加倍，直到赢了为止。似乎在所有人看来，只要本金足够，似乎此法无往而不利。那么这种做法可不可取？背后潜藏的玄机又在哪呢？我们从一个简单的游戏说起。所用的概率论知识，一个正常上过大学的理工科的人都应该会懂。游戏简单介绍在我们的家乡广泛流行着一个“机会游戏”，叫做“六合彩”，其游戏规则如下：人的年龄（虚岁）对原创 2021-06-13 11:01:25 · 4466 阅读 · 11 评论