声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1397|回复: 2

[FFT] 求助:频谱和功率谱的幅值差别为什么这么大?

[复制链接]
发表于 2013-9-6 20:39 | 显示全部楼层 |阅读模式

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

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

x
小弟近来很疑惑,作出的频谱和功率谱在幅值上差距很大,不知道是什么,肯定高手解惑 频谱.bmp

功率谱

功率谱

下面是程序:
y0=[0;0];Omega=0.8883;T1=2*pi/Omega;T=round(T1);tt=T1/100;max=500*T1;
[t,x]=ode45('vdp',[0:tt:max],y0);
y=x(:,1);dy=x(:,2);

h=x(round(((max/tt)+1)/2):end,1);nfft=2^nextpow2(length(h));hh=2*(abs(fft(h,nfft)))/length(h);fs=2*pi/tt;f=(0:(nfft-1))/nfft*fs;%求频率分布
figure;plot(f,hh);xlim([0,6]);xlabel('频率');ylabel('幅值');title('频谱');
h=x(round(((max/tt)+1)/2):end,1);nfft=2^nextpow2(length(h));hh=2*(abs(fft(h,nfft)).^2)/length(h);fs=2*pi/tt;f=(0:(nfft-1))/nfft*fs;%求频率分布
figure;plot(f,hh);xlim([0,6]);xlabel('频率');ylabel('幅值');title('功率谱')

function fy=vdp(t,x);
a=0.2;b=0.888523;c=0.394737;Omega=0.8883;g=0.2;
fy=[x(2);-a*x(2)-b^2*x(1)-3*c*x(1)^2-c*x(1)^3-g*cos(Omega*t)*x(1)-g*cos(Omega*t)];
end

回复
分享到:

使用道具 举报

发表于 2013-9-7 07:28 | 显示全部楼层
频谱和功率谱是有不同的因次,数量级不能比较。
有的地方认为功率谱和功率谱密度是一回事,有人认为不一样,我只认功率谱密度。
另外你这种波形,适用频谱法分析,不适用功率谱密度法,因为它不是随机信号。功率谱密度是定义在随机信号,而不是正弦或周期信号的,后者功率谱密度可能达到无穷大。

评分

1

查看全部评分

 楼主| 发表于 2013-9-7 08:47 | 显示全部楼层

谢谢主任的解释,这下明白了。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-5-19 18:53 , Processed in 0.062312 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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