微分方程的求解程序问题
这是一个微分方程的求解程序，其方程是这样的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=;
执行文件：
clc;clear all;
close all;
t0=0;
tf=60;
x0=';
=ode23('danziyoudukuangjia',,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.');
[ 本帖最后由 eight 于 2007-6-25 10:18 编辑 ] 程序本身没有什么问题. 这个结果与我师兄的结果之间还有差距，我想知道如果我通过对v(:,1)进行求导，这样通过对速度曲线求导从而得到加速度曲线的方式是否可行，应该怎样通过程序实现，请指教，请帮我写出程序的框架好吗，还有就是程序我知道没有问题，求解过程您看有没有问题 这个自己看看公式不就明白了?
页:
[1]