10/30/15 5:16 AM D:\Document\Semester V\tugasakustik.m 1 of 1

% Program Koefisien Pantulan (R) clear allclose all r1 = 1000;r2 = 900; c1 = 1000;c2 = 900; rho1 = r1/c1;rho2 = r2/c2; n=c1/c2;m=rho1/rho2; % Sudut Datang (thetai)thetai=linspace(0,pi/2,100); % sin (thetai) = (si) dan Cos (thetai) = cisi=sin(thetai);ci=cos(thetai); % Sudut Transmisi (thetat)% sin (thetat) = (st) dan Cos (thetat) = ctst = 1/n*si;thetat = asin(st);ct = cos(thetat); % Z1/Z2 = Z1Z2% R = (1 - Z1/Z2)/(1+Z1/Z2) Z1Z2 = n*m*ct./ci; R = abs((1-Z1Z2)./(1+Z1Z2)); plot(thetai*180/pi, R)xlabel('Angle of Insidence')ylabel('|R|')