应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » matlab
北京石油化工学院2026年研究生招生接收调剂公告
131/2返回列表12下一页
查看: 1342  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

1135725495

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

[求助] matlab 已有1人参与

用matlab软件计算置信区间，算出来的结果过大怎么解决？？？@月只蓝
回复此楼

» 猜你喜欢
1楼 2017-04-11 09:38:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

无虫

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
1135725495: 金币+5, ★★★很有帮助 2017-04-11 12:31:43
做参数拟合，首先先看决定系数（或相关系数）的高低，再看其他的。
置信区间跟设置的置信度也有关系，95%的置信度是最常用的，置信度越高，置信区间也越大。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2017-04-11 10:02:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

主管区长

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
2楼: Originally posted by 月只蓝 at 2017-04-11 10:02:44
做参数拟合，首先先看决定系数（或相关系数）的高低，再看其他的。
置信区间跟设置的置信度也有关系，95%的置信度是最常用的，置信度越高，置信区间也越大。

你好！以下是用matlab计算的结果,使用的ci = nlparci(k,residual,jacobian)函数，置信度为95%。算出来的决定系数在0.9以上呢。算出来的置信区间好大啊，不知问题出在哪？
ci =

   1.0e+03 *

   -0.1666    0.1667
   -0.0000    0.0000
   -0.1562    0.1562
   -6.0471    6.0552



以fmincon()的结果为初值，使用函数lsqnonlin()估计得到的参数值为:
        k1 = 0.05881980024
        k2 = 0.00000000000
        k3 = 0.00396318725
        k4 = 4.05955359944
  The sum of the squares is: 6.0e-06


output =
firstorderopt: 6.1534e-07
       iterations: 1
        funcCount: 10
     cgiterations: 0
        algorithm: 'trust-region-reflective'
          message: [1x425 char]
赞一下 回复此楼
3楼2017-04-11 11:16:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

新虫
引用回帖:
3楼: Originally posted by 1135725495 at 2017-04-11 11:16:48
你好！以下是用matlab计算的结果,使用的ci = nlparci(k,residual,jacobian)函数，置信度为95%。算出来的决定系数在0.9以上呢。算出来的置信区间好大啊，不知问题出在哪？
ci =

   1.0e+03 *

   -0.1666     ...

给出具体数据和方程，并指明因变量、自变量、待估计参数
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2017-04-11 11:22:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
4楼: Originally posted by 月只蓝 at 2017-04-11 11:22:39
给出具体数据和方程，并指明因变量、自变量、待估计参数...

下面是程序，方程 r1 =-dc/dt=(k(1)*c^2-0.25*k(2)*(9.404-c)^2)/(1+k(3)*c+0.5*k(4)*(9.404-c))^2，数据是c关于t的方程，因变化量c,自变量是t,待估计参数k(1),k(2),k(3),k(4)
function parafit
%  
% r1 =(k(1)*c^2-0.25*k(2)*(9.404-c)^2)/(1+k(3)*c+0.5*k(4)*(9.404-c))^2;
%
% dCAdt = - r1;
clear all
clc
%        t/min   CA     / mol/L
  Kinetics=[0        8.8599786
3        7.81716904
6        7.3661532
9        7.062404
12        6.82862056
15        6.63696704
18        6.47380764
21        6.33133704
31        5.95301412
41        5.6663802
51        5.43447756
61        5.23953264];
k0 = [0.04539  0.008031  0.07300  3.0230];         % 参数初值
lb = [0 0  0  0];                  % 参数下限
ub =[1  1   1  10];      % 参数上限
x0 = [8.8599786];
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  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('Predictive value','Experimental data')
x


% ------------------------------------------------------------------
function f = ObjFunc7Fmincon(k,x0,yexp)
tspan = [0  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  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)*x^2-0.25*k(2)*(9.404-x)^2)/(1+k(3)*x+0.5*k(4)*(9.404-x))^2;
赞一下 回复此楼
5楼2017-04-11 12:31:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

禁虫

【答案】应助回帖

★ ★ ★ ★ ★
1135725495: 金币+5, ★★★很有帮助 2017-04-11 14:46:28
内容已删除
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2017-04-11 13:48:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

兑换贵宾

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
6楼: Originally posted by 月只蓝 at 2017-04-11 13:48:24
换一下初值试试：

function parafit
%  
% r1 =(k(1)*c^2-0.25*k(2)*(9.404-c)^2)/(1+k(3)*c+0.5*k(4)*(9.404-c))^2;
%
% dCAdt = - r1;
clear all
clc
%        t/min   CA     / mol/L
  Kinetics= ...

你好，谢谢您的帮助，请问这个初值的选择有什么规律吗？我发现初值差一点，都会造成置信区间的范围很大
赞一下 回复此楼
7楼2017-04-11 15:17:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

无虫
引用回帖:
7楼: Originally posted by 1135725495 at 2017-04-11 15:17:47
你好，谢谢您的帮助，请问这个初值的选择有什么规律吗？我发现初值差一点，都会造成置信区间的范围很大...

多试
回复此楼

» 本帖已获得的红花（最新10朵）

1135725495
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
8楼2017-04-11 15:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
送红花一朵
引用回帖:
8楼: Originally posted by 月只蓝 at 2017-04-11 15:19:29
多试...

你好！按照您说的，我试了一下午到现在还是不行，置信区间的数量级没有变化。还有好几个方程需要计算，也是这种问题，就是换一下c的值。你能不能帮我试一下，或者你有没有什么方法可以找到比较好的初值。谢谢！
赞一下 回复此楼
9楼2017-04-11 20:49:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

无虫
引用回帖:
9楼: Originally posted by 1135725495 at 2017-04-11 20:49:19
你好！按照您说的，我试了一下午到现在还是不行，置信区间的数量级没有变化。还有好几个方程需要计算，也是这种问题，就是换一下c的值。你能不能帮我试一下，或者你有没有什么方法可以找到比较好的初值。谢谢！...

待估计参数不超过4个，都可以用5.5版本的1stopt软件计算的结果当初值。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
10楼2017-04-11 21:05:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1135725495 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 348求调剂 +3 zzzzyk123 2026-04-01 3/150 2026-04-01 06:16 by 无际的草原
[考研] 294分080500材料科学与工程求调剂 +12 柳溪边 2026-03-26 13/650 2026-04-01 01:37 by 1018329917
[考研] 311（085601）求调剂 +12 liziyeyeye 2026-03-28 13/650 2026-04-01 00:34 by fmesaito
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +3 快乐储蓄罐 2026-03-31 3/150 2026-03-31 20:20 by plum
[考研] 一志愿北京科技大学材料学硕328分求调剂 +4 1段时间 2026-03-31 5/250 2026-03-31 19:56 by wxiongid
[考研] 一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂! +9 真名有冰 2026-03-29 10/500 2026-03-31 18:01 by xhai2011
[考研] 一志愿中国科学院大学265求调剂 +6 恬淡ye 2026-03-31 7/350 2026-03-31 11:28 by 恬淡ye
[考研] 08工科求调剂286 +5 tgs_001 2026-03-28 5/250 2026-03-31 08:18 by 一只好果子?
[考研] 11408总分309，一志愿东南大学求调剂，不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +5 贝呗钡钡 2026-03-30 5/250 2026-03-30 18:16 by 无际的草原
[考研] 环境科学与工程334分求调剂 +6 王一一依依 2026-03-30 8/400 2026-03-30 11:52 by yjolah
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 340求调剂 +5 jhx777 2026-03-27 5/250 2026-03-28 04:18 by fmesaito
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[考研] 085600，材料与化工321分，求调剂 +9 大馋小子 2026-03-27 9/450 2026-03-27 14:30 by mmm just
[考研] 292求调剂 +4 求求了收下我吧� 2026-03-26 4/200 2026-03-27 10:37 by zhshch
[考研] 0703化学338求调剂！ +6 Zuhui0306 2026-03-26 7/350 2026-03-27 10:35 by shangxh
[考研] 机械学硕310分，数一英一，一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
信息提示
关闭
请填处理意见
关闭