|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
是这样,我有一列很大的数据,大概27万,我想从其中去4096个点来分析,用y=fft(x,4096)来做,但是横坐标是点数了,怎么把横坐标改为时间呢?我知道我这个程序是错误的,我是初学,希望各位指教!
p=load('f:\学习\mlab\xiangying.txt');
sf=1024;
fmin=0;
fmax=80;
x=p(:,2);
n=4096;
t=(0:1/sf:(n-1)/sf)';
nfft=4096;
ni=round(fmin*nfft/sf+1);
na=round(fmax*nfft/sf+1);
y=fft(x,nfft);
a=zeros(1,nfft);
a(ni:na)=y(ni:na);
a(nfft-na+1:nfft-ni+1)=y(nfft-na+1:nfft-ni+1);
y=ifft(a,nfft);
y=(real(y(1:n)))';
subplot(2,1,1);
plot(t,x);
xlabel('时间(s)');
ylabel('加速度(g)');
grid on;
subplot(2,1,2);
plot(t,y);
xlabel('时间(s)');
ylabel('加速度(g)');
grid on;
我知道这样,x y向量长度是不同的,那应该怎样呢? |
|