- 博客(23)
- 收藏
- 关注
RSS订阅原创 Neumann边界条件下二维Legendre-Galerkin方法快速计算
最近考虑Neumann 边界问题的 Legendre-Galerkin 方法求解, 由于书中只考了 Dirichlet 边界条件下的 特征值方法处理，通过查阅资料，整理了Neumann 边界条件下 Legendre-Galerkin方法的快速处理方法。对于二维的方程, Legendre-Galerkin 方法的基函数空间为两个一维的基函数空间直接张量积得到。离散后的线性系统如下：
2022-12-06 23:13:19
886
1
原创 Python自定义 colormap
from matplotlib.colors import LinearSegmentedColormapcm_data = [[0.2081, 0.1663, 0.5292], [0.2116238095, 0.1897809524, 0.5776761905], [0.212252381, 0.2137714286, 0.6269714286], [0.2081, 0.2386, 0.6770857143], [0.1959047619, 0.2644571429, 0.7279], [0
2021-10-07 13:54:22
836
原创 解决 Ubuntu 下 W: 无法下载 http://archive.getdeb.net/ubuntu/dists/xenial-getdeb/InRelease 500 的报错
之前安装了 Ubuntu 的一个文件, 由于东西没有删干净导致运行sudo apt-get update 出现:W: 无法下载 http://archive.getdeb.net/ubuntu/dists/xenial-getdeb/InRelease 500 Internal Server Error [IP: 2606:4700:3030::6815:4f84 80]W: 部分索引文件下载失败。如果忽略它们，那将转而使用旧的索引文件。解决方案:进入目录:cd /etc/apt/sou
2021-09-18 20:20:37
1127
原创 Matplotlib 一些喜欢的绘图参数整理
# === # === figsize, linewith # === # 方图: plt.figure(figsize = (12, 10)) # 线宽lw = 3 # 喜好设置 lw = 3 左右ms = 10 # markersize 为 10 左右# ===# === useful colors# ===color = ['r', 'green', 'orange', 'dodgerblue', 'darkmeganta'] # 一些好的配色# === # === pl
2021-09-18 11:09:45
351
原创 Ubuntu 下安装 vim-tex 以及 markdown-preview插件
最近换 Ubuntu 系统, 打算使用 neovim, 这两天成功配置了vim-tex插件以及 markdown-preview 插件,道路比较坎坷，写篇博客将这两天安装的过程进行整理.构建 neovim 的配置文件切换到 ~/.config/nvim下查看有没有 init.vim配置文件，如果没有,请手动创建. 这个文件就是 neovim 的配置文件, neovim 的所有配置文件全部在这个文件中运行.安装 git后面装 vim-plug 插件管理器以及其他的插件都是从 github 上面下载的,
2021-09-15 14:43:33
950
原创 2021-09-12 Linux 下 Anaconda下载以及安装
Linux 下 Anaconda 下载以及安装软件包下载以及安装首先去清华镜像网站下载 Linux 平台下的 Anaconda 软件包， 链接如下：https://mirrors.tuna.tsinghua.edu.cn/anaconda/archive/或者直接在终端下使用命令wget https://mirrors.tuna.tsinghua.edu.cn/anaconda/archive/Anaconda3-5.3.1-Linux-x86_64.sh这里注意看好版本以及对应的平台。下载
2021-09-12 02:07:06
1804
原创 2021.6.9 Python部分常用函数整理
import numpy as npimport matplotlib.pyolot as plt np.ix_ ; plt.spy 函数np.ix_([0, 1], [1,2]) # 构造出一个索引矩阵, 用于索引其它矩阵的 (0,1), (0, 2), (1, 0), (1, 2) 元plt.spy(M, color='r') # 绘制出矩阵 M 的稀疏矩阵结构, 其中 color 参数用来表示绘图的颜色.Legendre 多项式相关的函数from scipy.special.ortho
2021-06-09 17:51:19
190
原创 2021.6.2 Python 中的三角形网格
Python 中的三角形网格函数最近写有限元程序, 有些地方需要用到三角形网格, 这里记录一下, 免的以后忘记, 需要用到的包:import numpy as npimport matplotlib.pyplot as pltfrom matplotlib.patches import Polygonimport matplotlib.tri as tri生成三角形数据 tri.Triangulation()方法可以将平面上一组点集用来生成三角形数据对象(默认使用Delaunay三角化),
2021-06-03 08:45:56
2710
原创 2020-5-29有限元从0开始第九天(有限元的等价性)
继续给出参考资料镇楼Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 2020-5-27有限元从0开始第八天(插值算子)上次内容简要回忆在前面博客中给出了有限元局部插值算子, 以及全局插值算子的定义, 例子以及一些相关的性质, 在本节主要接招有限元的等价性.3.4 有限元的等价性在全局插值的应用中, 找到插值I
2021-05-29 16:14:51
569
4
原创 2020-5-27有限元从0开始第八天(插值算子)
继续给出参考资料镇楼Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 2021-5-26 第七天 (有限元的基本概念)昨日内容简要回忆在前面昨天的文章中, 讲述了有限元的基本概念, 给出了有限元的基本定义, 尤其是给出引理, 给出验证有限元定义第三条的方法, 即怎么去验证所给出的 N\mathcal{N}N, 定义
2021-05-27 11:32:08
798
原创 2021-5-26 第七天 (有限元的基本概念)
附上参考资料镇楼Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 从第三章开始看了,中间的东西抽空每天打点, Sobolev 空间的东西不做题还是没用.首先给出有限元的定义: 定义: (有限元) 令:K⊂RnK \subset \mathbb{R}^{n}K⊂Rn 表示内点非空并且具有分片光滑边界的有界闭集(网
2021-05-26 23:54:19
737
原创 2021-5-25有限元从0开始第六天(Lp空间,弱导数)
前面几天都在写下面的这个编程文档, 这个文档写起来东西真的多, 目前程序还没有调试对, 最后测试误差阶的时候出现了点问题, 不过没事, 东西还是得继续往下看的.2021-5-20 有限元方法从0开始学习笔记 继续接第三天 (一维线性有限元, 二次元编程计算推导)继续附上参考资料镇楼Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media,
2021-05-26 00:11:15
601
原创 2021-5-20 有限元方法从0开始学习笔记 继续接第三天 (一维线性有限元, 二次元编程计算推导)
今天的笔记将解第三天的内容, 对一维线性元进行程序实现, 第三天内容的链接接在下面:2021-5-18 有限元方法从0开始学习笔记 第三天 (有限元方法计算处理简介)继续附上参考资料镇楼Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 计算公式推导...
2021-05-24 16:25:08
1350
原创 2021-5-22 有限元方法从0开始学习笔记 第五天 (适应性逼近简介)
周末了, 早上好, 继续往下看, 给出参考资料如下:Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 第十三页上一次内容回忆2021-5-21 有限元方法从0开始学习笔记 第四天 (局部误差)给出了逐点误差估计,但是看的不是很懂, 还有一个地方没有推过去, 打算先放着.0.8 适应性逼近本章对适应性逼近方法的
2021-05-22 13:50:52
324
原创 2021-5-21 有限元方法从0开始学习笔记 第四天 (局部误差)
因为昨天要备课, 今天要讲一天的讨论课, 所以昨天没有写, 今天早上也没时间写, 从晚上开始写吧, 虽然少看了点, 但还是慢慢来把, 慢慢来总能全看会的.老样子, 给出参考资料:Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 第十二页上一次内容回忆2021-5-19 有限元方法从0开始学习笔记 第三天 (有限元方
2021-05-22 08:48:11
306
原创 2021-5-19 有限元方法从0开始学习笔记 第三天 (有限元方法计算处理简介)
又是新的一天, 今天上午有课, 因此可能只能写一点点了, 继续给出参考资料Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 从第十页.昨日内容简单回忆2021-5-18 有限元方法从0开始学习笔记 第二天 (一维线性元简介, 与差分格式的联系)给出了第一个具体的有限元空间, 一维的线性分片多项式空间, 在该空间
2021-05-20 09:27:19
466
原创 2021-5-18 有限元方法从0开始学习笔记 第二天 (一维线性元简介, 与差分格式的联系)
又是新的一天, 继续开始有限元从0开始学习, 参考资料Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. 昨天内容回顾2021-5-17 有限元方法从0开始学习笔记 第一天利用模型问题, 构造了其变分形式, 并且验证了当 fff 和 uuu 的光滑性足够时, 变分问题能够得到边值问题对古典解.阐述了 Ritz-
2021-05-19 00:35:17
591
2
原创 2021-5-17 有限元方法从0开始学习笔记 第一天
第一天从今天开始从0学习有限元方法, 由于记纸质笔记太麻烦, 字也比较丑, 无法静下心来, 因此选择在博客上记录笔记, 每天能看多少记录多少, 虽然速度可能会慢一点, 但希望能够静下心来, 坚持学习.参考资料 Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007. (0.1 节)弱形式考虑边值问题:{−d2udx2
2021-05-17 23:55:45
646
1
原创 离散Poisson系统的常用处理方法以及比较(Python实现)
本文主要介绍 Possion 方程标准五点差分方法离散后的系统的常用处理方法,并对他们进行比较.
2021-03-24 23:02:16
821
原创 一维Helmholtz方程的Chebyshev - Galerkin谱方法以及Python实现
加权余量法简介给定算子方程以及边界条件:{Tu−f=0u∈VBu−g=0u∈S\begin{cases} Tu - f = 0 & u \in V \\ Bu - g = 0 & u \in S\end{cases}{Tu−f=0Bu−g=0u∈Vu∈S其中uuu为待求函数, TTT 为区域内部 VVV 上算子, BBB 为区域边界上的算子, f,gf, gf,g 分别为定义在区域内部和区域边界的函数.处理上述问题的加权余量方法如下:首先,设方程组的近似解为:un
2020-12-22 16:58:17
1595
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人