【求助】求助FFT或功率谱分析

[zhuyjabc]

小弟最近采集到一组数据，但是画不出FFT或功率谱图，我的采集频率为25000HZ，请哪位高手指教一下 数据在附件里
我写的程序如下：
load tt.txt;
y=fft(tt,1024);
p=y.*conj(y)/1024;
ff=10000*(0:511)/1024;
plot(ff,p(1:512))

[foxstar]

本人也是初学者哈,我觉得可能是你的数据本来就没有对吧.但是我这有个程序可以画出来我给的数据信息.

A=2; %幅值
fs=100; %采样频率
f_singal=10; %信号频率
t=-1:1/fs:1; %信号步进
x=A*sin(2*pi*f_singal*t); %产生数据
save sinfile.txt x -ASCII

%画图
>> close all;
load sinfile.txt;
fs=100;
N=length(sinfile); %N样点的个数
fx=fft(sinfile);
df=fs/N; %df频率分辨率
n=0:N/2; %
f=n*df; %频率轴上横坐标
subplot(121);
plot(sinfile);
subplot(122);
plot(f,abs(fx(n+1)));

[gzqfast]

程序如下,你可以试一试


load tt.txt
Fs=25000;
num=1024;%FFT点数
f=Fs*(0:num/2-1)/num;
Y=abs(fft(tt,num));
plot(f,Y(1:num/2));

[limingaini119]

受用了，，感谢

[klant]

受益匪浅啊，感动的流泪了