function =vbr_rk_ex(z,u,t)
% function for 2
% dx dx
% m --+ c -- +k x = f(t)
% 2 dt
% dt dx
% where m=1,k=1,c=.1, and z(1)=x, z(2)=--
% dt
zd=[z(2);
-.1*z(2)-z(1)+u(2)];
这一段还是没看懂,求解释。 t=0:.5:20; % Creates time vector
u=; % Creates force matrix
x0=; % Creates initial state vector.
x=vbr_rk('vbr_rk_ex',u,t,x0); % Runs analysis.
plot(t,x(1,:)); % Plots displacement versus time.
plot(t,x(2,:)); % Plots velocity versus time. 每天来这里的目的就是为了学习,期待着有一天也能为这个论坛做一点贡献。 for l1=1:(n-1);
z1=z(:,l1);
u1=u(:,l1);
u2=u(:,l1+1);
这里面的u2用到了么?作用是什么? 不错的帖子!{:{39}:}
页:
1
[2]