24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

aqua426

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 173
帖子: 9
在线: 7.9小时
虫号: 1522880
注册: 2011-12-04
专业: 工程热物理相关交叉领域

[求助] 求助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函数求解？

回复此楼

» 猜你喜欢

第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有74人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有16人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

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人回复

1楼2013-07-19 10:54:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

somomo91

专家顾问 (职业作家)

专家经验: +3
计算强帖: 1
应助: 144 (高中生)
金币: 19353
红花: 115
帖子: 3854
在线: 1053.2小时
虫号: 2399301
注册: 2013-04-03
专业: 动力学与控制
管辖: 计算模拟

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-23 11:33:33

引用回帖:

5楼: Originally posted by aqua426 at 2013-07-19 14:23:59
好像挺简单谢谢...

用标准的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