H=zeros(mmm,'f') tpi=2*pi alpha=linspace(0,tpi,mmm) ae=sqrt(8./3.) H=ae*0.5*(1.-cos(alpha)) return H ######################################################################## print " The input file must have two columns: time(sec) & amplitude " a,b,num =read_two_columns() sr,dt,mean,sd,rms,skew,kurtosis,dur=signal_stats(a,b,num) sr,dt=sample_rate_check(a,b,num,sr,dt) ######################################################################## print " " print " Remove mean: 1=yes 2=no " mr_choice = GetInteger2() print " " print " Select Window: 1=Rectangular 2=Hanning " h_choice = GetInteger2() ########################################################################
from scipy.stats.stats import pearsonr import numpy as np ######################################################################## print ("Each file must have two columns: time(sec) & accel(G)") a1,b1,num1 =read_two_columns() sr1,dt1,mean1,sd1,rms1,skew1,kurtosis1,dur1=signal_stats(a1,b1,num1) sr1,dt1=sample_rate_check(a1,b1,num1,sr1,dt1) ######################################################################## a2,b2,num2 =read_two_columns() sr2,dt2,mean2,sd2,rms2,skew2,kurtosis2,dur2=signal_stats(a2,b2,num2) sr2,dt2=sample_rate_check(a2,b2,num2,sr2,dt2) ######################################################################## a1=array(a1) b1=array(b1) #print ("a1 time: ", a1)