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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 求解微分方程组，并输出图像（试写了一下报错）
31/1返回列表
查看: 876  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

whqs8426212

铜虫 (正式写手)

[求助] 求解微分方程组，并输出图像（试写了一下报错）

function dCdt = kineticsEqs(t,C)
k=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ];
dC(1)= -k(1)*C(1)*C(3)+k(2)*C(2)*C(3)-k(6)*C(1)*C(4)-k(7)*C(1)*C(5)+k(8)*C(2)*C(5)-k(9)*C(1)*C(6)+k(10)*C(2)*C(6);
dC(2)= k(1)*C(1)*C(3)-k(2)*C(2)*C(3)+ k(6)*C(1)*C(4)+ k(7)*C(1)*C(5)-k(8)*C(2)*C(5)+ k(9)*C(1)*C(6)-k(10)*C(2)*C(6);
dC(3)= -k(1)*C(1)*C(3)-k(2)*C(2)*C(3)-k(3)*C(4)*C(3)+k(7)*C(5)*C(1)+k(9)*C(1)*C(6)+k(11)*C(4)*c(4)+k(12)*C(5)*C(5)-k(15)*C(5)*C(6);                                                      
dC(4)=k(1)*C(1)*C(3)-k(3)*C(4)*C(3)-k(6)*C(1)*C(4)-k(11)*C(4)*C(4)-k(13)*C(4)*C(5)-k(14)*C(4)*C(6);
dC(5)= k(2)*C(2)*C(3)-k(4)*C(5)+k(5)*0.01*C(6)-k(7)*C(1)*C(5)-k(8)*C(2)*C(5)-k(12)*C(5)*C(5)-k(13)*C(4)*C(5)- k(15)*C(5)*C(6);                                                         
dC(6)= k(4)*C(5)-k(5)*0.01*C(6)-k(9)*C(1)*C(6)-k(10)*C(2)*C(6)-k(14)*0.00000000001*C(6)-k(15)*C(5)*C(6);     
dCdt = [ dC(1); dC(2);dC(3);dC(4); dC(5);dC(6)];
%这是微分方程组

键入M命令解微分方程：
>>tspan = [0 1 2 3 4 5 6 7 8 9 ]*100;
>> x0=[2 0 2 0 0 0];
>> ts=0:1:max(tspan);
>> [ts ys]=ode45(@KineticsEqs,ts,x0)

然后错误(刚刚用matlab希望大家能多多和大家交流)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-08-03 01:06:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenqu

铁杆木虫 (职业作家)

酷酷文曲
内容已删除
回复此楼
不羁年华
2楼2014-08-06 15:25:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqs8426212

铜虫 (正式写手)
内容已删除
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : start2主程序.rar
    • 2014-08-08 01:19:26, 2.45 K
3楼2014-08-08 01:19:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whqs8426212 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭