帮忙看一下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、[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 广告
2024-09-19 广告
随着AI技术的飞速发展,如今市面上涌现了许多实用易操作的AI生成工具1、简介:AiPPT: 这款AI工具智能理解用户输入的主题,提供“AI智能生成”和“导入本地大纲”的选项,生成的PPT内容丰富多样,可自由编辑和添加元素,图表类型包括柱状图...
点击进入详情页
本回答由AiPPT提供
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询