加速度转换成位移的matlab代码及说明

加速度转换成位移的 matlab 代码及说明由测量的加速度离散数据数据转化成位移数据一般不直接在时域进行积分处理，而是由时域转换成频域在频域中进行二次积分再转化到时域中得到位移结果。相关 matlab 处理程序方法参考王济老师的 matlab 在振动信号处理中的应用中的程序如下：％频域积分％%％％%%％%%%%%％％%％%%%%%%％%％％clear; clc; close all hidden%％％%%%%％%％％%%％%％％%fni=input(’频域积分—输入数据文件名：',’s’);fid=fopen(fni，'r'）;sf=fscanf(fid,’％f',1）;%采样频率fmin=fscanf(fid，’%f’,1);%最小截止频率fmax=fscanf(fid,'%f’,1)；％最大截止频率c=fscanf(fid，'%f’,1)；％单位变换系数it=fscanf(fid，'%f',1）；%积分次数sx=fscanf(fid，'％s',1)；%横向坐标轴的标注sy1=fscanf（fid，'％s’,1);％纵向坐标轴输入单位的标注sy2=fscanf（fid，’％s'，1)；％纵向坐标轴输出单位的标注fno=fscanf(fid,’％s'，1)；％输出数据文件名x=fscanf（fid,'%f’，［1,inf])；％输入数据存成行向量status=fclose（fid）；n=length(x);%建立时间向量t=0：1/sf：（n—1）/sf；％大于并最接近 n 的 2 的幂次方为 FFT 长度nfft=2^nextpow2（n）;％FFT 变换y=fft（x,nfft）;%计算频率间隔（Hz/s）df=sf/nfft;％计算指定频带对应频率数组的下标ni=round(fmin/df+1)；na=round（fmax/df+1）；%计算圆频率间隔（rad/s）dw=2*pi*df；％建立正的离散圆频率向量w1=0：dw:2＊pi＊(0.5＊sf—df)；%建立负的离散圆频率向量w2=2*pi＊（0.5＊sf-df)：—dw：0；%将正负圆频率向量组合成一个向量w=[w1，w2];%以积分次数为指数，建立圆频率变量向量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; ％进行一次积分的相位变换enda=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）; %IFFT 变换%取逆变换的实部 n 个元素并乘以单位变换系数为积分结果y=real(y（1：n)）*c;subplot(2,1，1); plot（t,x）; xlabel（sx)； ylabel（sy1); grid on; %绘制几分钱的时程曲线图形subplot(2,1,2）； plot(t,y）; xlabel(sx）； ylabel(sy2); ...