【答案】应助回帖

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有80人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有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人回复

1楼2015-05-25 21:08:36

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68121.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2015-05-26 13:56:42

1135725495

铁杆木虫 (著名写手)

应助: 1 (幼儿园)
金币: 6868.7
散金: 18
红花: 1
帖子: 1492
在线: 209.6小时
虫号: 2756771
注册: 2013-10-27
性别: MM
专业: 有机化工

引用回帖:

2楼: Originally posted by 月只蓝 at 2015-05-26 13:56:42
两组参考结果：

Correlation Coef. (R): 0.986609345793065
R-Square: 0.97339800120622
Determination Coef. (DC): 0.973395608126228
F-Statistic: 73.3952829093829
Parameters Best Estimate
--------- ...

您好，您的计算结果我看到了，我想说的是我上面编的Matlab程序正确吗？为什么会出现那样的计算结果？还有就是您计算的方法是最小二乘吗？我现在需要用matlab计算，因为写论文需要。所以您能不能帮忙看一下程序是否正确？

3楼2015-05-26 15:38:21

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68121.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

3楼: Originally posted by 1135725495 at 2015-05-26 15:38:21

您好，您的计算结果我看到了，我想说的是我上面编的Matlab程序正确吗？为什么会出现那样的计算结果？还有就是您计算的方法是最小二乘吗？我现在需要用matlab计算，因为写论文需要。所以您能不能帮忙看一下程 ...

拟合曲线和散点差异明显，可能是因为参数初值选取得不好。

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

4楼2015-05-26 16:25:47

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

引用回帖:

3楼: Originally posted by 1135725495 at 2015-05-26 15:38:21

程序有错，最后面函数写错了

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];
k0 = [0.0671 1.253  0.02 48.028];       % 参数初值
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,'o-',t,yexp(:,2),'*');legend('ca-pr','ca-real')

% ------------------------------------------------------------------
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);
f =  sum((x-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);
f1 = x - yexp(:,2);
f = [f1];
% ------------------------------------------------------------------
function dxdt = KineticEqs(t,x,k)
dxdt=-(k(1)*k(2)*x^2-0.25*k(3)*k(4)*(9.404-x)^2)/(1+k(2)*x+0.5*k(4)*(9.404-x));

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

5楼2015-05-26 16:56:17

1135725495

铁杆木虫 (著名写手)

应助: 1 (幼儿园)
金币: 6868.7
散金: 18
红花: 1
帖子: 1492
在线: 209.6小时
虫号: 2756771
注册: 2013-10-27
性别: MM
专业: 有机化工

引用回帖:

4楼: Originally posted by 月只蓝 at 2015-05-26 16:25:47
拟合曲线和散点差异明显，可能是因为参数初值选取得不好。...

改了初值了，我用别的软件计算出来的结果，用matlab计算出来的图不是一条斜线就是一条直线，感觉程序有问题。但是问题不知道出在哪里？

6楼2015-05-26 17:03:59

1135725495

铁杆木虫 (著名写手)

应助: 1 (幼儿园)
金币: 6868.7
散金: 18
红花: 1
帖子: 1492
在线: 209.6小时
虫号: 2756771
注册: 2013-10-27
性别: MM
专业: 有机化工

引用回帖:

5楼: Originally posted by dbb627 at 2015-05-26 16:56:17
程序有错，最后面函数写错了

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 / m ...

谢谢啊！请问这是改好的吗？我可以直接用吗？

7楼2015-05-26 17:05:15