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

本贴最后更新于 755 天前,其中的信息可能已经东海扬尘
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

相关帖子

欢迎来到这里!

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

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

推荐标签 标签

  • 链滴

    链滴是一个记录生活的地方。

    记录生活,连接点滴

    153 引用 • 3783 回帖 • 1 关注
  • SendCloud

    SendCloud 由搜狐武汉研发中心孵化的项目,是致力于为开发者提供高质量的触发邮件服务的云端邮件发送平台,为开发者提供便利的 API 接口来调用服务,让邮件准确迅速到达用户收件箱并获得强大的追踪数据。

    2 引用 • 8 回帖 • 483 关注
  • CSS

    CSS(Cascading Style Sheet)“层叠样式表”是用于控制网页样式并允许将样式信息与网页内容分离的一种标记性语言。

    198 引用 • 550 回帖
  • Oracle

    Oracle(甲骨文)公司,全称甲骨文股份有限公司(甲骨文软件系统有限公司),是全球最大的企业级软件公司,总部位于美国加利福尼亚州的红木滩。1989 年正式进入中国市场。2013 年,甲骨文已超越 IBM,成为继 Microsoft 后全球第二大软件公司。

    105 引用 • 127 回帖 • 382 关注
  • 微服务

    微服务架构是一种架构模式,它提倡将单一应用划分成一组小的服务。服务之间互相协调,互相配合,为用户提供最终价值。每个服务运行在独立的进程中。服务于服务之间才用轻量级的通信机制互相沟通。每个服务都围绕着具体业务构建,能够被独立的部署。

    96 引用 • 155 回帖 • 1 关注
  • JRebel

    JRebel 是一款 Java 虚拟机插件,它使得 Java 程序员能在不进行重部署的情况下,即时看到代码的改变对一个应用程序带来的影响。

    26 引用 • 78 回帖 • 664 关注
  • 工具

    子曰:“工欲善其事,必先利其器。”

    286 引用 • 729 回帖
  • 正则表达式

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

    31 引用 • 94 回帖
  • LeetCode

    LeetCode(力扣)是一个全球极客挚爱的高质量技术成长平台,想要学习和提升专业能力从这里开始,充足技术干货等你来啃,轻松拿下 Dream Offer!

    209 引用 • 72 回帖
  • Bug

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

    75 引用 • 1737 回帖 • 3 关注
  • Scala

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

    13 引用 • 11 回帖 • 130 关注
  • 锤子科技

    锤子科技(Smartisan)成立于 2012 年 5 月,是一家制造移动互联网终端设备的公司,公司的使命是用完美主义的工匠精神,打造用户体验一流的数码消费类产品(智能手机为主),改善人们的生活质量。

    4 引用 • 31 回帖 • 4 关注
  • OpenStack

    OpenStack 是一个云操作系统,通过数据中心可控制大型的计算、存储、网络等资源池。所有的管理通过前端界面管理员就可以完成,同样也可以通过 Web 接口让最终用户部署资源。

    10 引用 • 4 关注
  • Latke

    Latke 是一款以 JSON 为主的 Java Web 框架。

    71 引用 • 535 回帖 • 787 关注
  • WebComponents

    Web Components 是 W3C 定义的标准,它给了前端开发者扩展浏览器标签的能力,可以方便地定制可复用组件,更好的进行模块化开发,解放了前端开发者的生产力。

    1 引用
  • HBase

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

    17 引用 • 6 回帖 • 73 关注
  • Android

    Android 是一种以 Linux 为基础的开放源码操作系统,主要使用于便携设备。2005 年由 Google 收购注资,并拉拢多家制造商组成开放手机联盟开发改良,逐渐扩展到到平板电脑及其他领域上。

    334 引用 • 323 回帖
  • CSDN

    CSDN (Chinese Software Developer Network) 创立于 1999 年,是中国的 IT 社区和服务平台,为中国的软件开发者和 IT 从业者提供知识传播、职业发展、软件开发等全生命周期服务,满足他们在职业发展中学习及共享知识和信息、建立职业发展社交圈、通过软件开发实现技术商业化等刚性需求。

    14 引用 • 155 回帖
  • 周末

    星期六到星期天晚,实行五天工作制后,指每周的最后两天。再过几年可能就是三天了。

    14 引用 • 297 回帖 • 1 关注
  • jsoup

    jsoup 是一款 Java 的 HTML 解析器,可直接解析某个 URL 地址、HTML 文本内容。它提供了一套非常省力的 API,可通过 DOM,CSS 以及类似于 jQuery 的操作方法来取出和操作数据。

    6 引用 • 1 回帖 • 477 关注
  • frp

    frp 是一个可用于内网穿透的高性能的反向代理应用,支持 TCP、UDP、 HTTP 和 HTTPS 协议。

    20 引用 • 7 回帖
  • JVM

    JVM(Java Virtual Machine)Java 虚拟机是一个微型操作系统,有自己的硬件构架体系,还有相应的指令系统。能够识别 Java 独特的 .class 文件(字节码),能够将这些文件中的信息读取出来,使得 Java 程序只需要生成 Java 虚拟机上的字节码后就能在不同操作系统平台上进行运行。

    180 引用 • 120 回帖
  • 导航

    各种网址链接、内容导航。

    40 引用 • 173 回帖
  • 微软

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

    8 引用 • 44 回帖 • 1 关注
  • Flutter

    Flutter 是谷歌的移动 UI 框架,可以快速在 iOS 和 Android 上构建高质量的原生用户界面。 Flutter 可以与现有的代码一起工作,它正在被越来越多的开发者和组织使用,并且 Flutter 是完全免费、开源的。

    39 引用 • 92 回帖
  • Python

    Python 是一种面向对象、直译式电脑编程语言,具有近二十年的发展历史,成熟且稳定。它包含了一组完善而且容易理解的标准库,能够轻松完成很多常见的任务。它的语法简捷和清晰,尽量使用无异义的英语单词,与其它大多数程序设计语言使用大括号不一样,它使用缩进来定义语句块。

    543 引用 • 672 回帖
  • 运维

    互联网运维工作,以服务为中心,以稳定、安全、高效为三个基本点,确保公司的互联网业务能够 7×24 小时为用户提供高质量的服务。

    149 引用 • 257 回帖