| 查看: 1335 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求Matlab高手指导 已有1人参与
|
|||
|
Matlab运算结果的置信区间很大，和什么有关系呢？ function KineticsEst6 clear all clc tspan = [0 662.25]; k0 = [0.4587 0.4971 10 12 6 8 205 0.5653]; lb = [0 0 0 0 0 0 0]; ub = [50 50 100 100 50 50 1000 50]; P0 =[0.04020 4.02010 0.18794 0 0; 0.02228 4.45682 0.10418 0 0; 0.01541 4.62428 0.07206 0 0; 0.01336 4.67446 0.06244 0 0; 0.04020 4.02010 0.18794 0 0; 0.04020 4.02010 0.18794 0 0; 0.04020 4.02010 0.18794 0 0 ]; % 初始分压,MPa Pi=[0.01303 4.00374 0.96404 0.007095 0.012097; 0.00393 4.44636 0.53626 0.074475 0.005998; 0.00127 4.61707 0.37150 0.070931 0.003074; 0.00090 4.66795 0.32206 0.066794 0.002414; 0.00243 4.00519 0.96961 0.158768 0.006883; 0.00681 4.00390 0.96766 0.117792 0.009855; 0.01303 4.00374 0.96404 0.070945 0.012097 ]; % 经过Wc/F0后，各物质分压,MPa % 使用函数lsqnonlin()进行参数估计 options=optimset('MaxFunEvals',1000000,'MaxIter',400000) [k,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@ObjFunc,k0,lb,ub,[],P0,Pi); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\tk1 = %.4f ± %.4f\n',k(1),ci(1,2)-k(1)) fprintf('\tk2 = %.4f ± %.4f\n',k(2),ci(2,2)-k(2)) fprintf('\tk3 = %.4f ± %.4f\n',k(3),ci(3,2)-k(3)) fprintf('\tk4 = %.4f ± %.4f\n',k(4),ci(4,2)-k(4)) fprintf('\tk5 = %.4f ± %.4f\n',k(5),ci(5,2)-k(5)) fprintf('\tk6 = %.4f ± %.4f\n',k(6),ci(6,2)-k(6)) fprintf('\tk7 = %.4f ± %.4f\n',k(7),ci(7,2)-k(7)) fprintf('\tk8 = %.4f ± %.4f\n',k(8),ci(8,2)-k(8)) % ------------------------------------------------------------------ function f = ObjFunc(k,P0,Pi) % 目标函数 [m,n] = size(P0); Pcal = zeros(m,n); tspan =[0 662.25]; % 即Wc/F0,g.h/mol for i = 1:m [t PP] = ode45(@Euqations,tspan,P0(i,  ,[],k);Pcal(i, = PP(end,;end f= Pcal-Pi; % ------------------------------------------------------------------ function dPdt = Euqations(t, P, k) % here t = Wc / F0 denom = 1+k(3)*P(1)+k(5)*P(4)+k(5)*P(3)+k(6)*P(5); % k(3) = KDMM, k(4) = KME ,k(5)=KHPM,k(6)=KPDO,k(7)=Kp1,k(8)=Kp2 theA =k(3)*P(1)*P(2)*(1-P(4)*P(3)/k(7)*P(1)*P(2)^2) / denom; theB =k(5)* P(4)*P(2)*(1-P(5)*P(3)/k(8)*P(4)*P(2)^2)/ denom; r1 = k(1)*theA; r2 = k(2)*theB; dPDMMdt = -r1; dPHdt = -2*r1-2*r2; dPMEdt = r1+r2; dPHPMdt = r1-r2; dPPDOdt = r2; dPdt = [dPDMMdt;dPHdt;dPMEdt;dPHPMdt;dPPDOdt]; 使用函数lsqnonlin()估计得到的参数值为: k1 = 0.4599 ± 89509.6910 k2 = 0.4971 ± 101889.0949 k3 = 10.0001 ± 1946900.8019 k4 = 12.0000 ± 8639183.4320 k5 = 5.9997 ± 420897.9438 k6 = 7.9999 ± 2837488.5903 k7 = 205.0000 ± 1472315356.6507 k8 = 0.5653 ± 3696268.6136 结果和我设的初值一样，是不是就没有计算程序。 |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有184人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大师指点（焦急等待指点）
已经有5人回复
- origin8.0 非线性自定义函数拟合，求指导
已经有5人回复
- 跪求高手帮忙看看这组数据如何进行拟合
已经有7人回复
- ELK简单计算失败，求高手解答！
已经有4人回复
- matlab里拟合后，知道y坐标求x坐标
已经有4人回复
- 求在matlab中进行分类的程序
已经有5人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- 求Origin或Matlab高手非线性拟合公式
已经有5人回复
- 求助高手！！！！！！！！！！！matlab最小约束的多变量寻优函数求解方程
已经有5人回复
- 寻求matlab程序高手
已经有3人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- matlab 求解非线性方程 解出虚数解
已经有7人回复
- 学习matlab算法编程，推荐一些有用的资料或者书籍吧
已经有21人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- matlab回归三参数 求高手解救
已经有7人回复
- 【请教】 2D Gaussian Fitting based on Matlab
已经有4人回复
- 求SWMM软件高手指导
已经有15人回复
- matlab 气泡边缘提取 求助高手！！
已经有5人回复
- matlab中dydt的编写要求
已经有12人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
7楼2014-10-25 09:01:06
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-10-24 16:25:38
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2014-10-24 16:27:49
4楼2014-10-24 21:14:08
