Project Code for DFT in addition to IDFT
Here is the projection code for DFT in addition to IDFT. Just re-create it, glue it yours to the editor.Note: There mightiness live on bugs in addition to fault inward below code hence delight if you lot constitute whatever employment inward below code delight produce comment to own it. Your proposition is valuable to us.
Check Below code for the to a higher house given formula. Where the sequence of betoken is 4.
| clc clear closed all x = [0 i 2 3]; north = length(x); X = zeros(4,1); for k = 0:N-1 for n = 0:N-1 % Here inward below formula nosotros substitute the value of Omega. X(k+1) = X(k+1) + x(n+1)*exp(-j*pi/2*n*k); end end Answer = [X] |
Answer =
6.0000 + 0.0000i
-2.0000 + 2.0000i
-2.0000 - 0.0000i
-2.0000 - 2.0000i
Example 2 In given below nosotros perform DFT in addition to IDFT equally Linear Transformation.
| %% % All write Reserved Telecom-academy.blogspot.com % You are allowed to edit, re-create this nether the next status % You must render amount credit to Author. % You are non allowed to take away this Credit. %% clc % It clear upwardly the ascendancy window clear % It clears upwardly workspace closed all % this ascendancy used to closed all previous opened upwardly figure. Influenza A virus subtype H5N1 = [1,1,1,1;1,-1i,-1,1i;1,-1,1,-1;1,1i,-1,-1i ]; % Influenza A virus subtype H5N1 exceptional matrix operate to honour dft of sequnce(Lenght 4) B = [0 i 2 three ]; % nosotros volition honour sequence of the DFT. % We are changing row matrix equally per in addition to equally good this volition allows the % multiplication of ii matrices ... % X = [0; 1; 2; 3;] X = [0; 1; 2; 3;]; % We shop sequence inward colums matrix to let multiplication. DFT_Matrix = A*X; % Formulat operate to honour the DFT of the sequecne. %% % disp ascendancy is used display output inward the ascendancy window. disp('*************************Inputs Sequence***************************') disp(' X(n) = (0 i 2 3)') disp('**************Linear transformation of DFT method 1****************') disp('To honour DFT nosotros volition Multiply bellow matrix amongst sequnce(Colums Matrix)') disp('Matrix = [1,1,1,1;1,-1i,-1,1i;1,-1,1,-1;1,1i,-1,-1i ];') disp('DFT = ') disp(DFT_Matrix) disp('****************Linear transformation of DFT method 2***************') %% exp = 2.7182818284590452353602874713527; %Exact Value of EXP north = 4; W1 = exp^-((0+1i*2*pi*1)/N); % Formula is uses is to honour the Omega value W2 = exp^-((0+1i*2*pi*2)/N); % in addition to hence nosotros multiply amongst colum matrix W3 = exp^-((0+1i*2*pi*3)/N); W4 = exp^-((0+1i*2*pi*2)/N); W5 = exp^-((0+1i*2*pi*0)/N); W6 = exp^-((0+1i*2*pi*2)/N); W7 = exp^-((0+1i*2*pi*3)/N); W8 = exp^-((0+1i*2*pi*2)/N); W9 = exp^-((0+1i*2*pi*1)/N); Matrix = [1 i 1 i ; i W1 W2 W3 ; i W4 W5 W6 ; i W7 W8 W9]; X2 = Matrix*X; disp('DFT =') disp(X2) disp('******************************IDFT*********************************') W10 = exp^-((0-1i*2*pi*1)/N); % Omega value to honour IDFT. W11 = exp^-((0-1i*2*pi*2)/N); W12 = exp^-((0-1i*2*pi*3)/N); W13 = exp^-((0-1i*2*pi*2)/N); W14 = exp^-((0-1i*2*pi*0)/N); W15 = exp^-((0-1i*2*pi*2)/N); W16 = exp^-((0-1i*2*pi*3)/N); W17 = exp^-((0-1i*2*pi*2)/N); W18 = exp^-((0-1i*2*pi*1)/N); due west = [ i 1 i 1; i W10 W11 W12; i W13 W14 W15; i W16 W17 W18]; IDFT = 1/N*(W*X2); disp('IDFT =') disp(IDFT) disp('The output of IDFT tin live on writen equally X(n)= (0 i 2 3) ') disp('*******************************END ********************************') |
