| 查看: 1038 | 回复: 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); >> |
回复此楼» 猜你喜欢
一志愿哈工大 085600 277 12材科基求调剂
已经有26人回复
-
材料工程281还有调剂机会吗
已经有25人回复
-
339求调剂
已经有8人回复
-
材料考研调剂
已经有28人回复
-
调剂求收留
已经有27人回复
-
0831一轮调剂失败求助
已经有9人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有7人回复
-
求调剂
已经有17人回复
-
322求调剂
已经有6人回复
-
一志愿浙大生物325分求调剂
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
