[討論] fft能量積分已回收
請問各位大大
下以是我的寫的程式
RR=rri;
a = RR;
forg=a(1:512);
l=length(forg);
Tmeas=sum(forg)/1000;
Tmean=Tmeas/l;
Fs=1/Tmean;
Fmean=mean(forg);
F=forg-Fmean;
f = Fs*(0:(l/2))/l;
df = max(f)/(l/2);
FF = fft(F);
PffF = FF.*conj(FF)/l;
ddf = Fs/l;
m1 = 0.01/ddf; <----這裡是小弟要問的問題,因為要算某一段頻率的積分
m2 = 0.04/ddf; 我加上了這一段去算某一區間下的能量和,但是算出來的結果
m3 = 0.15/ddf; 跟套裝軟體算的不一樣,不知是那一段出了問題,我也用sin波
m4 = 0.4/ddf; 去做測試,試出來是錯的@@,請大大幫幫忙,是不是我設定的方式
錯了
nbtp = PffF ( m1 : m4 )*df;
TP = sum ( PffF ( m1 : m4 )*df) ;
vlfp = PffF((m1+1):m2)*df;
VLFP = sum (PffF((m1+1):m2)*df);
lfp = PffF((m2+1):m3)*df;
LFP = sum (PffF((m2+1):m3)*df);
hfp = PffF((m3+1):m4)*df;
HFP = sum (PffF((m3+1):m4)*df);
謝謝大家
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.125.31.64