Runge-Kutta实现多自由度系统响应的MATLAB程序

tywww123

谢谢，看看，学学

bauyo

受教了哦

kukumalu180614

其中rkf为动力微分方程，形如
function [zd]=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)];

这一段还是没看懂，求解释。

feizhongxiu

1. t=0:.5:20;                             % Creates time vector
   
2. u=[zeros(1,length(t));sin(t*1.1)];     % Creates force matrix
   
3. x0=[1 ;0];                             % Creates initial state vector.
   
4. x=vbr_rk('vbr_rk_ex',u,t,x0);           % Runs analysis.
   
5. plot(t,x(1,:));                        % Plots displacement versus time.
   
6. plot(t,x(2,:));                        % Plots velocity versus time.

复制代码

yap988

每天来这里的目的就是为了学习，期待着有一天也能为这个论坛做一点贡献。

lihaitao123

for l1=1:(n-1);
   z1=z(:,l1);
   u1=u(:,l1);
   u2=u(:,l1+1);

这里面的u2用到了么？作用是什么?

402144999

不错的帖子！