- 博客(410)
- 资源 (17)
- 收藏
- 关注
原创 曲面上的微积分(二)
曲面上的微积分(二)从上一次讨论结束,回去之后,我花了一点时间走了一遍证明。心中依然遗留了两个问题。对于Laplace-Beltrami方程的弱形式,∫γφ~Δγv~=−∫γ∇γφ~⋅∇γv~\int_{\gamma} \tilde{\varphi} \Delta_{\gamma} \tilde{v}=-\int_{\gamma} \nabla_{\gamma} \tilde{\varph...
2019-09-18 21:46:07 534
原创 曲面上的微积分(一)
曲面上的微积分(一)参数化方法、迹方法和窄带方法是三种主要方法去求解Laplace-Beltrami问题。首先,我会介绍一些关于曲面上微积分的知识。参数化曲面定义曲面的参数化表示,χi:Vi→Ui∩γ⊂Rn+1\chi_{i} : \mathcal{V}_{i} \rightarrow \mathcal{U}_{i} \cap \gamma \subset \mathbb{R}^{n+1}...
2019-09-11 22:24:30 2303 2
原创 linux下.bashrc的使用和个人用户家目录下的软件安装
linux下.bashrc的使用和个人用户家目录下的软件安装.bashrc文件简介:环境配置.bashrc文件主要保存个人的一些个性化设置,如命令别名替代、设置系统的环境变量(PATH)等。在公用的集群环境下,个人用户登录到该账户下的家目录,此只对这个用户的个性化设置起效。它们是一些隐藏文件,ls -a可查看。可以看到有四个相关的文件:.bashrc:当你登入shell时以及每次打开新的...
2019-08-29 00:06:41 1789
原创 linux下.bashrc的使用和个人用户家目录下的软件安装
linux下.bashrc的使用和个人用户家目录下的软件安装.bashrc文件简介:环境配置.bashrc文件主要保存个人的一些个性化设置,如命令别名替代、设置系统的环境变量(PATH)等。在公用的集群环境下,个人用户登录到该账户下的家目录,此只对这个用户的个性化设置起效。它们是一些隐藏文件,ls -a可查看。可以看到有四个相关的文件:.bashrc:当你登入shell时以及每次打开新的...
2019-08-28 23:15:54 3705
原创 驱动方腔流SIMPLE方法
驱动方腔流SIMPLE方法问题描述:二维驱动方腔流请在矩形区域上数值地求解不可压NS方程,你的报告将包含但不仅限于以下的讨论:问题描述和数值格式格式的时空精度分析比较非线性项处理的守恒格式和非守恒格式的不同不同雷诺数的影响(分别取100、1000、3000)数值格式:SIMPLE格式问题的数学表达考虑二维空间的不可压NS方程:ut+px=−(u2)x−(uv...
2019-06-15 17:05:41 9818 14
原创 曲面上变分的基础知识(方程组)
曲面上变分的基础知识(方程组)定理:曲面上变分的格林公式(向量形式)定理描述如下,为了方便,我们省略了严格的数学书写:∫ΓdivΓ(GP)⋅gds=−∫Γtr(G∇Γg)ds\int_{\Gamma} \operatorname{div}_{\Gamma}(\mathbf{GP} ) \cdot \mathbf{g} d s=-\int_{\Gamma} \operatorname{t...
2019-06-02 20:51:45 1175
原创 ubuntu安装搜狗输入法
安装简单,做个记录。到搜狗官网https://pinyin.sogou.com/linux/?r=pinyin下载搜狗输入法linux版本,ubuntu18的版本下载64 bit的。打开ubuntu的应用商店,搜索fcitx,将搜到的fcitx程序也就是小企鹅图标的全部安装上。安装完成后,我们再双击刚下载的deb程序文件就可以安装了。搜狗输入法安装完毕后我们打开命令行,输入:sudo ap...
2019-05-18 09:33:36 500
原创 基于cmake文件的Code::Blocks的简单使用
基于cmake文件的Code::Blocks的简单使用motivation在linux环境下可以使用vi+make的方式来书写和编译运行C++的代码,但是对于一个极大的工程文件而言,我们需要一个IDE。Code::Blocks小巧而“专情”,是一个在linux环境下完成C++工程的一个不错的选择。给你一个cmake工程文件,通过cmake获得到了编译运行代码所需的Makefile文件,Cod...
2019-05-15 21:51:40 3239
原创 差分格式的傅里叶稳定性分析方法
差分格式的傅里叶稳定性分析方法以下格式的傅里叶稳定性分析ujn+1−ujnΔt=θνuj+1n+1−2ujn+1+uj−1n+1Δx2+(1−θ)νuj+1n−2ujn+uj−1nΔx2\frac{u_{j}^{n+1}-u_{j}^{n}}{\Delta t}=\theta \nu \frac{u_{j+1}^{n+1}-2 u_{j}^{n+1}+u_{j-1}^{n+1}}{\Delt...
2019-05-07 21:47:18 7972 2
原创 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 2243 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 882 4
原创 如何优雅地书写英文文章
如何优雅地书写英文文章文章主要内容节选自 Nick Higham 在京的一个Talk,外加自己的一些思考。为什么要好好写文章1、让你文章更容易被接收2、展现你的交流技巧3、好好地写作,写完之后看着自己的作品,颇有成就感4、其他怎么写学习和模仿你认为写得好的文章。人家用了什么样的短语和单词,你觉得好的提前记下来。搞清楚:人家如何组织语言和公式?别人是如何把文章写得如此有趣?模仿它...
2019-04-27 22:05:52 1468 2
原创 曲面上的高阶数值积分
曲面上的高阶数值积分简介基本思想就是将3D区域或者2D曲面上的重积分写成累次积分的方式,再用高斯积分方法进行计算。它需要找单变量非线性函数的根以及估计被积函数、指定点水平集函数及其积分。Ω\OmegaΩ是三维区域,Ωh\Omega_hΩh是它的一个四面体剖分,Γ\GammaΓ是水平集函数L(x)L(x)L(x)决定的一个曲面。TTT是选定的一个四面体。我们想要做的其实就是计算:I−=...
2019-04-26 21:47:25 1617
原创 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 2711 2
原创 斯托克斯第一问题
斯托克斯第一问题模型建立假设有一块无限大的平板浸没在无界的静止流体中,突然以速度UUU沿其自身所在的平面运动起来,并且一直保持速度的大小和方向不变。请求解平板起动后流体运动随时间的变化过程。如图所示建立直角坐标系:考虑NS方程:ρ∂V⃗∂t+ρV⃗⋅∇V⃗=−∇p+ρg⃗+μ∇2V⃗\rho \frac{\partial \vec{V}}{\partial t}+\rho \vec{...
2019-04-18 11:45:16 2639 1
原创 柱坐标系下的流体力学控制方程组的微分形式的推导
柱坐标系下的流体力学控制方程组的微分形式的推导直角坐标系下描述我们以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 25970 8
原创 雷诺输运定理的证明
雷诺输运定理的证明雷诺输运定理描述的是微分符号如何放到求导符号里面的问题,它的一个表述如下: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 6646 5
原创 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 729
原创 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 14714 4
原创 keras入门实例:非线性拟合求拟合系数
目标要做的事情非常简单,就是用Keras去拟合函数y=x(1−x)(a0+a1x+a2x2)y=x(1-x)(a_0+a_1x+a_2x^2)y=x(1−x)(a0+a1x+a2x2)并求其系数a0,a1,a2a_0,a_1,a_2a0,a1,a2。数据已经给你了,放在csv文件中,没有header,两列多行,很轻易就能读进来。乍一看,神经网络逼近能力强,随便搭个网络都行。但请注意...
2019-03-31 17:14:37 3897 6
原创 Keras天坑:想当然的对层的直接运算带来的问题
天坑keras如何操作某一层的值(如让某一层的值取反加1等)?keras如何将某一层的神经元拆分以便进一步操作(如取输入的向量的第一个元素乘别的层)?keras如何重用某一层的值(如输入层和输出层乘积作为最终输出)?这些问题都指向同一个答案,即使用Lambda层。另外,如果想要更加灵活地操作层的话,推荐使用函数式模型写法,而不是序列式。Keras当中,任何的操作都是以网络层为单位,操作的实...
2019-03-31 16:24:04 2708 2
原创 auto.js小记:蚂蚁森林定时自动收取能量脚本
简介Autojs是一个支持无障碍服务的Android平台上的JavaScript IDE,其发展目标是JsBox和Workflow。同时有VS Code 插件可提供基础的在桌面开发的功能。下载地址:酷安官方文档:https://hyb1996.github.io/AutoJs-Docs/特性由无障碍服务实现的简单易用的自动操作函数悬浮窗录制和运行更专业&强大的选择器API,...
2019-03-26 22:04:35 18551 4
原创 同时连接了无线和有线,如何查看当前连接的是哪个
不论是win10还是win7,一个比较通用的方式就是查看跃点数。输入route print,打印路由信息。简单的原则就是:跃点数越小,网络优先级越高。有线和无线对应的网关和接口怎么看?找到“更改适配器选项”,进入网络连接设置。右击相应网络,查看状态。点击详细信息,即有:我们也可以修改跃点数来给无线和有线网络的优先级排序,详见网上。这是一个比较通用的方法。...
2019-03-23 09:09:17 5092
原创 tensorflow入门实例
陆老师说,学习一门语言,或者说学习一个工具最好最快的方式就是阅读代码,偶尔查查手册,亦或是借助一个google或者bing,深以为然。下面是tensorflow入门的小例子,看看就差不多了了解了tensorflow的运作方式。这几个例子是:三维拟合、CNN和GAN。这些例子我是从别处抄的,自己看了看,加了些注释。tensorflow三维拟合这里例子是官方学习网站上抄的。# coding...
2019-03-22 17:20:32 2060
原创 水平集方法的一个基本框架
水平集方法框架水平集方法是现代图像处理中很重要的一个方法,为了说清楚这个东西,我们先介绍几个基本的概念。零水平集对于一个函数 ϕ(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⃗∣ϕ...
2019-03-11 21:25:06 5641 7
原创 关于Anaconda下TensorFlow安装的小注记(win10)
tensorflow刚出来的时候,只是在linux和mac系统下安装比较方便。后来,在windows的python 3.5下的集成安装包也出现了。再后来,任意版本的python也能极方便地安装tensorflow了。下面简单介绍一下,在windows下,基于anaconda的tensorflow快速安装。假设你已经安装好了anaconda,如果是默认根环境配置的是其他版本的python,也...
2019-02-17 19:26:10 308
原创 关于使用dct求解零诺依曼边界条件PDE的一点说明
关于使用dct求解零诺依曼边界条件PDE的一点说明DCT的第二型定义如我的这篇文章中所述,对于纽曼边界条件的偏微分方程,谱方法中,我们一般使用离散余弦变换,即dct来完成。对于fft和dst而言,似乎定义比较唯一。但是对于,对于dct,它的定义形式就有很多种,维基百科中就列出了它的八种形式。一般常用的形式有两种,即DCT-I:二维形式形如下面这样:求和符号的两个撇撇表示求和符号中当m...
2019-01-17 20:47:39 3613
原创 多尺度(有限元)降阶模型下的DNN方法
多尺度模型的DNN解法今天,想和大家分享的是Eric Chung的三篇文章。Deep Multiscale Model Learning、Deep Global Model Reduction Learning、Prediction of Discretization of GMsFEM using Deep Learning。Eric是港中文(The Chinese University o...
2019-01-15 17:54:49 5768 4
原创 神经网络和有限元方法
用神经网络来表示有限元函数今天,我想分享的东西是许进超老师的一篇文章,ReLU Deep Neural Networks and Linear Finite Elements。你可以理解为,这是一篇译文。主要内容它主要是关于神经元数与层数,在用DNN表示线性有限元的时候。它告诉我们,当维数大等于2的时候,我们至少需要两个隐藏层,在用ReLu DNN表示线性有限元函数的时候。因为在其它的文献中...
2019-01-03 09:01:58 6813 3
原创 关于傅里叶级数(一般、指数、广义)、傅里叶变换、离散时间傅里叶变换、傅里叶变换的一个补充说明
关于傅里叶级数(一般、指数、广义)、傅里叶变换、离散时间傅里叶变换、傅里叶变换的一个补充说明傅里叶级数傅里叶级数定义式为:这个定义是我们在大学里学的。你们还记得这个级数什么时候收敛吗?比如说,一个充分条件是f(x)连续的周期函数,且傅里叶系数满足:则傅里叶级数一致收敛(Uniform convergent)。一般来说啊,在现实生活中,我们一般要求:1、周期性。2、周期内绝对可积(∈...
2018-12-21 10:38:57 7883 2
原创 使用傅里叶变换(包括差分傅里叶和傅里叶谱方法)及切比雪夫谱方法(配点)法求解PDE
傅里叶变换及切比雪夫谱方法(配点)法求解PDE以下的内容介绍的是傅里叶谱方法求解PDE、切比雪夫谱方法求解PDE(这里指的是配点法)以及一种先进行差分离散,再对离散系统的每个变量使用离散傅里叶级数展开的求解PDE的方法。因为时间仓促,文中的很多公式,我没能自己手打出来,而用简单地截图来替代,因此也导致了本文在排版上比较混乱,图片大小不一。另外,也因为截图的原因,导致符号不一致(因为截自不同...
2018-12-06 11:20:20 13256 42
原创 Ubuntu下phg的安装
Ubuntu下phg的安装PHG (Parallel Hierarchical Grid) 是科学与工程计算国家重点实验室正在发展的专门为三维自适应有限元设计的并行程序开发平台,其核心是分布式的层次网格结构。目前,PHG 处理的网格对象是三维四面体协调网格。PHG 采用 C 语言开发,基于 MPI 消息传递通信实现并行。PHG 通过面向对象的数据结构以及用户接口实现了并行网格剖分、动态负载平衡...
2018-11-23 15:23:17 2066
原创 Ubuntu(64 bit)下Tetgen和Tetview的安装
Ubuntu(64 bit)下Tetgen和Tetview的安装Tetgen的安装相对容易,Tetview因为官网上只有32位版本的,所以它所需要的很多动态链接库依赖也需要是32位的,所以,当提醒库缺失时,我们要想办法把缺失的库安上,并且需要安上32位的。Tetgen的下载和安装我们可以在官网下载源码,下载后解压,cd进目录,然后使用make安装,编译完成之后,目录下就会生成一个名为 t...
2018-11-22 18:25:06 2100
原创 matlab的PDE工具箱的简单使用
matlab的PDE工具箱的简单使用问题选择在一个二维的有界区域Ω\OmegaΩ上,matlab的PDE工具箱主要使用有限元方法解决下面四类问题:椭圆形方程(elliptic)−∇⋅(c∇u)+au=f-\nabla\cdot(c\nabla u)+au=f−∇⋅(c∇u)+au=f抛物型方程(parabolic)d∂u∂t−∇(c⋅∇u)+au=fd\frac{\partial u...
2018-11-15 08:51:07 53481 103
原创 用深度学习求解高维偏微分方程
用深度学习求解高维偏微分方程简述:我们一般求解的传统的PDE维数也就二维三维,但是在其他领域中,比如说金融学,通过数学建模构建出来的PDE的维数及其之高,维数升高带来的“维数灾难”问题亟待解决。这篇文章,介绍了一种基于深度学习的方法,可以处理一般的高维抛物型方程。这篇文章,用反向随机微分方程构造PDE,并用神经网络近似未知解的梯度,来求解高维的偏微分方程。高维PDE的来源:1、量子多体问...
2018-11-08 20:18:50 21427 50
原创 桥牌规则
桥牌规则大学选修了一门桥牌课程,受益匪浅。下面,不考虑套路,简单说说桥牌的规则,按流程进行的顺序来讲。发牌和理牌四人分坐东南西北四个位置,这个位置可以认为标定,南北两人为一方,东西两人为一方。一副牌去掉大小王为52张,依次按张发给四个人,没人12张。理牌,习惯上从左到右按花色和大小把牌分开,花色大小顺序依次为黑桃(S)>红桃(H)>方片(D)>梅花(C),点数大小为A&g...
2018-11-03 21:07:23 15092
原创 matlab做三维线性拟合(多元线性回归,准确来说不叫插值)
matlab三维拟合(多元线性回归)问题描述今天同学问了我一个问题,大概意思是给了你三列输入数据,一列输出数据,想用一个线性超平面做一个最小二乘拟合(注意这里不能叫插值)。一点思考刚听到这个问题,同学说的是做插值,说想要做一个插值,这种说法不准确的,不想说回归的话,你可以说这是一个拟合。插值和拟合的区别在哪呢?插值要求你做出来的函数(比如说多项式插值)要经过所有给定的数据点,而拟合只要求尽...
2018-10-24 21:40:25 27404 3
原创 LaTeX幻灯片通用简洁模板(Beamer)
LaTeX幻灯片通用简洁模板(Beamer)提供一个Ctex套装下使用WinEdt编辑器的xelatex编译可以通过的一个LaTex幻灯Beamer模板,使用方法为将代码拷贝进WinEdt文档,使用xelatex编译(一般Beamer都用这个编译,不用多想和尝试),点击菜单栏中的pdf查看按钮,可以看见效果。然后,如果想使用某些功能,比如说插图,可以把代码滑到底部,找到注释中相关内容,复制黏贴...
2018-10-11 11:01:54 58941 18
原创 矩阵相乘的OpenMP实现
矩阵相乘的OpenMP实现OpenMP简介OpenMP是基于共享内存的编程,不同于MPI,因为是共享内存,所以它不需要将计算结果丢来丢去共享。事实上,我们可以用很少的编译指导语句实现并行,而不需要关心底层的操作。举个简单的例子,如果程序有for循环,我们只要在for循环前加了一句#pragma omp parallel for,在循环语句前后不相关的情况下,对for循环就实现了并行。如下两种...
2018-10-03 00:22:48 7892 1
原创 使用MPI消息传递实现hello world的顺序输出
MPI:hello world顺序输出在跑MPI并行程序时,我们不知道哪个核心前进得快,故而在没有其他控制的条件下,谁先执行同一条代码,是不可知的。比如说,对于一个输出“hello world”的程序,如果采用多个节点来跑,打印顺序是不定的。下面是一个hellow程序:/* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- ...
2018-10-02 23:35:28 4649 9
各种课程大作业、项目等 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关注的人