请教高手

猴王

clc;clear all;
close all;
t0=0;
   tf=60;
   x0=[0,0,0,0,0,0]';
   %tol=0.0002;
   %trace=0;
   [t,x]=ode23('jianxiehezai',[0:0.0001:60],x0);
   d1=x(:,1);
   v1=x(:,4);
    d2=x(:,2);
   v2=x(:,5);
    d3=x(:,3);
   v3=x(:,6);
t=0:0.0001:60;
%ft=linspace(0,60,60000);%将0到60秒，分成60000份
F=[20*sin((2/25)*pi.*t);20*sin((2/25)*pi.*t);20*sin((2/25)*pi.*t)];

  m1=31.2;
  m2=31.2;
  m3=31.2;
  k1=67.51;
  k2=89.856;
  k3=89.856;
  c1=0.01;
  c2=0.01;
  c3=0.01;
M=[ m1+m2,-m2,0;-m2,m2+m3,-m3;0,-m3,m3];
C=[ c1+c2,-c2,0;-c2,c2+c3,-c3;0,-c3,c3];
K=[k1+k2,-k2,0;-k2,k2+k3,-k3;0,-k3,k3];
c=inv(M)*(F-C*[v1,v2,v3].'-K*[d1,d2,d3].');
xlabel('Time-sec');
ylabel('jiasudu');
plot(t,c)
这个程序是没有问题的，我现在就想问问关于采样频率的定义，和怎么用，用到什么地方，请各位高手帮助

eight

原帖由 猴王 于 2007-5-15 09:53 发表

                               
登录/注册后可看大图

clc;clear all;
close all;
t0=0;
   tf=60;
   x0=[0,0,0,0,0,0]';
   %tol=0.0002;
   %trace=0;
   [t,x]=ode23('jianxiehezai',[0:0.0001:60],x0);
   d1=x(:,1);
   v1=x(:,4);
    d2=x(:,2);
...

请移步信号处理方法版块