Dynamic Blinkie Text Generator at TextSpace.net

Feedburner

I heart FeedBurner

Senin, 12 Oktober 2015

matlab




TEORI
Sinyal Sinus : 
y = A sin (2πft + φ)
A = Amplitudo
f = frekuensi
t = waktu (dalam detik)
φ = fase sinyal
MEMULAI MATLAB
1.         Jalankan MATLAB
2.         Pada command window ketik “edit   jajal” kemudian tekan ENTER
            >> edit  jajal    [ENTER]
3.         Maka akan muncul  MATLAB  Editor berikut
Disinilah tempat untuk menuliskan program MATLAB
4.         ketiklah program berikut
t = 0:0.01:1;
f = 2;
y = 2*sin(2*pi*f*t);
plot(t,y);
5.         jalankan program dengan menekan tombol: F5, lihat hasilnya
SINYAL SINUS TUNGGAL
1.         Membangkitkan sinyal sinus tunggal
A = 3,              f = 1 Hz,          φ = 0,               t = 0 s/d 1
clc;                % membersihkan layar
clear all;          % membersihkan isi variabel
f = 1;              % frekuensi sinyal
t=0:0.01:1;         % t = 0, 0.01, 0.02,0.03,…, 1
y=3*sin(2*pi*f*t);  % sinyal sinus dengan amplitudo 3
plot(t,y);          % menggambar kurva dilayar
axis([0 1 -10 10])  % batasan sumbu x dan sumbu y
grid on             % tampilan layar dibuat kotak-2
2.         Mengamati perbedaan frekuensi sinyal dalam domain waktu
Gantilah nilai f dengan : 2, 3, 4 atau 5. Lihat hasilnya digambar dan bandingkan!
3.         Mengamati perbedaan Amplitudo
Gantilah nilai A dengan : 1, 2, 3, 4 atau 5. Lihat hasilnya digambar dan bandingkan!
4.         Menambahkan offset DC pada sinyal
y = K + A sin (2πft + φ)         K = konstanta (offset DC)
clc;                % membersihkan layar
clear all;          % membersihkan semua isi variabel yang pernah dibuat
f = 1;              % frekuensi sinyal
t=0:0.01:1;
y=4+3*sin(2*pi*f*t);
plot(t,y);          % menggambar kurva dilayar
axis([0 1 -10 10])  % batasan sumbu x dan sumbu y
grid on
Gantilah nilai K dengan : −5, −4, 3, 4 atau 5. Lihat hasilnya digambar dan bandingkan!
5.         Membangkitkan sinyal sinus tunggal, menampilkan suara dan menyimpan ke hardisk
clc;
clear all;
f=800;                      %frekuensi sinyal
fs=16000;                   %frekuensi sampling
t=1/fs:1/fs:1;
y=sin(2*pi*f*t);
sound(y,fs)                 %menampilkan suara
wavwrite(y,fs,'sinus1.wav') %menyimpan dg nama: sinus1.wav
plot(t,y);axis([0 1  -1 1]);
title('Sinyal Sinus (f=800 Hz), sampling 16000 Hz')
untuk memperpanjang atau memperpendek suara gantilah nilai t dengan berikut:
t=1/fs:1/fs:2;    t=1/fs:1/fs:3;    t=1/fs:1/fs:4;    t=1/fs:1/fs:0.5;
            perhatikan perbedaannya dengan program berikut
clc;
clear all;
fs=16000;
f=800;
t=0:0.001:1.0;
y=sin(2*pi*f*t);
wavplay(y,fs)
plot(t,y);axis([0 0.2 -1 1]);
title('Sinyal Sinus (f=800 Hz), sampling 16000 Hz')
6.         Mengamati perbedaan frekuensi sinyal dalam domain frekuensi
clc;
clear all;
f=800;                      %frekuensi sinyal
fs=16000;                   %frekuensi sampling
t=1/fs:1/fs:1;
y=sin(2*pi*f*t);
sound(y,fs)
Y_f=(abs(fft(y)));
plot(Y_f)
axis([0 1500 0 100])
Gantilah nilai f dengan : 200, 700, 900 atau 1000. Lihat hasilnya digambar dan suara kemudian bandingkan!
Gantilah nilai fs dengan : 8000, 10000, 15000 atau 20000. Lihat hasilnya digambar dan suara kemudian bandingkan!
7.         Mengamati perbedaan Amplitudo
Gantilah nilai A dengan : 2, 3, 4 atau 5. Lihat hasilnya digambar dan suara kemudian bandingkan!
8.         Membangkitkan sinyal Sirine
fs=16000;
t=0:1/fs:4;
alfa=2*pi*(150 + 250*t + 550*t.*t);
y= 10*sin(alfa);
sound(y,fs);
SINYAL SINUS GANDA
1.         Membangkitkan 2 buah sinyal sinus
A1 = 10,                       f 1= 1 Hz,        φ1 = 0,                        
A2 = 1,                         f 2= 20 Hz,      φ2 = 0,                        
t = 0 s/d 1
clc;                        % membersihkan layar
clear all;                  % membersihkan semua isi variabel yang pernah dibuat
f1 = 1;                     % frekuensi sinyal
f2 = 20;
t=0:0.01:1;
y1=10*sin(2*pi*f1*t);
plot(t,y1);
hold on;
y2=sin(2*pi*f2*t);
plot(t,y2);
axis([0 1 -15 15])          % batasan sumbu x dan sumbu y
grid on
2.         menampilkan 2 buah grafik dalam 2 window
clc;                        % membersihkan layar
clear all;                  % membersihkan semua isi variabel yang pernah dibuat
f1 = 1;                     % frekuensi sinyal
f2 = 20;
t=0:0.01:1;
y1=10*sin(2*pi*f1*t);
subplot(2,1,1);plot(t,y1);  % 2 baris, 1 kolom, posisi ke-1
axis([0 1 -15 15])          % batasan sumbu x dan sumbu y
grid on
y2=sin(2*pi*f2*t);
subplot(2,1,2);plot(t,y2);  % 2 baris, 1 kolom, posisi ke-2
axis([0 1 -15 15])          % batasan sumbu x dan sumbu y
grid on
OPERASI 2 BUAH SINYAL
1.         Penjumlahan
clc;                        % membersihkan layar
clear all;                  % membersihkan semua isi variabel yang pernah dibuat
f1 = 1;                     % frekuensi sinyal ke-1
f2 = 20;                    % frekuensi sinyal ke-2
t=0:0.01:1;
y1=10*sin(2*pi*f1*t);
y2=sin(2*pi*f2*t);
y3= y1+y2;
subplot(3,1,1);plot(t,y1);  % 3 baris, 1 kolom, posisi ke-1
axis([0 1 -15 15]);
subplot(3,1,2);plot(t,y2);  % 3 baris, 1 kolom, posisi ke-2
axis([0 1 -15 15]);
subplot(3,1,3);plot(t,y3);  % 3 baris, 1 kolom, posisi ke-3
axis([0 1 -15 15])          % batasan sumbu x dan sumbu y
grid on
2.         Pengurangan
clc;                        % membersihkan layar
clear all;                  % membersihkan semua isi variabel yang pernah dibuat
f1 = 1;                     % frekuensi sinyal
f2 = 20;
t=0:0.01:1;
y1=10*sin(2*pi*f1*t);
y2=sin(2*pi*f2*t);
y3= y1-y2;
subplot(3,1,1);plot(t,y1);  % 3 baris, 1 kolom, posisi ke-1
axis([0 1 -15 15]);
subplot(3,1,2);plot(t,y2);  % 3 baris, 1 kolom, posisi ke-2
axis([0 1 -15 15]);
subplot(3,1,3);plot(t,y3);  % 3 baris, 1 kolom, posisi ke-3
axis([0 1 -15 15])          % batasan sumbu x dan sumbu y
grid on
3.         Perkalian
clc;                        % membersihkan layar
clear all;                  % membersihkan semua isi variabel yang pernah dibuat
f1 = 1;                     % frekuensi sinyal
f2 = 20;
t=0:0.01:1;
y1=10*sin(2*pi*f1*t);
y2=sin(2*pi*f2*t);
y3= y1.*y2;
subplot(3,1,1);plot(t,y1);  % 3 baris, 1 kolom, posisi ke-1
axis([0 1 -15 15]);
subplot(3,1,2);plot(t,y2);  % 3 baris, 1 kolom, posisi ke-2
axis([0 1 -15 15]);
subplot(3,1,3);plot(t,y3);  % 3 baris, 1 kolom, posisi ke-3
axis([0 1 -15 15])          % batasan sumbu x dan sumbu y
grid on
4.         Penguatan dan Pelemahan Sinyal
y = A sin (2πft + φ)
            jika  A > 1       maka terjadi penguatan sinyal
            jika  0 < A < 1 maka terjadi pelemahan sinyal
clc;                        % membersihkan layar
clear all;                  % membersihkan semua isi variabel yang pernah dibuat
f = 1;                      % frekuensi sinyal
t=0:0.01:1;        
y1=10*sin(2*pi*f*t);        % Sinyal asli dengan amplitudo 10
y2=15*sin(2*pi*f*t);        % Sinyal asli diperkuat menjadi  3/2 kali
y3=5*sin(2*pi*f*t);         % Sinyal asli diperlemah menjadi 1/2 kali
subplot(3,1,1);plot(t,y1);  % 3 baris, 1 kolom, posisi ke-1
axis([0 1 -15 15]);
subplot(3,1,2);plot(t,y2);  % 3 baris, 1 kolom, posisi ke-2
axis([0 1 -15 15]);
subplot(3,1,3);plot(t,y3);  % 3 baris, 1 kolom, posisi ke-3
axis([0 1 -15 15]);         % batasan sumbu x dan sumbu y
grid on
6.         Pengamatan Operasi dua buah sinyal melalui suara
clc;
clear all;
fs= 16000;
f1= 852;
f2= 1209;
t=0:0.001:1.0;
y1=sin(2*pi*f1*t);
wavplay(y1,fs)
plot(t,y1);axis([0 1 -1 1]);
title('Sinyal Sinus Y1:(f=852 Hz), sampling 16000 Hz')
pause
y2=sin(2*pi*f2*t);
wavplay(y2,fs)
plot(t,y2);axis([0 1 -1 1]);
title('Sinyal Sinus Y2:(f=1209 Hz), sampling 16000 Hz')
pause
y3 = y1 + y2;
wavplay(y3,fs)
plot(t,y3);axis([0 1 -2 2]);
title('Sinyal Sinus Y3 = Y1+Y2')
pause
y4 = y1 - y2;
wavplay(y4,fs)
plot(t,y4);axis([0 1 -1 1]);
title('Sinyal Sinus Y4 = Y1-Y2')
pause
y5 = y1 .* y2;
wavplay(y5,fs)
plot(t,y5);axis([0 1 -1 1]);
title('Sinyal Sinus Y5 = Y1*Y2')
pause
FREKUENSI SUARA MUSIK
Nada dasar dari suara musik didasarkan pada nada A dengan frekuensi 440 Hz. Untuk frekuensi nada-nada yang lain dihitung menggunakan rumus berikut
f  = 440 x 2n/12  Hz
Dimana :
f  = Frekuensi dari nada-nada yang lain
440 = Ketentuan nilai untuk nada dasar A
n = langkah (oktaf) dari nada-nada yang akan di cari.
12 = jumlah oktaf  ( C, C#, D, D#, E, F, F#, G, G#, A, A#, B)
Sebagai contoh,  untuk mencari frekuensi nada F, maka n = −4, untuk mencari frekuensi nada B, maka n = 2.
F = 440 x 2 −4/12 = 349.2282 Hz ≈ 349.2 Hz
B = 440 x 2 2/12 = 493.8833 Hz ≈ 493.9 Hz
Berdasarkan rumus diatas, dapat dicari nilai-nilai frekuensi untuk nada-nada yang lain yaitu :
Oktaf 1
Oktaf 2
Oktaf 3
C
130.8
261.6
523.3
C#
138.6
277.2
554.4
D
146.8
293.7
587.3
D#
155.6
311.1
622.3
E
164.8
329.6
659.3
F
174.6
349.2
698.5
F#
185.0
370.0
740.0
G
196.0
394.0
784.0
G#
207.7
415.3
830.6
A
220.0
440.0
880.0
A#
233.1
466.2
932.3
B
247.0
493.9
987.8
%Judul Lagu: Gundul-gundul-Pacul
Fs=16000;
t=0:1/Fs:0.25;
c=sin(2*pi*262*t);
d=sin(2*pi*294*t);
e=sin(2*pi*330*t);
f=sin(2*pi*249*t);
g=sin(2*pi*392*t);
a=sin(2*pi*440*t);
b=sin(2*pi*494*t);
c1=sin(2*pi*523*t);
nol = [zeros(size(t))];
nada1 = [c,e,c,e,f,g,g,nol,b,c1,b,c1,b,g,nol,nol];
nada2 = [c,e,c,e,f,g,g,nol,b,c1,b,c1,b,g,nol];
nada3 = [c,nol,e,nol,g,nol,f,f,g,f,e,c,f,e,c,nol];
nada4 = [c,nol,e,nol,g,nol,f,f,g,f,e,c,f,e,c];
lagu=[nada1,nada2,nada3,nada4];
sound(lagu,Fs);
%Judul Lagu: Si-Unyil
Fs=16000;
t=0:1/Fs:0.25;          %panjang not ¼ detik
c=sin(2*pi*262*t);
d=sin(2*pi*294*t);
e=sin(2*pi*330*t);
f=sin(2*pi*249*t);
g=sin(2*pi*392*t);
a=sin(2*pi*440*t);
b=sin(2*pi*494*t);
c1=sin(2*pi*523*t);
nol = [zeros(size(t))];
unyil = [b,nol,a,nol,g,nol,a,nol,b,nol,a,nol,a,nol,nol,nol,b,nol,a,nol,g,nol,e,nol,g,nol,e,nol,e];
sound(unyil,Fs);
wavwrite(unyil,'unyil.wav');
%Nada Pesawat telepon 
%===========================
% Freq : 1209 | 1336 | 1477 
%===========================
% 697  :  1   |   2  |  3  
%---------------------------
% 770  :  4   |   5  |  6
%---------------------------
% 852  :  7   |   8  |  9 
%---------------------------
% 941  :  *   |   0  |  #
%---------------------------
Fs=9000;
t=0:0.001:1.5;
y9=sin(2*pi*852*t)+sin(2*pi*1477*t);
y8=sin(2*pi*852*t)+sin(2*pi*1336*t);
y7=sin(2*pi*852*t)+sin(2*pi*1209*t);
y6=sin(2*pi*770*t)+sin(2*pi*1477*t);
y5=sin(2*pi*770*t)+sin(2*pi*1336*t);
y4=sin(2*pi*770*t)+sin(2*pi*1209*t);
y3=sin(2*pi*697*t)+sin(2*pi*1477*t);
y2=sin(2*pi*697*t)+sin(2*pi*1336*t);
y1=sin(2*pi*697*t)+sin(2*pi*1209*t);
y0=sin(2*pi*941*t)+sin(2*pi*1336*t);
y = [y0,y1,y2,y3,y4,y5,y6,y7,y8,y9];
wavplay(y,Fs)
SINYAL WICARA PADA MANUSIA
Sinyal wicara adalah sinyal yang berubah-ubah terhadap waktu, bila diamati dalam selang waktu yang sangat singkat masih tampak stasioner. Tetapi bila diamati dalam selang waktu yang agak lama (> 1/5 detik), sinyalnya mulai menampakkan ucapan dari pembicara.   
Untuk menyajikan sebuah sinyal wicara ada tiga kondisi dasar yang harus diketahui, yaitu silence (S) (keadaan tenang) yaitu keadaan dimana sinyal wicara tidak diproduksi, unvoice (U) (tidak bersuara) yaitu vocal cord tidak berfibrasi, dan  voiced (V) (bersuara) dimana  vocal cord bervibrasi secara periodik sehingga menghasilkan sinyal wicara.
Merekam Suara dan menyimpan menjadi  “File.wav”
clear all;
fs=8000;
y= wavrecord(5.0*fs, fs, 'double'); %merekam suara
wavwrite(y,fs,'aiueo.wav');         %simpan rekaman ke hardisk
wavplay(y,fs);                      %mainkan hasil rekaman
figure,plot(y);                     %sinyal hasil rekaman di plot
Membaca dan memainkan file.wav
clear all;
fs=8000;
y=wavread('aiueo.wav');
wavplay(y,fs)
figure,plot(y);
Memotong sinyal
clear all;
fs=8000;
y=wavread('aiueo.wav');
wavplay(y,fs)
figure,plot(y);
A1= min(y);
A2= max(y);
t=length(y);
x1=0; x2=t;
axis([x1  x2   A1  A2]);
% suara ‘U’ diperkirakan berada pd t=16750 s/d 19500
y1=y(16750:19500);      %sinyal y dicopy pd t=16750 s/d 19500
wavplay(y1,fs)
figure,plot(y1)
axis([x1  x2   A1  A2]);
clear all;
fs=8000;
y=wavread('aiueo.wav');
wavplay(y,fs)
figure,plot(y);
A1= min(y);
A2= max(y);
t=length(y);
x1=0; x2=t;
axis([x1  x2   A1  A2]);
y1=y(16750:19500);
wavplay(y1,fs)
figure,plot(y1)
axis([x1  x2   A1  A2]);
y2=[y1' y' y1'];
t=length(y2);
wavplay(y2,fs)

Tidak ada komentar:

Posting Komentar