一种根据实验室仪器导出的双列 ASCII 数据集进行同时绘图的方法

本贴最后更新于 872 天前,其中的信息可能已经东海扬尘
mypath = ""; filename = {"", "", ""}; reflectances = Drop[#, 90]& @ Import[mypath <> # <> ".asc", "Table"]& /@ filenames; maxt = Max[ #[[2]]& /@ Select[#, #[[1]] > 600 &] ]& /@ reflectances; (* maxt[[-1]] = 103; *) absorbances = Table[ {#[[1]], maxt[[i]] - #[[2]]}& /@ #[[i]], {i, Length[#]}]&[reflectances]; maxa = Max[ #[[2]]& /@ Select[#, 300 < #[[1]] < 600 &] ]& /@ absorbances; (* maxa[[1]] = 100; *) unified = Table[ {#[[1]], #[[2]] / maxa[[i]]}& /@ #[[i]], {i, Length[#]}]&[absorbances]; plot1 = ListLinePlot[unified, PlotRange -> All, PlotLegends -> Placed[ LineLegend[ filenames, LegendFunction -> "Frame", LegendLayout -> "Column"], {{0.75, 1}, {0.5, 1}}], AxesLabel -> { ToExpression["\\lambda (\\rm{nm})", TeXForm, HoldForm], "Absorbance" }] Export[mypath <> "diagram_uniform.svg", plot1];

主要是在这其中学习 Mathematica 的语法。

Mathematica 的配置:

  • 下载 Wolfram Script
  • 下载 VSCode,安装其中的 Wolfram Language Notebook 扩展。
  • 总之照着它的说明配置就是了。

先说它的运算符吧。在 Mathematica 里面,字符串的连接用的不是 + 而是 <> 。因此 mypath <> "diagram_uniform.svg" 就是把我主动填写的目录字符串(注意这里代码显示的是空字符串,但具体执行的时候要补全!)和具体的文件名连在一起。

非常推荐把作出来的图 plot1 用 SVG 或 EPS 格式导出为矢量图。默认的位图 dpi 太低,导出的图特模糊。

然后是箭头 -> 。它在 Mathematica 官方付费的记事本里面会被自动转换为 ,我是相当的不喜欢这个设计,它会让初见者为如何输入它感到费解。但是 VSCode 里面就会忠诚地显示为 -> 啦,令人安心。它的意思是规定那些【可选】参数的性质。比如说 PlotLegends -> filenames 就是用 filenames 序列的每一个字符串逐一给 ListLinePlot 的主参数 unified 所对应的每条曲线标注在图例上。理所当然的这整一句都可以去掉,不影响作图,只是图例没了。

然后是括号 (), [], {}, [[]] 。圆括号 () 只用来支撑运算顺序,单方括号 [] 构筑成了几乎所有的函数调用,花括号 {} 用来构筑序列,包括矩阵,总之就是搭建有序元组,而双方括号 [[]] 则用来从有序元组中取出元素,是从 1 开始计数的。

然后是 @, /@ 。函数调用不是必须用 [] 的,它还有非括号的表达方式用的是 @ 比如输入 Sin @ Pi 会得到 。而相应地,函数调用就可以在序列上进行,相应的用的是 /@。比如输入 Sin @ {Pi/6, Pi/4, Pi/3, Pi/2, Pi} 会得到

In[1]:= Sin @ (Pi/6) 1 Out[1]= - 2 In[2]:= Sin @ {Pi/6, Pi/5, Pi/4, Pi/3, Pi/2, Pi} 1 5 Sqrt[5] 1 Sqrt[3] Out[2]= {-, Sqrt[- - -------], -------, -------, 1, 0} 2 8 8 Sqrt[2] 2

值得注意的是,这个 /@ 算符(或者说 f /@ expr )的波兰式(把算符写在参与运算各项的前面)表达就是 Map[f, expr],它把函数 f 作用到序列 expr 第一层的每一个元素上,得到一个元素数目和 expr 相同的序列。也就是说,它就是 map。在我的数据处理中,map 实在是发挥了重大作用。

然后是 #, & 。在我之前的帖子氢原子的 d 轨道概率密度函数的三维密度图 中就显式地定义过含参函数,诸如 r[x_, y_, z_] := Sqrt[x^2 + y^2 + z^2]; ,这里的 x_ 以下划线结尾表示对傀儡变元的声明,而之后应用该变元时则使用它去掉末尾下划线的形式 x。而这里我们想把事情做得更加简单,匿名函数就是以 & 结尾的一串表达式(不过如果它包含加减乘除等运算,那么还是需要用圆括号框住运算结果,这和函数的非括号调用遥相呼应,比如 (#+1)& @(5/2)),如果它参数只有一个,那么用 # 表示,否则依序用 #1, #2, ... 表示。比如匿名函数 #[[2]]&,它提取出来唯一参数的第 2 个元素。看上去搁这儿搁着,但是想把一列数对 pts 的第 2 个元素一起提出来时,对它的使用 #[[2]]& /@ pts 就可以简捷达成目的。除此之外,我程序中 Drop[#, 90]& @ Import[mypath <> # <> ".asc", "Table"]& 的使用就是把函数名 Drop 和第二参数 90 放在一起,而不是被超长的第一参数隔开,由此达到方便阅读的功效。不要紧,Mathematica 解释器会知道两个 # 和两个 & 该如何配对,大不了用圆括号。

除此之外,由于使用了相当长的变量名,所以诸如 unified = Table[ {#[[1]], #[[2]] / maxa[[i]]}& /@ #[[i]], {i, Length[#]}]&[absorbances]; 的处理,能有效地减小代码长度,方便更改变量名时少改几处,以及把作用的变量提到行尾,更加明显。

好了,我觉得我能讲的就这么多。其他诸如 Select[] 函数的使用啊,PlotLegends -> Placed[ LineLegend[]] 的具体细节啊,我都是偷的,是 🇰🇷 的,读者自己去看 MAA 官方文档就行了。

附赠一个效果图。

diagramuniform.jpg

相关帖子

欢迎来到这里!

我们正在构建一个小众社区,大家在这里相互信任,以平等 • 自由 • 奔放的价值观进行分享交流。最终,希望大家能够找到与自己志同道合的伙伴,共同成长。

注册 关于
请输入回帖内容 ...

推荐标签 标签

  • 开源中国

    开源中国是目前中国最大的开源技术社区。传播开源的理念,推广开源项目,为 IT 开发者提供了一个发现、使用、并交流开源技术的平台。目前开源中国社区已收录超过两万款开源软件。

    7 引用 • 86 回帖
  • 小薇

    小薇是一个用 Java 写的 QQ 聊天机器人 Web 服务,可以用于社群互动。

    由于 Smart QQ 从 2019 年 1 月 1 日起停止服务,所以该项目也已经停止维护了!

    34 引用 • 467 回帖 • 759 关注
  • Redis

    Redis 是一个开源的使用 ANSI C 语言编写、支持网络、可基于内存亦可持久化的日志型、Key-Value 数据库,并提供多种语言的 API。从 2010 年 3 月 15 日起,Redis 的开发工作由 VMware 主持。从 2013 年 5 月开始,Redis 的开发由 Pivotal 赞助。

    286 引用 • 248 回帖 • 18 关注
  • 正则表达式

    正则表达式(Regular Expression)使用单个字符串来描述、匹配一系列遵循某个句法规则的字符串。

    31 引用 • 94 回帖 • 1 关注
  • 新人

    让我们欢迎这对新人。哦,不好意思说错了,让我们欢迎这位新人!
    新手上路,请谨慎驾驶!

    52 引用 • 228 回帖
  • 以太坊

    以太坊(Ethereum)并不是一个机构,而是一款能够在区块链上实现智能合约、开源的底层系统。以太坊是一个平台和一种编程语言 Solidity,使开发人员能够建立和发布下一代去中心化应用。 以太坊可以用来编程、分散、担保和交易任何事物:投票、域名、金融交易所、众筹、公司管理、合同和知识产权等等。

    34 引用 • 367 回帖 • 1 关注
  • Node.js

    Node.js 是一个基于 Chrome JavaScript 运行时建立的平台, 用于方便地搭建响应速度快、易于扩展的网络应用。Node.js 使用事件驱动, 非阻塞 I/O 模型而得以轻量和高效。

    139 引用 • 269 回帖 • 2 关注
  • Outlook
    1 引用 • 5 回帖 • 1 关注
  • 倾城之链
    23 引用 • 66 回帖 • 160 关注
  • 一些有用的避坑指南。

    69 引用 • 93 回帖
  • DevOps

    DevOps(Development 和 Operations 的组合词)是一组过程、方法与系统的统称,用于促进开发(应用程序/软件工程)、技术运营和质量保障(QA)部门之间的沟通、协作与整合。

    56 引用 • 25 回帖 • 5 关注
  • V2Ray
    1 引用 • 15 回帖 • 1 关注
  • SQLServer

    SQL Server 是由 [微软] 开发和推广的关系数据库管理系统(DBMS),它最初是由 微软、Sybase 和 Ashton-Tate 三家公司共同开发的,并于 1988 年推出了第一个 OS/2 版本。

    21 引用 • 31 回帖 • 3 关注
  • Scala

    Scala 是一门多范式的编程语言,集成面向对象编程和函数式编程的各种特性。

    13 引用 • 11 回帖 • 157 关注
  • AWS
    11 引用 • 28 回帖 • 11 关注
  • jsDelivr

    jsDelivr 是一个开源的 CDN 服务,可为 npm 包、GitHub 仓库提供免费、快速并且可靠的全球 CDN 加速服务。

    5 引用 • 31 回帖 • 92 关注
  • WordPress

    WordPress 是一个使用 PHP 语言开发的博客平台,用户可以在支持 PHP 和 MySQL 数据库的服务器上架设自己的博客。也可以把 WordPress 当作一个内容管理系统(CMS)来使用。WordPress 是一个免费的开源项目,在 GNU 通用公共许可证(GPLv2)下授权发布。

    66 引用 • 114 回帖 • 204 关注
  • Bug

    Bug 本意是指臭虫、缺陷、损坏、犯贫、窃听器、小虫等。现在人们把在程序中一些缺陷或问题统称为 bug(漏洞)。

    76 引用 • 1742 回帖 • 1 关注
  • 微软

    微软是一家美国跨国科技公司,也是世界 PC 软件开发的先导,由比尔·盖茨与保罗·艾伦创办于 1975 年,公司总部设立在华盛顿州的雷德蒙德(Redmond,邻近西雅图)。以研发、制造、授权和提供广泛的电脑软件服务业务为主。

    8 引用 • 44 回帖 • 3 关注
  • 小说

    小说是以刻画人物形象为中心,通过完整的故事情节和环境描写来反映社会生活的文学体裁。

    31 引用 • 108 回帖
  • HBase

    HBase 是一个分布式的、面向列的开源数据库,该技术来源于 Fay Chang 所撰写的 Google 论文 “Bigtable:一个结构化数据的分布式存储系统”。就像 Bigtable 利用了 Google 文件系统所提供的分布式数据存储一样,HBase 在 Hadoop 之上提供了类似于 Bigtable 的能力。

    17 引用 • 6 回帖 • 68 关注
  • Q&A

    提问之前请先看《提问的智慧》,好的问题比好的答案更有价值。

    9314 引用 • 42389 回帖 • 115 关注
  • 黑曜石

    黑曜石是一款强大的知识库工具,支持本地 Markdown 文件编辑,支持双向链接和关系图。

    A second brain, for you, forever.

    21 引用 • 204 回帖
  • Word
    13 引用 • 40 回帖
  • abitmean

    有点意思就行了

    29 关注
  • golang

    Go 语言是 Google 推出的一种全新的编程语言,可以在不损失应用程序性能的情况下降低代码的复杂性。谷歌首席软件工程师罗布派克(Rob Pike)说:我们之所以开发 Go,是因为过去 10 多年间软件开发的难度令人沮丧。Go 是谷歌 2009 发布的第二款编程语言。

    498 引用 • 1395 回帖 • 260 关注
  • 深度学习

    深度学习(Deep Learning)是机器学习的分支,是一种试图使用包含复杂结构或由多重非线性变换构成的多个处理层对数据进行高层抽象的算法。

    54 引用 • 40 回帖