大家来总结一下：我这种解微分方程程序对吗？
方程组见附件，解法如下：(1式)*(-a3)+(2式)*(a1) 就得y2的微分方程
同理 在得出y1的微分方程
再用ode45
令y1=x1; =x2; y2=x3; =x4则程序如下：
function dx=jiaofeng(t,x)dx=
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=;=ode45('jiaofeng',,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' => '=A*+B*f
其中A即为特徵矩阵 A=, 求eig(A)即可得固有频率？但我又见：wn=eig(K,M)那种正确呢？还有就是振动频率如何求？不好意思这么多问题，概念不清，望大家指点。
回复 楼主 likefeng 的帖子
结果不是你想要的是什么意思？最后面的请在matlab中输入 help eig，你就知道二者之间的差别了
页:
[1]