Unix * 3 PC 2 Linux, Mac 4 Windows Cygwin Cygwin gnuplot Cygwin unix emulator online gnuplot 5 matplotlib 6 SuperMongo 7 gnuplot gnuplot OS *8 Uni

Size: px

Start display at page:

Download "Unix * 3 PC 2 Linux, Mac *4 Windows Cygwin Cygwin gnuplot Cygwin unix emulator online gnuplot *5 matplotlib *6 SuperMongo *7 gnuplot gnuplot OS *8 Uni"

ひでかむらかわ
5 years ago
Views:

1 ( ) Unix 1 *1 Unix Unix Unix Perl, Python *2 Unix 2 PC gnuplot * = / = (30 6)(30 + 6) 100 = *2 Perl, Python Python 1

2 Unix * 3 PC 2 Linux, Mac *4 Windows Cygwin Cygwin gnuplot Cygwin unix emulator online gnuplot *5 matplotlib *6 SuperMongo *7 gnuplot gnuplot OS *8 Unix Unix Unix 3 Unix emacs * pwd *3 Unix *4 *5 / *6 1 Python Python *7 2 *8 *9 emacs, vi, gedit 2

3 * 10 $ pwd /home/ichikawa /home/ichikawa /home/ichikawa * 11 $ wget ls * 12 $ ls Autumn2015.tar.gz ls $ ls Autumn2015.tar.gz Autumn2015.tar.gz Autumn2015.tar.gz $ ls A $ ls Autumn2015.tar.gz Unix.tar.gz tar $ tar zxvf Autumn2015.tar.gz * 13 ls $ ls Autumn2015/ Autumn2015.tar.gz *10 print working directory *11 *12 list *13 tar zxvf A 3

4 wget wget which $ which wget /usr/bin/wget * 14 wget wget URL which $ which wget $ wget * 15 wget URL mv YourDownloadDirectory. mv $ mv YourDownloadDirectory/Autumn2015.tar.gz. $ tar zxvf Autumn2015.tar.gz. tar Autumn2015 cp 3.2 Swift/BAT 70 ( BAT ) X * 16 (E = kev) 70 * 17 BAT * 18 1 BAT 70-month Survey Products Catalog: ASCII ( 300 KB) * 19 ( ) (SWIFT *14 /usr/bin/wget wget *15 wget *16 X hard X-ray X hard ( ) X soft ( ) X soft/hard E = 2 kev soft hard *17 detect NGC 1068 Hα *18 *19 4

5 1 Swift/BAT 70 JXXXX.X-XXXX ) (RA) (DEC) PC Mac Cmd+s * 20 wget * 21 $wget BAT 70m catalog 20nov2012.txt ls $ls Autumn2015/ Autumn2015.tar.gz BAT_70m_catalog_20nov2012.dat BAT mkdir * 22 $ mkdir BAT70mon $ ls Autumn2015/ Autumn2015.tar.gz BAT70mon/ BAT_70m_catalog_20nov2012.dat *20 *21 wget Autumn2015/ cp *22 Make directory 5

6 mv * 23 $ mv BAT_70m_catalog_20nov2012.dat BAT70mon/ $ ls Autumn2015/ Autumn2015.tar.gz BAT70mon/ ls BAT BAT BAT70mon BAT70mon cd * 24 $ cd BAT70mon/ $ pwd /home/ichikawa/bat70mon $ ls BAT_70m_catalog_20nov2012.dat BAT70mon BAT 3.3 BAT less $ less BAT_70m_catalog_20nov2012.dat BAT_NAME RA DEC SNR COUNTERPART_NAME deg deg 1 SWIFT J MASX J SWIFT J Fairall SWIFT J NGC SWIFT J MASX J deg 6 6 awk $ which awk /usr/bin/awk PC awk awk 6 *23 move *24 change directory 6

7 $ awk -F" " {print $6} BAT_70m_catalog_20nov2012.txt... 2MASX J PKS IGR J MASX J H PKS $ 6 awk {print $6} awk -F BAT 70m catalog 20nov2012.txt 6 ( ) 3 (R.A. * 25 ) 4 (Dec * 26 ) 8 (ctpt RA) 9 (ctpt DEC) * awk $ awk -F" " {print $8,$9} BAT_70m_catalog_20nov2012.txt $6 $8, $9 8 9 > ra vs dec.txt $ awk -F" " {print $8,$9} BAT_70m_catalog_20nov2012.txt > ra_vs_dec.txt $ ls BAT_70m_catalog_20nov2012.dat $ less ra_vs_dec.txt ra_vs_dec.txt CTPT_RA CTPT_DEC *25 R.A. Right Ascension *26 Declination *27 7

8 deg deg ra vs dec.txt > ( ) (ra vs dec.txt) $ awk NR>3{print} ra_vs_dec.txt > ra_vs_dec_nohead.txt NR awk NR 1 NR 2 m NR m NR>3{print} NR 3 print awk print ra vs dec.txt 4 * 28 less head head 10 $ head ra_vs_dec_nohead.txt n n *28 {print} 8

9 $ head -n ra_vs_dec_nohead.txt 3.4 gnuplot gnuplot $ gnuplot G N U P L O T Version 4.6 patchlevel 1 last modified Build System: Darwin x86_64 Copyright (C) , 1998, 2004, Thomas Williams, Colin Kelley and many others gnuplot home: faq, bugs, etc: type "help FAQ" immediate help: type "help" (plot window: hit h ) Terminal type set to x11 gnuplot> gnuplot gnuplot> plot sin(x) sin(x) gnuplot y = sin(x) 9

10 gnuplot gnuplot gnuplot> plot "ra_vs_dec_nohead.txt" u 1:2 3 ra vs dec nohead.txt 1 2 u 1:2 u using gnuplot gnuplot set se gnuplot "ra_vs_dec_nohead.txt" BAT R.A. Dec 3 (180, 60) (300, 40) ra vs dec.txt 0 x y 90 gnuplot> se xrange[0:360] gnuplot> se yrange[-90:90] gnuplot> plot "ra_vs_dec_nohead.txt" u 1:2 notitle 4 (.png ).eps.pdf.eps gnuplot> se term eps 10

11 BAT R.A. Dec Terminal type set to epscairo Options are transparent fontscale 0.5 size 5.00in, 3.00in gnuplot> se output "radec.eps" gnuplot> replot gnuplot> exit eps (se term eps) (se output radec.eps ) (replot) gnuplot (exit) 4 radec.eps $ ls BAT_70m_catalog_20nov2012.dat ra_vs_dec.txt ra_vs_dec_nohead.txt radec.eps eps Mac Linux open $ open radec.eps open Linux 3.5 BAT 18 z * X * 30 log 2 *29, redshift z λ ini λ obs λ obs = (1 + z)λ ini z z z = z cosmo calc *30 : erg/s 11

12 awk $ awk -F" " {print $18,$19} BAT_70m_catalog_20nov2012.dat NULL z NULL z Unix wc $ wc -l BAT_70m_catalog_20nov2012.dat 1213 BAT_70m_catalog_20nov2012.dat $ awk -F" " $19~/NULL/{print $18,$19} BAT_70m_catalog_20nov2012.dat wc -l NULL * 31 wc -l * z z 432 $ awk -F" " NR>3 && $19!~/NULL/{print $18,$19} BAT_70m_catalog_20nov2012.dat > z_vs_lx.txt awk 4 (NR> 3) 19 NULL ($19! /NULL/) z X log L X z vs lx.txt z sort $ sort z_vs_lx.txt *31 *

13 z z = X X 2 sort 2 sort -k $ sort -k2 z_vs_lx.txt X log L X /erg s 1 = erg s X sort head $ sort -k2 z_vs_lx.txt head sort -r ( ) $ sort -k2 -r z_vs_lx.txt X log L X /erg s 1 = z vs lx.txt gnuplot $ gnuplot gnuplot> plot "z_vs_lx.txt" u 1:2 5 z < 0.5 z = 4.0 log x log gnuplot 13

14 49 48 "z_vs_lx.txt" BAT z (x ) X (y ) gnuplot> se logscale x gnuplot> plot "z_vs_lx.txt" u 1: "z_vs_lx.txt" BAT z (x ) X (y ) 6 z = (z 1) ( z 0.01) (log L X = 42 44) BAT flux (z log L X ) BAT (z log L X ) z, X flux flux (luminosity) flux 14

15 BAT id 21 R.A., Dec plot (Active Galactic Nuclei; AGN) id =4 8 $ awk -F" " $21>=4 && $21<=8 {print $8,$9} BAT_70m_catalog_20nov2012.dat > ra_vs_dec_agn.txt gnuplot plot ( ) "ra_vs_dec_agn.txt" BAT AGN R.A. Dec 7 AGN AGN AGN gnuplot awk X (X-ray binary) (cataclysmic variable; CV) id = 1 9 id 15 * ((80, 70) (20, 75) ) 3.7 Unix awk Unix awk *33 awk or ( ) 15

16 BAT R.A. Dec Python 16

17 4 1. Unix Unix 2. BAT BAT 70m catalog 20nov2012.dat Unix wc tail * 34 BAT 70m catalog nohead.dat wc 3. 3 # BAT 70m catalog cmthead.dat 4. BAT z > 1 BAT 70m catalog zgt1.dat 5. BAT 70m catalog zgt1.dat 2 (BAT source type) 6. BAT 70m catalog zgt1.dat 2 * 35 BAT 70m catalog zgt1 ty7.dat BAT 70m catalog zgt1 ty8.dat 7. BAT 70m catalog zgt1 ty7.dat BAT 70m catalog zgt1 ty8.dat ( R. A., Dec) gnuplot *34 tail: *35 : awk 17

18 Unix awk google Mac Unix Unix for Mac OS X * 36 awk AWK * 37 awk * 38 BAT NASA Extragalactic Database (NED) * 39 z ( Blazar ) SIMBAD Astronomical Database * 40 Simbad * 41 NED Simbad NED Simbad ( ) AGN NED AGN AGN (Ichikawa, Kohei) ( ) (AGN) *36 *37 *38 *39 *40 *41 AladinLite 18

19 $ wc BAT_70m_catalog_20nov2012.dat BAT_70m_catalog_20nov2012.dat $ tail BAT_70m_catalog_20nov2012.dat > BAT_70m_catalog_nohead.txt $ wc BAT_70m_catalog_nohead.txt BAT_70m_catalog_nohead.txt 3. $ awk NR<=3{print "#"$0}NR>3 BAT_70m_catalog_20nov2012.dat > BAT_70m_catalog_cmthead.txt 4. $ awk -F" " $18>1.0 BAT_70m_catalog_nohead.txt > BAT_70m_catalog_zgt1.txt $ wc -l BAT_70m_catalog_zgt1.txt 30 BAT_70m_catalog_zgt1.txt $ awk -F" " {print $21} BAT_70m_catalog_zgt1.txt sort uniq 7 8 source type 7: Blazar/BL Lac 8: QSO $ awk -F" " $21==7{print $21} BAT_70m_catalog_zgt1.txt wc -l 11 $ awk -F" " $21==8{print $21} BAT_70m_catalog_zgt1.txt wc -l $ awk -F" " $21==7 {print} BAT_70m_catalog_zgt1.txt > BAT_70m_catalog_zgt1_ty7.txt $ awk -F" " $21==8 {print} BAT_70m_catalog_zgt1.txt > BAT_70m_catalog_zgt1_ty8.txt 7. $ gnuplot gnuplot> se datafile sepa " " gnuplot> plot "BAT_70m_catalog_zgt1_ty7.txt" u 8:9, "BAT_70m_catalog_zgt1_ty8.txt" u 8:9 gnuplot> se term eps Terminal type set to epscairo Options are transparent fontscale 0.5 size 5.00in, 3.00in gnuplot> se output "radec_zgt1.eps" gnuplot> replot gnuplot> exit 19

20 "BAT_70m_catalog_zgt1_ty7.txt" u 8:9 "BAT_70m_catalog_zgt1_ty8.txt" u 8: BAT R.A. Dec 20

1 1 Gnuplot gnuplot Windows gnuplot gp443win32.zip gnuplot binary, contrib, demo, docs, license 5 BUGS, Chang

1 1 Gnuplot gnuplot Windows gnuplot gp443win32.zip gnuplot binary, contrib, demo, docs, license 5 BUGS, Chang Gnuplot で微分積分 2011 年度前期数学解析 I 講義資料 (2011.6.24) 矢崎成俊 ( 宮崎大学 ) 1 1 Gnuplot gnuplot http://www.gnuplot.info/ Windows gnuplot 2011 6 22 4.4.3 gp443win32.zip gnuplot binary, contrib, demo, docs, license 5