python實現(xiàn)信號時域統(tǒng)計特征提取代碼-創(chuàng)新互聯(lián)

1.實驗數(shù)據(jù)需求

創(chuàng)新互聯(lián)是一家集網(wǎng)站建設,武陟企業(yè)網(wǎng)站建設,武陟品牌網(wǎng)站建設,網(wǎng)站定制,武陟網(wǎng)站建設報價,網(wǎng)絡營銷,網(wǎng)絡優(yōu)化,武陟網(wǎng)站推廣為一體的創(chuàng)新建站企業(yè)，幫助傳統(tǒng)企業(yè)提升企業(yè)形象加強企業(yè)競爭力。可充分滿足這一群體相比中小企業(yè)更為豐富、高端、多元的互聯(lián)網(wǎng)需求。同時我們時刻保持專業(yè)、時尚、前沿，時刻以成就客戶成長自我，堅持不斷學習、思考、沉淀、凈化自己，讓我們?yōu)楦嗟钠髽I(yè)打造出實用型網(wǎng)站。

為了對采集的壓力實驗數(shù)據(jù)做特征工程，需要對信號進行時域的統(tǒng)計特征提取，包含了均值、均方根、偏度、峭度、波形因子、波峰因子、脈沖因子、峭度因子等，現(xiàn)用python對其進行實現(xiàn)。

2.python實現(xiàn)

其中的輸入?yún)?shù)含義：

① data:實驗數(shù)據(jù)的DataFrame

② p1:所截取實驗信號的起始采樣點位置

③ p2:所截取實驗信號的終止采樣點位置

from pandas import Series
import math
pstf_list=[]
def psfeatureTime(data,p1,p2):
 #均值
 df_mean=data[p1:p2].mean()
 #方差
 df_var=data[p1:p2].var()
 #標準差
 df_std=data[p1:p2].std()
 #均方根
 df_rms=math.sqrt(pow(df_mean,2) + pow(df_std,2))
 #偏度
 df_skew=data[p1:p2].skew()
 #峭度
 df_kurt=data[p1:p2].kurt()
 sum=0
 for i in range(p1,p2):
  sum+=math.sqrt(abs(data[i]))
 #波形因子
 df_boxing=df_rms / (abs(data[p1:p2]).mean())
 #峰值因子
 df_fengzhi=(max(data[p1:p2])) / df_rms
 #脈沖因子
 df_maichong=(max(data[p1:p2])) / (abs(data[p1:p2]).mean())
 #裕度因子
 df_yudu=(max(data[p1:p2])) / pow((sum/(p2-p1)),2)
 featuretime_list = [df_mean,df_rms,df_skew,df_kurt,df_boxing,df_fengzhi,df_maichong,df_yudu]
 return featuretime_list 

本文標題：python實現(xiàn)信號時域統(tǒng)計特征提取代碼-創(chuàng)新互聯(lián)
標題來源：http://muchs.cn/article28/dcjejp.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián)，為您提供網(wǎng)站營銷、企業(yè)建站、外貿網(wǎng)站建設、網(wǎng)站收錄、用戶體驗、關鍵詞優(yōu)化

廣告

聲明：本網(wǎng)站發(fā)布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經(jīng)允許不得轉載，或轉載時需注明來源： 創(chuàng)新互聯(lián)