声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1076|回复: 2

[其他] 新手求助!

[复制链接]
发表于 2009-5-12 18:15 | 显示全部楼层 |阅读模式

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

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

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向量长度是不同的,那应该怎样呢?
回复
分享到:

使用道具 举报

 楼主| 发表于 2009-5-12 18:18 | 显示全部楼层
没人指导,苦啊,希望热心人的指教!
发表于 2009-5-13 09:06 | 显示全部楼层
从程序中看,LZ主要是把信号进行低通滤波,原始数据又较长,在这种情况下可以用低通滤波器,经滤波后x  y向量长度是相同的。

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 21:33 , Processed in 0.059159 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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