我对一组数据做了个振动加速度转换位移的分析，大神给看看

卡卡卡帅

数据为发动机1900rpm时的发动机加速度数据，按程序做完后结果如下所示，大神们给看看结果正确么？
程序如下：

1. fni=A;
   
2. i=(-1)^0.5;
   
3. sf=20480;                          
   
4. fmin=0.5;                        
   
5. fmax=10000;                        
   
6. c=1;                             
   
7. it=2;                             
   
8. x=fni';                             
   
9. n=length(x);
   
10. 
    
11. t=0:1/sf:(n-1)/sf;                 
    
12. nfft=2^nextpow2(n);                 
    
13. y=fft(x,nfft);                     
    
14. df=sf/nfft;                        
    
15. ni=round(fmin/df+1);               
    
16. na=round(fmax/df+1);               
    
17. dw=2*pi*df;                        
    
18. w1=0:dw:2*pi*0.5*sf;           
    
19. w2=-2*pi*(0.5*sf-df):dw:-dw;         
    
20. w=[w1,w2];                          
    
21. w=w.^it;                           
    
22. a=zeros(1,nfft);
    
23. a(2:nfft-1)=y(2:nfft-1)./w(2:nfft-1);        
    
24. if it == 2
    
25.    y=-a;                              
    
26. else
    
27.    a1=imag(a); a2=real(a); y=a1-a2*i;
    
28. end
    
29. a=zeros(1,nfft);
    
30. a(ni:na)=y(ni:na);                     
    
31. 
    
32. a(nfft-na+2:nfft-ni+2)=y(nfft-na+2:nfft-ni+2);            
    
33. 
    
34. y=ifft(a,nfft);
    
35. y=real(y(1:n))*c;
    
36. B=abs(y);
    
37. a=y;
    
38. subplot(2,1,1);
    
39. plot(t,a);
    
40. subplot(2,1,2);
    
41. plot(t,B);

复制代码

我的数据文件1670kb，稍微大了一点点，该怎么上传上来呢。

补充内容 (2014-7-3 15:44):
数据见4楼附件

卡卡卡帅

发动机加速度数据，打错字了。

westrongmc

发动机1900rpm转速稳定的话，结构应该是比较稳态的强迫振动响应，
不应该出现这种位移波形，这个积分结果似有问题。

采了10几秒数据？不妨将数据截取一下，少发几秒的数据，例如1秒的数据看一下？

卡卡卡帅

westrongmc 发表于 2014-6-27 17:37
发动机1900rpm转速稳定的话，结构应该是比较稳态的强迫振动响应，
不应该出现这种位移波形，这个积分结果 ...

1秒数据见附件，您帮我看一下。。

westrongmc

卡卡卡帅 发表于 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进行二次积分的位移结果。

KSI科尚 - Adat二次积分结果

备注：数据的开始部分和结束部分有积分处理的端点效应。

westrongmc

对加速度信号进行二次积分处理之后，得到的位移结果如下所示

KSI科尚-加速度二次积分处理的位移结果

楼主提供的加速度数据转贴在下面，为压缩格式，解压后为.mat格式。
采样频率 = 20480Hz；单位 = m/s^2

y.rar (1.64 MB, 下载次数: 3)

2014-7-26 22:16 上传

点击文件名下载附件
数据贡献者为楼主，解压后为.mat格式

mywl

不错的帖子 谢谢

westrongmc

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 ...

大神，能否请教下我上面贴的程序哪里出现问题了呢？
这么久，我还是没把这程序给弄好啊。