加速度转换成位移的matlab代码及说明.doc
需积分: 27 187 浏览量
2022-10-19
23:34:15
上传
评论 1
收藏 47KB DOC 举报
加速度转换成位移的 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];
%以积分次数为指数，建立圆频率变量向量
评论0