MATLAB program to identify tool wear through the characteristics of the sound produced by the chisel and workpiece in the turning process

 

PROGRAM MATLAB

TIME DOMAIN, FREQUENCY DOMAIN DAN FAST FOURIER TRANSFORM (FFT)

LANGKAH-LANGKAH MENJALAKANNYA

BY

AHMAD SUKRI HARAHAP, S.T 

DOSEN 

Dr. Eng, MEIFAL RUSLI

ISMET HARI MULYADI, Ph.D

 

1. Program untuk mendapatkan Time Domain

a)      Buka matlab

b)      Membuat program seperti yang telah dilampirkan

c)      Mikrofon diletakkan di dekat pahat dan benda kerja

d)     Mikrofon disambungkan ke laptop

e)      Menjalankan pemotongan dan program matlab secara bersamaan

f)       Untuk menjalankan program matlab, klik Run pada menu Editor

g)      Ketika selesai merekam, simpan data dalam bentuk exsl dan simpan gambar dalam bentuk (.JPG)

h)      Gambar yang muncul adalah time domain

 

PROGRAM :

 

clear all;

clc;

Fs = 40000;

recObj = audiorecorder(Fs,24,1);

ztime = 5;

disp('PEREKAMAN SINYAL SUARA MELALUI SOUND CARD');

disp('-----------------------------------------');

disp('SILAHKAN TEKAN ENTER UNTUK MEMULAI MEREKAM...');

pause

disp('Sedang Merekam...');

disp('Silahkan tunggu...');

disp('...................');

 

recordblocking(recObj, ztime);

disp('Selesai Merekam.');

 

play(recObj);

y=(getaudiodata(recObj))';

t=1/Fs:1/Fs:ztime;

N = length(t);

z = zeros(N,2);

plot(t,y);

xlabel('time (s)');

ylabel('voltage (volt)');

FILENAME= input ('Nama File Penyimpanan : ','s');

z(:,1) = t;

z(:,2) = y;

%save(FILENAME,'t','y','-ascii')

 

xlswrite(FILENAME,z);

 

% fid=fopen(FILENAME,'w');

%

% Formatspec1='Nama file data :  %s      \r';

% fprintf(fid,Formatspec1,FILENAME);

% Formatspec2='Sampling rate :  %d     \r';

% fprintf(fid,Formatspec2,Fs);

% fprintf(fid,'%3.5f      %3.5f  \r',z');

% status=fclose(fid);

 

2. Program untuk mendapatkan frequency domain

a)      Buka matlab

b)      Membuat program seperti yang telah dilampirkan

c)      Pastikan semua data exsl yang disimpan terbaca dalam bentuk nilai bukan text

d)     Pada filename ketik nama data yang disimpan

e)      Klik Run pada menu Editor yang akan menampilkan figure 1, 2, dan 3

f)       figure 1 adalah Time Domain

g)      figure 2 adalah Frequency Domain

h)      figure 3 adalah Spectrogram yaitu gabungan Time Domain dan Frequency Domain

i)        Simpan gambar dalam bentuk (.JPG)

 

PROGRAM :

 

clear all;clc;close all;

 

filename = 'n108tanpamemotong.xlsx';

tsound = xlsread(filename);

Fs = 10000;

 

t  = tsound(1:100000,1);

y = tsound(1:100000,2);

 

N  = length(t);

n  = 1:1:length(t);

p  = hann(N);

dT = max(t)/N;

 

f  = (t./dT)*1/dT/N;

fY = fft(y.*p/N);

Lp = 20*log(abs(fY)/0.00002);

 

 

figure(1)

plot(t,y); grid;

xlabel('time'); ylabel('Sound Pressure (Pa)')

axis([0 10 -0.45 0.45]);

%

figure (2)

set(gca,'FontName','Times New Roman',...

      'FontSize',12);

xlabel('Frequency [Hz]',...

      'FontName','Times New Roman',...

      'FontSize',12);

ylabel('Sound Pressure [dB]',...

      'FontName','Times New Roman',...

      'FontSize',12);

plot(f(1:N/5),Lp(1:N/5));grid;

xlabel('frequency [Hz]'); ylabel('SPL (dB)')

 

figure (3)

spectrogram(y./1000,3000,2000,3000,10000);grid;

axis([0 1500 0 10]); 

Don't be plagiarize if you want a blessing.