Advertisement
Guest User

Untitled

a guest
Aug 17th, 2017
56
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.49 KB | None | 0 0
  1. function [cb p] = aufg2a(i)
  2.  
  3. N = 20; m =5;
  4. f1 = sin(4*linspace(0,2*pi,N));
  5. f2 = max(pi/2 - abs(linspace(0,2*pi,N)-pi/2),pi/2-abs(linspace(0,2*pi,N)-3*pi/2));
  6. f3 = round((sin(2*linspace(0,2*pi,N))+1)/2);
  7.  
  8. f =[f1 ; f2 ; f3];
  9.  
  10. fk = f(i,:);
  11.  
  12. N = length(fk);
  13. M = N*2^m;
  14. x = 2*pi/N*[0:N-1]';
  15. fig = plot(x,fk,'or');
  16.  
  17. c = fft(fk)/N;
  18. cb = c;
  19. c = [c(1:N/2) zeros(1,M-N) c(N/2+1:N)];
  20. p = real( M*ifft(c) );
  21. hold on;
  22. plot(linspace(0,2*pi,M),p);
  23. hold off;
  24. saveas(fig,'aufg2a.pdf');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement