24小时热门版块排行榜

返回列表

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

hunter0541

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 1
在线: 1.8小时
虫号: 419436
注册: 2007-07-05
专业: 电工理论

[求助] 粘过来的m文件怎么都运行不了？

我是初学matlab
刚在论坛上看了个程序，如下：
function lee
x0=;
tspan=;
=ode45(@odefun111,tspan,x0);
plot(t,X);
legend('S1','S2','X1','X2')

function F=odefun111(t,x)
S1=x(1);S2=x(2);X1=x(3);X2=x(4);
Q=0.25;Qr=0.25;Qm=0.44;
a=5.654e-6;k=1.6625e-5;
ks=201;Si=405;V10=3;V20=9;
dS1=1/(V10+Q*t)*(Q*Si+Qr*S2-Qm*S1)-a*X1*S1;
dS2=1/(V20+Q*t)*(Qm*S1-Qr*S2)-k*X2*S2*(Si-S1)/(ks+Si-S1);
dX1=1/(V10+Q*t)*(Qr*X2-Qm*X1-Q*X1);
dX2=1/(V20+Q*t)*(Qm*X1-Qr*X2-Q*X2);
怎么我直接粘到matlab editor中就运行不了呢。
一大堆错误。求教了。

[ Last edited by hunter0541 on 2012-5-24 at 22:43 ]

回复此楼

» 猜你喜欢

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

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

K文件怎么运行已经有6人回复
matlab m文件运行出错已经有4人回复
电脑没安装matlab，如何运行M文件已经有3人回复
e-tex文件运行不了已经有8人回复
【求助】弱问下Mathematica里怎么运行*.nb文件已经有6人回复
【求助】Fortran程序，如何让exe文件可以单独运行？已经有6人回复
【求助】xshell的exe文件为什么不能运行？已经有3人回复
【求助】安装ansys时，运行Keygen程序,但不会自动生成一个文件License.dat,为什么? 已经有12人回复
【求助】请问会不会出现在两台电脑上运行同一个模拟文件出现结果不同的情况？已经有6人回复
【求助】关于MBR运行中活性污泥颜色变黄且粘度过大的问题已经有18人回复

1楼2012-05-24 22:25:36

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
cenwanglai: 金币+5, 谢谢~ 2012-05-26 22:23:17

function lee
x0=[2 2 0 0];%需要赋值
tspan=[0 1];%需要赋值
[t,X]=ode45(@odefun111,tspan,x0);
plot(t,X);
legend('S1','S2','X1','X2')

function F=odefun111(t,x)
S1=x(1);S2=x(2);X1=x(3);X2=x(4);
Q=0.25;Qr=0.25;Qm=0.44;
a=5.654e-6;k=1.6625e-5;
ks=201;Si=405;V10=3;V20=9;
dS1=1/(V10+Q*t)*(Q*Si+Qr*S2-Qm*S1)-a*X1*S1;
dS2=1/(V20+Q*t)*(Qm*S1-Qr*S2)-k*X2*S2*(Si-S1)/(ks+Si-S1);
dX1=1/(V10+Q*t)*(Qr*X2-Qm*X1-Q*X1);
dX2=1/(V20+Q*t)*(Qm*X1-Qr*X2-Q*X2);
F=[dS1;dS2;dX1;dX2];