声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1380|回复: 3

[FFT] 急急急!!!

[复制链接]
发表于 2007-5-29 16:59 | 显示全部楼层 |阅读模式

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

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

x
通过FFT变换的功率谱怎样转换为幅值谱,下面是我的程序,转换后的幅值和理论值相差很大,而且相位也相差较大.下面是我的程序和数据,请帮忙指点,万分感激!!

close all
clear all
A=load('wave1_1.txt');
A=A(:,1);
ttt=mean(A);
for i=1:512
    A(i)=A(i)-ttt;
end
Y=A;
subplot(311);
plot(Y);
xlabel('t');
ylabel('振幅');
title('原始信号');
Y6=fft(Y,512);
Pyy=Y6.*conj(Y6)/512;
c=sqrt(2*0.16*Y6);
f=1000*(0:255)/160;
subplot(312);
plot(f,c(1:256));
%plot(f,Pyy(1:256));
xlabel('频率');
ylabel('幅值');
title('功率谱');
p=angle(Y6);
g=p(1:256);
subplot(313);
plot(f,180*g/pi);
xlabel('频率');
ylabel('角度');
title('相频图');
回复
分享到:

使用道具 举报

发表于 2007-5-29 18:28 | 显示全部楼层
幅值不是c=sqrt(2*0.16*Y6);
而是c=abs(Y6)/512*2;
发表于 2007-5-29 22:21 | 显示全部楼层
功率谱=(幅值谱)^2/2
 楼主| 发表于 2007-5-30 08:23 | 显示全部楼层
谢谢各位,我来试下
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-12 00:16 , Processed in 0.066633 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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