声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2060|回复: 7

[编程技巧] [求助]大家帮我看看这个程序(功率谱计算)

[复制链接]
发表于 2006-6-3 15:48 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我做了一个信号的功率谱分析,可是得出的图像表示和预期结果不一样,以至于后面没法分析,大家能帮我修改下么?(数据和应得的结果图已经上传)
程序:
  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;
复制代码
1.bmp

心血管功能实验测试数据报告.txt

3.81 KB, 下载次数: 9

回复
分享到:

使用道具 举报

 楼主| 发表于 2006-6-3 21:01 | 显示全部楼层
<P>我的这种功率谱估计算法有问题吗?请各位帮帮忙</P>
 楼主| 发表于 2006-6-3 23:13 | 显示全部楼层
所得图形中,后面的是若干谐波,前面有一个略低的波峰,
可是我用了几种功率谱分析都得不到接近的图形……
后面的分析就完全进行不下去了……
发表于 2006-6-3 23:17 | 显示全部楼层
<P>为什么求功率谱还用到自相关系数?</P>
 楼主| 发表于 2006-6-3 23:44 | 显示全部楼层
功率谱估计有很多种,这是经典功率谱估计的一种(BT法),还可以用其他的。
 楼主| 发表于 2006-6-4 12:02 | 显示全部楼层
拜托啦~在线等待啊!!
发表于 2006-6-30 16:19 | 显示全部楼层
倒数第四行k=0:floor(nn/2);是n呢?还是n*n?<BR>
发表于 2006-6-30 16:22 | 显示全部楼层
<P>我用n计算,出来的结果和你的也不一样啊?</P>
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-12 19:40 , Processed in 0.065691 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表