帮忙看看这个功率谱分析的程序对不对~

Michelle25

% let s is signal
len=256;% adijust accoding to nquist rules
F=(1:len/2)/len;
fs=fft(s,len);
p=fs.*conj(fs)
p=p(1:len/2);
plot(F,p)

周期法是不是这样??

signalgood

% let s is signal
len=256;% adijust accoding to nquist rules
F=(0:len/2-1)/len;
fs=fft(s,len);
p=fs.*conj(fs)
p=p(1:len/2)*2*2/N/N;
plot(F,p)


我认为应该是这样

Michelle25

% let s is signal
len=256;% adijust accoding to nquist rules
F=(1:len/2)/len;
fs=fft(s,len);
p=abs(fs);
p=p.^2/length(p);
plot(F,p)

定义是这样的??????

Michelle25

pwelch(s,[],[],[],Fs)
是不是就可以得到welch功率谱了,中间省略的[][][]分别是什么呢?

Michelle25

:@Q  ~~~~~~~~~~~~

songzy41

用help pwelch看看各参数代表什么
[Pxx,F] = PWELCH(X,WINDOW,NOVERLAP,NFFT,Fs)

tujian510

% let s is signal
len=256;% adijust accoding to nquist rules
F=(1:len/2)/len;   %是不是应该这样：F=fs*(1:len/2)/len,fs是采样频率，这样才能得到相应采样点的频率轴（横轴坐标）
fs=fft(s,len);
p=fs.*conj(fs)
p=p(1:len/2);
plot(F,p)

zh19841004

:@) 以正弦信号为例来说一下吧：
n=1024;fs=10;/要制定采样频率和采样点才行的
t=(1:n)/fs;
yn=sin(t);/得到离散点
pxx=fft(yn,n);
f=fs*(0:n/2-1)/n;
p=abs(pxx);
p=fs*p.^2/n;
plot(f,p(1:n/2));
另外这并不是品均周期图法，要用PSD或PWELCH函数才行。
你这里也没进加窗处理呀；
wh=hanning(n);
yn=yn.*wh;