function y = kotak(x) t = x(1); % Referensi waktu t [second] RPM = x(2); % Kecepatan Motor [RPM] RPS = RPM/60; % Kecepatan Motor [RPS] n = 1/80; % Perbandingan Roda Gigi 80 kali motor berputar, 1 kali conveyor belt berputar cawan = n*RPS; % Jumlah cawan yang lewat per second [Hertz], frekuensi cawan periode_cawan = 1/cawan; % Periode berlalunya cawan di bawah katup [second] fungsi_sinus = sin(2*pi*cawan*t); % fungsi sinus periodik durasi = (1/periode_cawan)*2*pi; % konversi dari durasi dalam second ke dalam radian teta = (pi - durasi)/2; % sudut untuk menentukan batas level [radian] batas_level = sin(teta); % batas level untuk membuat kotak % Membuat gelombang kotak dari gelombang sinus: if fungsi_sinus >= batas_level y = 1; else y = 0; end