transformada rapida de fourier y transformada discreta de fourier
fourier
-
Upload
dicky-hernandez -
Category
Documents
-
view
1 -
download
0
description
Transcript of fourier
-
Nama : Dicky Dwi Arifiansyah
NRP : 2413100118
Laporan Praktikum Deret Fourier
Pada laporan ini saya mengambil data dari US. Energy Information Administration mengenai harga Crude Oil Prices: West Texas Intermediate (WTI) - Cushing, Oklahoma, atau yang biasa disebut
harga minyak dunia berikut ini adalah data yang terdapat di excel
DATE VALUE (US Dollar Per Barrel)
2015-04-06 52,08
2015-04-07 53,95
2015-04-08 50,44
2015-04-09 50,79
2015-04-10 51,63
2015-04-13 51,95
2015-04-14 53,30
2015-04-15 56,25
2015-04-16 56,69
2015-04-17 55,71
2015-04-20 56,37
2015-04-21 55,58
2015-04-22 56,17
2015-04-23 56,59
2015-04-24 55,98
2015-04-27 55,56
2015-04-28 57,05
2015-04-29 58,55
2015-04-30 59,62
2015-05-01 59,10
2015-05-04 58,92
Setelah itu buka matlab dan buat Variabel Editor baru kemudian dimasukkan data diatas menjadi
seperti ini
-
Setelah itu variable editor di ubah nama menjadi crudeoil seperti berikut
Ketik perintah plotFS untuk mengetahui hasil plot dari data diatas
-
Dan berikut ini hasil dari plot data
Dari plot data diatas untuk grafik Fourier series magnitude vs frequency in hertz membuktikan bahwa
frekuensi tertinggi berada di anka nol oleh karena itu data ini termasuk dalam distribusi normal.Oleh
karena itu dapat dicari tahu frekuensi dari gelombang yang ditunjukkan pada grafik plotting Periodic
signal terhadap Time dengan cara mencari masing-masing hasil kali antara frekuensi dengan
magnitude kemudian dibagi jumlah frekuensinya, sebagai berikut:
Frequency Magnitude Frequency x Magnitude
1 1,3195 1,3195
2 1,1609 2,3218
3 0,30695 0,92085
4 0,32021 1,28084
5 0,069002 0,34501
6 0,028775 0,17265
7 0,32 2,25176
8 0,22115 1,7692
9 0,26587 2,39283
10 0,31662 3,1662
=15,94064
Kemudian untuk mencari nilai frekuensi dari gelombang tersebut dengan cara mencari niai rata-rata
pada kolom Frequncy x Magnitude didapat kan sebagai berikut
-
= ( )
1
=15,94064
10= 1,594064
=1
=
1
1,594064= 0,627327
Maka frekuensi gelombang tersebut adalah 1,594064 Hz dan Periode gelombangnya adalah
0,627327 s1
Kesimpulan:
-Fourier series atau deret fourier dapat digunakan untuk mempermudah mencari frekuansi serta
periode suatu gelombang.
-Dari data diatas membuktikan bahwa data itu bukan merupakan data periodic tetapi data acak
sehingga bukan merupakan deret fourier. Fungsi matab diatas hanya untuk mengetahui berapa
frekuensi dan periode gelombangnya
-
LAMPIRAN
Syntax fungsi plotFS
function plotFS(x,T,M) % plotFS(x,T,M) % Periodic signals. x represents one period, from % t=0 to t=T-ts, where ts=T/length(x)is the sampling period. % (1) plot two periods of x(t), from -T to T. % (2) plot magnitude of FS coefficients from -M to M. % (3) plot phase of FS coefficients from -M to M. figure X=Fseries(x);
N=length(X); if N < M+1 X=[X,zeros(1,M+1-N)]; end Y=[conj(fliplr(X([2:M+1]))),X([1:M+1])]; subplot(3,1,2) f0=1/T; % fundamental frequency f=[-M*f0:f0:M*f0]; % frequency axis values stem(f,abs(Y),'b') title('Fourier Series magnitude vs. frequency in Hertz') subplot(3,1,3) stem(f,(1/pi)*angle(Y),'b') hold on title('Fourier Series phase/pi vs. frequency in Hertz') v=axis;axis([v(1:2),-1.2 1.2]); plot([v(1:2)],[0 0],'k') % do time domain panel subplot(3,1,1) L=512;N=length(X); if L
-
plot(t,xx(1:M),'g') mx=1.2*max(abs(x));axis([0 M/fs -mx mx]); xlabel('time in seconds') title('time snapshot') subplot(2,1,2) spectrm(x,fs) title('spectrum') sc=max(x); sound(x/sc,fs)