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 ]

回复此楼

» 猜你喜欢

求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有248人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有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];

赞一下

回复此楼

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

4楼2012-05-26 10:35:40

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

相关版块跳转我要订阅楼主 hunter0541 的主题更新

返回列表