《GMT绘时间序列及相分析的方法》由会员分享，可在线阅读，更多相关《GMT绘时间序列及相分析的方法（3页珍藏版）》请在金锄头文库上搜索。
1、GMT 绘时间序列及相分析的方法具体步骤如下：1 首先是准备好规则数据。如果数据格式是 matlab 处理和存储的矩阵格式，使用 MATLAB 则需要提取出来时间序列，并保存为ASCII格式数据。注意此时的数据不带有时间信息，只 是简单的一列数据而已，我们将在在后面添加时间信息。2为sigmaO数据添加这个方法不是固定的，可以自己编程实现，也可以使用已有的简单程序。注意时间数据的格 式是固定的几种。如果数据不多的话，可以使用Windows自带的excel添加。本例子的最终数据格式如下：第一列是时间，第二列到第四列是sigmaO的属性值，（time， sigmaOKu， sigmaOC， sig
2、maOC-Ku）。3 GMT 绘图将三列 sigmaO 的数据放到一幅图上，即同时绘制三条时间序列，且使用同一个坐标系。横 坐标为年月，具体设置为：年为大刻度，月为小刻度，界限设置为1993.1.1-2004.12.1。竖 坐标为sigma0的值，范围可当已知，设置为0-15。本例子参考GMT Cookbook中的程序21：GMT中Time格式需要按照数据的时间格式设置。在GMT的帮助文件中可以找到，有格里 高利，ISO之类。还要区分输入的时间格式和输出的时间格式。我们选择例21中使用的时间标准dd-o-yy。（GMT有细致的时间设置参数，注意在这一步别 出错误。）注意：数据文件中时间序列的日
3、期坐标可以和标注的日期坐标格式不一样，但是时间上要一 致。rem 主要代码，rem 内容参考文后注释。gmtsetINPUT_DATE_FORMATdd-o-yy PLOT_DATE_FORMAT o ANNOT_FONT_SIZE_PRIMARY +10pgmtsetTIME_FORMAT_PRIMARY abbreviated CHAR_ENCODING ISOLatin1+ 坐标轴设置，前面的是时间起点和终点，后面的0/300是竖坐标的界限。Set R=-R1999-08-11T00:00:00/2006-11-01T00:00:00/0/300gmtconvert为数据提取程序，可以提
4、取多列数据文件中的某一列或者几列数据。Gmtconvert -F0,2 -f0T -Hi RHAT_price.csv RHAT.envpsxy为绘制曲线的程序，可以设置点和线的样式。psxy -R -J -Gyellow -O -K RHAT.env example_21.ps4最终的GMT bat文件代码如下：gmtset INPUT_DATE_FORMAT dd-o-yy PLOT_DATE_FORMAT o ANNOT_FONT_SIZE_PRIMARY +10pgmtset TIME_FORMAT_PRIMARY abbreviated CHAR_ENCODING ISOLatin1
5、+REM Create a suitable -R stringset R=-R1993-01-01T00:00:00/2004-12-31T00:00:00/0/15REM Lay down the basemap:psbasemap %R% -JX9i/4i -Gwhite -K -Bs1Y:t/a:/:,db:WSen -Bpf1o/2WSen:.Time series of sigma0 over Xinjiang China: example_21.psgmtset OUTPUT_DATE_FORMAT dd-o-yygmtconvert -F0,2 -f0T XJsigma0_TS
6、.txt RHAT.envgmtconvert -F0,3 -f0T XJsigma0_TS.txt RHAT2.envgmtconvert -F0,1 -f0T XJsigma0_TS.txt RHAT1.envrem red C blue Ku black C-Kupsxy -R -J -O -K -Wthin -Sc0.1c -Gblack RHAT2.env example_21.pspsxy -R -J -O -K -Wthinner RHAT2.env example_21.pspsxy -R -J -O -K -Wthin,blue RHAT.env example_21.psp
7、sxy -R -J -O -K -St0.2c -Gblue RHAT.env example_21.pspsxy -R -J RHAT1.env -Wthin,red -O -K example_21.pspsxy -R -J -O -K -Si0.2c -Gred RHAT1.env example_21.ps rem gmtset INPUT_DATE_FORMAT dd-o-yyREM Clean up after ourselves:del RHAT.*del .gmt*5 时间序列的结果如下：至此，时间序列完成了。PS：网格虚线，也可以自动化实现。修改default中的grid p
8、rimary size,在粗细后面加上“，-”然后在psbasemap中的-B属性中加g控制网格密度。下一步是使用GMT分析其频率域以及趋势分析。（建议自己尝试）2 添加趋势。如下图所示。趋势不必在所有的时间序列中添加，目前只在明显的添加。但是要给出趋势线的主要参数。 并且在时间序列图中显示。使用MATLAB可以做到，但是不能与第一步的图结合。因此先全部用GMT作图。trend 1d3 FFT 绘图。116个数据FFT的MATLAB程序还没调好。这个也尝试用GMT绘制。spectrum 1dgrdfft?Matlab的字符串数组不知道怎么设置。用 excel 却可以很好的完成时间的添加。法为：
9、先写出第一个时间起点，然后设置单元格样 式，设置为 然后托动鼠标下拉即可。 然后再修改时间格式为 dd-mmm-yy.使用matlab编程难以实现，却使用excel实现。输入格式输出格式 In addition, you may use a singleoinstead of mm (to plot month name)You may specify either Gregorian calendar format or ISO week calendar format.As examples, try yy/mm/dd, yyyy=jjj, dd-o-yyyy, dd-mm-yy,yy-mm, etc.Controls how primary month-, week-, and weekday-names are formatted. Choose amongfull, abbreviated, and character.If the leading f, a, or c are replaced with F, A, and C the entire annotation will be in upper case.设置标注使用的时间界限。后面的/0/300 则是 Y 轴的界限。 提取第 0， 2 列， -f0T 表示时间
