24小时热门版块排行榜

>论坛更新日志 (2977)
>虫友互识 (299)
>文献求助 (260)
>导师招生 (169)
>论文投稿 (101)
>休闲灌水 (96)
>考博 (83)
>仿真模拟 (63)
>基金申请 (51)
>硕博家园 (40)
>博后之家 (38)
>找工作 (31)
>招聘信息布告栏 (29)
>公派出国 (25)
>考研 (25)
>SciFinder/Reaxys (24)

返回列表

jmcookie

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 969.9
散金: 74
红花: 1
帖子: 266
在线: 113.8小时
虫号: 601766
注册: 2008-09-13
性别: GG
专业: 食品科学基础

[求助] matalb求解一阶微分方程程序已有1人参与

两个微分方程如下：
dG/dt= VG*180*X；
dX/dt=μ*X；
其中，VG=3.42, 和μ=0.0617是已知的，时间范围为1-60min,
想用Matlab中的ODE45解算两个微分方程，不知道在Matlab中如何进行求解，本人新手，还希望各位能帮忙写个程序代码，在下不胜感激！

回复此楼

» 猜你喜欢

中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有202人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复
大豆异黄酮分离已经有0人回复
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标已经有10人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名已经有1人回复
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有3人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复

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

非线性二阶微分方程组求解，matlab 已经有12人回复
偏微分方程Matlab求解已经有7人回复
matlab解方程组求助攻已经有4人回复
请问什么是分段龙格库塔的方法？需要用matlab软件采用这种方法求二阶微分方程的解已经有3人回复
Matlab数值求解二阶常微分方程已经有9人回复
求助用matlab或者（1st0pt)编程解方程组已经有3人回复
求助，matlab求解二元二阶的常微分方程组已经有9人回复
matalb如何求方程组的积分已经有7人回复
matlab求解常微分方程已经有6人回复
Matlab用龙格库塔法求二阶非线性方程已经有13人回复
微分方程组参数拟合的问题（多参数动力系统参数估计）已经有11人回复
matlab处理常微分方程作图问题已经有7人回复
Matlab求解二阶偏微分方程组，希望能给出相应的求解代码已经有10人回复
分数阶微分方程的求解（用差分法，matlab编程求解）已经有9人回复
求助matlab求解复杂微分方程已经有7人回复
【求助】这个微分方程组simulink求解到底错在哪里？已经有3人回复
【求助】积分微分方程matlab求解已经有6人回复
【求助】急请微分方程高手看下这个方程组！已经有14人回复

1楼 2014-10-18 15:36:10

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68144.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.6小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jmcookie: 金币+10, ★★★★★最佳答案, 非常感谢 2014-10-19 09:18:23

原问题没给出G和X初值，这里按t=0时，二者都等于1取了。如果要修改初值，请自行修改代码中 y0=[1 1] 部分。运行以下代码时，请完整复制，粘贴进一个新建的m文件。

CODE:

function solve_ode

clear all;clc

format long

global VG u

VG=3.42;

u=0.0617;

tspan=1:60;

y0=[1 1];

[t y]=ode45(@odes,tspan,y0);

figure(1),plot(t,y(:,1),'b'),legend('G')

figure(2),plot(t,y(:,2),'r'),legend('X')

function f=odes(t,y)

global VG u

G=y(1);

X=y(2);

f(1)=VG*180*X;

f(2)=u*X;

f=f';