N=128;n=0:N-1;fund=3/N;%0+.0167/2;
x=sin(2*pi*fund*n)+sin(2*3*pi*fund*n)/3+...
sin(2*pi*5*fund*n)/5+sin(2*pi*7*fund*n)/7;
%x=x.*(hamming(N))';
DFT3D_1(x)
x=sin(2*pi*fund*n)+sin(2*3*pi*fund*n)/3+...
sin(2*pi*5*fund*n)/5+sin(2*pi*7*fund*n)/7;
%x=x.*(hamming(N))';
DFT3D_1(x)
function DFT3D_1(x)
N=length(x);
n=0:N-1;
%x=x.*(hamming(N))';
X=fft(x)*2/N;
X1=abs(X')*ones(1,N);X1=X1';
N=length(x);
n=0:N-1;
%x=x.*(hamming(N))';
X=fft(x)*2/N;
X1=abs(X')*ones(1,N);X1=X1';
f=(0:N-1)/N; % Harmonics
[ff,nn]=meshgrid(f,n);
W=cos(2*pi*f'*n);
[ff,nn]=meshgrid(f,n);
W=cos(2*pi*f'*n);
X1=W'.*X1;
X1(:,1)=x;
X1(:,1)=x;
plot3(ff,nn,X1,'LineWidth',1);grid on
xlabel('frequency');ylabel('time seconds');
axis([0 1 0 N -1.3 1.3]);
view([37.5 30])


No comments:
Post a Comment