第七章微积分数值计算.pptVIP

X=[1 3 5 2 4]; Y= gradient(X) Y = 2.0000 2.0000 -0.5000 -0.5000 2.0000 Y=gradient(X,2) Y= 1.0000 1.0000 -0.2500 -0.2500 1.0000 即 两边用前向和后向差分，中间用中心差分 7.1.4 拉普拉斯算子4*del2 由于内部算法的原因： U=4*del2(v,h)，对1维向量v以步长h求拉普拉斯算子时，返回一相同维数的向量U，且 默认的步长为1。 U=4*del2(v,h1,h2)，对矩阵v，横向(x方向)以步长h1，纵向(y方向)以步长 h2计算拉普拉斯算子。 4*del2(U) ans = 4 4 4 4 4 4 4 4 4 4 4 4 [x,y]= meshgrid(1:4,1:3); U = x.*x+y.*y U = 2 5 10 17 5 8 13 20 10 13 18 25 首先编写vdp.m如下： function fy=vdp(x,y) %注意表示一个微分方程组的函数必