声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1005|回复: 1

[FFT] 大家请看看,提点意见。附代码

[复制链接]
发表于 2009-12-17 11:41 | 显示全部楼层 |阅读模式

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

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

x
clc;
%预定参数;
node=5;
filepath='D:\history\';
msgbox('正在计算,请稍候...','运行程序','help','replace');
for i=1:node  
   A=load([filepath 'X-' num2str(i) '.txt']);
   eval(['P',num2str(i),'=','A(:,2)',';']);%分别保存在P1,P2...
   eval(['p_his','(:,i)','=','P',num2str(i),';']);%%
end %这里只是将信号组放在一个矩阵里面

msgbox('正在计算!','运行程序','help','replace');
for i=1:5
Fs=1/0.18; %采样频率
n=0:1/Fs:2952;
nfft=16400;
window2=hanning(16400); %汉宁窗
noverlap=20; %数据无重叠
range='half'; %频率间隔为[0 Fs/2],只计算一半的频率
[sp,f]=pwelch(p_his(:,i),window2,noverlap,nfft,Fs,range);
spp(:,i)=sp;ff(:,i)=f;

plot_spp=log10(abs(spp(:,i)));
figure(1);
grid on;
subplot(3,2,i);
plot(log10(abs(ff(:,i))),plot_spp);
end

根据上面的方法我求的的功率谱打出的图如下,怎么有负值啊?曲线好像是合理的,但是我要的是正值啊,幅值也不知道对不对,没有把握。
请大侠答疑啊。谢谢。
sxx.jpg
回复
分享到:

使用道具 举报

发表于 2009-12-17 19:38 | 显示全部楼层
取了对数,数值在1以下便为负值。

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 18:25 , Processed in 0.084065 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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