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

葡萄

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

gawk " BEGIN{for(i=0;i<4*3.14;i=i+0.4)print i,cos(i)}" > test1.dat
set ofn=1.ps
gmtset FRAME_PEN 1.6p ANOT_FONT_SIZE 19
psxy -R0/12.56/-1.5/1.5 -Ba2f0.5:"x":/a1f0.2g4:"y":WSne -W2/125/125/01ta -JX25/10 test.dat -K -Y5 > %ofn%
rem 画sin(x)曲线
gawk "{print $1,$3}" test.dat | psxy -R -W/255/0/0 -JX -K  -O >> %ofn%
rem 画cos(x)曲线
gawk "{print $1,$4}" test.dat | psxy -R -W/0/255/0 -JX -K  -O >> %ofn%
rem 画sin(x)+cos(x)曲线
psxy -R -JX -K -Sa0.2c -W3/255/0/0 -G255/255/0 -O test1.dat >> %ofn%
rem 以0.4为步长画cos(x)曲线上的数据点（点类型为五角星），线粗3，颜色为红色，内部充填黄色
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*