应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » Matlab错误求助（输入参数的数目不足）
21/1返回列表
查看: 12454  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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];
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-04-06 08:49:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiaofenger01

金虫 (初入文坛)
求助啊，有木有高手来指点一下
赞一下 回复此楼
2楼2015-04-09 10:34:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaofenger01 的主题更新
21/1返回列表
信息提示
关闭
请填处理意见
关闭