MATLAB中用散点图分析岩石中的化学成份
1、数据获取
一般的环境数据，如气温、水位数据等都是时间序列，即按时间顺序组织的数据，许多数据却是另外一种类型，例如岩石样本数据。可以从一些网站上获得这些数据，图中数据表示SiO2、TiO2、Al2O3、FeO、MgO、CaO、Na2O等七种元素在采样到的岩石中的含量。可到相关网站，根据自己的需求，选择数据。
2、 格式转换
网络提供的数据，一般是.asv格式，或者是excel格式，可以用相关工具，转换成.txt格式。
3、加载数据及统计维数
D=load('rock_data.txt');%加载数据文件
Ns=size(D);
N=Ns(1);%采样点的数目
4、循环统计
这个嵌套的for循环图绘制了各种成份i和j的所有组合。我们可以将范围限制在j>i，因为j=i的情况对应于绘制一个成份与自己的对比，并且j<i图是多余的。注意，外部for循环变量i的范围从1到M-1，内部for循环变量j的范围在i+1到M的间隔内。连续绘图之间的暂停是通过ginput()命令实现的；单击该图形表明是显示下一个图形。
5、sprintf()函数
它创建一个字符串，其中包含文本和变量的值。函数使用以字符%开头的占位符来指示变量值应该放在字符串的什么位置。因此，
i=2;
sprintf(‘element %d’,i);
返回字符串 ‘element 2’. %d 整数的占位符，被‘2’代替.
也可以在同一格式字符串中使用几个占位符，例如
i=2;
j=4;
sprintf(‘row %d column %d’, i, j);
返回字符串 ‘row 2 column 4’.
6、运行结果及分析
左图表示成份1（如SiO2）和成份2（TiO2）的关系，我们看到成份1分布在40%-50%，而成份2分布在1%以下。右图可以看到，随着成份3的增加，成份5呈减少的趋势。
