function hw4Prob2

clc
load('handel.mat')
hallelujah=y(1:65536);
hal_power=abs(fft(hallelujah));
plot(hal_power)

x=5
R=zeros(x);
for i=1:x
   for j=1:x
      if (i>1 && j>1)
         if ((i-1)==x-(j-1))
            R(i,j)=1; 
         else 
            R(i,j)=0; 
         end;
      else if (i==j)
         R(i,j)=1;
      else 
         R(i,j)=0; 
      end;
   end;
end;

fft_x=fft(y(1:x));
%[r1,c1]=size(R)
%[r2,c2]=size(fft_x)
o=fft_x-(R*(conj(fft_x)))

end;

