Mathematica 入门

news/2024/10/4 3:21:30

前言

Wolfram Mathematica(简称 MMA),是由 Wolfram Research 开发的科学计算软件。本文我们将介绍 Mathematica 的界面、语法和基本应用。类似的软件还有 MATLAB 和 Maple 等。

MMA 官网:https://www.wolfram.com/mathematica/

MMA 的安装及激活:Mathematica安装激活极简教程 - 科研小飞 (知乎)

本文适合有一定编程基础的读者。当然,如果你不了解编程也没关系,直接跳过相关部分即可。屏幕截图全部来自 Windows 11,Mathematica 13.3。(最新版本 14.0 已经在 2023 年 12 月发布,语法向下兼容)

演示使用的 MMA 版本
演示使用的 MMA 版本

1. 界面

安装并正确激活 MMA 后,新建笔记本,应该能看到如下的界面:

MMA 基本界面

一张图带你快速熟悉 MMA 的界面:

MMA 界面简介

几个需要注意的点:

  • Enter 默认换行,执行代码需要用 Shift+Enter
  • 指令面板左侧的红三角形也可执行代码,灰色的正方形用于终止正在进行的运算。
  • MMA 带有 Notebook(笔记本)功能,由于其主要用于计算而不是笔记,本文不详细介绍,若需要使用可自行查阅相关资料。

对于任意输入或输出,按下Ctrl+Shift+I即可显示输入(代码)形式,按下Ctrl+Shift+N即可显示数学形式。后面会对两种显示形式进行详细讲解。

2. 语法

Mathematica 实质是一个 Wolfram Language 的解释器,所以使用它之前必须学习 WL 的语法。
这里介绍基本语法和常用的指令。

2.1 帮助文档

首先来了解帮助文档的使用。输入 ? Solve 来获取 Solve 函数的说明:

help-solve

? 后面可以跟任何函数名,MMA 会直接在笔记本中输出简化版的帮助文档。可以选择输出右上角的 i 来获取更详细的文档(优先打开离线文档,不存在则默认打开在线文档)。菜单栏中的 “帮助 -> Wolfram 参考资料” 可以打开完整版文档。

2.2 注释

(* 注释 *)

如上,MMA 的注释使用 (**) 来标注,用法类似于 C/C++ 中的 /**/。注释可以添加在代码的任何位置,它们会被自动忽略。例如:

1 + 1 (*Hello*)
1 + (*World*)1 

这两行代码在实际执行中与 1 + 1 无区别:

comment-test

2.3 括号

在学习 MMA 的函数之前,一定要先学会括号的使用!!
官方参考文档:正确使用括号和大括号

MMA 中共有四种括号,分别为 ()[]{}[[]],具体作用如下:

圆括号 ()

用作对表达式编组和确定运算的优先次序:

圆括号的使用

方括号 []

表示函数的调用和传参:

方括号的使用

大括号 {}

表示列表:

大括号的使用

列表的具体使用后续会介绍。

双方括号 [[]]

表示对列表元素的访问,实质是函数 Part 的简写形式:

双方括号的使用

同样会在后面具体介绍。

2.4 运算与表达式

MMA 支持基本的数学运算,见下表:

名称 符号形式 函数形式 数学形式
加法 a + b Plus[a, b] \(a+b\)
减法 a - b Subtract[a, b] \(a-b\)
乘法 a * b[1] Times[a, b] \(a\times b\)
除法 a / b Divide[a, b] \(\frac ab\)
a ^ b Power[a, b] \(a^b\)
模余 - Mod[a, b] \(a\bmod b\)

数字和符号都可以参与运算:

数字和符号的运算

表达式是一个或多个运算/函数调用的组合,使用()指定优先级。上图中 4*5a+a3 x y/y 都是合法的表达式。

需要注意的是,MMA 只会对输入的表达式进行约分、合并同类项,而不会自动执行需要展开的化简。必要时可以使用 SimplifyFullSimplify 函数:

化简示例

关于函数的使用,后面会详解。

布尔表达式

表达式也可以是等式或不等式:

等式和不等式

注意等式用两个等号 == 连接,单个等号表示变量赋值。
不等式可以用 !=(不等于)、>(大于)、>=(大于等于)、<(小于)、<=(小于等于)连接。!=>=<= 在输入时会自动转换为相应的数学符号。

这类返回值为 True(真)或 False(假)的表达式,我们统称为布尔表达式。可以用 &&|| 运算符来连接两个布尔表达式,分别表示“与”“或”,所得结果仍为布尔表达式。

绝大部分布尔表达式都不能自动化简,但可以通过 SimplifyFullSimplifySolveReduce 来化简或求解:

布尔表达式的化简与求解

2.5 变量与常数

MMA 中,使用 变量名 = 值 的形式来声明或赋值变量:

创建变量

由上方 10a + 3 的计算结果可知,表达式中所有已声明的变量都会被替换成变量值。变量的声明也可以包含其他变量和符号,此时仍符合变量替换法则:

替换

用于解方程的 Solve 函数没有返回 \(a\) 的解,而是返回 \(x\) 的解,说明表达式在解析时 \(a\) 被自动替换为了 \(x\)

除此之外,还可以用类似于 Python 的语法,同时赋值多个变量

{a, b, c} = {1, 2, 3}

这行代码可以把变量 \(a,b,c\) 分别赋值为 \(1,2,3\)

若要取消声明一个变量,可以使用 Clear[变量名]

Clear的使用

在使用 MMA 的过程中,我们会发现某些特定符号不能声明为变量,且有一个固定不变的值。这些是内置常量

常量 名称 近似值 数学形式
Pi 圆周率 \(3.141592654\) \(\pi\)
E 自然常数 \(2.718281828\) \(\mathrm e\)
I 虚数单位 - \(i\)
Infinity 无穷大 - \(\infty\)
Degree - \(\degree\)

2.6 函数

在前面的讲解中,我们多次提到函数。现在,让我们详细解读函数的使用。

调用

要调用一个函数,只需输入 函数名[参数1, 参数2, ...] 即可:

函数的调用1

对于只有一个参数的函数,可以用 参数 // 函数名 的形式调用。这种形式常用于 Simplify(化简)、FullSimplify(完全简化)和 N(数值运算)的调用:

函数调用2

这种写法的最大优点在于可以链式调用:

链式调用

定义

说完了函数的调用,再来说说新函数的定义。

MMA 中,定义新函数的语法为 函数名[参数1_, 参数2_, ...] := 返回值。注意每个参数名后面都要加一个下划线_)。

定义新函数推荐用 :=,但使用 = 也可行。

定义新函数

拓展 - 定义递归函数
参考下面定义斐波那契数列的格式:
斐波那契数列

3. 内置函数

MMA 中有将近 \(6000\) 个内置函数[2],它们都以大写字母开头。下面介绍几种常用的内置函数:

3.1 数值运算 N

数值运算函数 N

  • N[expr] 给出 \(\mathrm{expr}\) 的数值值
  • N[expr, n] 尝试给出具有 \(n\) 位精度的结果

使用示例:

N示例

3.2 三角函数

MMA 内置了各类三角函数,它们全部使用习惯名称且首字母大写(SinSinhArcSin),这里不一一阐释。
只说一个注意点,MMA 中三角函数的参数默认是弧度,若要使用角度必须加上 Degree 单位:

三角函数

3.3 解方程和不等式 Solve/Reduce

解方程 Solve

Solve解方程

解不等式用 Reduce

Reduce解不等式

Reduce 还有更高级的使用,可以约化各种表达式:

Reduce的高级应用

这句话说的是:求解 \(a,b\) 的范围,使得存在 \(x\) 满足等式 \(ax=b\)。实际上就是在求解使一元一次方程有解的参数值。注意 SolveReduce 不一定总返回最简形式的表达式,很多情况下需要手动调用 SimplifyFullSimplify 进行化简。感兴趣的读者可以自行测试上面的例子中不使用 Simplify 的返回结果。

3.4 解方程的其他方法 NSolve/FindInstance

NSolve 用法同 Solve,但会返回数值解而不是精确解。相当于 Solve[...] // N
FindInstance 用法同 Solve,但添加了一个参数表示至多返回解的个数(默认为 \(1\)):

FindInstance的使用

3.5 极点值 Maximize/Minimize

Minimize 返回函数(在限制条件下)的最小值以及取到最小值的(一种)变量值:

Minimize的使用

如上图中的示例,Minimize 支持多个变量,可以指定条件,也可指定求解域。
Maximize 用法完全相同,返回最大值,此处略过。

3.6 偏导 D

偏导 D 的两种最常用用法:

  • D[f, x] 给出 \(f\) 关于 \(x\) 的偏导。
  • D[f, {x, n}] 给出 \(f\) 关于 \(x\)\(n\) 阶偏导。

示例略。注意 f'[x] 可以直接求出 f[x] 的导数:

符号求导

3.7 积分 Integrate

积分 Integrate

  • Integrate[f, x] 给出不定积分 \(\int f~\mathrm{d}x\)
  • Integrate[f, {x, x_min, x_max}] 给出定积分 \(\int_{x_{min}}^{x_{max}} f~\mathrm{d}x\)
  • Integrate[f, {x, x_min, x_max}, {y, y_min, y_max}, ...] 给出多重积分 \(\int_{x_{min}}^{x_{max}} \mathrm{d}x \int_{y_{min}}^{y_{max}} \mathrm{d}y\dots f\)

一个简单的定积分示例:

Integrate求定积分

3.8 展开 Expand/ExpandAll

Expand 很好理解,Expand[expr] 会展开表达式 \(\mathrm{expr}\) 中的乘积和正整数幂。限于篇幅,这里不提供使用范例,可参考官方文档。
ExpandAllExpand 的基础上,会展开表达式中任意位置的乘积和整数幂。如表达式 Sin[(1 + x)^3]Expand 不会展开其中的 (1+x)^3,而 ExpandAll 会。

3.9 因式分解 Factor

Factor[poly]:在整数上对一个多项式分解因式。使用示例参考官方文档。

3.10 绘图 Plot/Plot3D

PlotPlot3D 的用法较为复杂,这里只介绍最基本的函数绘图:

Plot和Plot3D

两个语句分别绘制:

  • \(\sin(x)\) 的图像,其中 \(-2\pi\le x\le 2\pi\)
  • \(\sin(x+y^2)\) 的图像,其中 \(-3\le x\le 3\)\(-2\le y\le 2\)

4. 总结

我自己当初学习 MMA 时,被网上杂乱的教程搞得混乱不堪。因此就想写一篇教程,涵盖所有常用语法知识点,并让初学者避开我自己踩的一些坑。
初衷是写个简明的教程,结果一写就是八千多字…… 也感谢认真读到这里的读者们!
后续可能还会更一些用 Mathematica 解决数学和实际问题的文章,敬请期待!


  1. 乘法也可简写为 a b(中间必须有空格)。MMA 中,大部分空格可省略,但是乘法的空格不能省略(若写作没有空格的 ab 会被认为是一个变量)。 ↩︎

  2. https://www.wolfram.com/language/fast-introduction-for-programmers/en/built-in-functions/ ↩︎

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.ryyt.cn/news/58132.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈,一经查实,立即删除!

相关文章

【算法笔记】位运算详解

0. 前言 突然想到位运算是个好东西,就来水一波文章了…… 注意:我把能想到的有关位运算的所有内容都放进来了,所以篇幅较长,请谅解!若有写的不清楚或者不够详细的地方欢迎在评论区补充,谢谢支持! 本文中参考代码均使用C++编写。 废话不多说,下面步入正题。1. 基本运算 …

【算法笔记】【专题】RMQ 问题:ST表/树状数组/线段树

0. 前言 好久没更算法笔记专栏了,正好学了新算法来更新…… 这也是本专栏的第一个专题问题,涉及到三种数据结构,如果写得有问题请各位大佬多多指教,谢谢! 1. 关于 RMQ 问题 RMQ 的全称是 Range Minimum/Maximum Query,即区间最大/最小值问题。 本文中,我们的算法以求最大…

js的对象

1.js有var开辟数组使用var a=[1,2,3,4,5,6]或者var a=new Array(1,1,1,1,1); 2.Array支持push添加操作,和indexof 查找,foreach遍历操作pop删除操作

【算法笔记】树形DP算法总结详解

0. 定义 树形DP,又称树状DP,即在树上进行的DP,是DP(动态规划)算法中较为复杂的一种。 1. 基础 令\(f[u]=~\)与树上顶点\(u\)有关的某些数据,并按照拓扑序(从叶子节点向上到根节点的顺序)进行\(\text{DP}\),确保在更新一个顶点时其子节点的dp值已经被更新好,以更新当前…

Seay安装和初步使用

作者网站现在已经无法访问:http://www.cnseay.com/2951/, 可以使用这个GitHub - f1tz/cnseay: Seay源代码审计系统 下载完安装包之后,解压到自己想要的电脑路径即可,无需进行任何额外的配置 利用工具对sql注入进行分析 进入软件之后,点击新建项目,选择需要分析的文件(这里…

代码整洁之道--读书笔记(5)

代码整洁之道简介: 本书是编程大师“Bob 大叔”40余年编程生涯的心得体会的总结,讲解要成为真正专业的程序员需要具备什么样的态度,需要遵循什么样的原则,需要采取什么样的行动。作者以自己以及身边的同事走过的弯路、犯过的错误为例,意在为后来者引路,助其职业生涯迈上更…

02网络参考模型

02网络参考模型02网络参考模型常见网络模型因为 OSI协议栈比较复杂 ,且TCP和IP两大协议在业界被广泛使用,所以 TCP/IP参考模型 成为了互联网的主流参考模型。OIS网络模型层级作用7.应用层 应用层 对应用程序提供接口。6.表示层进行数据格式的转换,以确保一个系统生成的应用层…

《痞子衡嵌入式半月刊》 第 107 期

痞子衡嵌入式半月刊: 第 107 期这里分享嵌入式领域有用有趣的项目/工具以及一些热点新闻,农历年分二十四节气,希望在每个交节之日准时发布一期。 本期刊是开源项目(GitHub: JayHeng/pzh-mcu-bi-weekly),欢迎提交 issue,投稿或推荐你知道的嵌入式那些事儿。 上期回顾 :《…