File: ggtest.m

package info (click to toggle)
hf 0.7.3-4etch1
  • links: PTS
  • area: main
  • in suites: etch
  • size: 3,684 kB
  • ctags: 3,156
  • sloc: ansic: 26,447; cpp: 4,909; sh: 3,785; makefile: 309
file content (11 lines) | stat: -rw-r--r-- 272 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
 
sz=1048576;
fd=fopen('ggtest','r');
[d,c]=fread(fd,[1 sz],'float32');
fclose(fd);
pwrdev = sum(d.^2)/sz
pwrdevdb = 10*log10(pwrdev)
f=fft(d);
fdiv=128;
ff=zeros(1,sz/fdiv);
for t=1:sz/fdiv,ff(t)=sum(abs(f(t*fdiv-fdiv+1:t*fdiv)).^2);end;
plot(ff(1:round(sz/fdiv/80))/c^2);