example_1_5.rar_phasor

clc clear all format short %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Example 1.5 clc clear all N=24; f0=60; t=0:(1/(f0*N)):.1; x=cos(120.1*pi*t); plot(t,x) close xt= x(1:N); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%DFT for n=1:N sum=0; for k=1:N sum=sum+xt(1,k)*exp(j*2*pi*(k-1)*(n-1)/N); end DFT(n,1)=sum; Samplenumber(n,1)=n-1; xtt(n,1)=xt(1,n); end X=DFT/N; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% for n=1:((N/2)-1) a0=2*X(1,1); ak(n,1)=2*real(X((n+1),1)); bk(n,1)=2*imag(X((n+1),1)); end outexam5 (:,1)=Samplenumber; outexam5 (:,2)=xtt; outexam5 (:,3)=DFT; outexam5 (:,4)=X amplitude=sqrt(ak(1,1)^2+bk(1,1)^2) phase=atan(bk(1,1)/ak(1,1))*(180/pi)