GMT画时间序列（三角函数）
本帖最后由 vfory 于 2013-3-21 18:13 编辑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 2}" threemap.dat | psxy -R -W/255/0/0 -JX -K -O >> %ofn%
rem 画cos(x)曲线
gawk "{print 3}" threemap.dat | psxy -R -W/0/255/0 -JX -K -O >> %ofn%
rem 画sin(x)+cos(x)曲线
gawk "{print 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
页:
[1]