| 查看: 1008 | 回复: 0 | ||
[求助]
用ode45求解矩阵微分方程，出现如下问题，求助各位大神，感激不尽
|
|
m文件如下： clc;clear; M=diag([1,2,3]); K=[9,-4,0;-5,11,-6;0,-6,13]; C=diag([1,4,7]); dx = @(t,x)[x(4:6); inv(M)*(-K*x(1:3)-C*x(4:6)+[sin(t);0;0])]; x0 = zeros(6,1); [t,x]=ode45(dx,[0 10],x0); 这样写可以运行。 但是如果将变量t写在前面： clc;clear; syms t M=diag([1,2,3]); K=[9,-4,0;-5,11,-6;0,-6,13]; C=diag([1,4,7]); F=[sin(t);0;0]; dx = @(t,x)[x(4:6); inv(M)*(-K*x(1:3)-C*x(4:6)+F)]; x0 = zeros(6,1); [t,x]=ode45(dx,[0 10],x0); 就会报错： Error using odearguments (line 111) Inputs must be floats, namely single or double. Error in ode45 (line 114) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in Untitled4 (line 9) [t,x]=ode45(dx,[0 10],x0); >> |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
现在用MATLAB求解一个循环后画出曲线遇到问题，求大神指点！！
已经有7人回复- 关于根据一个系统模型建立复杂微分方程解的问题，一般都是设其解为某种形式。
已经有3人回复
- matlab求解，求大神
已经有1人回复
- 微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- 这种偏微分方程如何求解啊？求指点
已经有20人回复
- 关于二阶偏微分方程求解
已经有21人回复
- 请教一下各位大神怎么求解 Runge–Kutta method的微分方程的各个系数
已经有12人回复
- 病态方程组求解问题，望指导啊
已经有7人回复
- 关于打靶法解二阶微分方程的数值解的问题
已经有6人回复
- 偏微分方程组的求解
已经有3人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- 各位大神，请帮我看看题目解法对不对
已经有6人回复
- 求大神看一下我这个，感觉是软件设置的问题
已经有7人回复
- 微分方程求解（通解中含hypergeom函数）
已经有17人回复
- 求微分方程 最大值问题
已经有4人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 求助利用matlab的ODE45求解微分方程
已经有16人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 求助：控制器u随状态变化而切换时的ode45求解
已经有3人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人