#conv_fast = fft.convolution(data,gaus) #t1 = time.time() #conv_slow = fft.convolution_slow(data_slow,gaus_slow) #t2 = time.time() # #print("Fast: ", t2-t1, "Slow: ", t1-t0) # ##plt.show() # ## assigment 4 problem 3 N=20 N2=fft.pow2(N) G=np.zeros(N2) H=np.zeros(N2) for i in range(0,round(.1*len(G)),1): G[i]=1 print(G) for i in range(round(.4*len(H)),round(.6*len(H)),1): H[i]=1 print(H) corr = fft.correlation(G,H) print(corr)
t2[j] = j*dt2 f2[j] = j if t2[j] <= 0.1: g2[j] = 1 h2[j] = 0 elif t2[j] <= 0.4: h2[j] = 0 g2[j] = 0 elif t2[j] <= 0.6: g2[j] = 0 h2[j] = 1 else: h2[j] = 0 g2[j] = 0 print("with N equal to " + str(N1)) print(fft.correlation(g1,h1)) print("with N equal to " + str(N2)) print(fft.correlation(g2,h2)) #fft.plot_c(f1[:0.5*N],[:N]) #ys = np.zeros(2*N) #yf = np.zeros(2*N) #sigma = PI / 16. / 2 #dsigma2 = 1. / (sigma*sigma) #C = 1. / (np.sqrt(2.*PI)*sigma) #I = 0 #for i in range(N): # x[i] = i*dx # f[i] = i