|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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('相频图'); |
|