MATLAB中如何將離散加速度積分求得速度

2021-05-11 15:26:38 字數 1403 閱讀 5111

1樓:昂徠教育科技****

方法一:

e=load('c:\users\administrator\desktop\in.txt');

d=e-sum(e(:,1))/length(e);

figure(1)

plot(d)

title('濾波前加速度')

%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 積分%%%%%%%%%%%%

dt=0.1;

sx=0;

for i=1:1:length(d)

sx=d(i,1)*dt+sx ;%從最值位置疊加加速度

vx(i,1)=sx; %速度矩陣

endfigure(2)

subplot(2,1,1)

plot(vx)

ylabel('速度')

%%% 直流分量確定和積分基線的標定 %%%%%%%%%%%%%%%%%%%%%

vx=vx-sum(vx(:,1))/length(vx);

% sx2=dt*cumtrapz(vx-0.04) %x方向位移矩陣,其中0.04為通過整週期積分得出的直流分量值。

sx2=dt*cumtrapz(vx); %x方向位移矩陣,其中vx為通過標定後的速度矩陣。

subplot(2,1,2)

plot(sx2)

ylabel('位移')

方法二:

將下面內容儲存為simpson.m

function y= simpson(f,a,b,m)

% f被積函式;a積分下限;b積分上限;m子區間個數(將x分為多少個區間)

if m==0

y=0;

else

h=(b-a)/(2*m);

s1=0;

s2=0;

for i=1:m

x=a+(2*i-1)*h;

s1=s1+feval(f,x);

endfor j=1:(m-1)

x=a+2*j*h;

s2=s2+feval(f,x);

endy=h/3*(feval(f,a)+2*s2+4*s1+feval(f,b));

end在命令列輸入

f=@(t)sin(t);

i=;m=floor(2*pi/0.01);

for i=0:m

y=simpson(f,0,2*pi/m*i,i);

i=[i,y];

endt=0:0.01:2*pi;

plot(t,i);

matlab 加速度積分成速度及位移 請神人幫忙

2樓:匿名使用者

你給我公式,我不可能知道你們的公式吧。還有你給的財富積分太低了這個可不是一般的程式設計呀

如何判斷加速度方向的正負,加速度怎麼判斷正負

看前一秒與後一秒的速度哪個大,後一秒大的話是加速為正,後一秒小的話減速為負 加速度怎麼判斷正負 看初速度和末速度的大小,初速度大加速度為負,初速度小加速度為正 怎麼看簡諧運動的速度加速度的大小及正負 做簡諧振動的物體 1 當物體向平衡位置運動時,速度增大 加速度減小。位移是正值,加速度 速度是負值 ...

加速度的方向如何判斷,高中物理加速度方向和速度方向怎麼判斷的

加速度的方向 與速度變化量 v的方向相同。與合外力的方向相同。當加速度方向與速度方向相同時,物體做加速運動 當加速度方向與速度方向相反時,物體做減速運動。也就是 加速運動 v 0,0 或者v 0,0 減速運動 v 0,0 或者v 0,0 加速度 acceleration 是速度變化量與發生這一變化所...

如何在ansys中加角加速度,ansys中施加角速度

gui方式 du方zhi 式 daodomega,domgx,domgy,domgz ansys中施加角速度 請問樓主的問題解決了嗎?我現在遇到和你同樣的問題。轉速不會施加 使用omega 命令可以施加角速度 其格式是 omega,omegx,omegy,omegz,kspin 如果對引數 omeg...