| 查看: 1477 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
1135725495铁杆木虫 (著名写手)
|
[求助]
求大家帮我看看下面中matlab中的程序问题出在哪里？ 已有1人参与
|
|||
|
求动力学参数 方程为：r1 =(k(1)*k(2)*c^2-0.25*k(3)*k(4)*(9.404-c)^2)/(1+k(2)*c+0.5*k(4)*(9.404-c)); t c 0 8.769713879 3 5.909461709 6 5.736047766 9 4.756936544 12 4.505799773 15 4.354160734 18 4.124025932 21 3.994725485 31 3.732949397 41 3.680613046 51 3.537445782 61 3.283577643 求k1,k2,k3,,k4 下面是写的程序，大家帮忙看看，为什么结果是直线呢？ function parafit % % r1 =(k(1)*k(2)*c^2-0.25*k(3)*k(4)*(9.404-c)^2)/(1+k(2)*c+0.5*k(4)*(9.404-c)); % % dCAdt = - r1; clear all clc % t/min CA / mol/L Kinetics=[0 7.029559717 1 3.960706476 2 3.163862085 3 2.963057756 6 2.820434159 9 2.698109429 12 2.556212898 15 2.447100967 18 2.207246045 21 2.13921916 31 2.000403189 41 1.917172016 51 1.851299168 61 1.636287388]; c=Kinetics( :,1)'; k0 = [0.671 0.02253 0.641 0.04028]; % 参数初值 lb = [0 0 0 0]; % 参数下限 ub = [100 100 100 100]; % 参数上限 x0 = [7.029559717]; yexp = Kinetics; warning off % 使用函数fmincon()进行参数估计 [k,fval,flag] = fmincon(@ObjFunc7Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp); fprintf('\n使用函数fmincon()估计得到的参数值为:\n') fprintf('\tk1 = %.11f\n',k(1)) fprintf('\tk2 = %.11f\n',k(2)) fprintf('\tk3 = %.11f\n',k(3)) fprintf('\tk4 = %.11f\n',k(4)) fprintf(' The sum of the squares is: %.1e\n\n',fval) k_fm= k; warning off % 使用函数lsqnonlin()进行参数估计 [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\tk1 = %.11f\n',k(1)) fprintf('\tk2 = %.11f\n',k(2)) fprintf('\tk3 = %.11f\n',k(3)) fprintf('\tk4 = %.11f\n',k(4)) fprintf(' The sum of the squares is: %.1e\n\n',resnorm) k_ls = k; output warning off % 以函数fmincon()估计得到的结果为初值，使用函数lsqnonlin()进行参数估计 k0 = k_fm; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n以fmincon()的结果为初值，使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\tk1 = %.11f\n',k(1)) fprintf('\tk2 = %.11f\n',k(2)) fprintf('\tk3 = %.11f\n',k(3)) fprintf('\tk4 = %.11f\n',k(4)) fprintf(' The sum of the squares is: %.1e\n\n',resnorm) k_fmls = k; output tspan = [0 1 2 3 6 9 12 15 18 21 31 41 51 61]; [t x] = ode45(@KineticEqs,tspan,x0,[],k_fmls); figure; plot(t,x(:,1),t,yexp(:,2),'*');legend('ca-pr','ca-real') x(:,1) figure;plot(t,x); hold on % ------------------------------------------------------------------ function f = ObjFunc7Fmincon(k,x0,yexp) tspan = [0 1 2 3 6 9 12 15 18 21 31 41 51 61]; [t x] = ode45(@KineticEqs,tspan,x0,[],k); y(:,2) = x(:,1); f = sum((y(:,2)-yexp(:,2)).^2) ; % ------------------------------------------------------------------ function f = ObjFunc7LNL(k,x0,yexp) tspan = [0 1 2 3 6 9 12 15 18 21 31 41 51 61]; [t x] = ode45(@KineticEqs,tspan,x0,[],k); y(:,2) = x(:,1); f1 = y(:,2) - yexp(:,2); f = [f1]; % ------------------------------------------------------------------ function dxdt = KineticEqs(t,x,k) Kinetics=[0 7.029559717 1 3.960706476 2 3.163862085 3 2.963057756 6 2.820434159 9 2.698109429 12 2.556212898 15 2.447100967 18 2.207246045 21 2.13921916 31 2.000403189 41 1.917172016 51 1.851299168 61 1.636287388]; c=Kinetics(:,1)'; for i=1:14 c(i)=c(:,i); r1 =(k(1)*k(2)*c(i)^2-0.25*k(3)*k(4)*(9.404-c(i))^2)/(1+k(2)*c(i)+0.5*k(4)*(9.404-c(i))); dCAdt = - r1; dxdt = [dCAdt]; end 结果如下：  H8_DJFZ~6NZC}07J0EXL1]L.png |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有149人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助MATLAB中程序看不懂的简单问题！谢谢！
已经有3人回复
- 如何求这个matlab程序画出的legend曲线的法矢量？
已经有6人回复
- 高分求将我的MATLAB程序改写为FORTRAN程序（1000BB）
已经有13人回复
- 求在matlab中进行分类的程序
已经有5人回复
- 用matlab画一个函数，能帮我看看程序哪里有问题吗
已经有8人回复
- 求教如何用matlab绘制下面两图
已经有5人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- matlab一运行就出现的问题
已经有7人回复
- 有个matlab程序，但是运行不下来，求高手
已经有23人回复
- matlab程序中调用m文件的问题
已经有13人回复
- 帮忙详解示范一下MATLAB中ginput函数的应用吧
已经有3人回复
- 实在是看不出这段matlab程序里哪里出错了，请大家帮忙看看！
已经有10人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
- matlab中遇到的画图小问题，程序没问题。标题要长@@@@@@@@@@@@@@@@@@@@@@
已经有8人回复
- 大家爱帮我看看这段MATLAB程序为啥出错了
已经有5人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
1135725495
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6868.7
- 散金: 18
- 红花: 1
- 帖子: 1492
- 在线: 209.6小时
- 虫号: 2756771
- 注册: 2013-10-27
- 性别: MM
- 专业: 有机化工
送红花一朵 
10楼2015-05-26 17:13:33
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
1135725495: 金币+10, ★有帮助 2015-05-26 15:35:53
感谢参与，应助指数 +1
1135725495: 金币+10, ★有帮助 2015-05-26 15:35:53
|
两组参考结果： Correlation Coef. (R): 0.986609345793065 R-Square: 0.97339800120622 Determination Coef. (DC): 0.973395608126228 F-Statistic: 73.3952829093829 Parameters Best Estimate -------------------- ------------- k1 0.0633244071641639 k2 -0.123667552854961 k3 0.0140433904054541 k4 -0.681989563148558 ====== Output Results ===== File: Data file - 1 No Obs. c Cal. c 1 5.909461709 6.03887691479301 2 5.736047766 5.36051185683508 3 4.756936544 4.91507500618295 4 4.505799773 4.59634816701972 5 4.354160734 4.35722735701414 6 4.124025932 4.17236164824689 7 3.994725485 4.02646466838814 8 3.732949397 3.71546906582311 9 3.680613046 3.55315135651946 10 3.537445782 3.46435427377623 11 3.283577643 3.41460522937653 Correlation Coef. (R): 0.98660934579373 R-Square: 0.973398001207533 Determination Coef. (DC): 0.973395608126228 F-Statistic: 73.3952829093831 Parameters Best Estimate -------------------- ------------- k1 0.481614078106072 k2 1.61663483596038 k3 0.0205007668163867 k4 46.4476555663092 ====== Output Results ===== File: Data file - 1 No Obs. c Cal. c 1 5.909461709 6.03887691237705 2 5.736047766 5.36051185622609 3 4.756936544 4.91507500636805 4 4.505799773 4.5963481675709 5 4.354160734 4.35722735771124 6 4.124025932 4.17236164896553 7 3.994725485 4.02646466905497 8 3.732949397 3.71546906608907 9 3.680613046 3.55315135632873 10 3.537445782 3.46435427319561 11 3.283577643 3.41460522849613 |
2楼2015-05-26 13:56:42
1135725495
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6868.7
- 散金: 18
- 红花: 1
- 帖子: 1492
- 在线: 209.6小时
- 虫号: 2756771
- 注册: 2013-10-27
- 性别: MM
- 专业: 有机化工
您好，您的计算结果我看到了，我想说的是我上面编的Matlab程序正确吗？为什么会出现那样的计算结果？还有就是您计算的方法是最小二乘吗？我现在需要用matlab计算，因为写论文需要。所以您能不能帮忙看一下程序是否正确？ 3楼2015-05-26 15:38:21
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-05-26 16:25:47
100