- 博客(10)
- 收藏
- 关注
原创 fealpy网格编号问题,
这就是编号为 0 的边, 所在的第一个剖分和第二个剖分的编号, 这条边在第一个剖分边的编号, 这条边在第二个剖分边的编号。用 fealpy 对求解域进行剖分, 我写下魏华祎老师已经实现的基本的网格信息点的函数,以四边形网格为例。这就是编号为0的剖分, 紧邻着的剖分编号, 对于第0个剖分, 挨着它的就是 2 1 0 以逆时针排列。↑ 这就是第0个剖分的边所在的全局编号(边的信息是按逆时针排布给的)我们来详细说下后三个, 并以详细介绍每个例子的第一行含义。这样建立了如下的四边形网格。
2023-12-13 01:10:08
90
原创 python 画多组图,部分子图不显示刻度
用```matplotlib.pyplot```画了多组图, 由于每组子图太多, 有的子图不需要展示Y轴和刻度,论坛中对这种情况的教程比较少,这句是不行的,这只会让每组图的最后一个子图不显示坐标轴刻度值.给一个算例,后续可以用if把想关闭的轴关掉。以下是我给出的解决方法。
2023-11-08 10:37:26
598
1
原创 云服务器的linux或没有图形化界面的linux 用python画多张图,并保存
服务器的linux系统是没有图形化界面的,所以用Python的matplotlib库画图可能不会出图,有时我们需要画多张图,所以在linux保存多张图片。在matplotlib调用前多写两句。首先,我们先创建三张图。
2023-11-04 20:51:55
139
1
原创 FFT的note
这个式子告诉我们,若'''n'''有'''N-1'''个取值, 则'''k'''也有'''N-1'''个取值,是一一对应的,他们数值上的对应关系'''k=pi*N/L**2''' , L是求解区间'''[-L, L]'''我们看离散傅里叶变换的微分性质:对于离散傅里叶算子'''F( )''',作用在'''u(x)'''上,我们有。fftfreq是取k的命令,我们不需要手动赋予,2*pi*k才是我们一开始取的k。注:fftfreq第二个参数 取得的是每个x的差值。我们对'''x(n)'''进行傅里叶变换,
2023-08-22 21:20:01
86
2
原创 连续两天遇到的报错类型,value cannot be converted to type double without overflow
rt
2023-08-20 21:28:07
440
1
原创 python中,定义一个pde函数,传入的数据维度可以是不一样的
在写FT_PINNs时,遇到问题,因为对x进行FFT变换,所以pde中的x成为了参数,我做了一个测试,
2023-08-18 16:09:03
93
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人