fork(1) download
copy 
  1. n = 500;
  2. t = linspace(0, 10*pi, n);
  3. x = sin(2*pi*t);
  4. y = sin(2*pi*t + pi/2);
  5.  
  6. ex=sum(abs(x).^2);    
  7. ey=sum(abs(y).^2);
  8. fprintf("Time domain:\nex=%.1f, ey=%.1f\n", ex, ey)
  9. ex=sum(abs(fft(x)).^2)/n;
  10. ey=sum(abs(fft(y)).^2)/n;  
  11. fprintf("Frequency domain (rectwin):\nex=%.1f, ey=%.1f\n", ex, ey)
  12. w = hamming(n)';
  13. ex=sum(abs(fft(x.*w)).^2)/sum(w.^2);
  14. ey=sum(abs(fft(y.*w)).^2)/sum(w.^2);
  15. fprintf("Frequency domain (hamming):\nex=%.1f, ey=%.1f\n", ex, ey)
Success #stdin #stdout 0.17s 65200KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
Time domain:
ex=250.1, ey=249.9
Frequency domain (rectwin):
ex=250.1, ey=249.9
Frequency domain (hamming):
ex=250.0, ey=250.0
https://ideone.com/8a9RPA
language:
Octave (octave 4.4.1)
created:
9 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!