| 查看: 1028 | 回复: 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); >> |
回复此楼» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有18人回复
-
面上可以超过30页吧？
已经有7人回复
-
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有5人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有6人回复
-
版面费该交吗
已经有13人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有10人回复
-
什么是人一生最重要的？
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
现在用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人