Simulasi Gelombang Dengan Google Colabs

Untuk memulai menggunakan Google Colabs dengan melakukan langkah-langkah berikut:

1. Buka google search engine, kemudian ketik colabs
2. Pilih Colaboratory, atau ketik alamat berikut: https://colab.research.google.com/
3. Akan muncul window berikut, kemudian pilih New notebook

Kemudian akan muncul window seperti berikut ini:

Membuat Gelombang

Ketik coding berikut untuk membuat sebuah gelombang:

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

frequency = 2000
amplitudo = 3
time = np.linspace(0, 2, 1000)

gelombang = amplitudo*np.sin(2*np.pi*frequency*time)

plt.plot(time, gelombang)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude (V)')
plt.title('Gelombang Sinus')
plt.show()

Membuat Osilasi Teredam (Damped Oscilation)

from scipy.integrate import solve_ivp

t = np.linspace(0,15,1000)

y = [0,1]

gamma = 1

omega_sqr = 100

def sho(t,y):

    solution = (y[1],(-gamma*y[1]-omega_sqr*y[0]))

    return solution

solution = solve_ivp(sho, [0,1000], y0 = y, t_eval = t)

plt.plot(t,solution.y[0])

plt.ylabel(“Position”)

plt.xlabel(“Time”)

plt.title(‘Damped Oscillator’, fontsize = 20)

Membuat Dua Gelombang dan Menggabungkannya

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

frequency1 = 2000
frequency2 = 3000
amplitudo1 = 3

time = np.linspace(0, 2, 1000)

gelombang_1 = amplitudo1*np.sin(2*np.pi*frequency1*time)
gelombang_2 = amplitudo1*np.sin(2*np.pi*frequency2*time)

Sinyal_Gabungan = gelombang_1 + gelombang_2

plt.plot(time, gelombang_1,'r')
plt.xlabel('Time (s)')
plt.ylabel('Amplitude (V)')
plt.title('Gelombang Gabungan, freq. {:d} '.format(frequency1))
plt.show()

plt.plot(time, gelombang_2,'b')
plt.xlabel('Time (s)')
plt.ylabel('Amplitude (V)')
plt.title('Gelombang 2, freq. {:d} '.format(frequency2))
plt.show()

plt.plot(time, Sinyal_Gabungan)
plt.xlabel('Time (s)')
plt.ylabel('Amplitude (V)')
plt.title('Gelombang Gabungan, freq. {:d}'.format(frequency1))
plt.show()

Menggeser Sinyal

import numpy as np
import matplotlib.pyplot as plt


freq = 2000
amp = 3
time = np.linspace(0, 2, 1000)

#sudut geser

theta = 180

signal1 = amp*np.sin(2*np.pi*freq*time)
signal2 = amp*np.sin((2*np.pi*freq*time) + theta)

plt.plot(time, signal1,'r')
plt.plot(time, signal2,'b')
plt.xlabel('Time (s)')
plt.ylabel('Amplitude (V)')
plt.title('Gelombang Sinus')

plt.show()