99999久久久久久亚洲,欧美人与禽猛交狂配,高清日韩av在线影院,一个人在线高清免费观看,啦啦啦在线视频免费观看www

熱線(xiàn)電話(huà)：13121318867

登錄
首頁(yè)精彩閱讀分析引力波數(shù)據(jù)的Python包：GWPY！
分析引力波數(shù)據(jù)的Python包：GWPY！
2016-02-14
收藏

美國(guó)科學(xué)家11日宣布，他們?nèi)ツ?月首次探測(cè)到引力波。這一發(fā)現(xiàn)印證了物理學(xué)大師愛(ài)因斯坦100年前的預(yù)言。宣布這一發(fā)現(xiàn)的，是激光干涉引力波天文臺(tái)（LIGO）的負(fù)責(zé)人。


這個(gè)機(jī)構(gòu)誕生于上世紀(jì)90年代，進(jìn)行引力波觀測(cè)已經(jīng)有近30年。那么觀測(cè)到的引力波數(shù)據(jù)的量應(yīng)該很大，科學(xué)家如何對(duì)這些數(shù)據(jù)進(jìn)行分析？有沒(méi)有用到Python編程語(yǔ)言？

答案是肯定的。筆者在Github上發(fā)現(xiàn)了一個(gè)專(zhuān)門(mén)用于分析引力波數(shù)據(jù)的Python包：GWPY。據(jù)維護(hù)者介紹，GWPY的代碼來(lái)自L(fǎng)IGO和另一個(gè)名叫Virgo的機(jī)構(gòu)，維護(hù)者將這兩個(gè)機(jī)構(gòu)科學(xué)家的Python代碼整理，最終的產(chǎn)品就是GWPY這個(gè)用戶(hù)友好的Python包。

在具體介紹GWPY之前，先給和筆者一樣的小白簡(jiǎn)單科普一下引力波和LIGO的相關(guān)知識(shí)。

什么是引力波？

上圖是兩個(gè)黑洞所產(chǎn)生的引力波的3-D模擬圖（NASA）。

首先，什么是引力波？在物理學(xué)上，引力波是愛(ài)因斯坦廣義相對(duì)論所預(yù)言的一種以光速傳播的時(shí)空波動(dòng)，如同石頭丟進(jìn)水里產(chǎn)生的波紋一樣，引力波被視為宇宙中的“時(shí)空漣漪”。通常引力波的產(chǎn)生非常困難，地球圍繞太陽(yáng)以每秒30千米的速度前進(jìn)，發(fā)出的引力波功率僅為200瓦，還不如家用電飯煲功率大。宇宙中大質(zhì)量天體的加速、碰撞和合并等事件才可以形成強(qiáng)大的引力波，但能產(chǎn)生這種較強(qiáng)引力波的波源距離地球都十分遙遠(yuǎn)，傳播到地球時(shí)變得非常微弱。

下面分享兩個(gè)優(yōu)秀的視頻，很好地解釋了引力波及背后的原理。第一個(gè)來(lái)自L(fǎng)IGO，第二個(gè)則是比較通俗的漫畫(huà)式講解。

LIGO科學(xué)家的解釋?zhuān)?/span>http://v.qq.com/boke/page/g/0/0/g0184mxwie0.html

漫畫(huà)式通俗解釋?zhuān)?/span>http://v.qq.com/page/j/x/u/j0184qlilxu.html

LIGO是什么？

激光干涉引力波觀測(cè)站（ Laser Interferometer Gravitational-Wave Observatory）LIGO是加州理工學(xué)院（Caltech）和麻省理工學(xué)院（MIT）的合作實(shí)驗(yàn)室，現(xiàn)在也有其他的大學(xué)參與。實(shí)驗(yàn)資金來(lái)源于美國(guó)國(guó)家科學(xué)基金會(huì)。LIGO是用來(lái)尋找宇宙中的引力波，從而可以驗(yàn)證黑洞的存在和檢驗(yàn)廣義相對(duì)論。
LIGO主要有兩個(gè)觀測(cè)點(diǎn)，位于路易斯安那Livingston Parish的LIGO Livingston觀測(cè)點(diǎn)，和華盛頓 Hanford的LIGO Hanford觀測(cè)點(diǎn)。除此之外，在加州Passadena 的Caltech校園中還有LIGO 40m Prototype 。

GWPY：LIGO用它分析引力波數(shù)據(jù)？
接下來(lái)是本文的重頭戲。我們一起來(lái)學(xué)習(xí)如何GWPY分析引力波數(shù)據(jù)。下面的介紹及示例均來(lái)自GWPY的官方文檔。

安裝

很簡(jiǎn)單，pip install gwpy就可以完成安裝。

不過(guò)安裝的過(guò)程可能會(huì)比較長(zhǎng)，因?yàn)間wpy使用的依賴(lài)包比較多，包括numpy、 scipy、 cycler、matplotlib、astropy等。

面向?qū)ο缶幊?/strong>

GWPY是一個(gè)面向?qū)ο缶幊痰腜ython包，也就是說(shuō)，數(shù)據(jù)對(duì)象是這個(gè)包的核心關(guān)注點(diǎn)。每一個(gè)數(shù)據(jù)對(duì)象都體現(xiàn)為一個(gè)類(lèi)實(shí)例，包含了其屬性和包含的數(shù)據(jù)。

如果想創(chuàng)建一個(gè)新的類(lèi)實(shí)例，建議使用標(biāo)準(zhǔn)的構(gòu)建器（constructor）。舉個(gè)例子，我們可以使用一個(gè)數(shù)據(jù)數(shù)組，生成一個(gè)TimeSeries對(duì)象：

>>> from gwpy.timeseries import TimeSeries>>> mydata = TimeSeries([1,2,3,4,5,6,7,8,9,10], sample_rate=1, epoch=0)
或者從在線(xiàn)數(shù)據(jù)服務(wù)器上下載：

>>> from gwpy.timeseries import TimeSeries>>> mydata = TimeSeries.fetch('H1:LDAS-STRAIN', 964656015, 964656615)
核心數(shù)據(jù)對(duì)象

據(jù)介紹，GWPY提供了4種核心數(shù)據(jù)對(duì)象，分別代表引力波探測(cè)器所產(chǎn)生的四種標(biāo)準(zhǔn)數(shù)據(jù)：

TimeSeries（時(shí)間序列數(shù)據(jù)）
Spectrum（光譜數(shù)據(jù)）
Spectrogram（光譜圖）
DataQualityFlag
引力波數(shù)據(jù)可視化

我們知道，將引力波探測(cè)器收集的數(shù)據(jù)可視化，對(duì)于理解引力波的特性、研究引力波信號(hào)來(lái)說(shuō)非常有幫助。gwpy.plotter模塊中提供了一些plot類(lèi)，可以直觀地展示相應(yīng)的數(shù)據(jù)類(lèi)型。

GWPY的核心數(shù)據(jù)對(duì)象里，大部分都內(nèi)置有一個(gè)plot()方法，可以讓研究人員快速對(duì)某個(gè)數(shù)據(jù)集進(jìn)行可視化展示。舉個(gè)例子：

>>> from gwpy.timeseries import TimeSeries
>>> data = TimeSeries.fetch('H1:LDAS-STRAIN', 968654552, 968654562)
>>> plot = data.plot()
>>> plot.show()
GWPY：利用公開(kāi)的LIGO數(shù)據(jù)進(jìn)行繪圖

我們接下來(lái)利用LIGO公開(kāi)的一些引力波時(shí)間序列數(shù)據(jù)進(jìn)行繪圖。我們可以直接在線(xiàn)加載這些數(shù)據(jù)。首先導(dǎo)入我們需要的模塊：

>>> from urllib2 import urlopen>>> from numpy import asarray>>> from gwpy.timeseries import TimeSeries


然后，下載數(shù)據(jù)，保存為文本字符串：

>>> data = urlopen('http://www.ligo.org/science/GW100916/L-strain_hp30-968654552-10.txt').read()


現(xiàn)在，我們可以對(duì)文本進(jìn)行解析，補(bǔ)充必要的元數(shù)據(jù)之后，就可以生成一個(gè)TimeSeries：

>>> ts = TimeSeries(asarray(data.splitlines(), dtype=float),>>>                 epoch=968654552, sample_rate=16384, unit='strain')


最后，我們就可以繪圖了：

>>> plot = ts.plot()>>> plot.set_title('LIGO Livingston Observatory data for GW100916')>>> plot.set_ylabel('Gravitational-wave strain amplitude')>>> plot.show()


來(lái)源/編程派


數(shù)據(jù)分析咨詢(xún)請(qǐng)掃描二維碼

若不方便掃碼，搜微信號(hào)：CDAshujufenxi

數(shù)據(jù)分析師資訊
更多

OK
客服在線(xiàn)
立即咨詢(xún)
客服在線(xiàn)
立即咨詢(xún)
') } function initGt() { var handler = function (captchaObj) { captchaObj.appendTo('#captcha'); captchaObj.onReady(function () { $("#wait").hide(); }).onSuccess(function(){ $('.getcheckcode').removeClass('dis'); $('.getcheckcode').trigger('click'); }); window.captchaObj = captchaObj; }; $('#captcha').show(); $.ajax({ url: "/login/gtstart?t=" + (new Date()).getTime(), // 加隨機(jī)數(shù)防止緩存 type: "get", dataType: "json", success: function (data) { $('#text').hide(); $('#wait').show(); // 調(diào)用 initGeetest 進(jìn)行初始化 // 參數(shù)1：配置參數(shù) // 參數(shù)2：回調(diào)，回調(diào)的第一個(gè)參數(shù)驗(yàn)證碼對(duì)象，之后可以使用它調(diào)用相應(yīng)的接口 initGeetest({ // 以下 4 個(gè)配置參數(shù)為必須，不能缺少 gt: data.gt, challenge: data.challenge, offline: !data.success, // 表示用戶(hù)后臺(tái)檢測(cè)極驗(yàn)服務(wù)器是否宕機(jī) new_captcha: data.new_captcha, // 用于宕機(jī)時(shí)表示是新驗(yàn)證碼的宕機(jī) product: "float", // 產(chǎn)品形式，包括：float，popup width: "280px", https: true // 更多配置參數(shù)說(shuō)明請(qǐng)參見(jiàn)：http://docs.geetest.com/install/client/web-front/ }, handler); } }); } function codeCutdown() { if(_wait == 0){ //倒計(jì)時(shí)完成 $(".getcheckcode").removeClass('dis').html("重新獲取"); }else{ $(".getcheckcode").addClass('dis').html("重新獲取("+_wait+"s)"); _wait--; setTimeout(function () { codeCutdown(); },1000); } } function inputValidate(ele,telInput) { var oInput = ele; var inputVal = oInput.val(); var oType = ele.attr('data-type'); var oEtag = $('#etag').val(); var oErr = oInput.closest('.form_box').next('.err_txt'); var empTxt = '請(qǐng)輸入'+oInput.attr('placeholder')+'！'; var errTxt = '請(qǐng)輸入正確的'+oInput.attr('placeholder')+'！'; var pattern; if(inputVal==""){ if(!telInput){ errFun(oErr,empTxt); } return false; }else { switch (oType){ case 'login_mobile': pattern = /^1[3456789]\d{9}$/; if(inputVal.length==11) { $.ajax({ url: '/login/checkmobile', type: "post", dataType: "json", data: { mobile: inputVal, etag: oEtag, page_ur: window.location.href, page_referer: document.referrer }, success: function (data) { } }); } break; case 'login_yzm': pattern = /^\d{6}$/; break; } if(oType=='login_mobile'){ } if(!!validateFun(pattern,inputVal)){ errFun(oErr,'') if(telInput){ $('.getcheckcode').removeClass('dis'); } }else { if(!telInput) { errFun(oErr, errTxt); }else { $('.getcheckcode').addClass('dis'); } return false; } } return true; } function errFun(obj,msg) { obj.html(msg); if(msg==''){ $('.login_submit').removeClass('dis'); }else { $('.login_submit').addClass('dis'); } } function validateFun(pat,val) { return pat.test(val); }