| 查看: 4460 | 回复: 12 | |||
[求助]
求助Matlab求解含已知变量的微分方程组（变量随时间变化的离散值已知） 已有2人参与
|
|||
|
微分方程组如下： dy1/dt=-a*y1/y3^m-b*y3*x(t) dy2/dt=-a*y2/y3^m-2*y1*x(t) dy3/dt=(1-y3)/a/y3^m-c/y3^m*(y1*x(t))^0.5 a,b,m是已知常数，x是随时间变化的变量，已知x在tspan内所取时间点的值，能否用ode函数求解？ |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有57人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 求助matlab，这个微分方程能否解出
已经有3人回复
- matlab解偏微分方程求助
已经有12人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- matlab求解偏微分方程组
已经有3人回复
- 求助matlab求解复杂微分方程
已经有7人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68138.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2013-07-19 11:09:21
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-23 11:33:33
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-23 11:33:33
|
用标准的MATLAB 程序的话，如下， 要注意的是，因为我不知道你的参数 x a b m c，随便取了值，所以这个程序不能正确运行，如果你的参数正确的话，应该没有问题的 ================================= clear；clc global x a b m c a = 1; b = 1; c = 1; m = 1; y_ini = [1, 0, 0]; t0 = 0; tend = 1e3; tspan = t0:0.5:tend; x = sin(tspan); [t,y] = ode45(@equ_sys, tspan, y_ini); ============ function [ dy ] = equ_sys( t, y ) global x a b m c dy = zeros(3,1); y1 = y(1) ; y2 = y(2) ; y3 = y(3) ; dy(1) = -a*y1/y3^m-b*y3*x(t) ; dy(2) = -a*y2/y3^m-2*y1*x(t) ; dy(3) = (1-y3)/a/y3^m-c/y3^m*(y1*x(t))^0.5; end |
6楼2013-07-19 23:37:33
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68138.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2013-07-19 10:57:22
3楼2013-07-19 10:59:06
5楼2013-07-19 14:23:59
7楼2015-01-29 21:22:33
8楼2015-01-29 21:25:56
9楼2016-11-22 09:29:41
10楼2018-04-10 00:09:02