功率谱程序中的一些问题（红色部分）
本帖最后由 mzbys 于 2015-12-14 16:05 编辑matlab中 = periodogram(x,window,nfft,fs)periodogram是用来计算功率谱密度的，参数中，X：所求功率谱密度的信号；window：所使用的窗口，默认是boxcar，其长度必须与x的长度一致；nfft：采样点数；fs：采样频率。例如：Fs=1000;n=0:1/Fs:1;x=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));nfft=1024;window=boxcar(length(n));=periodogram(x,window,nfft,Fs);P=10*log10(Pxx);plot(f,P);hold on;Pxx_1=abs(fft(x,nfft)).^2/length(n);t=0:round(nfft/2-1);f=t*Fs/nfft;%这里的nfft是频域的取样点数，和时域的取样点数N有什么区别吗，可以互换吗P_1=10*log10(Pxx_1(t+1));%这里t+1是什么意思plot(f,P_1,'r');legend('periodogram','公式');title('周期图法功率谱估计');
同问，请问下谁有好用的功率谱分析程序呢？ nfft是频域点数，和N是时域的取样点数，不可以互换
t+1是为了将index调整重1开始(比较下t=0:round(nfft/2-1)) 楼主搞定了？
页:
[1]