Maxima 对各种微积分的运算提供了强有力的支持。
可以这么说,在基本微积分运算能力上,Maxima 不输给任何商业软件。
### 求极限
求极限是微积分中最基本的运算。在Maxima 中用 limit 函数完成求极限的工作。
limit (expr, x, val)
计算 x 趋近于 val 时 expr 的极限。
如果极限不存在,范围值可能为 und (极限不存在),ind (极限不存在但是有界),infinity(发散)
limit (expr, x, val, dir)
计算 x 从 dir 方向趋近于 val 时 expr 的极限。
dir 为 plus 时求的是右极限,为 minus 时求的是左极限。
'limit(sin(x)/x,x,0);
limit(sin(x)/x,x,0);
'limit((1+1/x)^x,x,inf);
limit((1+1/x)^x,x,inf);
![](https://box.kancloud.cn/2016-01-24_56a4233f1918d.png)
limit((2^(1/x)+1)/(2^(1/x)-1),x,0);
'limit((2^(1/x)+1)/(2^(1/x)-1),x,0,minus);
limit((2^(1/x)+1)/(2^(1/x)-1),x,0,minus);
'limit((2^(1/x)+1)/(2^(1/x)-1),x,0,plus);
limit((2^(1/x)+1)/(2^(1/x)-1),x,0,plus);
![](https://box.kancloud.cn/2016-01-24_56a4233f252f4.png)
'limit((2*x^2+x)/(x^2-2),x,inf);
limit((2*x^2+x)/(x^2-2),x,inf);
![](https://box.kancloud.cn/2016-01-24_56a4233f30af6.png)
### 微分运算
diff 函数可以实现微分运算,有四种基本形式,最基本的形式是:
diff (expr, x) 求表达式 expr 对 x 的微分。
diff(sin(x)*x^3,x);
![](https://box.kancloud.cn/2016-01-24_56a4233f3c8a3.png)
diff(u(x)*v(x),x);
diff(u(x)*v,x);
![](https://box.kancloud.cn/2016-01-24_56a4233f479dc.png)
从上面的例子可以看到,一个函数如果不明确表明自变量,diff 函数是不会对其求导的。
diff (expr, x, n)
求表达式 expr 对 x 的 n 次微分。
diff(y^3*exp(-y^2),y,2);
![](https://box.kancloud.cn/2016-01-24_56a4233f53216.png)
diff (expr, x_1, n_1, …, x_m, n_m) 求的是混合偏导数。
diff(f(x,y),x,2,y,1);
![](https://box.kancloud.cn/2016-01-24_56a4233f5e386.png)
diff (expr) 计算的是全微分。
diff(f(x,y));
![](https://box.kancloud.cn/2016-01-24_56a4233f6c745.png)
### taylor 级数展开
函数f(x)的在x = a附近的幂级数可以通过powerseries (f(x), x, a)获得。
powerseries (f(x), x, a);
powerseries(1/(1-x^2), x, 0);
![](https://box.kancloud.cn/2016-01-24_56a4233f7adf6.png)
上面得到的结果中的求和指数 i2 看起来显得不那么专业,可以用 niceindices 函数将其变的看起来更专业些。
niceindices(powerseries(1/(1-x^2), x, 0));
![](https://box.kancloud.cn/2016-01-24_56a4233f895c5.png)
很多时候我们无法得到级数的解析表示,这时候可以用 taylor (f(x), x, a, n)得到函数f(x)在x = a附近第 n 阶项((x - a)^n)以下各项的泰勒级数
taylor(sin(x), x, 0, 8);
![](https://box.kancloud.cn/2016-01-24_56a4233f97921.png)
同样,对多元函数也可以进行 taylor 展开。
taylor (sin (y + x), x, 0, 3, y, 0, 3);
![](https://box.kancloud.cn/2016-01-24_56a4233fa46ae.png)
利用 pade 近似可以将 taylor 级数转化为多项式函数。比如下面的例子
taylor(sin(x), x, 0, 8);
pade(%,5,5);
![](https://box.kancloud.cn/2016-01-24_56a4233fb06c8.png)
### 积分运算
integrate(expr, x)
计算不定积分。
integrate(expr, x, a, b)
计算定积分,积分上下限分别为 a b。
integrate(1/(1+x^2),x);
integrate(1/(1+x^2),x, minf, inf);
![](https://box.kancloud.cn/2016-01-24_56a4233fbe08d.png)
risch (expr, x)
采用 Risch 方法计算不定积分。
risch (x^2*erf(x), x),ratsimp;
![](https://box.kancloud.cn/2016-01-24_56a4233fcb5e0.png)
changevar 函数利用变量替换,对待积分函数进行换元。这个功能对学习换元法求积分还是蛮有用的。
'integrate(x^2/(1+x^2),x,0,2);
changevar(%,x=tan(theta),theta,x);
![](https://box.kancloud.cn/2016-01-24_56a4233fd8739.png)
### 求和与求积运算
sum (expr, i, i_0, i_1)
计算求和。
product (expr, i, i_0, i_1)
计算求积。
sum (i^2, i, 1, 7);
sum (a[i], i, 1, 7);
sum(1/k^2,k,0,inf);
sum (1/3^i, i, 1, inf), simpsum;
![](https://box.kancloud.cn/2016-01-24_56a4233fe8ead.png)
product (x + i*(i+1)/2, i, 1, 4);
product (a(i), i, 1, n);
product (k, k, 1, n), simpproduct;
![](https://box.kancloud.cn/2016-01-24_56a4234002964.png)
### laplace 变换
Maxima 中提供了两个函数用来计算 laplace 变换和反变换。
laplace (expr, t, s) 对表达式 expr 进行 laplace 变换,积分变量为 t,变换后的参数为 s.
ilt (expr, s, t) 对表达式进行反变换。
需要说明的是这两个函数进行的都是单边 laplace 变换,Maxima 中还没有对双边拉式变换的支持。
laplace (exp (2*t + a) * sin(t) * t, t, s);
laplace ('diff (f (x), x), x, s);
ilt(1/(s+1),s,t);
![](https://box.kancloud.cn/2016-01-24_56a4234010d74.png)
提到 laplace 变换就不能不提提部分分式展开。
在大多数介绍 laplace 变换的初级课本中都只介绍利用部分分式展开的方式进行反变换。
Maxima 提供了 partfrac 函数来完成部分分式展开的工作。下面用一个例子来说明:
x/(x^3+4*x^2+5*x+2);
partfrac (x/(x^3+4*x^2+5*x+2), x);
![](https://box.kancloud.cn/2016-01-24_56a423401defa.png)
- 前言
- gnuplot 中的数据平滑
- 利用 Maxima 求解常微分方程
- Maxima 的绘图功能 1
- Maxima 的绘图功能 2
- Maxima 的绘图功能 3
- Maxima 矩阵及矢量运算 1
- Maxima 矩阵及矢量运算 3
- Maxima 矩阵及矢量运算 4
- maxima 学习笔记(coeff 函数的应用)
- maxima 代数表达式变换
- maxima 学习笔记(基本微积分)
- maxima 代数方程求解
- Maxima 中的复数运算
- Maxima 的三角函数化简功能
- Maxima 的基本微积分操作
- 利用 gnuplot 绘制时间序列图
- gnuplot 入门教程 1
- gnuplot 入门教程 2
- gnuplot 入门教程 3
- gnuplot 入门教程 4
- gnuplot 读取逗号分隔的数据文件
- scilab 读取处理 wav 文件
- scilab 读取处理 wav 文件 (2)
- Scilab 处理声音数据(补充)
- 利用 SCILAB 设计 FIR 滤波器(窗函数法)
- 利用 SCILAB 设计 FIR 滤波器(Minimax法)
- 利用 SCILAB 设计 iir 滤波器设计(模拟滤波器双线性变换法)
- Scilab 的绘图函数(1)
- Scilab 的绘图函数(2)
- Scilab 的绘图函数(3)
- Scilab 的绘图函数(4)