2 Windows 10 *1 3 Linux 3.1 Windows Bash on Ubuntu on Windows cygwin MacOS Linux OS Ubuntu OS Linux OS 1 GUI Windows Explorer Mac Finder 1 GUI

Size: px

Start display at page:

Download "2 Windows 10 *1 3 Linux 3.1 Windows Bash on Ubuntu on Windows cygwin MacOS Linux OS Ubuntu OS Linux OS 1 GUI Windows Explorer Mac Finder 1 GUI"

えりかすみだ
5 years ago
Views:

1 September I ll remember- 1,2 Linux 6 Linux Linux 2 3 Linux 2 (OS) Windows MacOS OS MacOS Linux 3 Windows Windows ( ) 1. Bash on Ubuntu on Windows ( Windows 10 ) 2. cygwin ( ) 3. VM Ware Linux OS ( ) 4. Linux OS ( ) Windows

2 2 Windows 10 *1 3 Linux 3.1 Windows Bash on Ubuntu on Windows cygwin MacOS Linux OS Ubuntu OS Linux OS 1 GUI Windows Explorer Mac Finder 1 GUI Graphical User Interface GUI GUI 1 $ 1 pwd ls cd *1 (yuhei.takagi@nao.ac.jp) 3 2

3 1 GUI Mac Windows Linux GUI / 3

4 2 cd./documents/images Documents images./ ls../../bin bin../ ls /bin bin /./../ / 2 mkdir rmdir 3.3 Unix OS Emacs vi Emacs try emacs Emacs $ *2 $ mkdir try emacs $ cd try emacs $ emacs & *2 emacs emacs emacs emacs 4

5 Emacs C-x C-f Emacs C-x M-x Ctrl x Alt x Ctrl x f *3 test.txt C-x C-s C-x C-c Emacs Emacs -nw E macs cat less less q 3.4 ssh ssh -X subaru01@ana03.adc.nao.ac.jp ssh -X ana03 ssh ssh 3.5 exit exit OS *3 5

6 4 4.1 gnuplot SuperMongo matplotlib gnuplot gnuplot *4 gnuplot OS * wget $ wget ls homework.tar.gz $ tar xvzf homework.tar.gz txt 3 *4 matplotlib Python Python SuperMongo(SM) *5 which $ which gnuplot /usr/local/bin/gnuplot Bash on Ubuntu on Windows Ubuntu Debian LinuxOS $ sudo apt-get install gnuplot $ sudo apt-get install gnuplot-x11 gnuplot 6

4.3 M45 M67 NGC 188 3 *6 3 cat less 3.3 Emacs M67 Yadav+08.txt 3 *7 # VizieR (http://vizier.u-strasbg.

7 4.3 M45 M67 NGC *6 3 cat less 3.3 Emacs M67 Yadav+08.txt 3 *7 # VizieR ( # RAJ2000 DEJ2000 B(V, Ic)mag e B(V, Ic) RAJ2000 DEJ2000 *8 B mag V mag Ic mag (magnutude: ) ( ) B V Ic Johnson-Cousins ( 4) 3 Yadav+08.txt emacs *6 *7 emacs M6 Tab *8 7

8 図 4 様々な測光システムのフィルタ透過率 (Bessell 2005, ARA&A, 43, 293) 縦軸が透過率横軸は波長れ異なっているので自分が今どの測光システムを使用しているかを注意する必要がある*9 4.4 図を作ってみる数字の羅列を眺めているだけでは面白くないので簡単な図を作ってみよう図を作るときに使用するのが 4.1 章でインストールした gnuplot だではターミナルに *9 厳密にはたとえ同じ測光システムのフィルターを使っていたとしても個々のフィルターの透過率や CCD の感度特性の違いにより系統誤差が生じる可能性がある異なる観測装置で取得したデータを比較する際はこの系統誤差の補正をする必要がある 8

9 5 gnuplot 1 $ gnuplot * 10 G N U P L O T ( ) Terminal type set to x11 gnuplot> gnuplot gnuplot M67 Yadav+08.txt gnuplot> set datafile separator ",", ( ) gnuplot gnuplot gnuplot> plot "M67 Yadav+08.txt" using 3:4 5 using 3:4 3 4 * 11 *10 gnuplot *11 using u with w line li 9

10 6 gnuplot 2 gnuplot> set xlabel "RA (deg)" gnuplot> set ylabel "DEC (deg)" gnuplot plot gnuplot> replot gnuplot> unset key gnuplot> plot "M67 Yadav+08.txt" using 3:4 pt 5 lc rgb "black" pt lc pointtype linecolor gnuplot pointtype linecolor

1.6 DSS colored 1.4 11.9 1.2 1 DEC (deg) 11.8 0.8 0.6 0.4 B-V 0.2 0 11.7-0.2 132.9 RA (deg) 132.8-0.4 1 Powered by Aladin 17.03 x 18.93 E N 7 gnuplot 3( ) M67 ( ) ( ) 11.6 12.

11 1.6 DSS colored DEC (deg) B-V RA (deg) Powered by Aladin x E N 7 gnuplot 3( ) M67 ( ) ( ) M67 Yadav+08.txt 7 M67 ( 7 ) M67 * 12 gnuplot gnuplot> quit B V B V ( ) B V B V * 13 B V Blue Visible ( 4) B V B *12 gnuplot *13 B V = 0 A0V A0V 11

12 M67 B V awk( ) awk awk 4.3 # * 14 $ grep -v ^# M67 Yadav+08.txt grep # -v # ˆ # # $ grep -v ^# M67 Yadav+08.txt head head 10 * 15 # awk B V B V $ grep -v ^# M67_Yadav+08.txt awk -F"," {print $5" "$7} awk F {print $5 $7} B V awk print 5 7 awk $ grep -v ^# M67_Yadav+08.txt awk -F"," {print $5" "$7" "$5-$7} *14 # gnuplot awk print # *15 head -n n 12

13 $ ( ) awk -F"," {print $5" "$7" "$5-$7} > M67_BvsVvsB-V.txt > >> > gnuplot gnuplot gnuplot $ emacs M67_cmd.plt & gnuplot # set xlabel "{/Times-Italic B-V} (mag)" # set ylabel "{/Times-Italic V} (mag)" # set xrange [-1:3] set yrange [22:6] # # unset key # ps plot "M67_BvsVvsB-V.txt" using 3:2 pt 5 ps 0.5 lc rgb "black" # set terminal postscript eps enhanced color "Times" 10 size 14cm, 10cm set output "M67_cmd.eps" replot 3 gnuplot eps gnuplot 13

14 V (mag) B-V (mag) 8 M67 $ gnuplot M67_cmd.plt * 16 8 $ gv M67_cmd.eps eps eps jpg png imagemagick convert OS jpg $ convert -density 300 M67_cmd.eps M67_cmd.jpg ( ) *16 gv eps eps 14

15 M67 B V M67 B V (HR Hertzsprung-Russel Diagram) M67 Yadav+08.txt e B e V B V 0 CCD CCD ( ) V V B V B V δ B V B V (δ B δ V ) δ B V = (δ B ) 2 + (δ V ) 2 (1) awk 10 grep # ( ) awk -F"," {printf("%5.3f %4.3f %4.3f %4.3f\n", $7,$8,$5-$7,sqrt($6**2+$8**2))} awk print printf printf {printf("( )", ( )} 15

16 V (mag) 14 V (mag) B-V (mag) B-V (mag) 9 ( ) M67 ( ) 0.1 M \n printf M67 V ev B-V eb-v.txt 5.2 plot "M67_V_eV_B-V_eB-V.txt" using 3:1:4:2 with xyerrorbar pt 5 ps 0.5 lc rgb "black" using 4 with xyerrorbar x y using x y x y 9 B V 0.1 awk grep # awk $ grep... awk -F"," {if ($6 <= 0.1 && $8 <= 0.1) print $0} awk if print if && (and) (or) 16

17 V (mag) B-V (mag) 10 M67 ( 2 ) print $0 * HR 7 M67 M67 M67 M67 M67 Yadav+08.txt pmra pmde (proper motion) mas/yr (mas milli arc second ) M67 pmra pmde 0 M67 M67 awk 0±5 mas 10 *17 print $0 printf 17

18 V (mag) B-V (mag) 11 M67 ( ) V ( V ) 10 * 18 M m M (pc) r m M = 5 log r 5 (2) awk awk log log(900)/log(10) 11 *

19 NGC188 Stetson+04.txt NGC188 (B V vs V ) 0.1 (NGC ) 2. M45 Stauffer+07.txt M45 (B V vs V ) (M ) 3. M67 NGC188 M45 A4 2-3 yuhei.takagi@nao.ac.jp yuu.niino@nao.ac.jp 9/15( ) 17:00 awk 3 19

20 A A.1 M67 Yadav+08.txt Yadav, R. K. S., et al. 2008, A&A, 484, 609 MPG/ESO 2.2m telescope (La Silla, Chile) Wide-Field Imager (WFI) B V I A.2 NGC188 Stetson+04.txt Stetson, P. B., et al. 2004, PASP, 116, ( ) Plaskett 1.8 m telescope, Dominion Astrophysical Observatory (Canada) f/5 Newtonian camera B V R B, V Stetson et al A.3 M45 Stauffer+07.txt Stauffer, J. R., et al. 2007, ApJS, 172, B, V Stauffer et al

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

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 2015 8 1 ( ) Unix 1 *1 Unix Unix Unix Perl, Python *2 Unix 2 PC gnuplot *1 100 10 10 6 10 = 10 7 1 1/3 3 10 7 10 7.5 1 24 3600 = (30 6)(30 + 6) 100 = 86400 1 10 7.5 *2 Perl, Python Python 1 Unix * 3 PC