matlab nlparci,matlab中nlparci函数的调用问题

本文档探讨了在MATLAB中使用nlparci函数计算置信区间时遇到的问题。通过示例代码展示了如何使用fmincon和lsqnonlin进行参数估计,并在估计算法后应用nlparci。尽管得到了参数值,但发现计算出的置信区间过大。寻求对大范围置信区间的解释。
摘要由CSDN通过智能技术生成

请问使用下面的程序计算出的置信区间很大,怎么回事?

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(

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值