请教这个问题的求解是不是出了问题

猴王

这是一个微分方程的求解程序，其方程是这样的mv''+cv'+kv=f,我想求出这个方程的v''，请指教是不是求解出了什么问题：
function xdot=danziyoudukuangjia(t,x)
F=20*sin(25*t);
  M=31.2;
  C=0.01;
  K=67.51;
xdot=[x(2);(1/M)*(F-C*x(2)-K*x(1))];
执行文件：
clc;clear all;
close all;
t0=0;
   tf=60;
   x0=[0,0]';

   [t,x]=ode23('danziyoudukuangjia',[0:0.005:60],x0);
   d=x(:,1);
   v=x(:,2);
   
t=0:0.005:60;

F=20*sin(25*t);                                
M=31.2;
C=0.01;
K=67.51;
s1=(1/M)*(F-C*v.'-K*d.');

zhlong

楼主，这是你的程序，出了什么问题你要告诉大家呀。

猴王

我想请教我的解法有什么问题没有，步骤有问题吗，我的不到我师兄的结果，这让我很头疼，我想知道如果让你来解决这个问题，你会怎么解决，如果可以的话，我想请你也编译以下这个问题的求解，我想对比一下

zhlong

我运行的d，v ，s1结果如下

猴王

如果您是自行编译的这个程序，可以说咱们的结果是一模一样，十分感谢，如果您方便的话，请您加我的qq，我师兄的结果和咱们两个的都还是有差别的，有比较大的差别，我不知道为什么，请您给与知道，我的qq是284285603，或者您可以留下你的号码，十分感谢你的帮助

zhlong

这么简单的一个程序，我不知道为什么会和你师兄的结果有很大的差异。如果方便，把你师兄做的相关结果和资料贴上来看看。

猴王

我把师兄的论文中的加速度时程曲线放到了附件里面，请你看看，还有他的一些条件，我觉得和咱们两个的模拟的信号都有些差距，希望您能仔细的看看谢谢了

zhlong

我不能确定你的程序与你师兄文中的模型是否完全吻合。
但是只就你的程序来说，编程和求解一点问题没有。

后知后觉

我也没看出来程序有什么问题。
不过你可以对照你们两个的程序看看啊。

光看图是不好找问题的。

猴王

可以把你的程序贴出来比较一下吗

zhlong

我用的就是你贴的程序

猴王

呵呵，我的意思是，请您能不能根据我的这个条件编译一个程序，来求解mv''+cv'+kv=f这个方程，然后画出v''的图形，我的程序我明白没有问题，但现在是和师兄的有比较大的差别，我不知道问题出在什么地方，所以想请您帮忙看看

猴王

请你再看看我的那个附件，里面有这样的一句话，用ansys求得结构的前三阶频率，这是为什么，为什么要求前三阶频率阿

zhlong

不好意思，对你的这个问题我也不清楚，你可以到基础力学板块去问问。