Differentiation of the position vector.

2 次查看(过去 30 天)
I have written a code for postion of a vibrating system. I am unable to determine the velocity and acceleration from the expression u. I want to diffrentiate u(i) in the same loop and plot it with t. How this can be done? Please Help.
m1 = 40;
m = m1*1000/9.81;
k = 3500*1000;
u(1) = 0.01; v(1) = 0.1;
wn = sqrt(k/m);
i=1;
for t=0:.001:1
u(i) = u(1)*cos(wn*t) + (v(1)/wn)*sin(wn*t);
i=i+1;
end
t= 0:0.001:1;
plot(t,u);
grid on

采纳的回答

fred  ssemwogerere
Hello, you could make use of the symbolic math toolbox. The page below should get you started:

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Acoustics, Noise and Vibration 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by