cara dan hasil dari PRAKTIKUM TRANSFORMASI FOORIER

cara dan hasil  dari PRAKTIKUM TRANSFORMASI FOORIER

adapun cara menggunakan  sinyal fouries ada dengan cara  kita melakukan program pada matlab.

4.1. Fenomena Gibb

1. Bangkitkan sebuah sinyal sinus dengan cara seperti berikut

%File name: fen_Gibb.m

clc;clf;

t=-3:6/100:3;

N=input(‘Masukan Jumlah Sinyal Yang Dikehendaki = ’);

c0=0.5; w0=pi; Fs=100;

xN=c0*ones(1,length(t));

for n=1:2:N;

theta=((-1)^((n-1)/2)-1)*pi/2;

xN=xN+2/n/pi*cos(n*w0*t+theta);

end subplot(211) plot(t,xN)

title('Phenomena Gibb') xlabel('Waktu(s)') ylabel('X(t)')

%Transformasi

Jumlah Sinyal 3

Jumlah Sinyal =5

Jumlah Sinyal =7

4.2. Pengamatan Frekuensi Pada Sinyal Tunggal

1. Bangkitkan sinyal sinus yang memiliki frekuensi f = 15 Hz, dan amplitudo 5 Volt.

Fs=100; t=(1:100)/Fs; f=15; A=5 ;

s=A*sin(2*pi*f*t); subplot(2,1,1) plot(t,s) xlabel('Waktu (s)')

2.   Lanjutkan langkah ini dengan memanfaatkan fungsi fft untuk mentranformasi sinyal ke dalam domain frekuensi

S=fft(s,512); w=(0:255)/256*(Fs/2); subplot(2,1,2) plot(w,abs(S(1:256))) xlabel('Frekuensi (Hz)')

4.  Cobalah anda merubah nilai f=5, 20 dan 30 Apa yang anda lihat pada gambar sinyal anda?

5.  Cobalah merubah nilai amplitudo dari 5 volt menjadi 7, 15 atau 20. Apa yang terjadi pada sinyal anda?

F =10                      A = 2

F =20                      A = 4

F = 5                        A = 5

F =10                      A = 2

F =20                      A = 4

F = 5                        A = 5

4.3. Pengamatan Frekuensi Pada Kombinasi 2 Sinyal

1. Caranya adalah dengan mengetik program berikut ini

clc;clf; fs=100; t=(1:400)/fs; f1=10;

s1=5*sin(2*pi*f1*t); f2=30; s2=3*sin(2*pi*f2*t); s=s1+s2; subplot(211)

plot(t,s)

title('Dua Sinyal Sinus') xlabel('Waktu(s)') ylabel('x(t)') S=fft(s,512); w=(0:255)/256*fs/2; Sab=abs(S); subplot(212) plot(w,Sab(1:256))

title('Sinyal pada domain frekuensi')

xlabel('Frekuensi(Hz)')

ylabel('x(f)')

2. Ubahlah nilai f2 10,25,20sedangkan f1 dan amplitudo tetap. Apa yang anda dapatkan dari langkah ini? Plot semua hasil running program.

f2 =10

f2 =25

f2 =20

3. Coba ubah nilai amplitudo pada sinyal KEDUA menjadi 1, 5 DAN 10  Apa yang anda dapatkan dari langkah ini

A = 1

A = 5

A = 10

4.4. Pengamatan Frekuensi Pada Kombinasi 4 Sinyal

1.  Pada percobaan berikut ini anda coba bangkitkan 4 sinyal sinus dengan frekuensi f1, f2, f3, dan f4.

clc;clf; fs=100; t=(1:100)/fs; f1=5;

s1=50*sin(2*pi*f1*t); f2=15; s2=40*sin(2*pi*f2*t); f3=25; s3=30*sin(2*pi*f3*t); f4=35; s4=20*sin(2*pi*f4*t); s=s1+s2+s3+s4;; subplot(211)

plot(t,s)

title('Empat Sinyal Sinus')

xlabel('Waktu(s)')

ylabel('x(t)') S=fft(s,512); w=(0:255)/256*fs/2; Sab=abs(S); subplot(212) plot(w,Sab(1:256))

title('Sinyal pada domain frekuensi')

xlabel('Frekuensi(Hz)')

ylabel('x(f)')

2. Ubah nilai f2 =10, f3 = 30 dan f4 =30, sedangkan amplituda tetap. Apa yang anda dapatkan dari langkah ini?

4.5. Pengamatan Frekuensi Pada Kombinasi 6 Sinyal

Fs=100; t=(1:200)/Fs; f1=2;

s1=20*sin(2*pi*f1*t); f2=5; s2=15*sin(2*pi*f2*t); f3=15; s3=10*sin(2*pi*f3*t); f4=20; s4=7*sin(2*pi*f4*t); f5=35; s5=5*sin(2*pi*f5*t); f6=45; s6=3*sin(2*pi*f6*t); s=s1+s2+s3+s4+s5+s6; subplot(211)

plot(t,s)

title('Enam Sinyal Sinus') xlabel('Waktu(s)') ylabel('x(t)')

S=fft(s,512); w=(0:255)/256*fs/2; Sab=abs(S); subplot(212) plot(w,Sab(1:256))

title('Sinyal pada domain frekuensi')

xlabel('Frekuensi (Hz)')

ylabel

4.6. Pengamatan Frekuensi Pada Sinyal Audio

1.   Buat program pemanggil file audio *.wav sebagai berikut :

clc;clf; [y,Fs]=wavread('gundul_pacul.wav'); Fs=16000;

subplot(211)

plot(y(100:10000))

title('Sinyal Audio') xlabel('Waktu(s)') ylabel('x(t)') S=fft(y);

%w=(0:255)/256*fs/2; Sab=abs(S); subplot(212) plot(Sab(100:10000))

% plot(Sab)

title('Sinyal audio domain frekuensi')

xlabel('Frekuensi(Hz)')

ylabel('x(f)')