MATLAB实践篇：初探与应用实例（一）

下面给出两道基于火星探测器MAVEN实测数据的基础练习题：

(数据下载地址：https://pds-ppi.igpp.ucla.edu/；打开后点击Mars-MAVEN)

%Task 1-MAG数据

1. 读取'mvn_mag_l2_2014358ss_20141224_v01_r01.sts'数据，并重新写入新文件，命名为：'mvn_mag_l2_2014358ss_20141224_v01_r01.txt'（要求：读取中将相关文件名称中含有日期的'2014358'和'20141224'作为变量赋值处理，可能用到'num2str'及'strcat'函数）

2. load重写的'.txt'文件，并识别整个数据中的每一列意义，分别给出Time(时间)；MAG(磁场三分量及总磁场)；卫星位置信息SC_pos(位置三分量及距离火星中心的距离)。注意：对应数据为Time 1列；MAG 4列；SC_pos 4列

3. 求出当天5:00-5:20am的平均磁场强度及其三分量并绘图

4. 绘制卫星在CYL（柱坐标系）下的运动轨迹

A.横坐标为x，纵坐标为根号下y^2+z^2，单位均为一个火星半径，用黑线表示，线宽为1 

B.火星中心位于(0,0)处，近似球体，单位为1，用红线表示，线宽为1

% Task 2-SWIA数据

1. 必应搜索'CDF (V3.8.0 version) Patch for MATLAB'

2. 下载与自己电脑配置符合的相关文件并将其解压保存在常用路径里

3. 读取'mvn_swi_l2_onboardsvymom_20150126_v01_r00.cdf'文件并查看其各cell（数组元胞）的名称（读取时别忘记添加路径，可详细阅读下载网站的example）

4. 给出其时间分辨率，即每个观测间隔为几秒

5. 挑选出其中的离子速度对应的元胞（可能用到'cell2mat'函数）

6. 画出当日6:00-7:50的速度三分量（x,y,z）及速度合成量（tot）的变化

  （在一张图里用'蓝、绿、红和黑'表示x,y,z,tot，线宽为1）

经过这两题，你应该可以较为快捷地熟悉Matlab语言环境和变量设置，文件数据读取和基本画图地一些操作，这都是在Matlab-基础篇中提及过的内容，希望可以尽自己最快的速度掌握，这才是刚刚打开Matlab大门的钥匙。

此后，请再来看看以下的逻辑判断相关内容，这一模块是为了帮助你提高程序处理中判断结构体的认知，以及自己在未来编程习惯中需要改变的一些习惯。

图一. 逻辑结构试读训练

当C和D出现分别出现以下四种情景时，对应的(1)-(7)到底时真还是假呢？是不是感觉很像高中数学的命题练习，只是穿插了T和F的运算。

• C=1；D=1;
• C=0；D=0；
• C=1；D=0；
• C=0；D=1。

如何检验？用Matlab命令行自行检验自己的计算结果哦。