example8_1.m

来自「用dsp解压mp3程序的算法」· M 代码 · 共 11 行

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% example8_1.m - Program for example 8-1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

In_sig = [1 1 1 1 0 0 0 0];        % define x(n)
N = length(In_sig);                % length N
X_in = In_sig*dftmtx(N);           % DFT
In_recover = (1/N)*(X_in*conj(dftmtx(N)));  % IDFT
figure;
subplot(2,1,1); stem(In_sig); title('Original signal');
subplot(2,1,2); stem(real(In_recover)); title('DFT->IDFT');