GMT画时间序列（三角函数）

vfory

C:\GMT_Samples\ex05
gmtset GRID_CROSS_SIZE 0 ANOT_FONT_SIZE 15 FRAME_WIDTH 0.2c
gmtset GRID_PEN 0.25pta
gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF
gmtset FRAME_WIDTH 0.1c
gmtset LABEL_FONT_SIZE 12 ANOT_FONT 1 LABEL_FONT 1
rem gmtset FRAME_PEN 1.6p ANOT_FONT_SIZE 19

set ofn=sinAndCosMap.ps

gawk " BEGIN{for(i=0;i<4*3.14;i=i+0.01)print i,sin(i),cos(i),sin(i)+cos(i)}" > threemap.dat
gawk " BEGIN{for(i=0;i<4*3.14;i=i+0.4)print i,cos(i)}" > cosx.dat

psxy -R0/12.56/-1.5/1.5 -Ba2f0.5:"x":/a1f0.2g4:"y":WSne -W2/125/125/01ta -JX25/10 threemap.dat -K -Y5 > %ofn%

rem 画sin(x)曲线
gawk "{print [Math Processing Error] 2}" threemap.dat | psxy -R -W/255/0/0 -JX -K -O >> %ofn%

rem 画cos(x)曲线
gawk "{print [Math Processing Error] 3}" threemap.dat | psxy -R -W/0/255/0 -JX -K -O >> %ofn%

rem 画sin(x)+cos(x)曲线
gawk "{print [Math Processing Error] 4}" threemap.dat | psxy -R -W/0/255/0 -JX -K -O >> %ofn%

rem 以0.4为步长画cos(x)曲线上的数据点（点类型为五角星），线粗3，颜色为红色，内部充填黄色
psxy -R -JX -K -Sa0.2c -W3/255/0/0 -G255/255/0 -O cosx.dat >> %ofn%

rem 绘制标签
echo 4.5 -1.2 27 0 7 LM --- y=sin(x)+cos(x)| pstext -G0/255/0 -R -JX -K -O >> %ofn%
echo 9 0.4 27 30 7 LM --- y=sin(x) | pstext -G125/125/0 -R -JX -K -O >> %ofn%
echo 1.5 0.3 27 30 7 LM --- y=cos(x) | pstext -G255/0/0 -R -JX -O >> %ofn%

del .gmt*
pause