×
Ir ao conteúdo
  • Cadastre-se

Python Período de uma onda senoidal


Posts recomendados

import cv2
import numpy as np
from skimage.io import imsave
import matplotlib.pyplot as plt
from statsmodels.graphics.tsaplots import plot_pacf
import statsmodels as sm
from PIL import Image 
from PIL.Image import NEAREST
from PIL.Image import BICUBIC
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import seaborn as sns
from sys import argv
from math import pi, sin
from scipy import fftpack
from PIL import Image 

#fs é a frequência de amostragem
fs = 20.0

t = np.linspace(0,17,17); 
#t= np.array([0,1,2,3,4,5,6,78,9,10,11,12,13,14,15,16,17])#np.linspace(0,1 , int(1*fs), endpoint=False) #números com intervalos especificados.
#precisa especificar onde começar e onde parar.

#s is the sum of sine wave(1Hz) and cosine wave(10 Hz)
s = np.array([0, 0.707, 1, 0.707, 0, -0.707, -1, -0.707, 0, 0.707, 1, 0.707, 0, -0.707, -1, -0.707, 0]) #np.asarray #1.0*np.sin(8 * 2.0 * np.pi * t) #np.sin(np.pi*t)+ np.sin(np.pi*t*2)
plt.plot(t, s)
plt.xlim(0,17) #(0,10)
plt.xlabel("tempo")
plt.title('Original Signal in Time Domain')
plt.show()

# Compute the one-dimensional discrete Fourier Transform.
fft_wave = np.fft.fft(s)

# Compute the Discrete Fourier Transform sample frequencies.
fft_fre = np.fft.fftfreq(n=s.size, d=1/fs) 

plt.subplot(111)#(211) #linhas, colunas e o número do gráfico.
plt.plot(fft_fre, fft_wave.real, label="Real part")
plt.xlim(-5,5)
plt.ylim(-10,10)
#plt.legend(loc=1) 
plt.title("FFT in Frequency Domain")

plt.subplot(212) 
plt.plot(fft_fre, fft_wave.imag,label="Imaginary part")
#plt.legend(loc=1)
plt.xlim(-5,5)
plt.ylim(-10,10)
plt.xlabel("frequency (Hz)")

plt.ylabel("Amplitude")
plt.xlabel("Time [s]")
plt.bar(fft_fre, np.abs(fft_wave), width=0.5) #np.abs é uma matriz de entrada
plt.show()

frequencias = fft_fre #Retorne as frequências de amostra da Transformada Discreta de Fourier. 
period_time = 1 / frequencias
amplitudes = np.abs(fft_wave) #(0, 0.707, 1, 0.707, 0, -0.707, -1, -0.707) np.abs= calcula o valor absoluto por elemento

fpos=frequencias[frequencias>0] #mascara para só pegar valores positivos
apos=amplitudes[frequencias>0]

plt.ylabel("Amplitude")
plt.xlabel("Frequency [Hz]")
plt.bar(fpos,apos, width=0.1) #plt.scatter p gráfico ser de pontos
plt.show()

idx=np.argmax(apos) #Retorna os índices dos valores máximos ao longo de um eixo.
print("Máxima amplitude corresponde a=",apos[idx])
tpCount     = len(apos)
timePeriod  = tpCount/fs

#peiodo = ii_arr[apos_max[0]:apos_max[1]] #priodo é a distancia entre duas amplitudes maximas
#print("O período corresponde a=", tpCount/fs)

Quero encontrar o período desta onda seno.. alguém consegue me ajudar a organizar o código para que isso seja possível?

 

Link para o comentário
Compartilhar em outros sites

Crie uma conta ou entre para comentar

Você precisa ser um usuário para fazer um comentário

Criar uma conta

Crie uma nova conta em nossa comunidade. É fácil!

Crie uma nova conta

Entrar

Já tem uma conta? Faça o login.

Entrar agora

Sobre o Clube do Hardware

No ar desde 1996, o Clube do Hardware é uma das maiores, mais antigas e mais respeitadas comunidades sobre tecnologia do Brasil. Leia mais

Direitos autorais

Não permitimos a cópia ou reprodução do conteúdo do nosso site, fórum, newsletters e redes sociais, mesmo citando-se a fonte. Leia mais

×
×
  • Criar novo...

Aprenda a ler resistores e capacitores

EBOOK GRÁTIS!

CLIQUE AQUI E BAIXE AGORA MESMO!