帮忙看一下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)
展开
 我来答
tianxiawulang
2016-10-22 · TA获得超过2.7万个赞
知道大有可为答主
回答量:4732
采纳率:89%
帮助的人:2650万
展开全部

三处错误:

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提供
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式