matlab加速度转化为位移-matlab数值积分实现加速度、速度、位移的转换（时域频域积。。。.pdfVIP

matlab加速度转化为位移,matlab数值积分实现加速度、速

度、位移的转换（时域频域积。。。

对此可以⽤滤波的⽅法⼤的趋势项去掉。

测试的代码如下

%测试积分对正弦信号的作⽤

clc

clear

closeall

%%原始正弦信号

ts=0.001;

fs=1/ts;

t=0:ts:1000*ts;

f=50;

dis=sin(2*pi*f*t);%位移

vel=2*pi*f.*cos(2*pi*f*t);%速度

acc=-(2*pi*f).^2.*sin(2*pi*f*t);%加速度

%多个正弦波的测试

%f1=400;

%dis1=sin(2*pi*f1*t);%位移

%vel1=2*pi*f1.*cos(2*pi*f1*t);%速度

%acc1=-(2*pi*f1).^2.*sin(2*pi*f1*t);%加速度

%dis=dis+dis1;

%vel=vel+vel1;

%acc=acc+acc1;

%

结：频域积分正常恢复信号，时域积分恢复加⼊的⾼频信息有误差

%加噪声测试

acc=acc+(2*pi*f).^2*0.2*randn(size(acc));

%结：噪声会使积分结果产⽣⼤的趋势项

figure

ax(1)=subplot(311);

plot(t,dis),title(位移)

ax(2)=subplot(312);

plot(t,vel),title(速度)

ax(3)=subplot(313);

plot(t,acc),title(加速度)

linkaxes(ax,x);

%由加速度信号积分算位移

[disint,velint]=IntFcn(acc,t,ts,2);

axes(ax(2));holdon

plot(t,velint,r),legend({原始信号,

恢复信号})

axes(ax(1));hold

on

plot(t,disint,r),legend({原始信号,恢复信号})

%%测试积分算⼦的频率特性

n=30;

amp=zeros(n,1);

f=[5:3040:10:480];

figure

fori=1:length(f)

fi=f(i);

acc=-(2*pi*fi).^2.*sin(2*pi*fi*t);%加速度

[disint,velint]=IntFcn(acc,t,ts,2);%积分算位移

amp(i)=sqrt(sum(disint.^2))/sqrt(sum(dis.^2));

plot(t,disint)

drawnow

%pause

end

close

figure

plot(f,amp)

title(位移积分的频率特性曲线)

xlabel(f)

ylabel(单位正弦波的积分位移幅值)

以上代码中使⽤IntFcn函数实现积分，它是封装之后的函数，可以实现时域积分和频域积分，其代码如下

%

积分操作由加速度求位移，可选时域积分和频域积分

function[disint,velint]=IntFcn(acc,t,ts,flag)

ifflag==1

%时域积分

[disint,velint]=IntFcn_Time(t,acc);

velenergy=sqrt(sum(velint.^2));

velint=detrend(velint);

velreenergy=sqrt(sum(velint.^2));

velint=

velint/velreenergy*velenergy;disenergy=sqrt(sum(disint.^2));

disint=detrend(disint);

disreenergy=sqrt(sum(disint.^2));

disint=disint/disreenergy*disenergy;%

此操作是为了弥补去趋势时能量的损失

%去除位移中的⼆次项

p=polyfit(t,disint,2);

disint=disint-polyval(p,t);

else

%频域积分

velin