帮忙看一下MATLAB程序 那里错了 拜托大神
m1=1814;m2=2903;k1=52.5;k2=140.1;F0=360;w=28;C1=52.5;C2=700.5;%微分方程A=[k1-k1;-k1k1+k2]...
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) 展开
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) 展开
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)
追问
非常感谢啊大神 ，多谢帮忙。可以加好友吗？想以后不会的请教
追答
可以在我的名字上面悬停鼠标，然后按【求助知友】，向我定向求助。
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
AiPPT
2024-09-19 广告
随着AI技术的飞速发展，如今市面上涌现了许多实用易操作的AI生成工具1、简介：AiPPT: 这款AI工具智能理解用户输入的主题，提供“AI智能生成”和“导入本地大纲”的选项，生成的PPT内容丰富多样，可自由编辑和添加元素，图表类型包括柱状图...
点击进入详情页
本回答由AiPPT提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
