| 查看: 1104 | 回复: 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 ] |
回复此楼» 猜你喜欢
LDPE黑十字为什么那么大
已经有0人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有159人回复
-
有机高分子材料论文润色/翻译怎么收费?
已经有256人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有159人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有157人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有152人回复
-
博士申请
已经有15人回复
-
【中科院杭州医学研究所27级博士招生】化学分析/生物材料与器件方向
已经有15人回复
-
静电纺丝膜厚度问题
已经有94人回复
-
常规静电纺丝制备
已经有106人回复
-
PLLA静电纺丝辊筒出现空隙带和纺丝膜透光有圆点，求解答
已经有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
