2019年04月_LSEC小陆的博客_陆嵩

12月 11月 10月 09月 08月 07月 06月 05月 04月 03月 02月 01月

原创 MacCormack差分格式的全局误差分析

MacCormack差分格式的全局误差分析问题请分析MacCormack差分格式的全局误差：un+12=un+hun′u_{n+\frac{1}{2}} = u_n + h u'_nun+21=un+hun′un+1=12(un+un+12+hun+12′)u_{n+1} =\frac{1}{2} \left( u_n + u_{n+\frac{1}{2}} + ...

2019-04-30 22:07:13 948 3

原创曲面逼近和等参trace有限元

曲面逼近和等参trace有限元用水平集方法来表示曲面在处理演化曲面问题上用得比较多，即Γ\GammaΓ是水平集函数的一个水平集：Γ={x∈Ω:ϕ(x)=0}\Gamma=\{x \in \Omega : \phi(x)=0\}Γ={x∈Ω:ϕ(x)=0}如果用ϕh\phi_hϕh来表示ϕ\phiϕ的插值，即ϕh=Ik(ϕ)\phi_{h}=I^{k}(\phi)ϕh=Ik(ϕ)，我们...

2019-04-28 20:40:38 308

原创如何优雅地书写英文文章

如何优雅地书写英文文章文章主要内容节选自 Nick Higham 在京的一个Talk，外加自己的一些思考。为什么要好好写文章1、让你文章更容易被接收2、展现你的交流技巧3、好好地写作，写完之后看着自己的作品，颇有成就感4、其他怎么写学习和模仿你认为写得好的文章。人家用了什么样的短语和单词，你觉得好的提前记下来。搞清楚：人家如何组织语言和公式？别人是如何把文章写得如此有趣？模仿它...

2019-04-27 22:05:52 722 2

原创曲面上的高阶数值积分

曲面上的高阶数值积分简介基本思想就是将3D区域或者2D曲面上的重积分写成累次积分的方式，再用高斯积分方法进行计算。它需要找单变量非线性函数的根以及估计被积函数、指定点水平集函数及其积分。Ω\OmegaΩ是三维区域，Ωh\Omega_hΩh是它的一个四面体剖分，Γ\GammaΓ是水平集函数L(x)L(x)L(x)决定的一个曲面。TTT是选定的一个四面体。我们想要做的其实就是计算：I−=...

2019-04-26 21:47:25 515

原创 Onsager原理

Onsager原理基本原理定义势能函数（自由能）：A(a)A(a)A(a)定义能量耗散函数：Φ(a˙,a)=12∑i,jζij(a)a˙ia˙j\Phi(\dot{a}, a)=\frac{1}{2} \sum_{i, j} \zeta_{i j}(a) \dot{a}_{i} \dot{a}_{j}Φ(a˙,a)=21i,j∑ζij(a)a˙ia˙j那么系统随时间演化由...

2019-04-21 20:57:44 924 1

原创斯托克斯第一问题

斯托克斯第一问题模型建立假设有一块无限大的平板浸没在无界的静止流体中，突然以速度UUU沿其自身所在的平面运动起来，并且一直保持速度的大小和方向不变。请求解平板起动后流体运动随时间的变化过程。如图所示建立直角坐标系：考虑NS方程：ρ∂V⃗∂t+ρV⃗⋅∇V⃗=−∇p+ρg⃗+μ∇2V⃗\rho \frac{\partial \vec{V}}{\partial t}+\rho \vec{...

2019-04-18 11:45:16 906

原创函数型数据主成分分析

函数型数据主成分分析基本思想函数型主成分分析（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 4918 6

原创柱坐标系下的流体力学控制方程组的微分形式的推导

柱坐标系下的流体力学控制方程组的微分形式的推导直角坐标系下描述我们以NS方程（Navier-Stokes Equation）为例，来推导控制方程的柱坐标表示。我这里考虑不可压的，密度和粘性系数都为常数的情况。此时，直角坐标系下的NS方程的表达形式为，ρDV⃗Dt=−∇p+ρg⃗+μ∇2V⃗∇⋅V⃗=0\begin{array}{c}{\rho \frac{D \vec{V}}{D t}=...

2019-04-13 15:35:03 7799 4

雷诺输运定理的证明雷诺输运定理描述的是微分符号如何放到求导符号里面的问题，它的一个表述如下：ddt∫Ω(t)fdV=∫Ω(t)∂f∂tdV+∫∂Ω(t)(v⋅n)fdA\frac{\mathrm{d}}{\mathrm{d} t} \int_{\Omega(t)} \mathbf{f} \mathrm{d} \mathrm{V}=\int_{\Omega(t)} \frac{\partial...

2019-04-13 15:25:15 2321 2

原创 Surfice PDE常用的微分几何工具

曲面偏微分常用工具稳定界面曲面上梯度的定义曲面上的梯度，无非就是梯度在曲面切面上的一个投影，我们关心得是投影算子P\mathbf{P}P，容易想到，向量v\mathbf{v}v的投影为：Pv=v−(n⋅v)n=v−nnTv=(I−nnT)v\mathbf{Pv}=\mathbf{v}-(\mathbf{n} \cdot \mathbf{v})\mathbf{n}=\mathbf{v}-\m...

2019-04-09 22:18:04 208

原创 Auto.js的VSCode扩展

Auto.js的电脑编辑、手机同步桌面编辑器Visual Studio Code的插件。可以让Visual Studio Code支持Auto.js开发。安装在VS Code中菜单"查看"->“扩展”->输入"Auto.js"或"hyb1996"搜索，即可看到"Auto.js-VSCodeExt"插件，安装即可。插件的更新也可以在这里更新。使用Step 1按 Ctrl+Sh...

2019-04-08 13:50:04 9613 3

傅里叶谱方法相关文献和代码汇总.rar

近20篇关于傅里叶谱方法的文章书籍（包括基于MATLAB工程实现），以及相关的代码。包括Nnmerical Recies（中文版和英文版）、Matlab微分方程高效解法：谱方法原理与实现、Spectral methods in Matlab等等。讲解参考我的博客：https://blog.csdn.net/lusongno1/article/details/84259533、https://blog.csdn.net/lusongno1/article/details/84849545等。

2019-12-14

马建国数学分析课堂讲义手稿

这是马建国老师数学分析上课用的讲义，红蓝两色，收集不易，请勿外传。

2021-01-09

博弈论与机制设计.rar

博弈论与机制设计相关资料，中文，内拉哈里，Narahari，囚徒困境，游戏理论，机制设计，VCG，二价密封拍卖，内容详实，清晰，我写得很清楚。

2020-06-08

beamer简明大方英文模板.rar

这是一份 beamer 的模板，简明大方，别具一格，看起来舒心，用起来顺手。欢迎下载。是在 metro 主题上的一种改进。

2020-05-20

song_data.csv

我的博客 https://blog.csdn.net/lusongno1/article/details/81039022 中使用的数据。Python 音乐推荐的用户数据。包含 [user_id,song_id,listen_count,title,artist,song] 等标签。

2020-03-13

单位分解有限元方法(PUFEM)_PPT

单位分解有限元方法(PUFEM)介绍的PPT，参考文档见：https://blog.csdn.net/lusongno1/article/details/76599281。

2019-12-14

Fast Marching.rar

快速行进算法（fast marching）完整的运行部分和函数输入变量的说明。采用了三种Fast Marching方法，包括传统的一阶fast marching方法、二阶msfm方法、和matlab工具箱方法。文档说明可参考博客：https://blog.csdn.net/lusongno1/article/details/88409735。

2019-12-14

双曲守恒律的ENO格式和WENO格式.zip

双曲守恒律的ENO格式和WENO格式的MATLAB代码，求解的是Burgers方程的初值问题，分别使用有限体积法4阶ENO格式，有限体积法3阶和5阶WENO方法求解。时间方向用三阶TVD Runge-Kutta方法。研究格式在解光滑情况下和有间断的情况下的数值精度并作图。说明文档参考博客https://blog.csdn.net/lusongno1/article/details/79070049。

2019-12-14