首页 > 其他分享 >Mathematica 入门

Mathematica 入门

时间:2024-09-08 23:24:22浏览次数:11  
标签:入门 Mathematica Solve 使用 文档 MMA 表达式 函数

前言

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

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/ ↩︎

标签:入门,Mathematica,Solve,使用,文档,MMA,表达式,函数
From: https://www.cnblogs.com/stanleys/p/18403717/mathematica-intro

相关文章

  • SpringMVC入门概念
    SpringMVC是目前主流的实现MVC设计模式的框架,是Spring框架的一个分支产品,以Springloc容器为基础,并利用容器的特性来简化它的配置。什么意思呢?其实就是springmvc里需要用到各种各样的组件,然后这些组件它其实就是对象,那这些对象怎么去使用呢?就交给IOC容器来管理。也就是说......
  • Redis 入门 - 安装最全讲解(Windows、Linux、Docker)
    经过上一章节的介绍,相信大家对Redis已经有了大致的认知,今天主要给大家详细讲解Redis在Windows、Linux、Docker下的安装过程。01Windows下面给大家介绍三种在Windows环境下安装Redis的方式:官方建议方式、可执行文件方式、脚本方式。1、官方建议方式Redis官方是不支持直接......
  • Redis 入门 - 图形化管理工具如何选择,最全分类
    工欲善其事必先利其器,上一章Redis服务环境已经搭建完成,现在就需要一个趁手的工具,有个好工具可以做到事半功倍。Redis图形化管理工具五花八门,可供选择的很多,大家可以根据自己的需求应用场景进行选择。就是因为选择过多所以才写了这篇文章,和大家分享下到底有哪些工具可供选择,好给......
  • 第 18 场 小白入门赛(蓝桥)水浒专场!!
    第二题:情报传递1代码1(c++):#include<iostream>usingnamespacestd;intmain(){intt;cin>>t;while(t--){inta,b,c;cin>>a>>b>>c;intsteps=0;while(a<b){......
  • JavaWeb学习笔记,关于HTML的入门标签及属性
    一.HTML入门结构标签以及特点   以上标签即为HTML的入门标签,包括了HTML的基本框架结构标签以及部分常用标签,需要注意的是,HTML的语法松散,但我们更要严格要求自己,使用正常符合要求的代码格式,以免后期出现错误而无法及时发现问题,值得提起的还有,<h1>到<h6>是HTML中预定义好......
  • Java教程:入门基础【十万字详解】
    ✨博客主页:https://blog.csdn.net/m0_63815035?type=blog......
  • C++入门基础(内容太干,噎住了)
     文章目录1.缺省参数2.函数重载2.1重载条件:1.参数类型不同2.参数个数不同3.参数类型顺序不同 2.2不做重载条件情况:1.返回值相同时2.当全缺省遇见无参数3.引用3.1引用特性:3.2引用的使用1.缺省参数1.缺省参数是声明或定义函数时为函数的参数指定⼀个缺省值。......
  • redis基础——SpringDataRedis入门
    redis基础——SpringDataRedis入门最近在学习redis,学到了redis的java客户端,其中最常用的是Jedis和lettuce,而SpringDataRedis是spring整合了Jedis和lettuce的产物,它提供了RedisTemplate工具类,封装各种对redis的操作,将不同数据类型的API封装到了不同类型中,避免了代码臃肿。S......
  • 入门pytorch
    ###卷积神经网络模型 卷积神经网络(简称CNN)是一种专为图像输入而设计的网络。它最明显的特征就是具有三个层次,卷积层,池化层,全连接层。 借用一张图,下图很好的表示了什么是卷积(提取特征),什么是池化(减少数据量),而全连接层就是一个简单普通的神经网络。  如下代码,该代码定......
  • CTF入门教程(非常详细)从零基础入门到竞赛,看这一篇就够了!
       一、CTF简介CTF(CaptureTheFlag)中文一般译作夺旗赛,在网络安全领域中指的是网络安全技术人员之间进行技术竞技的一种比赛形式。CTF起源于1996年DEFCON全球黑客大会,以代替之前黑客们通过互相发起真实攻击进行技术比拼的方式。发展至今,已经成为全球范围网络安全圈流行的......