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

本贴最后更新于 739 天前,其中的信息可能已经东海扬尘
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 会得到 0。而相应地,函数调用就可以在序列上进行,相应的用的是 /@。比如输入 Sin @ {Pi/6, Pi/4, Pi/3, Pi/2, Pi} 会得到 \{\frac12, \frac{1}{\sqrt{2}}, \frac{\sqrt{3}}{2}, 1, 0\}

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

相关帖子

欢迎来到这里!

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

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

推荐标签 标签

  • Wide

    Wide 是一款基于 Web 的 Go 语言 IDE。通过浏览器就可以进行 Go 开发,并有代码自动完成、查看表达式、编译反馈、Lint、实时结果输出等功能。

    欢迎访问我们运维的实例: https://wide.b3log.org

    30 引用 • 218 回帖 • 625 关注
  • CentOS

    CentOS(Community Enterprise Operating System)是 Linux 发行版之一,它是来自于 Red Hat Enterprise Linux 依照开放源代码规定释出的源代码所编译而成。由于出自同样的源代码,因此有些要求高度稳定的服务器以 CentOS 替代商业版的 Red Hat Enterprise Linux 使用。两者的不同在于 CentOS 并不包含封闭源代码软件。

    238 引用 • 224 回帖
  • ReactiveX

    ReactiveX 是一个专注于异步编程与控制可观察数据(或者事件)流的 API。它组合了观察者模式,迭代器模式和函数式编程的优秀思想。

    1 引用 • 2 回帖 • 154 关注
  • etcd

    etcd 是一个分布式、高可用的 key-value 数据存储,专门用于在分布式系统中保存关键数据。

    5 引用 • 26 回帖 • 526 关注
  • 爬虫

    网络爬虫(Spider、Crawler),是一种按照一定的规则,自动地抓取万维网信息的程序。

    106 引用 • 275 回帖
  • 正则表达式

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

    31 引用 • 94 回帖 • 1 关注
  • SVN

    SVN 是 Subversion 的简称,是一个开放源代码的版本控制系统,相较于 RCS、CVS,它采用了分支管理系统,它的设计目标就是取代 CVS。

    29 引用 • 98 回帖 • 683 关注
  • Jenkins

    Jenkins 是一套开源的持续集成工具。它提供了非常丰富的插件,让构建、部署、自动化集成项目变得简单易用。

    53 引用 • 37 回帖
  • Spark

    Spark 是 UC Berkeley AMP lab 所开源的类 Hadoop MapReduce 的通用并行框架。Spark 拥有 Hadoop MapReduce 所具有的优点;但不同于 MapReduce 的是 Job 中间输出结果可以保存在内存中,从而不再需要读写 HDFS,因此 Spark 能更好地适用于数据挖掘与机器学习等需要迭代的 MapReduce 的算法。

    74 引用 • 46 回帖 • 561 关注
  • API

    应用程序编程接口(Application Programming Interface)是一些预先定义的函数,目的是提供应用程序与开发人员基于某软件或硬件得以访问一组例程的能力,而又无需访问源码,或理解内部工作机制的细节。

    77 引用 • 429 回帖 • 2 关注
  • 新人

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

    52 引用 • 228 回帖
  • 宕机

    宕机,多指一些网站、游戏、网络应用等服务器一种区别于正常运行的状态,也叫“Down 机”、“当机”或“死机”。宕机状态不仅仅是指服务器“挂掉了”、“死机了”状态,也包括服务器假死、停用、关闭等一些原因而导致出现的不能够正常运行的状态。

    13 引用 • 82 回帖 • 53 关注
  • 安装

    你若安好,便是晴天。

    132 引用 • 1184 回帖
  • iOS

    iOS 是由苹果公司开发的移动操作系统,最早于 2007 年 1 月 9 日的 Macworld 大会上公布这个系统,最初是设计给 iPhone 使用的,后来陆续套用到 iPod touch、iPad 以及 Apple TV 等产品上。iOS 与苹果的 Mac OS X 操作系统一样,属于类 Unix 的商业操作系统。

    84 引用 • 139 回帖
  • Redis

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

    286 引用 • 248 回帖 • 74 关注
  • Sublime

    Sublime Text 是一款可以用来写代码、写文章的文本编辑器。支持代码高亮、自动完成,还支持通过插件进行扩展。

    10 引用 • 5 回帖 • 1 关注
  • PWA

    PWA(Progressive Web App)是 Google 在 2015 年提出、2016 年 6 月开始推广的项目。它结合了一系列现代 Web 技术,在网页应用中实现和原生应用相近的用户体验。

    14 引用 • 69 回帖 • 150 关注
  • 创业

    你比 99% 的人都优秀么?

    84 引用 • 1399 回帖 • 1 关注
  • ActiveMQ

    ActiveMQ 是 Apache 旗下的一款开源消息总线系统,它完整实现了 JMS 规范,是一个企业级的消息中间件。

    19 引用 • 13 回帖 • 668 关注
  • ZeroNet

    ZeroNet 是一个基于比特币加密技术和 BT 网络技术的去中心化的、开放开源的网络和交流系统。

    1 引用 • 21 回帖 • 637 关注
  • JWT

    JWT(JSON Web Token)是一种用于双方之间传递信息的简洁的、安全的表述性声明规范。JWT 作为一个开放的标准(RFC 7519),定义了一种简洁的,自包含的方法用于通信双方之间以 JSON 的形式安全的传递信息。

    20 引用 • 15 回帖 • 6 关注
  • Vim

    Vim 是类 UNIX 系统文本编辑器 Vi 的加强版本,加入了更多特性来帮助编辑源代码。Vim 的部分增强功能包括文件比较(vimdiff)、语法高亮、全面的帮助系统、本地脚本(Vimscript)和便于选择的可视化模式。

    29 引用 • 66 回帖 • 5 关注
  • 资讯

    资讯是用户因为及时地获得它并利用它而能够在相对短的时间内给自己带来价值的信息,资讯有时效性和地域性。

    55 引用 • 85 回帖
  • 房星科技

    房星网,我们不和没有钱的程序员谈理想,我们要让程序员又有理想又有钱。我们有雄厚的房地产行业线下资源,遍布昆明全城的 100 家门店、四千地产经纪人是我们坚实的后盾。

    6 引用 • 141 回帖 • 584 关注
  • Telegram

    Telegram 是一个非盈利性、基于云端的即时消息服务。它提供了支持各大操作系统平台的开源的客户端,也提供了很多强大的 APIs 给开发者创建自己的客户端和机器人。

    5 引用 • 35 回帖 • 1 关注
  • jQuery

    jQuery 是一套跨浏览器的 JavaScript 库,强化 HTML 与 JavaScript 之间的操作。由 John Resig 在 2006 年 1 月的 BarCamp NYC 上释出第一个版本。全球约有 28% 的网站使用 jQuery,是非常受欢迎的 JavaScript 库。

    63 引用 • 134 回帖 • 726 关注
  • MongoDB

    MongoDB(来自于英文单词“Humongous”,中文含义为“庞大”)是一个基于分布式文件存储的数据库,由 C++ 语言编写。旨在为应用提供可扩展的高性能数据存储解决方案。MongoDB 是一个介于关系数据库和非关系数据库之间的产品,是非关系数据库当中功能最丰富,最像关系数据库的。它支持的数据结构非常松散,是类似 JSON 的 BSON 格式,因此可以存储比较复杂的数据类型。

    90 引用 • 59 回帖 • 4 关注