Question

帮忙看一下MATLAB程序 那里错了 拜托大神

m1 = 1814; m2 = 2903;
k1 = 52.5; k2 = 140.1;
F0= 360; w = 28;
C1 = 52.5; C2 = 700.5;
% 微分方程
A=[k1 -k1;-k1 k1+k2];
dx = @(t,x)[x(3:4); (-A*x(1:2)-x(3:4).*[C1 -C1;-C1;C1+C2]+[0;F0*sin

(w*t)])./[m1;m2]];
x0=zeros(4,1)
ode45(dx,[0 5],x0)

Answer 1

三处错误：

1、[C1 -C1;-C1;C1+C2] 不满足矩阵拼接的要求；

2、上述C矩阵与x_dot的关系不应该是点乘，而应该是矩阵右乘；

3、sin函数后面不应该换行。

 

改后的代码：

m1 = 1814; m2 = 2903; 
k1 = 52.5; k2 = 140.1;
F0= 360; w = 28;
C1 = 52.5; C2 = 700.5;
% 微分方程
A=[k1 -k1;-k1 k1+k2];
dx = @(t,x)[x(3:4); (-A*x(1:2)-[C1 -C1;-C1 C1+C2]*x(3:4)+[0;F0*sin(w*t)])./[m1;m2]];
x0=zeros(4,1);
ode45(dx,[0 5],x0)

追问

非常感谢啊大神 ，多谢帮忙。可以加好友吗？想以后不会的请教

追答

可以在我的名字上面悬停鼠标，然后按【求助知友】，向我定向求助。