前言

本文主要介绍了数值微分、数值积分的方法，以及MATLAB利用数值微积分方法进行微积分运算的方法。

如果读者只是为了快速了解MATLAB数值运算方法则可直接阅读微、积分运算章节，但为了更加熟练的使用微积分运算，还是建立读者先本文顺序先学习数值微积分的方法。在学习过程可以发现，MATLAB微积分运算章节与数值微积分章节都是一一对应的。

数值微分

导数的差分公式

我们在高等数学中学习过，函数在附近某点取值的近似公式——泰勒公式为：f(x+h)=f(x)+hf′(x)1!+h2f′′(x)2!+⋯+hnf(n)(x)n!+o(hn)f(x+h)=f\left(x\right)+h\frac{f^{\prime}\left(x\right)}{1 !}+h^2\frac{f''(x)}{2!}+\cdots+h^{n}\frac{f^{(n)}\left(x\right)}{n !}+o\left(h^{n}\right)f(x+h)=f(x)+h1!f′(x)+h22!f′′(x)+⋯+hnn!f(n)(x)+o(hn)

由泰勒展开式 f(x+h)=f(x)+hf′(x)+o(h2)f(x+h)=f(x)+hf'(x)+o(h^{2})f(x+h)=f(x)+hf′(x)+o(h2)，f(x−h)=f(x)−hf′(x)+o(h2)f(x-h)=f(x)-hf'(x)+o(h^{2})f(x−h)=f(x)−hf′(x)+o(h2)，得一阶导数的差分公式：

向前差分	向后差分	中心差分
f′(x)=f(x+h)−f(x)h+o(h)f'(x)=\frac{f(x+h)-f(x)}{h}+o(h)f′(x)=hf(x+h)−f(x)+o(h)	f′(x)=f(x)−f(x−h)h+o(h)f'(x)=\frac{f(x)-f(x-h)}{h}+o(h)f′(x)=hf(x)−f(x−h)+o(h)	f′(x)=f(x+h)−f(x−h)2h+o(h2)f'(x)=\frac{f(x+h)-f(x-h)}{2h}+o(h^{2})f′(x)=2hf(x+h)−f(x−h)+o(h2)

无穷小量的阶数代表了算法的精度，当公式形式变为f(k)x=[]+o(hn)f^{(k)}x=[\ ]+o(h^n)f(k)x=[ ]+o(hn)时，o(hn)o(h^n)o(hn)表示公式为n阶精度的差分公式。

一阶的4阶精度的差分公式：f′(x)=8[f(x+h)−f(x−h)]−[f(x+2h)−f(x−2h)]12h+o(h4)f'(x)=\frac{8\left[ f(x+h)-f(x-h)\right]-\left[f(x+2h)-f(x-2h)\right]}{12h}+o(h^4)f′(x)=12h8[f(x+h)−f(x−h)]−[f(x+2h)−f(x−2h)]+o(h4)；

这里前一个阶数指所要表示的 f(n)f^{(n)}f(n) 的阶数，第二个阶数指公式对应的精度。

二阶的2阶精度的差分公式：f′′(x)=f(x+h)−2f(x)+f(x−h)h2+o(h2)f''(x)=\frac{f(x+h)-2f(x)+f(x-h)}{h^2}+o(h^2)f′′(x)=h2f(x+h)−2f(x)+f(x−h)+o(h2)；

三阶的2阶精度的差分公式：f′′′(x)=[f(x+2h)−f(x−2h)]−2[f(x+h)−f(x−h)]2h3+o(h2)f'''(x)=\frac{\left[f(x+2h)-f(x-2h)\right]-2\left[f(x+h)-f(x-h)\right]}{2h^3}+o(h^2)f′′′(x)=2h3[f(x+2h)−f(x−2h)]−2[f(x+h)−f(x−h)]+o(h2)。

偏微分的差分公式

nabla算子的运算

∇f=(∂xf,∂yf,∂f)\nabla f=(\partial_x f,\partial_y f, \partial f)∇f=(∂xf,∂yf,∂f)， ∇⋅F⃗=∂xFx+∂yFy+∂zFz\nabla\cdot \vec{F}=\partial_xF_x+\partial_yF_y+\partial_zF_z∇⋅F=∂xFx+∂yFy+∂zFz， ∇2f=(∂x2+∂y2+∂z2)f\nabla^2 f=(\partial_x^2+\partial_y^2+\partial_z^2)f∇2f=(∂x2+∂y2+∂z2)f。

∇×F⃗=(∂yFz−∂zFy,∂zFx−∂xFz,∂xFy−∂yFx)\nabla\times \vec{F}=\left(\partial_y F_z-\partial_zF_y,\ \partial_zF_x-\partial_xF_z,\ \partial_xF_y-\partial_yF_x\right)∇×F=(∂yFz−∂zFy, ∂zFx−∂xFz, ∂xFy−∂yFx)

nabla算子对应的数值差分公式

∇f=(∂xf,∂yf,∂zf)=(fi+1,j,k−fi−1,j,k2h,fi,j+1,k−fi,j−1,k2h,fi,j,k+1−fi,j,k−12h)∇2f=(∂x2f+∂y2f+∂z2f)=(fi+1,j,k+fi−1,j,k−2fi,j,kh2+fi,j+1,k+fi,j−1,k−2fi,j,kh2+fi,j,k+1+fi,j,k−1−2fi,j,kh2)F=(Fx(r),Fy(r),Fz(r));r=(x,y,z)∇⋅F=∂xFx(r)+∂yFy(r)+∂zFz(r)=Fxi+1,j,k−Fxi−1,j,k2h+Fyi,j+1,k−Fyi,j+1,k2h+Fzi,j,k+1−Fzi,j,k−12h∂yFz(r)−∂zFy(r)=Fzi,j+1,k−Fzi,j−1,k2h−Fyi,j,k+1−Fyi,j,k−12h\begin{aligned} &\nabla f=\left(\partial_{x} f, \partial_{y} f, \partial_{z} f\right)=\left(\frac{f_{i+1, j, k}-f_{i-1, j, k}}{2 h}, \frac{f_{i, j+1, k}-f_{i, j-1, k}}{2 h}, \frac{f_{i, j, k+1}-f_{i, j, k-1}}{2 h}\right) \\ &\nabla^{2} f=\left(\partial_{x}^{2} f+\partial_{y}^{2} f+\partial_{z}^{2} f\right) \\ &=\left(\frac{f_{i+1, j, k}+f_{i-1, j, k}-2 f_{i, j, k}}{h^{2}}+\frac{f_{i, j+1, k}+f_{i, j-1, k}-2 f_{i, j, k}}{h^{2}}+\frac{f_{i, j, k+1}+f_{i, j, k-1}-2 f_{i, j, k}}{h^{2}}\right) \\ &\mathbf{F}=(F x(\mathbf{r}), F y(\mathbf{r}), F z(\mathbf{r})) ; \quad \mathbf{r}=(x, y, z) \\ &\nabla \cdot \mathbf{F}=\partial_{x} F_{x}(\mathbf{r})+\partial_{y} F_{y}(\mathbf{r})+\partial_{z} F_{z}(\mathbf{r}) \\ &=\frac{F x_{i+1, j, k}-F x_{i-1, j, k}}{2 h}+\frac{F y_{i, j+1, k}-F y_{i, j+1, k}}{2 h}+\frac{F z_{i, j, k+1}-F z_{i, j, k-1}}{2 h} \\ &\partial_{y} F z(\mathbf{r})-\partial_{z} F y(\mathbf{r})=\frac{F z_{i, j+1, k}-F z_{i, j-1, k}}{2 h}-\frac{F y_{i, j, k+1}-F y_{i, j, k-1}}{2 h} \end{aligned}∇f=(∂xf,∂yf,∂zf)=(2hfi+1,j,k−fi−1,j,k,2hfi,j+1,k−fi,j−1,k,2hfi,j,k+1−fi,j,k−1)∇2f=(∂x2f+∂y2f+∂z2f)=(h2fi+1,j,k+fi−1,j,k−2fi,j,k+h2fi,j+1,k+fi,j−1,k−2fi,j,k+h2fi,j,k+1+fi,j,k−1−2fi,j,k)F=(Fx(r),Fy(r),Fz(r));r=(x,y,z)∇⋅F=∂xFx(r)+∂yFy(r)+∂zFz(r)=2hFxi+1,j,k−Fxi−1,j,k+2hFyi,j+1,k−Fyi,j+1,k+2hFzi,j,k+1−Fzi,j,k−1∂yFz(r)−∂zFy(r)=2hFzi,j+1,k−Fzi,j−1,k−2hFyi,j,k+1−Fyi,j,k−1

数值积分

我们先推导数值积分的近似计算公式。

一维定积分 I=∫abf(x)dxI=\int_{a}^{b} f(x) \mathrm{d} xI=∫abf(x)dx 的几何意义是曲线下面所包围的面积，按照积分的定义，它是许多小曲边梯形面积之和的极限 I=∫abf(x)dx=lim⁡n→∞∑k=1nf(ξk)b−anI=\int_{a}^{b} f(x) \mathrm{d} x=\lim _{n \rightarrow \infty} \sum_{k=1}^{n} f\left(\xi_{k}\right) \frac{b-a}{n}I=∫abf(x)dx=limn→∞∑k=1nf(ξk)nb−a。

令 h=b−anh=\frac{b-a}{n}h=nb−a，将微分的数值计算过程倒转可得到“向前积分”和“向后积分”公式，它们都叫矩形公式，都是用矩形来替代曲边梯形。Ln=h∑k=0n−1fk,Rn=h∑k=1nfkL_{n}=h \sum_{k=0}^{n-1} f_{k},\ R_{n}=h \sum_{k=1}^{n} f_{k}Ln=h∑k=0n−1fk, Rn=h∑k=1nfk。

将两式平均得出近似计算定积分的梯形公式：Tn=h∑k=1n−1fk+h2(f0+fn)T_{n}=h \sum_{k=1}^{n-1} f_{k}+\frac{h}{2}\left(f_{0}+f_{n}\right)Tn=h∑k=1n−1fk+2h(f0+fn)；

将中心差分的过程倒转得到了中点法(积分)公式：Mn=h∑k=0n−1f(xk+12)M_n=h\sum_{k=0}^{n-1} f(x_{k+\frac{1}{2}})Mn=h∑k=0n−1f(xk+21)。

在高等数学中介绍过，辛普森(Simpson)公式：Sn=h3(f0+f2n+2∑k=1n−1f2k+4∑k=0n−1f2k+1)S_n=\frac{h}{3}\left(f_0+f_{2n}+2\sum_{k=1}^{n-1}f_{2k}+4\sum_{k=0}^{n-1}f_{2k+1} \right)Sn=3h(f0+f2n+2∑k=1n−1f2k+4∑k=0n−1f2k+1)。

在a,b两点间，辛普森公式为：S=h6(f(a)+4f(a+b2)+f(b))S=\frac{h}{6}(f(a)+4f(\frac{a+b}{2})+f(b))S=6h(f(a)+4f(2a+b)+f(b))

在a,b区间内，辛普森公式使用了f(a),f(b)f(a),f(b)f(a),f(b)以及中点f(a+b2)f(\frac{a+b}{2})f(2a+b)，更进一步，还可以再多使用两个点f(3a+b4)f(\frac{3a+b}{4})f(43a+b)和f(a+3b4)f(\frac{a+3b}{4})f(4a+3b)得到复合辛普森公式：S=h12(f(a)+4f(3a+b4)+a+b2+4f(a+3b4)+f(b))S=\frac{h}{12}(f(a)+4f(\frac{3a+b}{4})+\frac{a+b}{2}+4f(\frac{a+3b}{4})+f(b))S=12h(f(a)+4f(43a+b)+2a+b+4f(4a+3b)+f(b))。

对数值积分的方法做个简单整理

矩形公式：Ln=h∑k=0n−1fk,Rn=h∑k=1nfkL_{n}=h \sum_{k=0}^{n-1} f_{k},\ R_{n}=h \sum_{k=1}^{n} f_{k}Ln=h∑k=0n−1fk, Rn=h∑k=1nfk。

梯形公式：Tn=h∑k=1n−1fk+h2(f0+fn)T_{n}=h \sum_{k=1}^{n-1} f_{k}+\frac{h}{2}\left(f_{0}+f_{n}\right)Tn=h∑k=1n−1fk+2h(f0+fn)；

中点法：Mn=h∑k=0n−1f(xk+12)M_n=h\sum_{k=0}^{n-1} f(x_{k+\frac{1}{2}})Mn=h∑k=0n−1f(xk+21)。

在a,b两点间，辛普森公式为：S=h6(f(a)+4f(a+b2)+f(b))S=\frac{h}{6}(f(a)+4f(\frac{a+b}{2})+f(b))S=6h(f(a)+4f(2a+b)+f(b))

在a,b区间内，复合辛普森公式：S=h12(f(a)+4f(3a+b4)+a+b2+4f(a+3b4)+f(b))S=\frac{h}{12}(f(a)+4f(\frac{3a+b}{4})+\frac{a+b}{2}+4f(\frac{a+3b}{4})+f(b))S=12h(f(a)+4f(43a+b)+2a+b+4f(4a+3b)+f(b))。

六阶牛顿-柯林特法则

牛顿-柯林特法则是在辛普森公式等四阶积分法则上进行的修正。SSS 和 S2S_{2}S2 是同一个积分的近似值 (S2S_2S2的步长是SSS的一半)，因此它们之间的差可以用于估计误差 E=S2−SE=S_{2}-SE=S2−S。而且他们还可以结合起来得到一个更准确的近似值QQQ，而由步长和原来积分法则是4阶可知，S2S_2S2的准确值大约是SSS的242^424倍，所以 Q−S=16(Q−S2)Q-S=16(Q-S_2)Q−S=16(Q−S2)，因此 Q=S2+(S2−S)/15Q=S_2+(S_2-S)/15Q=S2+(S2−S)/15。

数值积分的阶数

数值积分的阶数指使用该法则不能精确计算的多项式的最低次数。

以中点法为例，中点法只能对 f(x)=Constf(x)=Constf(x)=Const 以及 f(x)=kx+Constf(x)=kx+Constf(x)=kx+Const 等函数进行精确计算，对于f(x)=kx2f(x)=kx^2f(x)=kx2则无法精确计算，因此中点法阶数为2阶。

类似可得，辛普森和复合辛普森公式阶数为4。

微分运算

微分运算与数值微分的方法相对应，先做个简单的复习。

数值微分的常有方法有向前差分 f′(x)=f(x+h)−f(x)h+o(h)f'(x)=\frac{f(x+h)-f(x)}{h}+o(h)f′(x)=hf(x+h)−f(x)+o(h)、向后差分 f′(x)=f(x)−f(x−h)h+o(h)f'(x)=\frac{f(x)-f(x-h)}{h}+o(h)f′(x)=hf(x)−f(x−h)+o(h)、中心差分 f′(x)=f(x+h)−f(x−h)2h+o(h2)f'(x)=\frac{f(x+h)-f(x-h)}{2h}+o(h^{2})f′(x)=2hf(x+h)−f(x−h)+o(h2) 与二阶的差分公式：f′′(x)=f(x+h)−2f(x)+f(x−h)h2+o(h2)f''(x)=\frac{f(x+h)-2f(x)+f(x-h)}{h^2}+o(h^2)f′′(x)=h2f(x+h)−2f(x)+f(x−h)+o(h2)。

差分运算 diff

%% diff
% 向量差分
a = [1 4 5 7 8 9];
diff(a)     % 类似于向前差分, ans=[3 1 2 1 1]
diff(a,2)   % diff(a,2)=diff(diff(a)), ans=[-2 1 -1 0]% 矩阵差分
b = [1 3 4; 4 6 7; 2 1 7];
diff(b)     % 矩阵差分是竖直方向操作的, ans=[3 3 3; -2 -5 0]
diff(b,3)   % diff(b,3)=diff(diff(diff(b))), ans=[-3 5]

梯度运算 gradient

%% gradient
% 向量梯度
a = [1 4 5 7 8 9];
gradient(a)     % 对向量各点取中心差分, 端点处则取向前/向后差分
% ans = [-3.0000 2.0000 1.5000 1.5000 1.0000 1.0000]% 矩阵梯度
b = [1 3 4; 4 6 7; 2 1 7];
gradient(b)     % 矩阵梯度又变为了横向的中心差分操作, 同样在端点处变为横向向前/向后差分
% ans = [2.0000 1.5000 1.0000; 2.0000 1.5000 1.0000 -1.0000 2.5000 6.0000]
[px, py]=gradient(b)   % 若使用两个矩阵承接矩阵梯度的结果, 则除px为gradient(b)的结果外, 还可以得到py为纵向的差分运算的结果。

离散拉普拉斯算符 del2

%% del2
x = [1 2 3 4 5 6 7 8];
fx=x.^2;
4*del2(fx)  % 4*del2相当于f''(需要乘4为del2函数的性质), ans=[2 2 2 2 2 2 2 2]
gx=x.^3;
4*del2(gx)  % ans=[6 12 18 24 30 36 42 48]

不要忘记del2也是数值计算，它在计算二阶导数时必然后对应公式，由于fx,gx阶数较低才可准确计算出导数。del2对应的公式为二阶的精度为2阶的差分公式 f′′(x)=f(x+h)−2f(x)+f(x−h)h2+o(h2)f''(x)=\frac{f(x+h)-2f(x)+f(x-h)}{h^2}+o(h^2)f′′(x)=h2f(x+h)−2f(x)+f(x−h)+o(h2)，h默认为1。

一旦del2要处理的函数阶数超过二阶差分公式的精度，那么del2就不可用f’'(x)进行计算。

a = [1 4 5 7 8 9];
4*del2(a)   % ans=[-5 -2 1 -1 0 1]
% 对于左端点 ans(1)=2*a(1)-5*a(2)+4*a(3)-a(4), 右端点同理% 矩阵拉普拉斯算子
b = [1 3 4; 4 6 7; 2 1 7];
4*del2(b)   % 对列表del2公式扩展 I(i,j)=[f(i+1,j)+f(i,j+1)+f(i-1,j)+f(i,j-1)-4f(i,j)]/4, ans=[-6 -9 -4; -6 -9 -4; 2, -1 4]
% 对于左下角的点, I(i,j)=(u(i+2,j)+u(i,j+2)+2*u(i,j))/4-(u(i+1,j)+u(i,j+1))/2
% 对于左边界点, I(i,j)=(u(i+1,j)+u(i-1,j)+u(i,j+1)+u(i,j))/4-(u(i,j)+u(i,j+1))/2

对于左下角的点 I(i,j)=u(i+2,j)+u(i,j+2)+2∗u(i,j)4−u(i+1,j)+u(i,j+1)2I(i,j)=\frac{u(i+2,j)+u(i,j+2)+2*u(i,j)}{4}-\frac{u(i+1,j)+u(i,j+1)}{2}I(i,j)=4u(i+2,j)+u(i,j+2)+2∗u(i,j)−2u(i+1,j)+u(i,j+1)
对于左边界点 I(i,j)=u(i+1,j)+u(i−1,j)+u(i,j+2)+u(i,j)4−u(i,j)+u(i,j+1)2I(i,j)=\frac{u(i+1,j)+u(i-1,j)+u(i,j+2)+u(i,j)}{4}-\frac{u(i,j)+u(i,j+1)}{2}I(i,j)=4u(i+1,j)+u(i−1,j)+u(i,j+2)+u(i,j)−2u(i,j)+u(i,j+1)

积分运算

数值积分内容回顾

矩形公式：Ln=h∑k=0n−1fk,Rn=h∑k=1nfkL_{n}=h \sum_{k=0}^{n-1} f_{k},\ R_{n}=h \sum_{k=1}^{n} f_{k}Ln=h∑k=0n−1fk, Rn=h∑k=1nfk；
梯形公式：Tn=h∑k=1n−1fk+h2(f0+fn)T_{n}=h \sum_{k=1}^{n-1} f_{k}+\frac{h}{2}\left(f_{0}+f_{n}\right)Tn=h∑k=1n−1fk+2h(f0+fn)；
中点法：Mn=h∑k=0n−1f(xk+12)M_n=h\sum_{k=0}^{n-1} f(x_{k+\frac{1}{2}})Mn=h∑k=0n−1f(xk+21)；
在a,b两点间，辛普森公式：S=h6(f(a)+4f(a+b2)+f(b))S=\frac{h}{6}(f(a)+4f(\frac{a+b}{2})+f(b))S=6h(f(a)+4f(2a+b)+f(b))；
在a,b区间内，复合辛普森公式：S=h12(f(a)+4f(3a+b4)+a+b2+4f(a+3b4)+f(b))S=\frac{h}{12}(f(a)+4f(\frac{3a+b}{4})+\frac{a+b}{2}+4f(\frac{a+3b}{4})+f(b))S=12h(f(a)+4f(43a+b)+2a+b+4f(4a+3b)+f(b))。

梯形积分 trapz T=h∑k=1n−1fk+h2(f0+fn)T=h \sum_{k=1}^{n-1} f_{k}+\frac{h}{2}\left(f_{0}+f_{n}\right)T=h∑k=1n−1fk+2h(f0+fn)

%% 梯形积分 trapz
% 正如我们在数值积分篇讲过的梯形积分, 对于列表，我们同样可以将公式带入求其积分
a = [1 4 5 7 8 9];
trapz(a)    % ans = 29
x = [0.3: 0.01: 0.4];
y=x.^2;
trapz(y)    % ans = 1.2335% 对于矩阵, trapz对列向量积分, trapz([],2)可对行向量积分
b = [1 3 4; 4 6 7; 2 1 7];
trapz(b)    % ans = [5.5000 8.0000 12.5000]%% 累计梯形积分 cumtrapz
% 累计梯形积分的第n个元素是原矢量前n个元素的梯形积分
cumtrapz(b) % ans = [0 0 0; 2.5000 4.5000 5.5000; 5.5000 8.0000 12.5000]

函数积分 quad, quadl

quad函数使用自适应辛普森算法，可以自动调整步长。

辛普森算法为 ∫abf(x)dx≈b−a6[f(a)+4f(a+b2)+f(b)]\int_a^b f(x)\mathrm{d}x \approx \frac{b-a}{6}\left[f(a)+4f(\frac{a+b}{2})+f(b)\right]∫abf(x)dx≈6b−a[f(a)+4f(2a+b)+f(b)]，而我们可以知道积分区间越小，其值越准确，但相应的也会增加计算量，因此我们决定均衡二者，当在区间[a,b]内的积分近似等于区间[a,c]与区间[c,b]的积分和即∣∫abf(x)dx−[∫acf(x)dx+∫cbf(x)dx]∣<ϵ\mid\int_a^b f(x)\mathrm{d}x-[\int_a^c f(x)\mathrm{d}x+\int_c^b f(x)\mathrm{d}x]\mid<\epsilon∣∫abf(x)dx−[∫acf(x)dx+∫cbf(x)dx]∣<ϵ时，该区间积分完成，否则对拆分后的小区间继续重复上述过程。这样的算法可自动适应被积函数，在被积函数变化剧烈处，将区间划分较密，在被积函数变化缓慢处，区间便比较稀疏了。

quadl函数则使用精度更高的自适应洛巴托(Lobatto)法，也称高斯积分法。

%% quad, quadl
clc;clear all;close all
% f=inline('exp(x.^(-3))','x');        % 建立积分函数
f=@(x)exp(x.^(-3));       % 使用@定义的匿名函数可以生成更高效的程序
format long
quadl(f,0.5,0.6)%% 二重积分 dblquad (double->dbl)
ff=@(x,y)(y.*sin(x)+x.*cos(y));
dblquad(ff,pi,2*pi,0,pi)%% 三重积分 triplequad
fff=@(x,y,z)(y*sin(x)+z*cos(y));
triplequad(fff,0,pi,0,1,-1,1)

符号工具箱

MATLAB也能使用符号运算工具箱实现求导函数和推导不定积分公式。

diff(differential) 符号变量微分，taylor 泰勒展开，int(integral) 符号变量积分，symsum 求级数和，limit 求极限，simple 将表达式化简。

符号工具箱内函数用法见：MATLAB解方程

总结

正确理解数值微分和数值积分的计算方法对数值计算极其重要，它们是数值计算的逻辑，如果对这部分概念不清楚则很可能在后续的编程中出现数值计算失控的情况。

数值计算与MATLAB微积分相关推荐

Matlab—微积分运算
Matlab-微积分运算 1. 符号方法 1.1 极限 1.2 求导 1.3 偏导 1.4 积分 2. 数值方法 2.1 梯度 2.2 定积分 2.3 高精度数值积分 2.4 重积分 2.5 计算椭圆 ...
matlab 画xos函数,振荡积分的数值计算与Matlab实现
第 29卷第 12期(下) 2013年 12月赤峰学院学报 (自然科学版 ) Journal of Chifeng University(Natural Science Edit ...
matlab微积分计算
matlab微积分计算 syms x定义一个变量求极限 limit(f, a)计算表达式极限时的值第2个参数为变量第3个参数为趋近的值,默认为0 第4个参数表示左极限或右极限 'left' 'r ...
matlab计算下列极限,MATLAB微积分计算极限，又快又好
有时候我们在使用MATLAB的时候,想微积分计算极限,怎么操作呢,下面来分享一下方法工具/材料电脑 MATLAB 操作方法 01 第一步我们使用MATLAB计算极限的limit函数,去计算函数f( ...
微分法MATLAB语言程序,matlab微积分运算命令与例题pdf
一.用Matlab进行微积分运算直接在word计算可以用notebook syms x n;s1=symsum(sin(x)/n^2, n,1, inf) s1 = (pi^2*sin(x))/6 ...
matlab计算微积分极限,Matlab微积分
MATLAB提供了解决微分和积分微积分的各种方法,求解任何程度的微分方程和极限计算.可以轻松绘制复杂功能的图形,并通过求解原始功能以及其衍生来检查图形上的最大值,最小值和其他固定点. 本章将介绍微积分 ...
matlab的数值计算功能,MATlAB数值计算功能
三．逆矩阵及行列式(Revers and determinant of matrix) 1．方阵的逆和行列式(Revers and determinant of square matrix) 若a ...
数值计算与matlab笔试样卷,数值分析试卷及其答案.doc
数值分析期末考试一. 设,若要确保其近似数的相对误差限为0.1%,则它的近似数至少取几位有效数字?(4分) 解:设有位有效数字. 因为,所以可得的第一位有效数字为8(1分) 又因为,令,可知至少具有 ...
matlab中复合中点式程序,《现代数值计算》Matlab程序整理（23页）-原创力文档
第2 章线性方程组的直接解法 1.1 写一个求解线性代数方程x b 的列主元素高斯消去法程序,在程序入口输入n 的 1 值,其中ij ,bj ln j ,1 i , j  n i  j  ...

数值计算与MATLAB微积分

目录

前言