大家来总结一下：我这种解微分方程程序对吗？

likefeng

方程组见附件，解法如下：

(1式)*(-a3)+(2式)*(a1) 就得y2的微分方程
同理 在得出y1的微分方程
再用ode45
令y1=x1; =x2; y2=x3; =x4则程序如下：

function dx=jiaofeng(t,x)

dx=[x(2);

0.366663*(0.36972*x(4)-1.524*x(2)+4.009*x(3)-12.1949*x(1));

x(4);

-0.73322*(-0.05152*x()-10533069*x(3)+157.876*x(1))];

x0=[0;0.2;0;0.3];

[t,x]=ode45('jiaofeng',[0,40],x0);

subplot(2,2,1);

plot(t,x(:,1));

hold on;

subplot(2,2,2);

plot(t,x(:,2));

hold on;

subplot(2,2,3);

plot(t,x(:,3));

hold on;

subplot(2,2,4);

plot(t,x(:,4))

但结果并不是我想要的。

另外，

运动方程Mx"+Cx'+Kx=f, 设q1=x, q2=x' => [q1; q2]'=A*[q1;q2]+B*f
其中A即为特徵矩阵 A=[zeros(n,n) ones(n,n); -inv(M)*K, -inv(M)*C], 求eig(A)即可得固有频率？

但我又见：wn=eig(K,M)

那种正确呢？

还有就是振动频率如何求？

不好意思这么多问题，概念不清，望大家指点。

无水1324

结果不是你想要的是什么意思？
最后面的请在matlab中输入 help eig，你就知道二者之间的差别了