一种根据实验室仪器导出的双列 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

相关帖子

欢迎来到这里!

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

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

推荐标签 标签

  • 心情

    心是产生任何想法的源泉,心本体会陷入到对自己本体不能理解的状态中,因为心能产生任何想法,不能分出对错,不能分出自己。

    59 引用 • 369 回帖
  • GraphQL

    GraphQL 是一个用于 API 的查询语言,是一个使用基于类型系统来执行查询的服务端运行时(类型系统由你的数据定义)。GraphQL 并没有和任何特定数据库或者存储引擎绑定,而是依靠你现有的代码和数据支撑。

    4 引用 • 3 回帖
  • ZooKeeper

    ZooKeeper 是一个分布式的,开放源码的分布式应用程序协调服务,是 Google 的 Chubby 一个开源的实现,是 Hadoop 和 HBase 的重要组件。它是一个为分布式应用提供一致性服务的软件,提供的功能包括:配置维护、域名服务、分布式同步、组服务等。

    59 引用 • 29 回帖 • 1 关注
  • OpenCV
    15 引用 • 36 回帖 • 2 关注
  • 代码片段

    代码片段分为 CSS 与 JS 两种代码,添加在 [设置 - 外观 - 代码片段] 中,这些代码会在思源笔记加载时自动执行,用于改善笔记的样式或功能。

    用户在该标签下分享代码片段时需在帖子标题前添加 [css] [js] 用于区分代码片段类型。

    129 引用 • 862 回帖
  • RYMCU

    RYMCU 致力于打造一个即严谨又活泼、专业又不失有趣,为数百万人服务的开源嵌入式知识学习交流平台。

    4 引用 • 6 回帖 • 53 关注
  • 分享

    有什么新发现就分享给大家吧!

    247 引用 • 1794 回帖 • 1 关注
  • Ant-Design

    Ant Design 是服务于企业级产品的设计体系,基于确定和自然的设计价值观上的模块化解决方案,让设计者和开发者专注于更好的用户体验。

    17 引用 • 23 回帖
  • WiFiDog

    WiFiDog 是一套开源的无线热点认证管理工具,主要功能包括:位置相关的内容递送;用户认证和授权;集中式网络监控。

    1 引用 • 7 回帖 • 603 关注
  • 链滴

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

    记录生活,连接点滴

    168 引用 • 3834 回帖
  • 大数据

    大数据(big data)是指无法在一定时间范围内用常规软件工具进行捕捉、管理和处理的数据集合,是需要新处理模式才能具有更强的决策力、洞察发现力和流程优化能力的海量、高增长率和多样化的信息资产。

    93 引用 • 113 回帖
  • 钉钉

    钉钉,专为中国企业打造的免费沟通协同多端平台, 阿里巴巴出品。

    15 引用 • 67 回帖 • 292 关注
  • 锤子科技

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

    4 引用 • 31 回帖 • 7 关注
  • GitLab

    GitLab 是利用 Ruby 一个开源的版本管理系统,实现一个自托管的 Git 项目仓库,可通过 Web 界面操作公开或私有项目。

    46 引用 • 72 回帖 • 1 关注
  • CSS

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

    198 引用 • 541 回帖 • 1 关注
  • 微软

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

    8 引用 • 44 回帖
  • 房星科技

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

    6 引用 • 141 回帖 • 589 关注
  • Vditor

    Vditor 是一款浏览器端的 Markdown 编辑器,支持所见即所得、即时渲染(类似 Typora)和分屏预览模式。它使用 TypeScript 实现,支持原生 JavaScript、Vue、React 和 Angular。

    364 引用 • 1840 回帖 • 2 关注
  • MyBatis

    MyBatis 本是 Apache 软件基金会 的一个开源项目 iBatis,2010 年这个项目由 Apache 软件基金会迁移到了 google code,并且改名为 MyBatis ,2013 年 11 月再次迁移到了 GitHub。

    171 引用 • 414 回帖 • 375 关注
  • 知乎

    知乎是网络问答社区,连接各行各业的用户。用户分享着彼此的知识、经验和见解,为中文互联网源源不断地提供多种多样的信息。

    10 引用 • 66 回帖
  • Facebook

    Facebook 是一个联系朋友的社交工具。大家可以通过它和朋友、同事、同学以及周围的人保持互动交流,分享无限上传的图片,发布链接和视频,更可以增进对朋友的了解。

    4 引用 • 15 回帖 • 439 关注
  • Webswing

    Webswing 是一个能将任何 Swing 应用通过纯 HTML5 运行在浏览器中的 Web 服务器,详细介绍请看 将 Java Swing 应用变成 Web 应用

    1 引用 • 15 回帖 • 635 关注
  • jsDelivr

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

    5 引用 • 31 回帖 • 92 关注
  • 服务

    提供一个服务绝不仅仅是简单的把硬件和软件累加在一起,它包括了服务的可靠性、服务的标准化、以及对服务的监控、维护、技术支持等。

    41 引用 • 24 回帖
  • SQLServer

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

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

    FreeMarker 是一款好用且功能强大的 Java 模版引擎。

    23 引用 • 20 回帖 • 462 关注
  • Outlook
    1 引用 • 5 回帖 • 2 关注