MATLAB科学计算科学计算03(微积分与泰勒傅里叶级数)(微积分与泰勒傅里叶级数)
文章目录文章目录1 微积分解析解1.1 单变量函数的极限1.2 多变量函数的极限2 级数2.1 Taylor泰勒 幂级数展开2.2 Fourier傅里叶级数
展开2.3 级数求和3 数值微分4 数值积分4.1 梯形法:4.2 单变量数值积分4.3 integral 数值积分4.4 quadgk()函数(自适应高斯-
勒让德积分法)4.5 双重积分的数值解4.6 三重定积分的数值求解
1 微积分解析解微积分解析解
1.1 单变量函数的极限单变量函数的极限
求解析解:要先将变量符号化;求解析解:要先将变量符号化;
syms x y a b
x0处的极限
格式1: L= limit( fun, x, x0)
x0处左右极限
格式2: L= limit( fun, x, x0, ‘left’ 或 ‘right’)
无穷处
L=limit(fun,x,inf)
1.2 多变量函数的极限多变量函数的极限
格式: L1=limit(limit(f,x,x0),y,y0)
或 L1=limit(limit(f,y,y0), x,x0)
函数的导数和高阶导数
格式: y=diff(fun,x) %求导数
y= diff(fun,x,n) %求n阶导数
多元函数的偏导:
格式: f=diff(diff(f,x,m),y,n)
或 f=diff(diff(f,y,n),x,m)
简化,化简函数
simplify(fun)
美化函数
pretty(fun)
隐函数的偏导数:
格式:F= - diff(f,xj)/diff(f,xi)
不定积分:
格式: F=int(fun,x)
定积分与无穷积分计算:
格式: I=int(f,x,a,b)
格式: I=int(f,x,a,inf)
例题,多重积分解析解,matlab求解
syms x y z
int(int(int(4xzexp(-x2*y-z2),x,0,1),y,0,pi),z,0,pi)
ans =
-(exp(-pi^2) – 1)(eulergamma + log(pi) – ei(-pi))
Ei(n,z)为指数积分,无解析解,但可求其数值解:
vpa(ans,60)