2016年8月7日

高鐵上的氣壓變化-新竹至台中

繼上回的《高鐵上的氣壓變化-新竹至台北
今天搭高鐵前往台中,我又做了一次高鐵上的觀測氣壓

原來這整段有31段隧道,顯見竹苗一帶的丘陵地形,這些隧道多到我沒辦法在圖上標示了。

編號11的新城隧道一直到編號41的神岡隧道

1.觀察氣壓變化,我認為新苗隧道、頭份一號隧道和頭份二號隧道至少都有一個豎井。氣壓在隧道內都會有一次大的變化。

2.途中氣壓變化最大的不在隧道段,是在兩個隧道之間(苑裡隧道-神岡隧道)

3.上回推測某些隧道表現出來的氣壓變化可能是跟明挖覆蓋有關,但是觀察神岡隧道(最右邊那個)卻沒看到這個現象。


高鐵氣壓偵測-新竹台中

今天回程的時候也進行了一次觀測,但是時間是六點多,無法藉由外界天光與隧道的光線差異去偵測隧道所在,於是想了應該也可以透過偵測聲音改變。


新竹到苗栗段的聲音是怎麼變化的呢?偵測平均音量,結果是這樣,可以觀測到隧道內的確有比較吵。我還想作的是去作FFT,把特定一個頻率段抓出來,其實在隧道裡會聽到一個聲音,幾乎是一聽就知道進隧道了。中間的低谷(8:23左右)是苗栗站,過了苗栗站,整個高鐵就吵起來了,不過不是乘客的吵,是列車的整體聲音。


sound


回程(台中-新竹)的高鐵倒是一致性的都很吵,不過這班沒有停苗栗
sound2


作圖的時候,為了要把隧道段自動highlight,特別去找了資料來參考。上回因為隧道少,自己手動用inkscape畫一畫,這次多達31個,根本沒辦法手動。

參考
Plotting the recession - some basic data analysis with gnuplot

基本就是設定這個函數
g(x) = (x<1000?GPVAL_Y_MAX:GPVAL_Y_MIN)

然後先跑一次圖
plot "light" using 2:4
獲得GPVAL_Y_MAX和GPVAL_Y_MIN的數值

然後用下面這指令畫圖的時候,就會自動把小於1000的亮度值Highlight,用個灰色方塊畫出來(其實是個被填滿的曲線)
plot  "light" using 2:(g($4)) with filledcurve lc rgb "#eeeeee" notitle,\

-----------------------------------
gnuplot的命令稿
-----------------------------------
set xdata time
set datafile separator ","
set timefmt "%H:%M:%S"
set format x "%H:%M"
set xrange ["08:12:00":"08:42:00"]
set xlabel "Time"
set ylabel "light"
set y2label "pressure"
set ytics nomirror
set y2tics
g(x) = (x<1000?GPVAL_Y_MAX:GPVAL_Y_MIN)
plot "light" using 2:4
plot  "light" using 2:(g($4)) with filledcurve lc rgb "#eeeeee" notitle,\
      "light" using 2:4 with lines lc rgbcolor "#f4c2c2" title "light" axis x1y1,\
      "pressure" using 2:4 with lines title "pressure" lc rgbcolor "blue"  axis x1y2,\
pause -1