| 查看: 998 | 回复: 2 | |||
[交流]
用Matlab求反应速率常数已有1人参与
|
|
[kopt fopt]=fminsearch(@calmin,[0,0,0],[],inicon1,t,aobs1); ---------------------------------------- function s=calmin(k,inicon1,t,aobs1) [T,D]=ode45(@kinect,t,inicon1,[],k); s=((D(:,1)-aobs1).^2); --------------------------------------- %动力学方程目的就是求k1、k2、k3% function z=kinect(t,d,k) A=d(1);B=d(2);C=d(3);D=d(4);E=d(5); k1=k(1);k2=k(2);k3=k(3); z=[-k1*A*B-k2*A*C-k3*A*D; -k1*A*B; k1*A*B-k2*A*C; k2*A*C-k3*A*D; k3*A*D;]; ------------------------------- 代码的意义是：通过已知条件，初始浓度，若干个时间内物质A的浓度和反应的方程来求得反应速率常数k1、k2、k3，我运行了一下 报错如下： ??? Subscripted assignment dimension mismatch. Error in ==> fminsearch at 195 fv(:,1) = funfcn(x,varargin{:}); Error in ==> zqm1 at 1 [kopt fopt]=fminsearch(@calmin,[0,0,0],[],inicon1,t,aobs1); [ Last edited by sddtc888 on 2014-4-10 at 13:31 ] |
回复此楼» 猜你喜欢
二苯甲酮亚胺的合成
已经有1人回复
-
聚左旋乳酸的作用与机制：从胶原再生到组织修复
已经有0人回复
-
有机高分子材料论文润色/翻译怎么收费?
已经有267人回复
-
可生物降解聚酯聚己内酯的性能分析与发展趋势
已经有0人回复
-
带膜支架的膜和支架之间的结合力如何测试
已经有2人回复
-
西交利物浦大学奖学金博士招生（生物传感或机器学习方向）
已经有1人回复
-
南方科技大学招收金属材料方向博士生
已经有23人回复
-
可降解聚酯材料在医疗器械中的应用趋势与创新方向
已经有0人回复
-
可降解微球如何提升药物精准治疗效果
已经有0人回复
-
静电纺丝膜分层问题
已经有0人回复
-
什么脱膜剂可以完全清洗干净啊？
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求高版本1stopt跑一个程序
已经有7人回复- 求助用matlab三次样条法算一个反应的反应速率
已经有13人回复
- 求1stopt拟合动力学参数
已经有7人回复
- 符号方程组 4个方程 3个参数 求哥哥姐姐帮忙看一下哈 急~
已经有7人回复
- 未完成的matlab
已经有5人回复
- 求助循环反应反应速率-线性微分方程组求解
已经有16人回复
- 非线性方程组 输出结果太太太长，如何简化显示完全？
已经有7人回复
- 拜求！matlab拟合动力学参数
已经有14人回复
- 求助 平衡常数的问题
已经有8人回复
- 关于酶的Km及IC50
已经有14人回复
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- matlab中什么命令是可以算正交实验数据的k1,k2,k3值的？
已经有1人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- 考博物化5题求解
已经有21人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】催化反应动力学
已经有5人回复
2楼2014-04-10 14:23:25
3楼2014-04-10 14:25:38