Sinyal dan FFT Pada Python

Seri Belajar Pyhton

1. Download WinPython : https://www.youtube.com/watch?v=GDCdaGmnrVw
2. Install Package Python: https://youtu.be/c_kohGVBH_I
3. Coding Python Membuat Sinyal: https://youtu.be/YXO-Cq3tkAM
4. Membuat Multi Sinyal Dengan Python : https://youtu.be/URUeI6gbdvM
5. Menampilkan Sinyal dalam FFT Python : https://youtu.be/cqG_zlwwR68

-------------------------
Link untuk mendownload WinPython
https://winpython.github.io/
WinPython 3.7: https://sourceforge.net/projects/winpython/files/WinPython_3.7/3.7.7.1/
-------------------------


----------- Script 1 -------
import numpy as np
from scipy import signal as sg
import matplotlib.pyplot as plt

----------- Script 3 -------
import numpy as np
from scipy import signal as sg
import matplotlib.pyplot as plt

freq = 20
amp = 1
time = np.linspace(0,2,1000)

plt.plot(time,signal1)
plt.xlabel('Time (s)')
plt.ylabel('Amplitudo (V)')

----------- Script 4 -------
import numpy as np
from scipy import signal as sg
import matplotlib.pyplot as plt
from numpy.fft import fft, ifft

#sampling rate
sr = 2000

#sampling interval
ts = 1.0/sr
t = np.arange(0,1,ts)

# gelombang pertama
freq1 = 1
amp1 = 3
signal1 = amp1*np.sin(2*np.pi*freq1*t)

# gelombang kedua
freq2 = 4
amp2 = 0.5
signal2 = amp2*np.sin(2*np.pi*freq2*t)

# gelombang ketiga
freq3 = 7
amp3 = 1
signal3 = amp3*np.sin(2*np.pi*freq3*t)

# gelombang campuran/jumlah
signalAll = signal1 + signal2 + signal3

plt.figure(figsize = (32, 6))
plt.subplot(1,3,1)
plt.plot(t, signal1)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')

plt.subplot(1,3,2)
plt.plot(t, signal2)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')

plt.subplot(1,3,3)
plt.plot(t, signal3)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')

---- script 5 ----------------
import numpy as np
from scipy import signal as sg
import matplotlib.pyplot as plt
from numpy.fft import fft, ifft

#sampling rate
sr = 2000

#sampling interval
ts = 1.0/sr
t = np.arange(0,1,ts)

# gelombang pertama
freq1 = 1
amp1 = 3
signal1 = amp1*np.sin(2*np.pi*freq1*t)

# gelombang kedua
freq2 = 4
amp2 = 0.5
signal2 = amp2*np.sin(2*np.pi*freq2*t)

# gelombang ketiga
freq3 = 7
amp3 = 1
signal3 = amp3*np.sin(2*np.pi*freq3*t)

# gelombang campuran/jumlah
signalAll = signal1 + signal2 + signal3

plt.figure(figsize = (32, 6))
plt.subplot(1,3,1)
plt.plot(t, signal1)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')

plt.subplot(1,3,2)
plt.plot(t, signal2)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')

plt.subplot(1,3,3)
plt.plot(t, signal3)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')

plt.figure(figsize = (12, 6))
plt.subplot(1,2,1)
plt.plot(t, signalAll)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')

X = fft(signalAll)
N = len(X)
n = np.arange(N)
T = N/sr
freqFFT = n/T

plt.subplot(1,2,2)
plt.stem(freqFFT, np.abs(X),'b', markerfmt=" ", basefmt="-b")
plt.xlabel('Freq (Hz)')
plt.ylabel('FFT Amplitude |X(freq)|')
plt.xlim(0,10)