Maxima介绍

Size: px

Start display at page:

Download "Maxima介绍"

召允左
5 years ago
Views:

1 Maxima 介绍李东风北京大学数学科学学院目录 1 Maxima 初步认识 1 2 多项式和分式计算 6 3 三角函数化简 9 4 函数和微积分 12 5 方程 16 6 图形 19 Maxima 是什么 1 Maxima 初步认识计算机代数软件可以进行多项式化简分解因式分式化简通分积分微分等公式推导这种运算又称为符号计算, 与数值计算相对著名的计算机代数系统有 Maple Mathematica Maxima Maxima 是一个免费的计算机代数软件, 可安装在各种不同的操作系统中 1

2 1 MAXIMA 初步认识 2 具有计算机代数系统必须的各种功能用 Lisp 语言编制, 是一个有悠久历史的软件, 开发历史可追溯至 1968 年到 1982 年之间在 MIT 开发的 DOE Macsyma, 后来转由 Texas 大学 William F. Schelter 教授维护, 在 1998 年转为 GNU GPL 版权协议安装与运行在 MS Windows 系统中, 可以下载安装 wxmaxima 软件 ( 比如, 文件 maxima-clisp-sbcl win64.exe wxmaxima 运行界面是一个窗口形式, 在窗口中逐行输入命令, 命令以分号结尾, 用 Shift+RETURN 键运行每行命令命令用 $ 结尾则不显示命令的结果可以用 % 访问上一个公式的输出 ; 用 %i1, %o1 等访问历史的输入与输出用 %th(k) 访问向上数第 k 个结果用冒号定义一个变量运行界面示例精确计算四则运算符号用 + - * /, 乘方用 ˆ 表示 Maxima 中不带小数点的四则计算结果保存为有理式, 不作近似如 12/5 + 6*(2^3 + 1/11);

3 1 MAXIMA 初步认识 3 近似计算为了求以上结果的近似值, 用 numer 函数开关, 或 float 函数, 如 float(%); 或 12/5 + 6*(2^3 + 1/11), numer; 这样,Maxima 可以当作一个高级计算器使用

4 1 MAXIMA 初步认识 4 数值类型 Maxima 的数值类型有整数有理数普通浮点数高精度的浮点数普通浮点数如 1.2, 1.2E-3 高精度浮点数如 1.2B0, 1.2B-3, 有效位数用变量 fpprec 控制函数 float(x) 把 x 转为浮点型, 函数 bfloat(x)x 转为高精度浮点型有高精度浮点数参加运算的算式高精度浮点数 ; 没有高精度浮点数但有浮点数算式浮点数特殊常数用 %e %pi 表示 e 和 π, 用 %i 表示虚数单位 inf, minf, 分别表示 +,, infinity 表示复数无穷大复数用 %i 表示虚数单位, 如此可以输入复数, 如 z1 : 5 + 3*%i; 用 realpart 和 imagpart 返回复数的实部和虚部用 conjugate 返回复数的复共轭用 abs 求复数模, 用 carg 和求辐角

5 1 MAXIMA 初步认识 5 rectform 把复数表为直角坐标形式, 如 rectform(z1); 显示为 3 i + 5 polarform 把复数表为直角坐标形式, 如 polarform(z1); 显示为 34 e i atan( 3 5) 列表用方括号把逗号分隔的数值组合在一起作为一个列表, 如 [2,3,4] + [7,8,9] [9,11,13]

6 2 多项式和分式计算 6 多项式合并 2 多项式和分式计算考虑如下的一个推导问题甲乙比赛, 每局比赛甲赢的概率 p (0.5, 1) 问 : 三局两胜还是五局三胜对甲有利? 三局两胜中甲赢的概率为 P ( 三局两胜中甲赢 ) = P ( 甲赢 2 局 ) + P ( 甲赢 3 局 ) =C3p 2 2 (1 p) + p 3 = 3p 2 2p 3, 五局三胜中甲赢的概率为 P ( 五局三胜中甲赢 ) =P ( 甲赢 3 局 ) + P ( 甲赢 4 局 ) + P ( 甲赢 5 局 ) =C5p 3 3 (1 p) 2 + C5p 4 4 (1 p) + p 5 = 10p 3 15p 4 + 6p 5, 在 wxmaxima 中用如下两个命令定义了这两个概率 : P1 : 3*p^2*(1-p) + p^3; P2 : 10*p^3*(1-p)^2 + 5*p^4*(1-p) + p^5; 用 expand 函数可以把多项式合并同类项化简 : expand(p1); expand(p2); 如下命令计算两个多项式的差并合并同类项 : P3 : expand(p2-p1); 结果得到两个概率的差为 6p 5 15p p 3 3p 2

7 2 多项式和分式计算 7 多项式因式分解以上的两种赛制甲获胜概率之差简单地化为 3p 2 (2p 3 5p 2 + 4p 1), 为判断其是否在 p (0.5, 1) 总为正值, 可以分解因式 Maxima 用 factor 函数分解因式, 如 : factor(p3); 在 p (0.5, 1) 为正值 3 (p 1) 2 p 2 (2 p 1). 推导结果导出 wxmaxima 的结果显示为数学公式形式, 选中结果, 可以选编辑复制为 LaTeX 菜单, 把结果导出到 L A TEX 文件中也可以把公式复制为纯文本图片分式通分合并考虑如下分式 : R1 : (x+1)^2/(x-1) + 1/(x+1); (x + 1) 2 x x + 1

8 2 多项式和分式计算 8 expand 函数可以分别展开分子分母中的多项式, 但不能通分合并 : expand(r1); 1 x x2 x x x x 1 函数 ratexpand 可以把分式的分子分母中多项式合并同类项, 然后通分, 结果表示为分子每个不同幂次单独一个分式的加法 : ratexpand(r1); x 3 x x2 x x x 2 1 函数 ratsimp 把分式之和通分合并变成一个分式 : ratsimp(r1); x x x x 2 1 factor 函数可以对分式的分子分母分别进行因式分解, 如 factor(%); x (x x + 4) (x 1) (x + 1)

9 3 三角函数化简 9 三角函数化简 3 三角函数化简 Maxima 用如下函数进行三角函数化简 : trigexpand 和差化积 ; trigreduce 积化和差 ; trigsimp 利用 sin 2 x + cos 2 x = 1 之类化简 ; trigrat 把三角有理分式化为 sin 和 cos 的线性函数和差化积考虑如下三角函数有理式 : T1 : sin(2*x)/cos(x) + cos(2*x); 结果显示为 sin (2 x) cos (x) + cos (2 x) trigexpand 把 sin(2x) 展开为 2 sin x cos x, 把 cos(2x) 展开为 cos 2 x sin 2 x: T2 : trigexpand(t1); 结果显示为 sin (x) sin (x) + cos (x) 2

10 3 三角函数化简 10 上述结果中同时有 cos 2 x 和 sin 2 x, 用 trigsimp 化简 trigsimp(%); 2 sin (x) + 2 cos (x) 2 1 积化和差如下程序中 T3 : trigreduce(t2); sin 2 x, cos 2 x 被化为一次式, 结果显示为 cos (2 x) 用 expand 简单地合并同类项 : expand(%); cos (2 x) sin (x) 1 2 cos (2 x) + 2 sin (x) 用 trigsimp 或 trigrat 也可以合并同类项

11 3 三角函数化简 11 三角函数有理式化简对 T1 sin (2 x) cos (x) + cos (2 x) 这样的三角函数有理分式, 用 ratsimp 可以直接化简为三角函数线性式, 或分子分母都是三角函数线性项的有理式 : trigrat(t1); cos (2 x) + 2 sin (x) 又如, T4 : sin(x)^2 / cos(x) + cos(2*x); 为 cos (2 x) + sin (x)2 cos (x) 用 trigrat 作用 : T4 : sin(x)^2 / cos(x) + cos(2*x); cos (3 x) cos (2 x) + cos (x) cos (x) 三角函数有理分式 T4 又可以用 trigexpand 和 trigsimp 作用 : trigexpand(%);

12 4 函数和微积分 12 3 cos (x) sin (x) 2 + sin (x) 2 + cos (x) 3 cos (x) 2 + cos (x) cos (x) trigsimp(%); 2 cos (x) 3 cos (x) 2 cos (x) + 1 cos (x) 4 函数和微积分常用初等函数 sqrt() abs() max() min() sign() exp() log() sin() cos() tan() cot() sec() csc() asin() acos() atan() acot() asec() acsc() sinh() cosh() tanh() coth() sech() csch() asinh() acosh() atanh() acoth() asech() acsch() 组合函数用 x! 计算阶乘, 如 5!; 120 用 binomial(n, k) 计算 C k n, 如 binomial(5,2); 10

13 4 函数和微积分 13 自定义变量和自定义函数用冒号定义变量, 可以是数值, 也可以是代数式, 如 : g : 9.8; P1 : 3*p^2*(1-p) + p^3; 用 := 定义函数, 如 f(x,y) := 1/2*g*x^2*y; 结果显示为 f (x, y) := 1 2 g x2 y 调用如 f(2,10); 代换带有未知数的式子可以用逗号格式把未知数代换为数值或其他代数式, 如 P1 : 3*p^2*(1-p) + p^3; P1, p=0.5; P1, p=1; 后面两式的结果分别为 0.5 和 1 用 a 替换 p, 如

14 4 函数和微积分 14 expand(p1), p=1/2 + a; 2 a a 级数和用 sum 函数表示级数和, 如 sum(i^2, i, 1, 10); 而表示 10 i=1 i2, 385 sum(i^2, i, 1, n); 显示为 10 i 2 i=1 级数和化简用 simpsum 进行级数化简, 如 %, simpsum; 2 n n 2 + n 6

15 4 函数和微积分 15 连乘积 n i=m i2 可以表示为 product(iˆ2, i, m, n) 用 simpproduct 化简连乘积, 如 product(i, i, 1, n), simpproduct; n! 微分 diff(f(x), x, n) 求 d n f(x) dx n 如 diff(sin(x)*exp(2*x), x, 1); 2 e 2 x sin (x) + e 2 x cos (x) 又 diff(sin(x)*exp(2*x), x, 2); 3 e 2 x sin (x) + 4 e 2 x cos (x)

16 5 方程 16 不定积分 integrate(f(x), x) 计算不定积分如 integrate(x / (1+x)^2, x); log (x + 1) + 1 x + 1 定积分 integrate(f(x), x, a, b) 计算定积分如 integrate(x*sin(x)*exp(-x), x, 0, %pi); b a f(x) dx 又如 (π + 1) e π integrate(exp(-1/2*x^2), x, minf, inf); 2 π 5 方程

17 5 方程 17 一元方程用 solve(eqn, x) 求解关于未知数 x 的一元方程, 如 eq1 : a*x^2 + b*x + c = 0; 显示为 a x 2 + b x + c = 0 求解是符号运算 : eq1 : a*x^2 + b*x + c = 0; b2 4 a c + b [x =, x = 2 a b2 4 a c b ] 2 a 以上一般解代入具体参数计算, 如 %, a=1, b=2, c=1; [x = 1, x = 1] 又如, 以下的三次方程 ( 注意可以只给出方程左端 ) solve(x^3 + 2*x^2 + 3*x + 4, x);

18 5 方程 18 [ 5 x = 9 ( 5 2 x = ( 5 2 x = ( 3 ) i ( ) ) 1 3 ( 3 i ) ( ( ) 1 ( 35 3 ) 3 i , ) ( ) i ( ] ) ) , 二元方程组二元方程求解也使用 solve, 其中方程和未知数都用方括号表示的列表表示如 eq1 : y + 2*c = 0; eq2 : 2*x*y - c*y = 5; 对应于 y + 2 c = 0 2 x y c y = 5 求解如 solve([eq1, eq2], [x, y]);

19 6 图形 19 [[x = 2 c2 5, y = 2 c]] 4 c 结果表示为解的列表, 每个解又是包含 x, y 两部分的列表一个二元二次方程组的例子 : solve([ (x-1)^2 + y^2 = 4, x*y=1], [x, y]); [ [x = , y = ], [x = i , y = i ], [x = i , y = i ], [x = , y = ] ] 6 图形曲线图函数的曲线图使用 plot2d 函数, 用如 plot2d(sin(x), [x, 0, 2*%pi]); 可以指定纵坐标范围, 如

20 6 图形 20 plot2d(sin(x), [x, 0, 2*%pi], [y, -1.2, 1.2]); 图形文件作图结果可以复制到 Windows 剪贴板, 然后粘贴到 Word 中, 或粘贴到画图程序中再保存为图形文件也可以直接保存为图形文件, 如 plot2d(sin(x), [x, 0, 2*%pi], [gnuplot_term, png], [gnuplot_out_file, "test-save.png"]); 多条曲线多条曲线只要指定多个函数, 如 plot2d([sin(x), cos(x)], [x, 0, 2*%pi]);

21 6 图形 21 对数坐标轴 plot2d 中加选项 [logy] 使得 y 轴为对数轴, 加选项 [logx, logy] 使得 x, y 轴都为对数轴, 如 plot2d(2^x, [x, 0, 3]); plot2d(2^x, [x, 0, 3], [logy]); 参数方程作图在 plot2d 中把要绘图的函数用 [parametrix, y(t), x(t), [t, tmin, tmax]] 表示, 可以进行参数方程作图例如, 如下参数方程表示的曲线 x(t) = cos(t), y(t) = sin(3t), t [0, 2π] 可以作图如下 :

22 6 图形 22 plot2d([parametric, cos(t), sin(3*t), [t, 0, 2*%pi]], [x, -1.2, 1.2], [y, -1.2, 1.2], [nticks, 200]); 极坐标参数曲线可以化为直角坐标后参数作图, 如 [ z(t) = e cos t 2 cos 4t sin 5 程序如 t 12 ] e it, t [ 8π, 8π], r : (exp(cos(t)) - 2*cos(4*t) - sin(t/12)^5); plot2d([parametric, r*sin(t), r*cos(t), [t, -8*%pi, 8*%pi]], [nticks, 2000]); 三维曲面图

23 6 图形 23 三维图用 plot3d 如 plot3d(sin(sqrt(x^2+y^2))/sqrt(x^2+y^2), [x, -12, 12], [y, -12, 12])$ 三维参数曲面图

24 6 图形 24 三维图用 plot3d 和参数方程, 有两个参数如 plot3d([cos(x)*(3 + y*cos(x/2)), sin(x)*(3 + y*cos(x/2)), y*sin(x/2)], [x, -%pi, %pi], [y, -1, 1], ['grid, 50, 15])$

1 2 5.2.................................... 13 5.3....................................... 14 5.4............................... 15 6 Maxima 17 6.1....

1 2 5.2.................................... 13 5.3....................................... 14 5.4............................... 15 6 Maxima 17 6.1.... Maxima dbzhang 1 2 2 3 2.1 Maxima..................................... 3 2.2......................................... 4 2.3 Maxima........................................ 4 2.4........................................