| 查看: 12605 | 回复: 1 | ||
xiaofenger01金虫 (初入文坛)
|
[求助]
Matlab错误求助（输入参数的数目不足）
|
|
下面这段程序错在哪里？应该怎么改？ function Cha5demo4 T0 = 180; R = 8.315; k0 = [5.78052E+10 3.92317E+12 1.64254E+4 6.264E+8]; Ea = [124670 150386 77954 111528]; C0 = [1 0 0 0 0]; tspan = [0:1:300]; opt=odeset('reltol',1e-4) [t,T] = ode23(@Ramp,tspan,T0,opt) [t,C] = ode45(@MassEquations,tspan,C0,opt,k0,Ea,R,T) plot(t,C(:,1),'r-',t,C(:,2),'k:',t,C(:,3),'b-.',t,C(:,4),'k--'); xlabel('Time (s)'); ylabel('Concentration (kmol/m^3)'); legend('A','B','C','D') function dTdt = Ramp(t,T) dTdt = 1/6; function dCdt = MassEquations(t,tspan,C,k0,Ea,R,T) for i=1:4 k(i) = k0(i)*exp(-Ea(i)./(R.*(T+273.15))); end rA = -(k(1)+k(2))*C(1); rB = k(1)*C(1)-k(3)*C(2); rC = k(2)*C(1)-k(4)*C(3); rD = k(3)*C(2)-k(4)*C(4); rE = k(4)*C(3); dCdt = [rA; rB; rC; rD; rE]; |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有245人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助一个matlab仿真问题
已经有3人回复
- 用Matlab求一个二阶微分方程，出错了，求指导（求助）
已经有13人回复
- 求助，有关于matlab作图的问题
已经有7人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- matlab “Inner matrix dimensions must agree.”一直不知道哪里出错了 求助
已经有6人回复
- matlab程序问题求解答
已经有10人回复
- matlab解决重心法选址问题，运行不了啊，求助各位大神
已经有10人回复
- 如何在matlab里写.m函数，老是出错 有谁知道具体格式不
已经有4人回复
- matlab 图片叠加问题求助
已经有4人回复
- 求助MATLAB做三维图问题
已经有4人回复
- matlab语句问题
已经有5人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 大家爱帮我看看这段MATLAB程序为啥出错了
已经有5人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- 求助 matlab交点问题
已经有19人回复
- 求助Matlab解方程问题
已经有6人回复
- 求助matlab-simulink运行时报错问题
已经有3人回复
- 【求助】MATLAB导入数据精度问题。。。
已经有4人回复
- 【求助】MATLAB错误
已经有6人回复
- 【求助】程序出错
已经有4人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】请教Matlab高手根据积分画图问题
已经有4人回复
xiaofenger01
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 709.2
- 帖子: 44
- 在线: 93.6小时
- 虫号: 2318750
- 注册: 2013-03-05
- 性别: MM
- 专业: 化学反应工程
2楼2015-04-09 10:34:42
10