我试图写一个函数来取任何一般函数/数列的导数。具体来说,我使用的是。问题是,我不能计算导数的边界点,因为中心差分公式使用超出界限的指数。我的代码在下面
import numpy as np
n = 20000 # number of points in array
xs = np.linspace(start=-2*np.pi, stop=2*np.pi, num=n) # x values
y = np.array([np.sin(i) for i in xs]) # our function, sine
def deriv(f, h):
"""
C
我的目标是在这两个电路曲线形状之间找到一条平滑的最佳拟合线。有没有比我的算法更好的算法,可以找到两条直线之间的一组点(或一条曲线),比如这个例子?
到目前为止,我的算法取内部部分,并为每个点找到最接近的点,但是这不起作用,因为(看第一个角)。
(红色是内部,绿色是外部,蓝色是我找到的优化后的点)
这是我的jsfiddle:
这就是算法:
for (i = 0; i < coords[0].length; i++) {
var currentI = coords[0][i];
j = 0;
var currentJ = coords[0][j];
current
下面的程序解决了一维热方程,它代表的杆,其两端保持在零温度与初始条件10*np.sin(np.pi*x)。
在这个计算中,Dirichlet边界条件(两端温度为零)是如何工作的?我被告知矩阵A的上、下行包含两个非零元素,而缺失的第三个元素是Dirichlet条件。但我不明白这个条件是怎样影响计算的。如果A中缺少元素,u_{0}或u_{n}如何为零?
下面的有限差分法使用的是。
import numpy as np
import scipy.linalg
# Number of internal points
N = 200
# Calculate Spatial Step-Size
h =