我对一组数据做了个振动加速度转换位移的分析，大神给看看
本帖最后由 牛小贱 于 2014-7-3 20:29 编辑数据为发动机1900rpm时的发动机加速度数据，按程序做完后结果如下所示，大神们给看看结果正确么？
程序如下：
fni=A;
i=(-1)^0.5;
sf=20480;
fmin=0.5;
fmax=10000;
c=1;
it=2;
x=fni';
n=length(x);
t=0:1/sf:(n-1)/sf;
nfft=2^nextpow2(n);
y=fft(x,nfft);
df=sf/nfft;
ni=round(fmin/df+1);
na=round(fmax/df+1);
dw=2*pi*df;
w1=0:dw:2*pi*0.5*sf;
w2=-2*pi*(0.5*sf-df):dw:-dw;
w=;
w=w.^it;
a=zeros(1,nfft);
a(2:nfft-1)=y(2:nfft-1)./w(2:nfft-1);
if it == 2
y=-a;
else
a1=imag(a); a2=real(a); y=a1-a2*i;
end
a=zeros(1,nfft);
a(ni:na)=y(ni:na);
a(nfft-na+2:nfft-ni+2)=y(nfft-na+2:nfft-ni+2);
y=ifft(a,nfft);
y=real(y(1:n))*c;
B=abs(y);
a=y;
subplot(2,1,1);
plot(t,a);
subplot(2,1,2);
plot(t,B);我的数据文件1670kb，稍微大了一点点，该怎么上传上来呢。
补充内容 (2014-7-3 15:44):
数据见4楼附件 发动机加速度数据，打错字了。 本帖最后由 westrongmc 于 2014-6-27 20:51 编辑
发动机1900rpm转速稳定的话，结构应该是比较稳态的强迫振动响应，
不应该出现这种位移波形，这个积分结果似有问题。
采了10几秒数据？不妨将数据截取一下，少发几秒的数据，例如1秒的数据看一下？
本帖最后由 卡卡卡帅 于 2014-7-3 15:43 编辑
westrongmc 发表于 2014-6-27 17:37
发动机1900rpm转速稳定的话，结构应该是比较稳态的强迫振动响应，
不应该出现这种位移波形，这个积分结果 ...
1秒数据见附件，您帮我看一下。。 本帖最后由 westrongmc 于 2014-7-3 23:43 编辑
卡卡卡帅 发表于 2014-7-3 15:40
1秒数据见附件，您帮我看一下。。
几个问题：
1. 采样频率是多少？看程序里面好像是20480Hz。
如果附件是1秒的数据，为何只有10240个点？
2. 加速度的单位是什么？ 是m/s^2吗？
3. 转速对吗？
如果按照20480Hz的采样率看，其主频在38Hz左右，也就是转速不是1900rpm，而是在2300rpm左右，故有此疑问。
请将全程的数据发给我，ksi@chinaksi.com谢谢！
----------------------------------------------------------------------
下面是按照20480Hz采样，加速度单位为m/s^2，对Adat进行二次积分的位移结果。
备注：数据的开始部分和结束部分有积分处理的端点效应。
本帖最后由 westrongmc 于 2014-7-26 22:33 编辑
对加速度信号进行二次积分处理之后，得到的位移结果如下所示
楼主提供的加速度数据转贴在下面，为压缩格式，解压后为.mat格式。
采样频率 = 20480Hz；单位 = m/s^2
不错的帖子 谢谢 {:{03}:} 本帖最后由 westrongmc 于 2014-7-27 21:22 编辑
westrongmc 发表于 2014-7-26 22:08
对加速度信号进行二次积分处理之后，得到的位移结果如下所示
楼主提供的加速度数据下载链接，无需积分。
采样频率=20480Hz，单位=m/s^2
http://www.chinaksi.com/download/5417/
westrongmc 发表于 2014-7-27 21:20
楼主提供的加速度数据下载链接，无需积分。
采样频率=20480Hz，单位=m/s^2
http://www.chinaksi.com/do ...
大神，能否请教下我上面贴的程序哪里出现问题了呢？
这么久，我还是没把这程序给弄好啊。
页:
[1]