| 查看: 948 | 回复: 4 | ||||
[求助]
粘过来的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 ] |
回复此楼» 猜你喜欢
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有247人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有20人回复
-
考博自荐
已经有1人回复
-
山东大学第二批博士研究生招生
已经有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人回复
2楼2012-05-25 00:03:11
877427909
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1146.5
- 帖子: 73
- 在线: 40.8小时
- 虫号: 1161829
- 注册: 2010-12-03
- 性别: GG
- 专业: 机构学与机器人
3楼2012-05-25 22:17:07
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
感谢参与，应助指数 +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]; |
4楼2012-05-26 10:35:40
5楼2012-05-29 16:18:55
5