[求助]大家帮我看看这个程序（功率谱计算）

flyingbone

我做了一个信号的功率谱分析，可是得出的图像表示和预期结果不一样，以至于后面没法分析，大家能帮我修改下么？（数据和应得的结果图已经上传）
程序：

1. fidin=fopen('8.txt');
   
2. fidout=fopen('mkmatlab.txt','w');
   
3. while ~feof(fidin)
   
4. tline=fgetl(fidin);
   
5. if double(tline(1))>=48&&double(tline(1))<=57
   
6. fprintf(fidout,'%s\n\n',tline);
   
7. continue
   
8. end
   
9. end
   
10. fclose(fidout);
    
11. a=importdata('MKMATLAB.txt');
    
12. %%
    
13. nfft=512;
    
14. fs=181.8;
    
15. %%
    
16. x=255-a;
    
17. x=x-mean(x);
    
18. r=xcorr(x,'biased');
    
19. n=length(r);
    
20. py=abs(fft(r)/n);
    
21. k=0:floor(nn/2);
    
22. py=10*log10(py(k+1));
    
23. plot(k/n,py);
    
24. grid;

复制代码

flyingbone

<P>我的这种功率谱估计算法有问题吗？请各位帮帮忙</P>

flyingbone

所得图形中，后面的是若干谐波，前面有一个略低的波峰，
可是我用了几种功率谱分析都得不到接近的图形……
后面的分析就完全进行不下去了……

zhao197642

<P>为什么求功率谱还用到自相关系数?</P>

flyingbone

功率谱估计有很多种，这是经典功率谱估计的一种（BT法），还可以用其他的。

flyingbone

拜托啦～在线等待啊！！

bettylong

倒数第四行k=0:floor(nn/2);是n呢？还是n*n?<BR>

bettylong

<P>我用n计算，出来的结果和你的也不一样啊？</P>